| 1 | /* Position of planets mercury to neptune; from:
 | 
|---|
| 2 | ftp://ftp.bdl.fr/pub/ephem/planets/vsop87/
 | 
|---|
| 3 | from README:
 | 
|---|
| 4 | 
 | 
|---|
| 5 | ==========================                         ===========================
 | 
|---|
| 6 |                               BUREAU DES LONGITUDES
 | 
|---|
| 7 |                             PLANETARY SOLUTION VSOP87
 | 
|---|
| 8 |                                   1996, January
 | 
|---|
| 9 | ==========================                         ===========================
 | 
|---|
| 10 | 
 | 
|---|
| 11 | These files and programs are associated to :
 | 
|---|
| 12 | 
 | 
|---|
| 13 | Planetary Theories in rectangular and spherical variables: VSOP87 solution.
 | 
|---|
| 14 |     Bretagnon P., Francou G.
 | 
|---|
| 15 |     Astron. Astrophys. 202, 309 (1988).
 | 
|---|
| 16 | 
 | 
|---|
| 17 | Theorie du mouvement de l'ensemble des planetes (VSOP82).
 | 
|---|
| 18 |     Bretagnon P.
 | 
|---|
| 19 |     Astron. Astrophys. 114, 278 (1982).
 | 
|---|
| 20 | 
 | 
|---|
| 21 | ==============================================================================
 | 
|---|
| 22 | 
 | 
|---|
| 23 | Description:
 | 
|---|
| 24 |     The Planetary solutions VSOP87 (Variations Seculaires des Orbites
 | 
|---|
| 25 |     Planetaires) are analytical solutions of the motion of the planets in
 | 
|---|
| 26 |     different versions. The main version VSOP87 consists of the series in
 | 
|---|
| 27 |     elliptic elements as in the case of VSOP82 solution and the other
 | 
|---|
| 28 |     versions VSOP87 (A-B-C-D-E) are built in rectangular and spherical
 | 
|---|
| 29 |     variables.
 | 
|---|
| 30 | 
 | 
|---|
| 31 | Authors' Address:
 | 
|---|
| 32 |     P. Bretagnon, G. Francou
 | 
|---|
| 33 |     Bureau des Longitudes, CNRS URA 707
 | 
|---|
| 34 |     77, Avenue Denfert-Rochereau
 | 
|---|
| 35 |     75014, Paris, France
 | 
|---|
| 36 |     Tel    : (33) 1 40 51 22 69  (33) 1 40 51 22 60
 | 
|---|
| 37 |     Fax    : (33) 1 46 33 28 34
 | 
|---|
| 38 |     E-mail : pierre@bdl.fr  francou@bdl.fr
 | 
|---|
| 39 | 
 | 
|---|
| 40 | Contents:
 | 
|---|
| 41 |     The main version of VSOP87 is similar to the previous theory VSOP82.
 | 
|---|
| 42 |     In the both cases the constants of integration have been determined by
 | 
|---|
| 43 |     fitting to the numerical integration DE200 of the Jet Propulsion
 | 
|---|
| 44 |     Laboratory. The various versions of VSOP87 are different from one to
 | 
|---|
| 45 |     another in the type of coordinates and the reference frame.
 | 
|---|
| 46 |     VSOP87  : heliocentric elliptic    variables; equinox and ecliptic J2000.
 | 
|---|
| 47 |     VSOP87A : heliocentric rectangular variables; equinox and ecliptic J2000.
 | 
|---|
| 48 |     VSOP87B : heliocentric spherical   variables; equinox and ecliptic J2000.
 | 
|---|
| 49 |     VSOP87C : heliocentric rectangular variables; equinox and ecliptic of date.
 | 
|---|
| 50 |     VSOP87D : heliocentric spherical   variables; equinox and ecliptic of date.
 | 
|---|
| 51 |     VSOP87E : barycentric  rectangular variables; equinox and ecliptic J2000.
 | 
|---|
| 52 | ...
 | 
|---|
| 53 | ==============================================================================
 | 
|---|
| 54 | User feed-back is encouraged. Unless otherwise specified, send comments and bug
 | 
|---|
| 55 | reports to:                    E-mail     : comments@bdl.fr
 | 
|---|
| 56 |                                Fax        : (33) 1 46 33 28 34
 | 
|---|
| 57 |                                Postal mail: Bureau des longitudes
 | 
|---|
| 58 |                                             77 avenue Denfert Rochereau
 | 
|---|
| 59 |                                             F-75014 PARIS
 | 
|---|
| 60 | ==============================================================================
 | 
|---|
| 61 |   implemented for C: stern
 | 
|---|
| 62 | */
 | 
|---|
| 63 | 
 | 
|---|
| 64 | #define VSOP_ASCALE     1e8     /* amplitude factor as stored */
 | 
|---|
| 65 | 
 | 
|---|
| 66 | /* coding flags */
 | 
|---|
| 67 | #define VSOP_SPHERICAL  1       /* version in data.c uses spherical coords */
 | 
|---|
| 68 | #define VSOP_GETRATE    0       /* calculate time derivatives of coordinates */
 | 
|---|
| 69 | 
 | 
|---|
| 70 | /* data tables */
 | 
|---|
| 71 | extern double vx_mercury[][3];
 | 
|---|
| 72 | extern int vn_mercury[][3];
 | 
|---|
| 73 | extern double vx_venus[][3];
 | 
|---|
| 74 | extern int vn_venus[][3];
 | 
|---|
| 75 | extern double vx_earth[][3];
 | 
|---|
| 76 | extern int vn_earth[][3];
 | 
|---|
| 77 | extern double vx_mars[][3];
 | 
|---|
| 78 | extern int vn_mars[][3];
 | 
|---|
| 79 | extern double vx_jupiter[][3];
 | 
|---|
| 80 | extern int vn_jupiter[][3];
 | 
|---|
| 81 | extern double vx_saturn[][3];
 | 
|---|
| 82 | extern int vn_saturn[][3];
 | 
|---|
| 83 | extern double vx_uranus[][3];
 | 
|---|
| 84 | extern int vn_uranus[][3];
 | 
|---|
| 85 | extern double vx_neptune[][3];
 | 
|---|
| 86 | extern int vn_neptune[][3];
 | 
|---|
| 87 | 
 | 
|---|
| 88 | extern int vsop87 P_((double mjd, int obj, double prec, double *ret));
 | 
|---|
| 89 | 
 | 
|---|
| 90 | 
 | 
|---|
| 91 | /* For RCS Only -- Do Not Edit
 | 
|---|
| 92 |  * @(#) $RCSfile: vsop87.h,v $ $Date: 2001-10-22 12:08:28 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $
 | 
|---|
| 93 |  */
 | 
|---|