Computation

The program VSOP2013.f (Fortran) computes the planetary elliptic elements between 1890 and 2000 using the files VSOP2013.

This program computes also the planetary rectangular coordinates referred to the dynamical ecliptic frame of J2000 and referred to the equatorial frame of ICRS.

Starting date : 1890 June 26 12h (DJ 2411545.0)

Number of dates : 11

Step : 4000 days.

The results of this computation are given in the file: VSOP2013.outThe file VSOP2013.ctl is given for checking these results.

The mean elements of the elliptic variables of the series VSOP2013 are given in the file: VSOP2013-secular.dat.

The astronomical constants consistent with the solution VSOP2013 are the constants used in the numerical integration INPOP10a, especially:

AU	0.1495978706910000D+09	km
CLIGHT	0.2997924580000000D+06	m/s
GM_Mer	0.4912547451450812D-10	au3/day2
GM_Ven	0.7243452486162703D-09	au3/day2
GM_EMB	0.8997011603631609D-09	au3/day2
GM_Mar	0.9549535105779258D-10	au3/day2
GM_Ves	0.3939673413269574D-13	au3/day2
GM_Iri	0.2299798032187260D-14	au3/day2
GM_Bam	0.1388563508297703D-14	au3/day2
GM_Cer	0.1408056343979966D-12	au3/day2
GM_Pal	0.3296275038741825D-13	au3/day2
GM_Jup	0.2825345842083778D-06	au3/day2
GM_Sat	0.8459715185680659D-07	au3/day2
GM_Ura	0.1292024916781969D-07	au3/day2
GM_Nep	0.1524358900784276D-07	au3/day2
GM_Plu	0.2188699765425970D-11	au3/day2
GM_Sun	0.2959122083684144D-03	au3/day2
EMRAT	0.8130056999999999D+02
RSUN	0.6960000000000000D+06	km
RMOON	0.1738000000000000D+04	km
REARTH	0.6378136600000000D+04	km
J2SUN	0.2400000000000000D-06