Changeset 39457 for trunk/Ohana/src/skycalc
- Timestamp:
- Mar 11, 2016, 10:23:42 PM (10 years ago)
- Location:
- trunk/Ohana
- Files:
-
- 2 edited
-
. (modified) (1 prop)
-
src/skycalc/src/astro.c (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/Ohana
-
Property svn:mergeinfo
set to
/branches/eam_branches/ohana.20160226 merged eligible
-
Property svn:mergeinfo
set to
-
trunk/Ohana/src/skycalc/src/astro.c
r27491 r39457 1 1 # include <skycalc_internal.h> 2 # define UNUSED_PARAM(x)(void)(x) 2 3 3 4 /* Takes a coordinate pair and precesses it using matrix procedures … … 104 105 equatorial (all assumed of date.) */ 105 106 void SC_eclrot(double jd, double *x, double *y, double *z) { 107 UNUSED_PARAM(x); 106 108 107 109 double incl; … … 112 114 113 115 incl = (23.439291 + T * (-0.0130042 - 0.00000016 * T))/DEG_IN_RADIAN; 114 /* 1992 Astron Almanac, p. B18, dropping the115 cubic term, which is 2 milli-arcsec! */116 /* 1992 Astron Almanac, p. B18, dropping the 117 cubic term, which is 2 milli-arcsec! */ 116 118 ypr = cos(incl) * *y - sin(incl) * *z; 117 119 zpr = sin(incl) * *y + cos(incl) * *z; 118 120 *y = ypr; 119 121 *z = zpr; 120 /* x remains the same. */122 /* x remains the same. */ 121 123 } 122 124 123 125 /* Given a julian date in 1900-2100, returns the correction 124 delta t which is:125 TDT - UT (after 1983 and before 1993)126 ET - UT (before 1983)127 an extrapolated guess (after 1993).126 delta t which is: 127 TDT - UT (after 1983 and before 1993) 128 ET - UT (before 1983) 129 an extrapolated guess (after 1993). 128 130 129 For dates in the past (<= 1993) the value is linearly130 interpolated on 5-year intervals; for dates after the present,131 an extrapolation is used, because the true value of delta t132 cannot be predicted precisely. Note that TDT is essentially the133 modern version of ephemeris time with a slightly cleaner134 definition.131 For dates in the past (<= 1993) the value is linearly 132 interpolated on 5-year intervals; for dates after the present, 133 an extrapolation is used, because the true value of delta t 134 cannot be predicted precisely. Note that TDT is essentially the 135 modern version of ephemeris time with a slightly cleaner 136 definition. 135 137 136 Where the algorithm shifts there is an approximately 0.1 second137 discontinuity. Also, the 5-year linear interpolation scheme can138 lead to errors as large as 0.5 seconds in some cases, though139 usually rather smaller. */138 Where the algorithm shifts there is an approximately 0.1 second 139 discontinuity. Also, the 5-year linear interpolation scheme can 140 lead to errors as large as 0.5 seconds in some cases, though 141 usually rather smaller. */ 140 142 141 143 double SC_etcorr (double jd) {
Note:
See TracChangeset
for help on using the changeset viewer.
