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/ |