gust86.h   gust86.h 
skipping to change at line 60 skipping to change at line 60
#ifdef __cplusplus #ifdef __cplusplus
extern "C" { extern "C" {
#endif #endif
#define GUST86_MIRANDA 0 #define GUST86_MIRANDA 0
#define GUST86_ARIEL 1 #define GUST86_ARIEL 1
#define GUST86_UMBRIEL 2 #define GUST86_UMBRIEL 2
#define GUST86_TITANIA 3 #define GUST86_TITANIA 3
#define GUST86_OBERON 4 #define GUST86_OBERON 4
void GetGust86Coor(double jd,int body,double *xyz); void GetGust86Coor(const double jd, const int body, double *xyz);
/* Return the rectangular coordinates of the given satellite /* Return the rectangular coordinates of the given satellite
and the given julian date jd expressed in dynamical time (TAI+32.184s) . and the given julian date jd expressed in dynamical time (TAI+32.184s) .
The origin of the xyz-coordinates is the center of the planet. The origin of the xyz-coordinates is the center of the planet.
The reference frame is "dynamical equinox and ecliptic J2000", The reference frame is "dynamical equinox and ecliptic J2000",
which is the reference frame in VSOP87 and VSOP87A. which is the reference frame in VSOP87 and VSOP87A.
According to vsop87.doc VSOP87 coordinates can be transformed to According to vsop87.doc VSOP87 coordinates can be transformed to
FK5 (=J2000=ICRF for our accuracy requirements) by FK5 (=J2000=ICRF for our accuracy requirements) by
X cos(psi) -sin(psi) 0 1 0 0 X X cos(psi) -sin(psi) 0 1 0 0 X
Y = sin(psi) cos(psi) 0 * 0 cos(eps) -sin(eps) * Y Y = sin(psi) cos(psi) 0 * 0 cos(eps) -sin(eps) * Y
skipping to change at line 89 skipping to change at line 89
ignored in the frame definition. The system is aligned with the ignored in the frame definition. The system is aligned with the
IAU-sponsored J2000 frame of the Radio Source Catalog of the IAU-sponsored J2000 frame of the Radio Source Catalog of the
International Earth Rotational Service (ICRF). International Earth Rotational Service (ICRF).
The ICRF is thought to differ from FK5 by at most 0.01 arcsec. The ICRF is thought to differ from FK5 by at most 0.01 arcsec.
From this I conclude that in the context of stellarium From this I conclude that in the context of stellarium
ICRF, J2000 and FK5 are the same, while the transformation ICRF, J2000 and FK5 are the same, while the transformation
ICRF <-> VSOP87 must be done with the matrix given above. ICRF <-> VSOP87 must be done with the matrix given above.
*/ */
void GetGust86OsculatingCoor(const double jd0, const double jd, const int b void GetGust86OsculatingCoor(const double jd0, const double jd, const int b
ody,double *xyz); ody, double *xyz);
/* The oculating orbit of epoch jd0, evatuated at jd, is returned. /* The oculating orbit of epoch jd0, evaluated at jd, is returned.
*/ */
#ifdef __cplusplus #ifdef __cplusplus
} }
#endif #endif
#endif #endif
 End of changes. 2 change blocks. 
4 lines changed or deleted 4 lines changed or added

This html diff was produced by rfcdiff 1.41. The latest version is available from http://tools.ietf.org/tools/rfcdiff/