mirror
/
XEphem
mirror of https://github.com/XEphem/XEphem.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

Revert "Astronomical Almanac 2020 Pluto elements (#2 #22)" 

This reverts commit 26e5d1abc3.  It broke
Pluto’s motion, which users noticed as a big jump at the dates where we
switch from the short-term solution to these long-term orbital elements.
(Reported to the mailing list on groups.io: thank you, Maxime Gommeaux!)
This commit is contained in:
Brandon Rhodes 2022-01-27 14:27:21 -05:00
parent 5317eae6eb
commit 139b2ea2f1
1 changed files with 9 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
libastro/plans.c
View File

@ -45,17 +45,17 @@ pluto_ell (
double mj, /* epoch */ double mj, /* epoch */
double *ret) /* ecliptic coordinates {l,b,r} at equinox of date */ double *ret) /* ecliptic coordinates {l,b,r} at equinox of date */
{ {
/* Mean orbital elements of Pluto. /* mean orbital elements of Pluto.
* Astronomical Almanac 2020 page G3. * The origin of these is somewhat obscure.
*/ */
double a = 39.789, /* semimajor axis, au */ double a = 39.543, /* semimajor axis, au */
e = 0.252, /* excentricity */ e = 0.2490, /* excentricity */
inc0 = 17.097, /* inclination, deg */ inc0 = 17.140, /* inclination, deg */
Om0 = 110.297, /* long asc node, deg */ Om0 = 110.307, /* long asc node, deg */
omeg0 = 115.058, /* arg of perihel, deg */ omeg0 = 113.768, /* arg of perihel, deg */
mjp = 2459000.5 - MJD0, /* epoch of perihel */ mjp = 2448045.539 - MJD0, /* epoch of perihel */
mjeq = J2000, /* equinox of elements */ mjeq = J2000, /* equinox of elements */
n = 0.0039; /* daily motion, deg */ n = 144.9600/36525.; /* daily motion, deg */
double inc, Om, omeg; /* orbital elements at epoch of date */ double inc, Om, omeg; /* orbital elements at epoch of date */
double ma, ea, nu; /* mean, excentric and true anomaly */ double ma, ea, nu; /* mean, excentric and true anomaly */