+ * EPOCH = 2455493.5 ! 2010-Oct-24.0000000 (TDB) + * 2456805.500000000 = A.D. 2014-May-28 00:00:00.0000 TDB + * EC= 6.410446888795743E-01 QR= 1.243139403466938E+00 IN= 7.040739880944700E+00 + * OM= 5.015120446775977E+01 W = 1.277401708802635E+01 Tp= 2457247.537876891904 + * N = 1.529272576444755E-01 MA= 2.924003597119001E+02 TA= 2.205859235916598E+02 + * A = 3.463214960064703E+00 AD= 5.683290516662468E+00 PR= 2.354060391489698E+03 + */ + private static final double axisCGAU = 3.463214960064703E+00; // Semi-major axis [au] + private static final double eccentricityCG = 6.410446888795743E-01; // Eccentricity [-] + private static final double inclinationCG = 7.040739880944700E+00; // Inclination [degrees] + private static final double argPeriapsisCG = 1.277401708802635E+01; // Arg perifocus [degrees] + private static final double longNodeCG = 5.015120446775977E+01; // Long asc node [degrees] + private static final double periapsisPassageCG = 2457247.537876891904; // Time of periapsis [JD] + private static final double meanMotionCG = 1.529272576444755E-01; // Mean motion [degrees/day] private static final double[] CGORBITPARS = new double[] - {axisCG, eccentricityCG, inclinationCG, argPerihelionCG, longNodeCG, - perihelionPassageCG, meanMotionCG}; - - + {axisCGAU, eccentricityCG, inclinationCG, argPeriapsisCG, longNodeCG, + periapsisPassageCG, meanMotionCG}; + /** * Keplerian orbital parameters for comet D/1993 F2-A (Shoemaker-Levy 9). * https://ssd.jpl.nasa.gov/sbdb.cgi?sstr=1993%20F2-A * https://en.wikipedia.org/wiki/Comet_Shoemaker–Levy_9 - * Comet Shoemaker–Levy 9 (formally designated D/1993 F2) was a comet that - * broke apart in July 1992 and collided with Jupiter in July 1994, providing - * the first direct observation of an extraterrestrial collision of + * Comet Shoemaker–Levy 9 (formally designated D/1993 F2) was a comet that + * broke apart in July 1992 and collided with Jupiter in July 1994, providing + * the first direct observation of an extraterrestrial collision of * Solar System objects. * Orbital parameters are not valid before 1992-Jul-15 00:00 UT * Orbital paramerers are not valid after 1994-Jul-16 20:11 UT * REMARK: Parameters mass and mu = G*M are unknown. + * Orbital Elements at Epoch 2449480.5 (1994-May-08.0) */ - private static final double axisSL9 = 6.86479462772464; // [au] - private static final double eccentricitySL9 = 0.216209166902718; // [-] - private static final double inclinationSL9 = 6.00329387351007; // [degrees] + private static final double axisSL9 = 6.864794627724644; // [au] + private static final double eccentricitySL9 = 0.2162091669027183; // [-] + private static final double inclinationSL9 = 6.003293873510072; // [degrees] private static final double argPerihelionSL9 = 354.8935191875186; // [degrees] private static final double longNodeSL9 = 220.5376550079234; // [degrees] private static final double perihelionPassageSL9 = 2449435.603196492293; // [JED] private static final double meanMotionSL9 = 0.05479775297461272; // [degrees/day] private static final double orbitalPeriodSL9 = 6569.612446823952; // [days] private static final double[] SL9ORBITPARS = new double[] - {axisSL9, eccentricitySL9, inclinationSL9, argPerihelionSL9, longNodeSL9, - perihelionPassageSL9, meanMotionSL9}; - + {axisSL9, eccentricitySL9, inclinationSL9, argPerihelionSL9, longNodeSL9, + perihelionPassageSL9, meanMotionSL9}; + /** * Keplerian orbital parameters for comet C/1995 O1 (Hale-Bopp). * https://ssd.jpl.nasa.gov/sbdb.cgi?sstr=1995%20O1 @@ -670,7 +1665,7 @@ the JPL Solar System Dynamics web site (http://ssd.jpl.nasa.gov/). * of certainty is difficult, Hale–Bopp met or exceeded most predictions * when it passed perihelion on April 1, 1997. It was visible to the naked * eye for a record 18 months, twice as long as the previous record holder, - * the Great Comet of 1811. Accordingly, Hale–Bopp was dubbed the Great + * the Great Comet of 1811. Accordingly, Hale–Bopp was dubbed the Great * Comet of 1997. * REMARK: Parameters mass and mu = G*M are unknown. */ @@ -683,28 +1678,28 @@ the JPL Solar System Dynamics web site (http://ssd.jpl.nasa.gov/). private static final double meanMotionHB = 0.0003733635782842797; // [degrees/day] private static final double orbitalPeriodHB = 964207.6006832551; // [days] private static final double[] HBORBITPARS = new double[] - {axisHB, eccentricityHB, inclinationHB, argPerihelionHB, longNodeHB, - perihelionPassageHB, meanMotionHB}; - + {axisHB, eccentricityHB, inclinationHB, argPerihelionHB, longNodeHB, + perihelionPassageHB, meanMotionHB}; + /** * Keplerian orbital parameters for asteroid 3122 Florence. * https://ssd.jpl.nasa.gov/sbdb.cgi?sstr=3122 * https://echo.jpl.nasa.gov/asteroids/Florence/Florence_planning.html * https://www.scientias.nl/forse-aardscheerder-schiet-op-1-september-planeet/ * https://en.wikipedia.org/wiki/3122_Florence - * 3122 Florence, provisional designation 1981 ET3, is a stony asteroid of + * 3122 Florence, provisional designation 1981 ET3, is a stony asteroid of * the Amor group, classified as near-Earth object and potentially hazardous - * asteroid, approximately 5 kilometers in diameter. It was discovered on + * asteroid, approximately 5 kilometers in diameter. It was discovered on * 2 March 1981 by American astronomer Schelte Bus at Siding Spring Observatory. - * Florence orbits the Sun at a distance of 1.0–2.5 AU once every 2 years and - * 4 months (859 days). Its orbit has an eccentricity of 0.42 and an + * Florence orbits the Sun at a distance of 1.0–2.5 AU once every 2 years and + * 4 months (859 days). Its orbit has an eccentricity of 0.42 and an * inclination of 22° with respect to the ecliptic. - * Florence is classified as a potentially hazardous asteroid (PHA), due to - * both its absolute magnitude (H ≤ 22) and its minimum orbit intersection - * distance (MOID ≤ 0.05 AU). - * On 2017-Sep-01 it will pass 0.04723 AU (7,066,000 km; 4,390,000 mi) - * from Earth, brightening to apparent magnitude 8.5, when it will be - * visible in small telescopes for several nights as it moves through the + * Florence is classified as a potentially hazardous asteroid (PHA), due to + * both its absolute magnitude (H ≤ 22) and its minimum orbit intersection + * distance (MOID ≤ 0.05 AU). + * On 2017-Sep-01 it will pass 0.04723 AU (7,066,000 km; 4,390,000 mi) + * from Earth, brightening to apparent magnitude 8.5, when it will be + * visible in small telescopes for several nights as it moves through the * constellations Piscis Austrinus, Capricornus, Aquarius and Delphinus. * Naming citation was published on 6 April 1993 (M.P.C. 21955). * REMARK: Parameters mass and mu = G*M are unknown. @@ -718,13 +1713,114 @@ the JPL Solar System Dynamics web site (http://ssd.jpl.nasa.gov/). private static final double meanMotionFlorence = 0.418854854065512; // [degrees/day] private static final double orbitalPeriodFlorence = 859.4862790910698; // [days] private static final double[] FLORENCEORBITPARS = new double[] - {axisFlorence, eccentricityFlorence, inclinationFlorence, argPerihelionFlorence, - longNodeFlorence, perihelionPassageFlorence, meanMotionFlorence}; - + {axisFlorence, eccentricityFlorence, inclinationFlorence, argPerihelionFlorence, + longNodeFlorence, perihelionPassageFlorence, meanMotionFlorence}; + + /** + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS + * Ephemeris Type [change] : ELEMENTS + * Target Body [change] : Ultima Thule (3713011) [2486958] + * Observer Location [change] : Sun (body center) [500@10] + * Time Span [change] : Start=2019-01-01, Stop=2019-01-02, Step=30 m + * Table Settings [change] : defaults + * Display/Output [change] : default (formatted HTML) + * + * This pre-computed trajectory is consistent with the New Horizons spacecraft + * Kuiper-Belt extended mission, with the planned 3500 km flyby of Ultima Thule + * on 2019-Jan-1 @ 05:33 UTC. + * + * 2458484.729166667 = A.D. 2019-Jan-01 05:30:00.0000 TDB + * EC= 4.098363091038795E-02 QR= 4.272551661427943E+01 IN= 2.451767154639498E+00 + * OM= 1.589787485075718E+02 W = 1.746745157307338E+02 Tp= 2471795.170112812892 + * N = 3.314454945917162E-03 MA= 3.158831431737062E+02 TA= 3.124912422291994E+02 + * A = 4.455139452399387E+01 AD= 4.637727243370831E+01 PR= 1.086151436282029E+05 + * + private static final double axisUTAU = 4.455139452399387E+01; // Semi-major axis [au] + private static final double eccentricityUT = 4.098363091038795E-02; // Eccentricity [-] + private static final double inclinationUT = 2.451767154639498E+00; // Inclination [degrees] + private static final double argPeriapsisUT = 1.746745157307338E+02; // Arg perifocus [degrees] + private static final double longNodeUT = 1.589787485075718E+02; // Long asc node [degrees] + private static final double periapsisPassageUT = 2471795.170112812892; // Time of periapsis [JD] + private static final double meanMotionUT = 3.314454945917162E-03; // Mean motion [degrees/day] + private static final double[] UTORBITPARS = new double[] + {axisUTAU, eccentricityUT, inclinationUT, argPeriapsisUT, longNodeUT, + periapsisPassageUT, meanMotionUT}; + */ + + /** + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS + * Ephemeris Type [change] : ELEMENTS + * Target Body [change] : Ultima Thule (3713011) [2486958] + * Observer Location [change] : Sun (body center) [500@10] + * Time Span [change] : Start=2006-01-19, Stop=2006-01-20, Step=30 m + * Table Settings [change] : defaults + * Display/Output [change] : default (formatted HTML) + * + * This pre-computed trajectory is consistent with the New Horizons spacecraft + * Kuiper-Belt extended mission, with the planned 3500 km flyby of Ultima Thule + * on 2019-Jan-1 @ 05:33 UTC. + * + * 2453755.291666667 = A.D. 2006-Jan-19 19:00:00.0000 TDB + * EC= 3.697176831510071E-02 QR= 4.274830978442073E+01 IN= 2.451474757830545E+00 + * OM= 1.589482983942412E+02 W = 1.775741605683169E+02 Tp= 2472638.909251091536 + * N = 3.332607599482501E-03 MA= 2.970683125324245E+02 TA= 2.932176608369780E+02 + * A = 4.438946686913731E+01 AD= 4.603062395385389E+01 PR= 1.080235188973050E+05 + * + private static final double axisUTAU = 4.438946686913731E+01; // Semi-major axis [au] + private static final double eccentricityUT = 3.697176831510071E-02; // Eccentricity [-] + private static final double inclinationUT = 2.451474757830545E+00; // Inclination [degrees] + private static final double argPeriapsisUT = 1.775741605683169E+02; // Arg perifocus [degrees] + private static final double longNodeUT = 1.589482983942412E+02; // Long asc node [degrees] + private static final double periapsisPassageUT = 2472638.909251091536; // Time of periapsis [JD] + private static final double meanMotionUT = 3.332607599482501E-03; // Mean motion [degrees/day] + private static final double[] UTORBITPARS = new double[] + {axisUTAU, eccentricityUT, inclinationUT, argPeriapsisUT, longNodeUT, + periapsisPassageUT, meanMotionUT}; + */ + + /** + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS + * Ephemeris Type [change] : ELEMENTS + * Target Body [change] : Ultima Thule (3713011) [2486958] + * Observer Location [change] : Sun (body center) [500@10] + * Time Span [change] : Start=2018-12-01, Stop=2019-01-01, Step=1 MO + * Table Settings [change] : defaults + * Display/Output [change] : default (formatted HTML) + *
+ * This pre-computed trajectory is consistent with the New Horizons spacecraft + * Kuiper-Belt extended mission, with the reconstructed 3537.7 km flyby of + * 2014 MU69 on 2019-Jan-1 @ 05:34:31 TDB (05:33:22 UTC). + *
+ * 2458453.500000000 = A.D. 2018-Dec-01 00:00:00.0000 TDB + * EC= 4.081683015901894E-02 QR= 4.272578650883825E+01 IN= 2.451830205640409E+00 + * OM= 1.589767972590836E+02 W = 1.747632634476068E+02 Tp= 2471820.637855731416 + * N = 3.315288289463411E-03 MA= 3.156840844032503E+02 TA= 3.122948602673823E+02 + * A = 4.454392847189091E+01 AD= 4.636207043494358E+01 PR= 1.085878417102203E+05 + * 2458484.500000000 = A.D. 2019-Jan-01 00:00:00.0000 TDB + * EC= 4.097873475827028E-02 QR= 4.272563856314676E+01 IN= 2.451769758298457E+00 + * OM= 1.589786644650622E+02 W = 1.746720281547497E+02 Tp= 2471794.512226176914 + * N = 3.314466137893298E-03 MA= 3.158844151813905E+02 TA= 3.124930091663442E+02 + * A = 4.455129423264393E+01 AD= 4.637694990214109E+01 PR= 1.086147768668468E+05 + *
+ * Orbital parameter of January 1, 2019 are used to define the orbit of Ultima Thule + */ + private static final double axisUTAU = 4.455129423264393E+01; // Semi-major axis [au] + private static final double eccentricityUT = 4.097873475827028E-02; // Eccentricity [-] + private static final double inclinationUT = 2.451769758298457E+00; // Inclination [degrees] + private static final double argPeriapsisUT = 1.746720281547497E+02; // Arg perifocus [degrees] + private static final double longNodeUT = 1.589786644650622E+02; // Long asc node [degrees] + private static final double periapsisPassageUT = 2471794.512226176914; // Time of periapsis [JD] + private static final double meanMotionUT = 3.314466137893298E-03; // Mean motion [degrees/day] + private static final double[] UTORBITPARS = new double[] + {axisUTAU, eccentricityUT, inclinationUT, argPeriapsisUT, longNodeUT, + periapsisPassageUT, meanMotionUT}; + /** - * Define some of the orbit parameters for moon using data from HORIZONS - * web interface https://ssd.jpl.nasa.gov/horizons.cgi#results. - * Settings: + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS * Ephemeris type : OBSERVER * Target Body : Moon [Luna] [301] * Observer Location : Geocentric [500] @@ -733,153 +1829,1692 @@ the JPL Solar System Dynamics web site (http://ssd.jpl.nasa.gov/). * Display/Output : default (formatted HTML) * NOTE: ORBIT PARAMETERS ARE NOT CORRECTED FOR DATE */ - private static double axisMoonMeter = 3.844E08; // Semi-major axis [m] - private static double eccentricityMoon = 0.05490; // Eccentricity [-] - private static double inclinationMoon = 5.145; // Inclination [degrees] - private static double orbitPeriodMoon = 27.321582; // Orbit period [days] - private static double meanMotionMoon = 360.0/orbitPeriodMoon; // Mean motion [degrees/day] - private static double axisMoonAU = axisMoonMeter / ASTRONOMICALUNIT; + private static final double axisMoonAU = 2.548289534512777E-03; // Semi-major axis [au] + private static final double eccentricityMoon = 6.476694137484437E-02; // Eccentricity [-] + private static final double inclinationMoon = 5.240010960708354E+00; // Inclination [degrees] + private static final double argPeriapsisMoon = 3.081359025079810E+02; // Arg perifocus [degrees] + private static final double longNodeMoon = 1.239837037681769E+02; // Long asc node [degrees] + private static final double periapsisPassageMoon = 2451533.965359302238; // Time of periapsis [JD] + private static final double meanMotionMoon = 1.335975862260855E+01; // Mean motion [degrees/day] private static final double[] MOONORBITPARS = new double[] - {axisMoonAU, eccentricityMoon, inclinationMoon, 0.0, 0.0, 0.0, meanMotionMoon}; - + {axisMoonAU, eccentricityMoon, inclinationMoon, argPeriapsisMoon, longNodeMoon, + periapsisPassageMoon, meanMotionMoon}; + + /** + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS + * Ephemeris Type [change] : ELEMENTS + * Target Body [change] : Io (JI) [501] + * Observer Location [change] : Jupiter System Barycenter [500@5] + * Time Span [change] : Start=2000-01-01, Stop=2000-01-02, Step=1 d + * Table Settings [change] : defaults + * Display/Output [change] : default (formatted HTML) + */ + private static final double axisIoAU = 2.821786546733507E-03; // Semi-major axis [au] + private static final double eccentricityIo = 3.654784965339888E-03; // Eccentricity [-] + private static final double inclinationIo = 2.212609179741271E+00; // Inclination [degrees] + private static final double argPeriapsisIo = 6.218469675691234E+01; // Arg perifocus [degrees] + private static final double longNodeIo = 3.368501231726219E+02; // Long asc node [degrees] + private static final double periapsisPassageIo = 2451545.103514090180; // Time of periapsis [JD] + private static final double meanMotionIo = 2.031615704411821E+02; // Mean motion [degrees/day] + private static final double[] IOORBITPARS = new double[] + {axisIoAU, eccentricityIo, inclinationIo, argPeriapsisIo, longNodeIo, + periapsisPassageIo, meanMotionIo}; + + /** + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS + * Ephemeris Type [change] : ELEMENTS + * Target Body [change] : Europa (JII) [502] + * Observer Location [change] : Jupiter System Barycenter [500@5] + * Time Span [change] : Start=2000-01-01, Stop=2000-01-02, Step=1 d + * Table Settings [change] : defaults + * Display/Output [change] : default (formatted HTML) + */ + private static final double axisEuropaAU = 4.484929379399280E-03; // Semi-major axis [au] + private static final double eccentricityEuropa = 9.470425146083724E-03; // Eccentricity [-] + private static final double inclinationEuropa = 1.790857714257787E+00; // Inclination [degrees] + private static final double argPeriapsisEuropa = 2.557899602714836E+02; // Arg perifocus [degrees] + private static final double longNodeEuropa = 3.326257958798038E+02; // Long asc node [degrees] + private static final double periapsisPassageEuropa = 2451545.154986763373; // Time of periapsis [JD] + private static final double meanMotionEuropa = 1.013931372961153E+02; // Mean motion [degrees/day] + private static final double[] EUROPAORBITPARS = new double[] + {axisEuropaAU, eccentricityEuropa, inclinationEuropa, argPeriapsisEuropa, + longNodeEuropa, periapsisPassageEuropa, meanMotionEuropa}; + + /** + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS + * Ephemeris Type [change] : ELEMENTS + * Target Body [change] : Ganymede (JIII) [503] + * Observer Location [change] : Jupiter System Barycenter [500@5] + * Time Span [change] : Start=2000-01-01, Stop=2000-01-02, Step=1 d + * Table Settings [change] : defaults + * Display/Output [change] : default (formatted HTML) + */ + private static final double axisGanymedeAU = 7.156339844320714E-03; // Semi-major axis [au] + private static final double eccentricityGanymede = 1.318103012416448E-03; // Eccentricity [-] + private static final double inclinationGanymede = 2.214135822185767E+00; // Inclination [degrees] + private static final double argPeriapsisGanymede = 3.167413036642092E+02; // Arg perifocus [degrees] + private static final double longNodeGanymede = 3.431712649776430E+02; // Long asc node [degrees] + private static final double periapsisPassageGanymede = 2451546.588401503861; // Time of periapsis [JD] + private static final double meanMotionGanymede = 5.030036883436198E+01; // Mean motion [degrees/day] + private static final double[] GANYMEDEORBITPARS = new double[] + {axisGanymedeAU, eccentricityGanymede, inclinationGanymede, argPeriapsisGanymede, + longNodeGanymede, periapsisPassageGanymede, meanMotionGanymede}; + + /** + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS + * Ephemeris Type [change] : ELEMENTS + * Target Body [change] : Callisto (JIV) [504] + * Observer Location [change] : Jupiter System Barycenter [500@5] + * Time Span [change] : Start=2000-01-01, Stop=2000-01-02, Step=1 d + * Table Settings [change] : defaults + * Display/Output [change] : default (formatted HTML) + */ + private static final double axisCallistoAU = 1.258560648085115E-02; // Semi-major axis [au] + private static final double eccentricityCallisto = 7.432943295907821E-03; // Eccentricity [-] + private static final double inclinationCallisto = 2.016903900389733E+00; // Inclination [degrees] + private static final double argPeriapsisCallisto = 1.632112921781330E+01; // Arg perifocus [degrees] + private static final double longNodeCallisto = 3.379421848030697E+02; // Long asc node [degrees] + private static final double periapsisPassageCallisto = 2451541.062475862447; // Time of periapsis [JD] + private static final double meanMotionCallisto = 2.156802147671815E+01; // Mean motion [degrees/day] + private static final double[] CALLISTOORBITPARS = new double[] + {axisCallistoAU, eccentricityCallisto, inclinationCallisto, argPeriapsisCallisto, + longNodeCallisto, periapsisPassageCallisto, meanMotionCallisto}; + + /** + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS + * Ephemeris Type [change] : ELEMENTS + * Target Body [change] : Io (JI) [501] + * Observer Location [change] : Jupiter (body center) [500@599] + * Time Span [change] : Start=1979-07-09, Stop=1977-07-10, Step=30 m + * Table Settings [change] : defaults + * Display/Output [change] : default (formatted HTML) + * * + * 2444064.437500000 = A.D. 1979-Jul-09 22:30:00.0000 TDB + * EC= 4.464069846305494E-03 QR= 2.808511485227889E-03 IN= 2.179905298844204E+00 + * OM= 3.378976714091401E+02 W = 2.090801050866877E+02 Tp= 2444064.145622990560 + * N = 2.032332690733699E+02 MA= 5.931911877792535E+01 TA= 5.976031049336378E+01 + * A = 2.821105095417602E-03 AD= 2.833698705607314E-03 PR= 1.771363525476900E+00 + * + private static final double axisIoAU = 2.821105095417602E-03; // Semi-major axis [au] + private static final double eccentricityIo = 4.464069846305494E-03; // Eccentricity [-] + private static final double inclinationIo = 2.179905298844204E+00; // Inclination [degrees] + private static final double argPeriapsisIo = 2.090801050866877E+02; // Arg perifocus [degrees] + private static final double longNodeIo = 3.378976714091401E+02; // Long asc node [degrees] + private static final double periapsisPassageIo = 2444064.145622990560; // Time of periapsis [JD] + private static final double meanMotionIo = 2.032332690733699E+02; // Mean motion [degrees/day] + private static final double[] IOORBITPARS = new double[] + {axisIoAU, eccentricityIo, inclinationIo, argPeriapsisIo, longNodeIo, + periapsisPassageIo, meanMotionIo}; + */ + + /** + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS + * Ephemeris Type [change] : ELEMENTS + * Target Body [change] : Europa (JII) [502] + * Observer Location [change] : Jupiter (body center) [500@599] + * Time Span [change] : Start=1979-07-09, Stop=1977-07-10, Step=30 m + * Table Settings [change] : defaults + * Display/Output [change] : default (formatted HTML) + * 2444064.437500000 = A.D. 1979-Jul-09 22:30:00.0000 TDB + * EC= 9.387957322704407E-03 QR= 4.444762203181402E-03 IN= 2.253032329893990E+00 + * OM= 3.503598726198932E+02 W = 7.745870708442311E+00 Tp= 2444064.095360492356 + * N = 1.013214937274340E+02 MA= 3.466608598681508E+01 TA= 3.528393136496152E+01 + * A = 4.486884887012564E-03 AD= 4.529007570843726E-03 PR= 3.553046710586796E+00 + * + private static final double axisEuropaAU = 4.486884887012564E-03; // Semi-major axis [au] + private static final double eccentricityEuropa = 9.387957322704407E-03; // Eccentricity [-] + private static final double inclinationEuropa = 2.253032329893990E+00; // Inclination [degrees] + private static final double argPeriapsisEuropa = 7.745870708442311E+00; // Arg perifocus [degrees] + private static final double longNodeEuropa = 3.503598726198932E+02; // Long asc node [degrees] + private static final double periapsisPassageEuropa = 2444064.095360492356; // Time of periapsis [JD] + private static final double meanMotionEuropa = 1.013214937274340E+02; // Mean motion [degrees/day] + private static final double[] EUROPAORBITPARS = new double[] + {axisEuropaAU, eccentricityEuropa, inclinationEuropa, argPeriapsisEuropa, + longNodeEuropa, periapsisPassageEuropa, meanMotionEuropa}; + */ + + /** + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS + * Ephemeris Type [change] : ELEMENTS + * Target Body [change] : Ganymede (JIII) [503] + * Observer Location [change] : Jupiter (body center) [500@599] + * Time Span [change] : Start=1979-07-09, Stop=1979-07-10, Step=30 m + * Table Settings [change] : defaults + * Display/Output [change] : default (formatted HTML) + * 2444064.437500000 = A.D. 1979-Jul-09 22:30:00.0000 TDB + * EC= 1.247036336893339E-03 QR= 7.147635725762605E-03 IN= 2.040028856429598E+00 + * OM= 3.417668823947965E+02 W = 2.384131483164773E+02 Tp= 2444061.339921466075 + * N = 5.030068952448489E+01 MA= 1.558103361030715E+02 TA= 1.558688073671426E+02 + * A = 7.156560216399621E-03 AD= 7.165484707036636E-03 PR= 7.156959544754600E+00 + * + private static final double axisGanymedeAU = 7.156560216399621E-03; // Semi-major axis [au] + private static final double eccentricityGanymede = 1.247036336893339E-03; // Eccentricity [-] + private static final double inclinationGanymede = 2.040028856429598E+00; // Inclination [degrees] + private static final double argPeriapsisGanymede = 2.384131483164773E+02; // Arg perifocus [degrees] + private static final double longNodeGanymede = 3.417668823947965E+02; // Long asc node [degrees] + private static final double periapsisPassageGanymede = 2444061.339921466075; // Time of periapsis [JD] + private static final double meanMotionGanymede = 5.030068952448489E+01; // Mean motion [degrees/day] + private static final double[] GANYMEDEORBITPARS = new double[] + {axisGanymedeAU, eccentricityGanymede, inclinationGanymede, argPeriapsisGanymede, + longNodeGanymede, periapsisPassageGanymede, meanMotionGanymede}; + */ + + /** + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS + * Ephemeris Type [change] : ELEMENTS + * Target Body [change] : Callisto (JIV) [504] + * Observer Location [change] : Jupiter (body center) [500@599] + * Time Span [change] : Start=1979-07-09, Stop=1979-07-10, Step=30 m + * Table Settings [change] : defaults + * Display/Output [change] : default (formatted HTML) + * 2444064.437500000 = A.D. 1979-Jul-09 22:30:00.0000 TDB + * EC= 7.835734158106549E-03 QR= 1.249249486072064E-02 IN= 2.072797100534001E+00 + * OM= 3.394822939858797E+02 W = 3.598648006728227E+02 Tp= 2444063.691855364013 + * N = 2.155397593895479E+01 MA= 1.607160654440537E+01 TA= 1.632254239586581E+01 + * A = 1.259115581039419E-02 AD= 1.268981676006773E-02 PR= 1.670225488882388E+01 + * + private static final double axisCallistoAU = 1.259115581039419E-02; // Semi-major axis [au] + private static final double eccentricityCallisto = 7.835734158106549E-03; // Eccentricity [-] + private static final double inclinationCallisto = 2.072797100534001E+00; // Inclination [degrees] + private static final double argPeriapsisCallisto = 3.598648006728227E+02; // Arg perifocus [degrees] + private static final double longNodeCallisto = 3.394822939858797E+02; // Long asc node [degrees] + private static final double periapsisPassageCallisto = 2444063.691855364013; // Time of periapsis [JD] + private static final double meanMotionCallisto = 2.155397593895479E+01; // Mean motion [degrees/day] + private static final double[] CALLISTOORBITPARS = new double[] + {axisCallistoAU, eccentricityCallisto, inclinationCallisto, argPeriapsisCallisto, + longNodeCallisto, periapsisPassageCallisto, meanMotionCallisto}; + */ + + /** FOR TEST PURPOSES + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS + * Ephemeris Type [change] : ELEMENTS + * Target Body [change] : Callisto (JIV) [504] + * Observer Location [change] : Sun (body center) [500@10] + * Time Span [change] : Start=2000-01-01, Stop=2000-01-02, Step=30 m + * Table Settings [change] : defaults + * Display/Output [change] : default (formatted HTML) + * 2451544.500000000 = A.D. 2000-Jan-01 00:00:00.0000 TDB + * EC= 1.463402256336346E+00 QR= 4.810961611782971E+00 IN= 1.173099259407824E+00 + * OM= 1.207096686230750E+02 W = 2.948438791339357E+02 Tp= 2451678.653757016174 + * N = 2.946411930880803E-02 MA=-3.952722302454301E+00 TA= 3.407729998288843E+02 + * A =-1.038182603990403E+01 AD= 6.684586453809735E+91 PR= 1.157407291666667E+95 + * + private static final double axisCallistoAUT = -1.038182603990403E+01; // Semi-major axis [au] + private static final double eccentricityCallistoT = 1.463402256336346E+00; // Eccentricity [-] + private static final double inclinationCallistoT = 1.173099259407824E+00; // Inclination [degrees] + private static final double argPeriapsisCallistoT = 2.948438791339357E+02; // Arg perifocus [degrees] + private static final double longNodeCallistoT = 1.207096686230750E+02; // Long asc node [degrees] + private static final double periapsisPassageCallistoT = 2451678.653757016174; // Time of periapsis [JD] + private static final double meanMotionCallistoT = 2.946411930880803E-02; // Mean motion [degrees/day] + private static final double[] CALLISTOORBITPARSTEST = new double[] + {axisCallistoAUT, eccentricityCallistoT, inclinationCallistoT, argPeriapsisCallistoT, + longNodeCallistoT, periapsisPassageCallistoT, meanMotionCallistoT}; + */ + + /** + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS + * Ephemeris Type [change] : ELEMENTS + * Target Body [change] : Mimas (SI) [606] + * Observer Location [change] : Saturn System Barycenter [500@6] + * Time Span [change] : Start=1980-11-12, Stop=1980-11-12, Step=1 d + * Table Settings [change] : defaults + * Display/Output [change] : default (formatted HTML) + *
+ * Nov 12, 1980 is date of flyby of Voyager 1. + *
+ * 2444555.500000000 = A.D. 1980-Nov-12 00:00:00.0000 TDB + * EC= 2.294076748949336E-02 QR= 1.216152212442891E-03 IN= 2.833548787776371E+01 + * OM= 1.662380994388547E+02 W = 3.110660983844796E+02 Tp= 2444555.468582638539 + * N = 3.794894129614744E+02 MA= 1.192255598685642E+01 TA= 1.248129166805152E+01 + * A = 1.244706740366238E-03 AD= 1.273261268289585E-03 PR= 9.486430654036375E-01 + */ + private static final double axisMimasAU = 1.244706740366238E-03; // Semi-major axis [au] + private static final double eccentricityMimas = 2.294076748949336E-02; // Eccentricity [-] + private static final double inclinationMimas = 2.833548787776371E+01; // Inclination [degrees] + private static final double argPeriapsisMimas = 3.110660983844796E+02; // Arg perifocus [degrees] + private static final double longNodeMimas = 1.662380994388547E+02; // Long asc node [degrees] + private static final double periapsisPassageMimas = 2444555.468582638539; // Time of periapsis [JD] + private static final double meanMotionMimas = 3.794894129614744E+02; // Mean motion [degrees/day] + private static final double[] MIMASORBITPARS = new double[] + {axisMimasAU, eccentricityMimas, inclinationMimas, argPeriapsisMimas, + longNodeMimas, periapsisPassageMimas, meanMotionMimas}; + + /** + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS + * Ephemeris Type [change] : ELEMENTS + * Target Body [change] : Enceladus (SII) [606] + * Observer Location [change] : Saturn System Barycenter [500@6] + * Time Span [change] : Start=1980-11-12, Stop=1980-11-12, Step=1 d + * Table Settings [change] : defaults + * Display/Output [change] : default (formatted HTML) + *
+ * Nov 12, 1980 is date of flyby of Voyager 1. + *
+ * 2444555.500000000 = A.D. 1980-Nov-12 00:00:00.0000 TDB + * EC= 4.483389662009040E-03 QR= 1.584906509048897E-03 IN= 2.805397904088017E+01 + * OM= 1.695244324305456E+02 W = 2.726068154542818E+02 Tp= 2444556.014053585008 + * N = 2.623417662340597E+02 MA= 2.251422746212319E+02 TA= 2.247795280296231E+02 + * A = 1.592044263842871E-03 AD= 1.599182018636845E-03 PR= 1.372255760749931E+00 + */ + private static final double axisEnceladusAU = 1.592044263842871E-03; // Semi-major axis [au] + private static final double eccentricityEnceladus = 4.483389662009040E-03; // Eccentricity [-] + private static final double inclinationEnceladus = 2.805397904088017E+01; // Inclination [degrees] + private static final double argPeriapsisEnceladus = 2.726068154542818E+02; // Arg perifocus [degrees] + private static final double longNodeEnceladus = 1.695244324305456E+02; // Long asc node [degrees] + private static final double periapsisPassageEnceladus = 2444556.014053585008; // Time of periapsis [JD] + private static final double meanMotionEnceladus = 2.623417662340597E+02; // Mean motion [degrees/day] + private static final double[] ENCELADUSORBITPARS = new double[] + {axisEnceladusAU, eccentricityEnceladus, inclinationEnceladus, argPeriapsisEnceladus, + longNodeEnceladus, periapsisPassageEnceladus, meanMotionEnceladus}; + + /** + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS + * Ephemeris Type [change] : ELEMENTS + * Target Body [change] : Tethys (SIII) [606] + * Observer Location [change] : Saturn System Barycenter [500@6] + * Time Span [change] : Start=1980-11-12, Stop=1980-11-12, Step=1 d + * Table Settings [change] : defaults + * Display/Output [change] : default (formatted HTML) + *
+ * Nov 12, 1980 is date of flyby of Voyager 1. + *
+ * 2444555.500000000 = A.D. 1980-Nov-12 00:00:00.0000 TDB + * EC= 1.765221092417035E-03 QR= 1.970712342111551E-03 IN= 2.700778748710897E+01 + * OM= 1.702432284406243E+02 W = 3.229877171493513E+02 Tp= 2444555.393274717033 + * N = 1.899825673459838E+02 MA= 2.027594327456301E+01 TA= 2.034618692144220E+01 + * A = 1.974197236714391E-03 AD= 1.977682131317230E-03 PR= 1.894910701698181E+00 + */ + private static final double axisTethysAU = 1.974197236714391E-03; // Semi-major axis [au] + private static final double eccentricityTethys = 1.765221092417035E-03; // Eccentricity [-] + private static final double inclinationTethys = 2.700778748710897E+01; // Inclination [degrees] + private static final double argPeriapsisTethys = 3.229877171493513E+02; // Arg perifocus [degrees] + private static final double longNodeTethys = 1.702432284406243E+02; // Long asc node [degrees] + private static final double periapsisPassageTethys = 2444555.393274717033; // Time of periapsis [JD] + private static final double meanMotionTethys = 1.899825673459838E+02; // Mean motion [degrees/day] + private static final double[] TETHYSORBITPARS = new double[] + {axisTethysAU, eccentricityTethys, inclinationTethys, argPeriapsisTethys, + longNodeTethys, periapsisPassageTethys, meanMotionTethys}; + + /** + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS + * Ephemeris Type [change] : ELEMENTS + * Target Body [change] : Dione (SIV) [606] + * Observer Location [change] : Saturn System Barycenter [500@6] + * Time Span [change] : Start=1980-11-12, Stop=1980-11-12, Step=1 d + * Table Settings [change] : defaults + * Display/Output [change] : default (formatted HTML) + *
+ * Nov 12, 1980 is date of flyby of Voyager 1. + *
+ * 2444555.500000000 = A.D. 1980-Nov-12 00:00:00.0000 TDB + * EC= 2.175557127959605E-03 QR= 2.514064366570496E-03 IN= 2.803416128061381E+01 + * OM= 1.695809000369315E+02 W = 3.193174176587208E+02 Tp= 2444556.303023209330 + * N = 1.317697228848951E+02 MA= 2.541858542513983E+02 TA= 2.539461679426581E+02 + * A = 2.519545782356523E-03 AD= 2.525027198142550E-03 PR= 2.732038833491903E+00 + */ + private static final double axisDioneAU = 2.519545782356523E-03; // Semi-major axis [au] + private static final double eccentricityDione = 2.175557127959605E-03; // Eccentricity [-] + private static final double inclinationDione = 2.803416128061381E+01; // Inclination [degrees] + private static final double argPeriapsisDione = 3.193174176587208E+02; // Arg perifocus [degrees] + private static final double longNodeDione = 1.695809000369315E+02; // Long asc node [degrees] + private static final double periapsisPassageDione = 2444556.303023209330; // Time of periapsis [JD] + private static final double meanMotionDione = 1.317697228848951E+02; // Mean motion [degrees/day] + private static final double[] DIONEORBITPARS = new double[] + {axisDioneAU, eccentricityDione, inclinationDione, argPeriapsisDione, + longNodeDione, periapsisPassageDione, meanMotionDione}; + + /** + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS + * Ephemeris Type [change] : ELEMENTS + * Target Body [change] : Rhea (SV) [606] + * Observer Location [change] : Saturn System Barycenter [500@6] + * Time Span [change] : Start=1980-11-12, Stop=1980-11-12, Step=1 d + * Table Settings [change] : defaults + * Display/Output [change] : default (formatted HTML) + *
+ * Nov 12, 1980 is date of flyby of Voyager 1. + *
+ * 2444555.500000000 = A.D. 1980-Nov-12 00:00:00.0000 TDB + * EC= 1.077873250919286E-03 QR= 3.519529004414233E-03 IN= 2.776323154892035E+01 + * OM= 1.699392472553282E+02 W = 1.616205502165373E+02 Tp= 2444553.785297347233 + * N = 7.968347273234254E+01 MA= 1.366334620792505E+02 TA= 1.367181923615887E+02 + * A = 3.523326704022749E-03 AD= 3.527124403631264E-03 PR= 4.517875384387965E+00 + */ + private static final double axisRheaAU = 3.523326704022749E-03; // Semi-major axis [au] + private static final double eccentricityRhea = 1.077873250919286E-03; // Eccentricity [-] + private static final double inclinationRhea = 2.776323154892035E+01; // Inclination [degrees] + private static final double argPeriapsisRhea = 1.616205502165373E+02; // Arg perifocus [degrees] + private static final double longNodeRhea = 1.699392472553282E+02; // Long asc node [degrees] + private static final double periapsisPassageRhea = 2444553.785297347233; // Time of periapsis [JD] + private static final double meanMotionRhea = 7.968347273234254E+01; // Mean motion [degrees/day] + private static final double[] RHEAORBITPARS = new double[] + {axisRheaAU, eccentricityRhea, inclinationRhea, argPeriapsisRhea, + longNodeRhea, periapsisPassageRhea, meanMotionRhea}; + + /** + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS + * Ephemeris Type [change] : ELEMENTS + * Target Body [change] : Titan (SVI) [606] + * Observer Location [change] : Saturn System Barycenter [500@6] + * Time Span [change] : Start=1980-11-12, Stop=1980-11-12, Step=1 d + * Table Settings [change] : defaults + * Display/Output [change] : default (formatted HTML) + *
+ * Nov 12, 1980 is date of flyby of Voyager 1. + *
+ * 2444555.500000000 = A.D. 1980-Nov-12 00:00:00.0000 TDB + * EC= 2.883100757260863E-02 QR= 7.930970353173764E-03 IN= 2.774693775242916E+01 + * OM= 1.693222701040926E+02 W = 1.546395634824284E+02 Tp= 2444553.225263091270 + * N = 2.257349567824485E+01 MA= 5.134876377910118E+01 TA= 5.398732093364459E+01 + * A = 8.166416365241102E-03 AD= 8.401862377308443E-03 PR= 1.594790656845183E+01 + */ + private static final double axisTitanAU = 8.166416365241102E-03; // Semi-major axis [au] + private static final double eccentricityTitan = 2.883100757260863E-02; // Eccentricity [-] + private static final double inclinationTitan = 2.774693775242916E+01; // Inclination [degrees] + private static final double argPeriapsisTitan = 1.546395634824284E+02; // Arg perifocus [degrees] + private static final double longNodeTitan = 1.693222701040926E+02; // Long asc node [degrees] + private static final double periapsisPassageTitan = 2444553.225263091270; // Time of periapsis [JD] + private static final double meanMotionTitan = 2.257349567824485E+01; // Mean motion [degrees/day] + private static final double[] TITANORBITPARS = new double[] + {axisTitanAU, eccentricityTitan, inclinationTitan, argPeriapsisTitan, + longNodeTitan, periapsisPassageTitan, meanMotionTitan}; + + /** + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS + * Ephemeris Type [change] : ELEMENTS + * Target Body [change] : Hyperion (SVII) [606] + * Observer Location [change] : Saturn System Barycenter [500@6] + * Time Span [change] : Start=1980-11-12, Stop=1980-11-12, Step=1 d + * Table Settings [change] : defaults + * Display/Output [change] : default (formatted HTML) + *
+ * Nov 12, 1980 is date of flyby of Voyager 1. + *
+ * 2444555.500000000 = A.D. 1980-Nov-12 00:00:00.0000 TDB + * EC= 1.272431943004429E-01 QR= 8.666112626930228E-03 IN= 2.766098902279464E+01 + * OM= 1.679386964197648E+02 W = 1.854300936975505E+02 Tp= 2444562.071466295980 + * N = 1.684238675061948E+01 MA= 2.493208231217351E+02 TA= 2.365127508721484E+02 + * A = 9.929584702560890E-03 AD= 1.119305677819155E-02 PR= 2.137464275879776E+01 + */ + private static final double axisHyperionAU = 9.929584702560890E-03; // Semi-major axis [au] + private static final double eccentricityHyperion = 1.272431943004429E-01; // Eccentricity [-] + private static final double inclinationHyperion = 2.766098902279464E+01; // Inclination [degrees] + private static final double argPeriapsisHyperion = 1.854300936975505E+02; // Arg perifocus [degrees] + private static final double longNodeHyperion = 1.679386964197648E+02; // Long asc node [degrees] + private static final double periapsisPassageHyperion = 2444562.071466295980; // Time of periapsis [JD] + private static final double meanMotionHyperion = 1.684238675061948E+01; // Mean motion [degrees/day] + private static final double[] HYPERIONORBITPARS = new double[] + {axisHyperionAU, eccentricityHyperion, inclinationHyperion, argPeriapsisHyperion, + longNodeHyperion, periapsisPassageHyperion, meanMotionHyperion}; + + /** + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS + * Ephemeris Type [change] : ELEMENTS + * Target Body [change] : Iapetus (SVIII) [606] + * Observer Location [change] : Saturn System Barycenter [500@6] + * Time Span [change] : Start=1980-11-12, Stop=1980-11-12, Step=1 d + * Table Settings [change] : defaults + * Display/Output [change] : default (formatted HTML) + *
+ * Nov 12, 1980 is date of flyby of Voyager 1. + *
+ * 2444555.500000000 = A.D. 1980-Nov-12 00:00:00.0000 TDB + * EC= 2.856225061595786E-02 QR= 2.312483296093841E-02 IN= 1.749827907600319E+01 + * OM= 1.403973722109327E+02 W = 2.261903994839276E+02 Tp= 2444517.435074096080 + * N = 4.537354975447288E+00 MA= 1.727140809408426E+02 TA= 1.731149457095495E+02 + * A = 2.380475020205992E-02 AD= 2.448466744318143E-02 PR= 7.934137883151000E+01 + */ + private static final double axisIapetusAU = 2.380475020205992E-02; // Semi-major axis [au] + private static final double eccentricityIapetus = 2.856225061595786E-02; // Eccentricity [-] + private static final double inclinationIapetus = 1.749827907600319E+01; // Inclination [degrees] + private static final double argPeriapsisIapetus = 2.261903994839276E+02; // Arg perifocus [degrees] + private static final double longNodeIapetus = 1.403973722109327E+02; // Long asc node [degrees] + private static final double periapsisPassageIapetus = 2444517.435074096080; // Time of periapsis [JD] + private static final double meanMotionIapetus = 4.537354975447288E+00; // Mean motion [degrees/day] + private static final double[] IAPETUSORBITPARS = new double[] + {axisIapetusAU, eccentricityIapetus, inclinationIapetus, argPeriapsisIapetus, + longNodeIapetus, periapsisPassageIapetus, meanMotionIapetus}; + + /******************************************************************************* + https://ssd.jpl.nasa.gov/horizons/app.html#/ + Ephemeris / WWW_USER Tue Dec 28 02:57:07 2021 Pasadena, USA / Horizons + ******************************************************************************* + Target body name: Phoebe (609) {source: sat427l_merged_DE438} + Center body name: Saturn Barycenter (6) {source: sat427l_merged_DE438} + Center-site name: BODY CENTER + ******************************************************************************* + Start time : A.D. 1980-Nov-12 00:00:00.0000 TDB + Stop time : A.D. 1980-Nov-13 00:00:00.0000 TDB + Step-size : 1440 minutes + NOTE: Nov 12, 1980 is date of flyby of Voyager 1 + ******************************************************************************* + Center geodetic : 0.00000000,0.00000000,0.0000000 {E-lon(deg),Lat(deg),Alt(km)} + Center cylindric: 0.00000000,0.00000000,0.0000000 {E-lon(deg),Dxy(km),Dz(km)} + Center radii : (undefined) + Keplerian GM : 8.4597056200480170E-08 au^3/d^2 + Output units : AU-D, deg, Julian Day Number (Tp) + Output type : GEOMETRIC osculating elements + Output format : 10 + Reference frame : Ecliptic of J2000.0 + ******************************************************************************* + JDTDB + EC QR IN + OM W Tp + N MA TA + A AD PR + ******************************************************************************* + 2444555.500000000 = A.D. 1980-Nov-12 00:00:00.0000 TDB + EC= 1.855076035617937E-01 QR= 7.010169845245748E-02 IN= 1.733587164130361E+02 + OM= 2.574241853583012E+02 W = 3.474806387916817E+02 Tp= 2444303.988814169541 + N = 6.599909779942061E-01 MA= 1.659951135126032E+02 TA= 1.701863482979023E+02 + A = 8.606795933149750E-02 AD= 1.020342202105375E-01 PR= 5.454620017596063E+02 + *******************************************************************************/ + private static final double axisPhoebeAU = 8.606795933149750E-02; // Semi-major axis [au] + private static final double eccentricityPhoebe = 1.855076035617937E-01; // Eccentricity [-] + private static final double inclinationPhoebe = 1.733587164130361E+02; // Inclination [degrees] + private static final double argPeriapsisPhoebe = 3.474806387916817E+02; // Arg perifocus [degrees] + private static final double longNodePhoebe = 2.574241853583012E+02; // Long asc node [degrees] + private static final double periapsisPassagePhoebe = 2444303.988814169541; // Time of periapsis [JD] + private static final double meanMotionPhoebe = 6.599909779942061E-01; // Mean motion [degrees/day] + private static final double[] PHOEBEORBITPARS = new double[] + {axisPhoebeAU, eccentricityPhoebe, inclinationPhoebe, argPeriapsisPhoebe, + longNodePhoebe, periapsisPassagePhoebe, meanMotionPhoebe}; + + /** + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS + * Ephemeris Type [change] : ELEMENTS + * Target Body [change] : Miranda (UV) [705] + * Center : Uranus System Barycenter [500@7] + * Time Span [change] : Start=1986-01-24, Stop=1986-01-24, Step=6 h + * Table Settings [change] : defaults + * Display/Output [change] : default (formatted HTML) + *
+ * Jan 24, 1986 is date of flyby of Voyager 2. + *
+ * 2446455.250000000 = A.D. 1986-Jan-24 18:00:00.0000 TDB + * EC= 1.334113739191030E-03 QR= 8.664871813807257E-04 IN= 1.019201605356678E+02 + * OM= 1.690199433071529E+02 W = 3.223465505057740E+02 Tp= 2446455.060737770516 + * N = 2.548265562097606E+02 MA= 4.822904210236982E+01 TA= 4.834318747087165E+01 + * A = 8.676447181199060E-04 AD= 8.688022548590863E-04 PR= 1.412725601893964E+00 + */ + private static final double axisMirandaAU = 8.676447181199060E-04; // Semi-major axis [au] + private static final double eccentricityMiranda = 1.334113739191030E-03; // Eccentricity [-] + private static final double inclinationMiranda = 1.019201605356678E+02; // Inclination [degrees] + private static final double argPeriapsisMiranda = 3.223465505057740E+02; // Arg perifocus [degrees] + private static final double longNodeMiranda = 1.690199433071529E+02; // Long asc node [degrees] + private static final double periapsisPassageMiranda = 2446455.060737770516; // Time of periapsis [JD] + private static final double meanMotionMiranda = 2.548265562097606E+02; // Mean motion [degrees/day] + private static final double[] MIRANDAORBITPARS = new double[] + {axisMirandaAU, eccentricityMiranda, inclinationMiranda, argPeriapsisMiranda, + longNodeMiranda, periapsisPassageMiranda, meanMotionMiranda}; + + /** + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS + * Ephemeris Type [change] : ELEMENTS + * Target Body [change] : Ariel (UI) [701] + * Center : Uranus System Barycenter [500@7] + * Time Span [change] : Start=1986-01-24, Stop=1986-01-24, Step=6 h + * Table Settings [change] : defaults + * Display/Output [change] : default (formatted HTML) + *
+ * Jan 24, 1986 is date of flyby of Voyager 2. + *
+ * 2446455.250000000 = A.D. 1986-Jan-24 18:00:00.0000 TDB + * EC= 1.828146549847673E-03 QR= 1.273608250361064E-03 IN= 9.771964547613906E+01 + * OM= 1.676311983381491E+02 W = 3.268039389316039E+02 Tp= 2446454.735053625889 + * N = 1.428903374678187E+02 MA= 7.358086117257147E+01 TA= 7.378193771139888E+01 + * A = 1.275940857237032E-03 AD= 1.278273464113000E-03 PR= 2.519414583096482E+00 + */ + private static final double axisArielAU = 1.275940857237032E-03; // Semi-major axis [au] + private static final double eccentricityAriel = 1.828146549847673E-03; // Eccentricity [-] + private static final double inclinationAriel = 9.771964547613906E+01; // Inclination [degrees] + private static final double argPeriapsisAriel = 3.268039389316039E+02; // Arg perifocus [degrees] + private static final double longNodeAriel = 1.676311983381491E+02; // Long asc node [degrees] + private static final double periapsisPassageAriel = 2446454.735053625889; // Time of periapsis [JD] + private static final double meanMotionAriel = 1.428903374678187E+02; // Mean motion [degrees/day] + private static final double[] ARIELORBITPARS = new double[] + {axisArielAU, eccentricityAriel, inclinationAriel, argPeriapsisAriel, + longNodeAriel, periapsisPassageAriel, meanMotionAriel}; + + /** + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS + * Ephemeris Type [change] : ELEMENTS + * Target Body [change] : Umbriel (UII) [702] + * Center : Uranus System Barycenter [500@7] + * Time Span [change] : Start=1986-01-24, Stop=1986-01-24, Step=6 h + * Table Settings [change] : defaults + * Display/Output [change] : default (formatted HTML) + *
+ * Jan 24, 1986 is date of flyby of Voyager 2. + *
+ * 2446455.250000000 = A.D. 1986-Jan-24 18:00:00.0000 TDB + * EC= 4.264444554427566E-03 QR= 1.770865367553330E-03 IN= 9.769328000824999E+01 + * OM= 1.675902215804727E+02 W = 2.949950690755990E+02 Tp= 2446456.513212617952 + * N = 8.683327893081464E+01 MA= 2.503111063928964E+02 TA= 2.498518364737965E+02 + * A = 1.778449466696910E-03 AD= 1.786033565840491E-03 PR= 4.145875917997222E+00 + */ + private static final double axisUmbrielAU = 1.778449466696910E-03; // Semi-major axis [au] + private static final double eccentricityUmbriel = 4.264444554427566E-03; // Eccentricity [-] + private static final double inclinationUmbriel = 9.769328000824999E+01; // Inclination [degrees] + private static final double argPeriapsisUmbriel = 2.949950690755990E+02; // Arg perifocus [degrees] + private static final double longNodeUmbriel = 1.675902215804727E+02; // Long asc node [degrees] + private static final double periapsisPassageUmbriel = 2446456.513212617952; // Time of periapsis [JD] + private static final double meanMotionUmbriel = 8.683327893081464E+01; // Mean motion [degrees/day] + private static final double[] UMBRIELORBITPARS = new double[] + {axisUmbrielAU, eccentricityUmbriel, inclinationUmbriel, argPeriapsisUmbriel, + longNodeUmbriel, periapsisPassageUmbriel, meanMotionUmbriel}; + + /** + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS + * Ephemeris Type [change] : ELEMENTS + * Target Body [change] : Titania (UIII) [703] + * Center : Uranus System Barycenter [500@7] + * Time Span [change] : Start=1986-01-24, Stop=1986-01-24, Step=6 h + * Table Settings [change] : defaults + * Display/Output [change] : default (formatted HTML) + *
+ * Jan 24, 1986 is date of flyby of Voyager 2. + *
+ * 2446455.250000000 = A.D. 1986-Jan-24 18:00:00.0000 TDB + * EC= 2.699684405758504E-03 QR= 2.907991492514086E-03 IN= 9.784186655153967E+01 + * OM= 1.676200525421978E+02 W = 1.913287392441945E+02 Tp= 2446458.712360218167 + * N = 4.136011948310419E+01 MA= 2.167963676760602E+02 TA= 2.166115685681959E+02 + * A = 2.915863403473766E-03 AD= 2.923735314433446E-03 PR= 8.704036750838251E+00 + */ + private static final double axisTitaniaAU = 2.915863403473766E-03; // Semi-major axis [au] + private static final double eccentricityTitania = 2.699684405758504E-03; // Eccentricity [-] + private static final double inclinationTitania = 9.784186655153967E+01; // Inclination [degrees] + private static final double argPeriapsisTitania = 1.913287392441945E+02; // Arg perifocus [degrees] + private static final double longNodeTitania = 1.676200525421978E+02; // Long asc node [degrees] + private static final double periapsisPassageTitania = 2446458.712360218167; // Time of periapsis [JD] + private static final double meanMotionTitania = 4.136011948310419E+01; // Mean motion [degrees/day] + private static final double[] TITANIAORBITPARS = new double[] + {axisTitaniaAU, eccentricityTitania, inclinationTitania, argPeriapsisTitania, + longNodeTitania, periapsisPassageTitania, meanMotionTitania}; + + /** + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS + * Ephemeris Type [change] : ELEMENTS + * Target Body [change] : Oberon (UIV) [704] + * Center : Uranus System Barycenter [500@7] + * Time Span [change] : Start=1986-01-24, Stop=1986-01-24, Step=6 h + * Table Settings [change] : defaults + * Display/Output [change] : default (formatted HTML) + *
+ * Jan 24, 1986 is date of flyby of Voyager 2. + *
+ * 2446455.250000000 = A.D. 1986-Jan-24 18:00:00.0000 TDB + * EC= 3.479588686913168E-04 QR= 3.899090046075020E-03 IN= 9.783394065511513E+01 + * OM= 1.677619229560692E+02 W = 3.089743224447537E+02 Tp= 2446454.454100156669 + * N = 2.673396445911760E+01 MA= 2.127755813031945E+01 TA= 2.129203341544224E+01 + * A = 3.900447241284487E-03 AD= 3.901804436493955E-03 PR= 1.346601625623177E+01 + */ + private static final double axisOberonAU = 3.900447241284487E-03; // Semi-major axis [au] + private static final double eccentricityOberon = 3.479588686913168E-04; // Eccentricity [-] + private static final double inclinationOberon = 9.783394065511513E+01; // Inclination [degrees] + private static final double argPeriapsisOberon = 3.089743224447537E+02; // Arg perifocus [degrees] + private static final double longNodeOberon = 1.677619229560692E+02; // Long asc node [degrees] + private static final double periapsisPassageOberon = 2446454.454100156669; // Time of periapsis [JD] + private static final double meanMotionOberon = 2.673396445911760E+01; // Mean motion [degrees/day] + private static final double[] OBERONORBITPARS = new double[] + {axisOberonAU, eccentricityOberon, inclinationOberon, argPeriapsisOberon, + longNodeOberon, periapsisPassageOberon, meanMotionOberon}; + + /** + * https://ssd.jpl.nasa.gov/horizons.cgi#results + * Results from HORIZONS + * Ephemeris Type [change] : ELEMENTS + * Target Body [change] : Triton (NI) [801] + * Center : Neptune System Barycenter [500@8] + * Time Span [change] : Start=1989-08-25, Stop=1989-08-26, Step=6 h + * Table Settings [change] : defaults + * Display/Output [change] : default (formatted HTML) + *
+ * Aug 25, 1989 is date of flyby of Voyager 2. + *
+ * 2447763.500000000 = A.D. 1989-Aug-25 00:00:00.0000 TDB
+ * EC= 1.217220754327991E-05 QR= 2.370934201605266E-03 IN= 1.309092040930801E+02
+ * OM= 2.131742023027621E+02 W = 2.781135619426180E+02 Tp= 2447763.663316546939
+ * N = 6.125517670911833E+01 MA= 3.499960160500586E+02 TA= 3.499957737408385E+02
+ * A = 2.370963061459727E-03 AD= 2.370991921314189E-03 PR= 5.877054305296145E+00
+ */
+ private static final double axisTritonAU = 2.370963061459727E-03; // Semi-major axis [au]
+ private static final double eccentricityTriton = 1.217220754327991E-05; // Eccentricity [-]
+ private static final double inclinationTriton = 1.309092040930801E+02; // Inclination [degrees]
+ private static final double argPeriapsisTriton = 2.781135619426180E+02; // Arg perifocus [degrees]
+ private static final double longNodeTriton = 2.131742023027621E+02; // Long asc node [degrees]
+ private static final double periapsisPassageTriton = 2447763.663316546939; // Time of periapsis [JD]
+ private static final double meanMotionTriton = 6.125517670911833E+01; // Mean motion [degrees/day]
+ private static final double[] TRITONORBITPARS = new double[]
+ {axisTritonAU, eccentricityTriton, inclinationTriton, argPeriapsisTriton,
+ longNodeTriton, periapsisPassageTriton, meanMotionTriton};
+
+ /*******************************************************************************
+ https://ssd.jpl.nasa.gov/horizons/app.html#/
+ Ephemeris / WWW_USER Thu Dec 30 07:57:55 2021 Pasadena, USA / Horizons
+ *******************************************************************************
+ Target body name: Nereid (802) {source: nep097}
+ Center body name: Neptune Barycenter (8) {source: DE441}
+ Center-site name: BODY CENTER
+ *******************************************************************************
+ Start time : A.D. 1989-Aug-25 00:00:00.0000 TDB
+ Stop time : A.D. 1989-Aug-26 00:00:00.0000 TDB
+ Step-size : 360 minutes
+ *******************************************************************************
+ Center geodetic : 0.00000000,0.00000000,0.0000000 {E-lon(deg),Lat(deg),Alt(km)}
+ Center cylindric: 0.00000000,0.00000000,0.0000000 {E-lon(deg),Dxy(km),Dz(km)}
+ Center radii : (undefined)
+ Keplerian GM : 1.5243573478851939E-08 au^3/d^2
+ Output units : AU-D, deg, Julian Day Number (Tp)
+ Output type : GEOMETRIC osculating elements
+ Output format : 10
+ Reference frame : Ecliptic of J2000.0
+ *******************************************************************************
+ JDTDB
+ EC QR IN
+ OM W Tp
+ N MA TA
+ A AD PR
+ *******************************************************************************
+ 2447763.500000000 = A.D. 1989-Aug-25 00:00:00.0000 TDB
+ EC= 7.542482788502861E-01 QR= 9.057091227970883E-03 IN= 5.054626152808997E+00
+ OM= 3.198523853673589E+02 W = 2.965172967734192E+02 Tp= 2447727.526164759882
+ N = 9.998336879159416E-01 MA= 3.596785235641266E+01 TA= 1.306064118236043E+02
+ A = 3.685464006355112E-02 AD= 6.465218889913137E-02 PR= 3.600598823094126E+02
+ *******************************************************************************/
+ private static final double axisNereidAU = 3.685464006355112E-02; // Semi-major axis [au]
+ private static final double eccentricityNereid = 7.542482788502861E-01; // Eccentricity [-]
+ private static final double inclinationNereid = 5.054626152808997E+00; // Inclination [degrees]
+ private static final double argPeriapsisNereid = 2.965172967734192E+02; // Arg perifocus [degrees]
+ private static final double longNodeNereid = 3.198523853673589E+02; // Long asc node [degrees]
+ private static final double periapsisPassageNereid = 2447727.526164759882; // Time of periapsis [JD]
+ private static final double meanMotionNereid = 9.998336879159416E-01; // Mean motion [degrees/day]
+ private static final double[] NEREIDORBITPARS = new double[]
+ {axisNereidAU, eccentricityNereid, inclinationNereid, argPeriapsisNereid,
+ longNodeNereid, periapsisPassageNereid, meanMotionNereid};
+
+ /*******************************************************************************
+ https://ssd.jpl.nasa.gov/horizons/app.html#/
+ Ephemeris / WWW_USER Thu Dec 30 08:05:02 2021 Pasadena, USA / Horizons
+ *******************************************************************************
+ Target body name: Proteus (808) {source: nep097}
+ Center body name: Neptune Barycenter (8) {source: DE441}
+ Center-site name: BODY CENTER
+ *******************************************************************************
+ Start time : A.D. 1989-Aug-25 00:00:00.0000 TDB
+ Stop time : A.D. 1989-Aug-26 00:00:00.0000 TDB
+ Step-size : 360 minutes
+ *******************************************************************************
+ Center geodetic : 0.00000000,0.00000000,0.0000000 {E-lon(deg),Lat(deg),Alt(km)}
+ Center cylindric: 0.00000000,0.00000000,0.0000000 {E-lon(deg),Dxy(km),Dz(km)}
+ Center radii : (undefined)
+ Keplerian GM : 1.5243556197896240E-08 au^3/d^2
+ Output units : AU-D, deg, Julian Day Number (Tp)
+ Output type : GEOMETRIC osculating elements
+ Output format : 10
+ Reference frame : Ecliptic of J2000.0
+ *******************************************************************************
+ JDTDB
+ EC QR IN
+ OM W Tp
+ N MA TA
+ A AD PR
+ *******************************************************************************
+ 2447763.500000000 = A.D. 1989-Aug-25 00:00:00.0000 TDB
+ EC= 6.179682141605470E-04 QR= 7.865777304444289E-04 IN= 2.892163404392493E+01
+ OM= 4.814562822915969E+01 W = 1.162927241864745E+02 Tp= 2447763.503025105689
+ N = 3.203690961820102E+02 MA= 3.590308496370286E+02 TA= 3.590296509601386E+02
+ A = 7.870641110475628E-04 AD= 7.875504916506969E-04 PR= 1.123703891200150E+00
+ *******************************************************************************/
+ private static final double axisProteusAU = 3.685464006355112E-02; // Semi-major axis [au]
+ private static final double eccentricityProteus = 7.542482788502861E-01; // Eccentricity [-]
+ private static final double inclinationProteus = 5.054626152808997E+00; // Inclination [degrees]
+ private static final double argPeriapsisProteus = 2.965172967734192E+02; // Arg perifocus [degrees]
+ private static final double longNodeProteus = 3.198523853673589E+02; // Long asc node [degrees]
+ private static final double periapsisPassageProteus = 2447727.526164759882; // Time of periapsis [JD]
+ private static final double meanMotionProteus = 9.998336879159416E-01; // Mean motion [degrees/day]
+ private static final double[] PROTEUSORBITPARS = new double[]
+ {axisProteusAU, eccentricityProteus, inclinationProteus, argPeriapsisProteus,
+ longNodeProteus, periapsisPassageProteus, meanMotionProteus};
+
+ /*******************************************************************************
+ Revised: Aug 19, 2015 134340 Pluto 999
+
+ Pre-computed solution PLU043/DE430 reversion; fit to data through mid-2015.
+
+ PHYSICAL DATA (updated 2021-Jun-07; Mc= Charon mass, radius is IAU 2009):
+ Mass x10^22 (kg) = 1.307+-0.018 Volume, 10^10 km^3 = 0.697
+ GM (planet) km^3/s^2 = 869.96 Density (R=1195 km) = 1.86 g/cm^3
+ GM 1-sigma, km^3/s^2 = 0.08 Surface gravity = 0.611 m/s^2
+ Vol. mean radius (km) = 1188.3+-1.6 Mass ratio (Mc/Mp) = 0.117
+ Sidereal rot. period = 153.29335198 h Sid. rot. rat, rad/s = 0.0000113856
+ Mean solar day, h = 153.2820 Mean orbit velocity = 4.67 km/s
+ Sidereal orbit period = 249.58932 yr Escape speed, km/s = 1.21
+ Perihelion Aphelion Mean
+ Solar Constant (W/m^2) 1.56 0.56 0.88
+ Maximum Planetary IR (W/m^2) 0.8 0.3 0.5
+ Minimum Planetary IR (W/m^2) 0.8 0.3 0.5
+ *******************************************************************************
+ *******************************************************************************
+ Ephemeris / WWW_USER Sun Nov 21 04:48:35 2021 Pasadena, USA / Horizons
+ *******************************************************************************
+ Target body name: Pluto (999) {source: plu043_merged}
+ Center body name: Pluto Barycenter (9) {source: plu043_merged}
+ Center-site name: BODY CENTER
+ *******************************************************************************
+ Start time : A.D. 2015-Jul-01 00:00:00.0000 TDB
+ Stop time : A.D. 2015-Jul-02 00:00:00.0000 TDB
+ Step-size : 60 minutes
+ *******************************************************************************
+ Center geodetic : 0.00000000,0.00000000,0.0000000 {E-lon(deg),Lat(deg),Alt(km)}
+ Center cylindric: 0.00000000,0.00000000,0.0000000 {E-lon(deg),Dxy(km),Dz(km)}
+ Center radii : (undefined)
+ Keplerian GM : 2.8034366453558377E-15 au^3/d^2
+ Output units : AU-D, deg, Julian Day Number (Tp)
+ Output type : GEOMETRIC osculating elements
+ Output format : 10
+ Reference frame : Ecliptic of J2000.0
+ *******************************************************************************
+ JDTDB
+ EC QR IN
+ OM W Tp
+ N MA TA
+ A AD PR
+ *******************************************************************************
+ $$SOE
+ 2457204.500000000 = A.D. 2015-Jul-01 00:00:00.0000 TDB
+ EC= 7.744741847030238E-03 QR= 1.400195619029954E-05 IN= 1.128726302861868E+02
+ OM= 2.274086440915950E+02 W = 3.447560503405338E+02 Tp= 2457207.641122277826
+ N = 5.722951259150261E+01 MA= 1.802351030427898E+02 TA= 1.802314963573664E+02
+ A = 1.411124413325199E-05 AD= 1.422053207620445E-05 PR= 6.290460702848140E+00
+ *******************************************************************************/
+ private static final double axisPlutoAU = 1.411124413325199E-05; // Semi-major axis [au]
+ private static final double eccentricityPluto = 7.744741847030238E-03; // Eccentricity [-]
+ private static final double inclinationPluto = 1.128726302861868E+02; // Inclination [degrees]
+ private static final double argPeriapsisPluto = 3.447560503405338E+02; // Arg perifocus [degrees]
+ private static final double longNodePluto = 2.274086440915950E+02; // Long asc node [degrees]
+ private static final double periapsisPassagePluto = 2457207.641122277826; // Time of periapsis [JD]
+ private static final double meanMotionPluto = 5.722951259150261E+01; // Mean motion [degrees/day]
+ private static final double[] PLUTOORBITPARSRELATIVE = new double[]
+ {axisPlutoAU, eccentricityPluto, inclinationPluto, argPeriapsisPluto,
+ longNodePluto, periapsisPassagePluto, meanMotionPluto};
+
+ /*******************************************************************************
+ https://ssd.jpl.nasa.gov/horizons/app.html#/
+ Ephemeris / WWW_USER Sat Nov 20 05:44:40 2021 Pasadena, USA / Horizons
+ *******************************************************************************
+ Target body name: Charon (901) {source: plu043_merged}
+ Center body name: Pluto Barycenter (9) {source: plu043_merged}
+ Center-site name: BODY CENTER
+ *******************************************************************************
+ Start time : A.D. 2015-Jul-01 00:00:00.0000 TDB
+ Stop time : A.D. 2015-Jul-02 00:00:00.0000 TDB
+ Step-size : 60 minutes
+ *******************************************************************************
+ Center geodetic : 0.00000000,0.00000000,0.0000000 {E-lon(deg),Lat(deg),Alt(km)}
+ Center cylindric: 0.00000000,0.00000000,0.0000000 {E-lon(deg),Dxy(km),Dz(km)}
+ Center radii : (undefined)
+ Keplerian GM : 1.5389673638556205E-12 au^3/d^2
+ Output units : AU-D, deg, Julian Day Number (Tp)
+ Output type : GEOMETRIC osculating elements
+ Output format : 10
+ Reference frame : Ecliptic of J2000.0
+ *******************************************************************************
+ JDTDB
+ EC QR IN
+ OM W Tp
+ N MA TA
+ A AD PR
+ *******************************************************************************
+ 2457204.500000000 = A.D. 2015-Jul-01 00:00:00.0000 TDB
+ EC= 1.209835207171253E-03 QR= 1.167769278684351E-04 IN= 1.128726338127190E+02
+ OM= 2.274086243380147E+02 W = 3.430150233953275E+02 Tp= 2457204.464913493488
+ N = 5.622288855037032E+01 MA= 1.972664746750097E+00 TA= 1.977444225564479E+00
+ A = 1.169183798407318E-04 AD= 1.170598318130285E-04 PR= 6.403086167966530E+00
+ *******************************************************************************/
+ private static final double axisCharonAU = 1.169183798407318E-04; // Semi-major axis [au]
+ private static final double eccentricityCharon = 1.209835207171253E-03; // Eccentricity [-]
+ private static final double inclinationCharon = 1.128726338127190E+02; // Inclination [degrees]
+ private static final double argPeriapsisCharon = 3.430150233953275E+02; // Arg perifocus [degrees]
+ private static final double longNodeCharon = 2.274086243380147E+02; // Long asc node [degrees]
+ private static final double periapsisPassageCharon = 2457204.464913493488; // Time of periapsis [JD]
+ private static final double meanMotionCharon = 5.622288855037032E+01; // Mean motion [degrees/day]
+ private static final double[] CHARONORBITPARS = new double[]
+ {axisCharonAU, eccentricityCharon, inclinationCharon, argPeriapsisCharon,
+ longNodeCharon, periapsisPassageCharon, meanMotionCharon};
+
+ /*******************************************************************************
+ https://ssd.jpl.nasa.gov/horizons/app.html#/
+ Ephemeris / WWW_USER Mon Dec 27 03:00:39 2021 Pasadena, USA / Horizons
+ *******************************************************************************
+ Target body name: Nix (902) {source: plu043_merged}
+ Center body name: Pluto Barycenter (9) {source: plu043_merged}
+ Center-site name: BODY CENTER
+ *******************************************************************************
+ Start time : A.D. 2015-Jul-01 00:00:00.0000 TDB
+ Stop time : A.D. 2015-Jul-02 00:00:00.0000 TDB
+ Step-size : 1440 minutes
+ *******************************************************************************
+ Center geodetic : 0.00000000,0.00000000,0.0000000 {E-lon(deg),Lat(deg),Alt(km)}
+ Center cylindric: 0.00000000,0.00000000,0.0000000 {E-lon(deg),Dxy(km),Dz(km)}
+ Center radii : (undefined)
+ Keplerian GM : 2.1750820795064077E-12 au^3/d^2
+ Output units : AU-D, deg, Julian Day Number (Tp)
+ Output type : GEOMETRIC osculating elements
+ Output format : 10
+ Reference frame : Ecliptic of J2000.0
+ *******************************************************************************
+ JDTDB
+ EC QR IN
+ OM W Tp
+ N MA TA
+ A AD PR
+ *******************************************************************************
+ 2457204.500000000 = A.D. 2015-Jul-01 00:00:00.0000 TDB
+ EC= 6.907136142536396E-03 QR= 3.242328766008259E-04 IN= 1.128878636501151E+02
+ OM= 2.274275512295376E+02 W = 2.629010752961010E+02 Tp= 2457197.616501887795
+ N = 1.432384705760017E+01 MA= 9.859817417836248E+01 TA= 9.937974495630276E+01
+ A = 3.264879734825708E-04 AD= 3.287430703643157E-04 PR= 2.513291286568057E+01
+ *******************************************************************************/
+ private static final double axisNixAU = 3.264879734825708E-04; // Semi-major axis [au]
+ private static final double eccentricityNix = 6.907136142536396E-03; // Eccentricity [-]
+ private static final double inclinationNix = 1.128878636501151E+02; // Inclination [degrees]
+ private static final double argPeriapsisNix = 2.629010752961010E+02; // Arg perifocus [degrees]
+ private static final double longNodeNix = 2.274275512295376E+02; // Long asc node [degrees]
+ private static final double periapsisPassageNix = 2457197.616501887795; // Time of periapsis [JD]
+ private static final double meanMotionNix = 1.432384705760017E+01; // Mean motion [degrees/day]
+ private static final double[] NIXORBITPARS = new double[]
+ {axisNixAU, eccentricityNix, inclinationNix, argPeriapsisNix,
+ longNodeNix, periapsisPassageNix, meanMotionNix};
+
+ /*******************************************************************************
+ https://ssd.jpl.nasa.gov/horizons/app.html#/
+ Ephemeris / WWW_USER Mon Dec 27 03:07:42 2021 Pasadena, USA / Horizons
+ *******************************************************************************
+ Target body name: Hydra (903) {source: plu043_merged}
+ Center body name: Pluto Barycenter (9) {source: plu043_merged}
+ Center-site name: BODY CENTER
+ *******************************************************************************
+ Start time : A.D. 2015-Jul-01 00:00:00.0000 TDB
+ Stop time : A.D. 2015-Jul-02 00:00:00.0000 TDB
+ Step-size : 1440 minutes
+ *******************************************************************************
+ Center geodetic : 0.00000000,0.00000000,0.0000000 {E-lon(deg),Lat(deg),Alt(km)}
+ Center cylindric: 0.00000000,0.00000000,0.0000000 {E-lon(deg),Dxy(km),Dz(km)}
+ Center radii : (undefined)
+ Keplerian GM : 2.1750719564057657E-12 au^3/d^2
+ Output units : AU-D, deg, Julian Day Number (Tp)
+ Output type : GEOMETRIC osculating elements
+ Output format : 10
+ Reference frame : Ecliptic of J2000.0
+ *******************************************************************************
+ JDTDB
+ EC QR IN
+ OM W Tp
+ N MA TA
+ A AD PR
+ *******************************************************************************
+ 2457204.500000000 = A.D. 2015-Jul-01 00:00:00.0000 TDB
+ EC= 1.311241066181010E-02 QR= 4.300411732442002E-04 IN= 1.128180809725054E+02
+ OM= 2.277384657437538E+02 W = 3.225272550766830E+02 Tp= 2457206.570122104604
+ N = 9.289576619947059E+00 MA= 3.407694420955177E+02 TA= 3.402667716739093E+02
+ A = 4.357549713768184E-04 AD= 4.414687695094366E-04 PR= 3.875311165710064E+01
+ *******************************************************************************/
+ private static final double axisHydraAU = 4.357549713768184E-04; // Semi-major axis [au]
+ private static final double eccentricityHydra = 1.311241066181010E-02; // Eccentricity [-]
+ private static final double inclinationHydra = 1.128180809725054E+02; // Inclination [degrees]
+ private static final double argPeriapsisHydra = 3.225272550766830E+02; // Arg perifocus [degrees]
+ private static final double longNodeHydra = 2.277384657437538E+02; // Long asc node [degrees]
+ private static final double periapsisPassageHydra = 2457206.570122104604; // Time of periapsis [JD]
+ private static final double meanMotionHydra = 9.289576619947059E+00; // Mean motion [degrees/day]
+ private static final double[] HYDRAORBITPARS = new double[]
+ {axisHydraAU, eccentricityHydra, inclinationHydra, argPeriapsisHydra,
+ longNodeHydra, periapsisPassageHydra, meanMotionHydra};
+
+ /*******************************************************************************
+ https://ssd.jpl.nasa.gov/horizons/app.html#/
+ Ephemeris / WWW_USER Mon Dec 27 03:12:09 2021 Pasadena, USA / Horizons
+ *******************************************************************************
+ Target body name: Kerberos (904) {source: plu043_merged}
+ Center body name: Pluto Barycenter (9) {source: plu043_merged}
+ Center-site name: BODY CENTER
+ *******************************************************************************
+ Start time : A.D. 2015-Jul-01 00:00:00.0000 TDB
+ Stop time : A.D. 2015-Jul-02 00:00:00.0000 TDB
+ Step-size : 1440 minutes
+ *******************************************************************************
+ Center geodetic : 0.00000000,0.00000000,0.0000000 {E-lon(deg),Lat(deg),Alt(km)}
+ Center cylindric: 0.00000000,0.00000000,0.0000000 {E-lon(deg),Dxy(km),Dz(km)}
+ Center radii : (undefined)
+ Keplerian GM : 2.1750934275185082E-12 au^3/d^2
+ Output units : AU-D, deg, Julian Day Number (Tp)
+ Output type : GEOMETRIC osculating elements
+ Output format : 10
+ Reference frame : Ecliptic of J2000.0
+ *******************************************************************************
+ JDTDB
+ EC QR IN
+ OM W Tp
+ N MA TA
+ A AD PR
+ *******************************************************************************
+ $$SOE
+ 2457204.500000000 = A.D. 2015-Jul-01 00:00:00.0000 TDB
+ EC= 1.400088452161128E-02 QR= 3.857110920101890E-04 IN= 1.129962142318857E+02
+ OM= 2.278774257927614E+02 W = 2.828746358131037E+02 Tp= 2457206.379643523134
+ N = 1.092152253453862E+01 MA= 3.394714309036033E+02 TA= 3.388994559448614E+02
+ A = 3.911880710187546E-04 AD= 3.966650500273201E-04 PR= 3.296243713836812E+01
+ *******************************************************************************/
+ private static final double axisKerberosAU = 3.911880710187546E-04; // Semi-major axis [au]
+ private static final double eccentricityKerberos = 1.400088452161128E-02; // Eccentricity [-]
+ private static final double inclinationKerberos = 1.129962142318857E+02; // Inclination [degrees]
+ private static final double argPeriapsisKerberos = 2.828746358131037E+02; // Arg perifocus [degrees]
+ private static final double longNodeKerberos = 2.278774257927614E+02; // Long asc node [degrees]
+ private static final double periapsisPassageKerberos = 2457206.379643523134; // Time of periapsis [JD]
+ private static final double meanMotionKerberos = 1.092152253453862E+01; // Mean motion [degrees/day]
+ private static final double[] KERBEROSORBITPARS = new double[]
+ {axisKerberosAU, eccentricityKerberos, inclinationKerberos, argPeriapsisKerberos,
+ longNodeKerberos, periapsisPassageKerberos, meanMotionKerberos};
+
+ /*******************************************************************************
+ https://ssd.jpl.nasa.gov/horizons/app.html#/
+ Ephemeris / WWW_USER Mon Dec 27 03:15:37 2021 Pasadena, USA / Horizons
+ *******************************************************************************
+ Target body name: Styx (905) {source: plu043_merged}
+ Center body name: Pluto Barycenter (9) {source: plu043_merged}
+ Center-site name: BODY CENTER
+ *******************************************************************************
+ Start time : A.D. 2015-Jul-01 00:00:00.0000 TDB
+ Stop time : A.D. 2015-Jul-02 00:00:00.0000 TDB
+ Step-size : 1440 minutes
+ *******************************************************************************
+ Center geodetic : 0.00000000,0.00000000,0.0000000 {E-lon(deg),Lat(deg),Alt(km)}
+ Center cylindric: 0.00000000,0.00000000,0.0000000 {E-lon(deg),Dxy(km),Dz(km)}
+ Center radii : (undefined)
+ Keplerian GM : 2.1750964648933585E-12 au^3/d^2
+ Output units : AU-D, deg, Julian Day Number (Tp)
+ Output type : GEOMETRIC osculating elements
+ Output format : 10
+ Reference frame : Ecliptic of J2000.0
+ *******************************************************************************
+ JDTDB
+ EC QR IN
+ OM W Tp
+ N MA TA
+ A AD PR
+ *******************************************************************************
+ $$SOE
+ 2457204.500000000 = A.D. 2015-Jul-01 00:00:00.0000 TDB
+ EC= 2.810881721111728E-02 QR= 2.823335773580959E-04 IN= 1.128512646171002E+02
+ OM= 2.271145505212035E+02 W = 1.093419872582783E+02 Tp= 2457204.978474052157
+ N = 1.706650374625700E+01 MA= 3.518341207991095E+02 TA= 3.513601516816546E+02
+ A = 2.904991652953653E-04 AD= 2.986647532326348E-04 PR= 2.109395136534363E+01
+ *******************************************************************************/
+ private static final double axisStyxAU = 2.904991652953653E-04; // Semi-major axis [au]
+ private static final double eccentricityStyx = 2.810881721111728E-02; // Eccentricity [-]
+ private static final double inclinationStyx = 1.128512646171002E+02; // Inclination [degrees]
+ private static final double argPeriapsisStyx = 1.093419872582783E+02; // Arg perifocus [degrees]
+ private static final double longNodeStyx = 2.271145505212035E+02; // Long asc node [degrees]
+ private static final double periapsisPassageStyx = 2457204.978474052157; // Time of periapsis [JD]
+ private static final double meanMotionStyx = 1.706650374625700E+01; // Mean motion [degrees/day]
+ private static final double[] STYXORBITPARS = new double[]
+ {axisStyxAU, eccentricityStyx, inclinationStyx, argPeriapsisStyx,
+ longNodeStyx, periapsisPassageStyx, meanMotionStyx};
+
+ /*******************************************************************************
+ https://ssd.jpl.nasa.gov/horizons/app.html#/
+ Ephemeris / WWW_USER Tue Dec 28 01:37:30 2021 Pasadena, USA / Horizons
+ *******************************************************************************
+ Target body name: Phobos (401) {source: mar097}
+ Center body name: Mars Barycenter (4) {source: mar097}
+ Center-site name: BODY CENTER
+ *******************************************************************************
+ Start time : A.D. 2000-Jan-01 00:00:00.0000 TDB
+ Stop time : A.D. 2000-Jan-02 00:00:00.0000 TDB
+ Step-size : 1440 minutes
+ *******************************************************************************
+ Center geodetic : 0.00000000,0.00000000,0.0000000 {E-lon(deg),Lat(deg),Alt(km)}
+ Center cylindric: 0.00000000,0.00000000,0.0000000 {E-lon(deg),Dxy(km),Dz(km)}
+ Center radii : (undefined)
+ Keplerian GM : 9.5495483556275204E-11 au^3/d^2
+ Output units : AU-D, deg, Julian Day Number (Tp)
+ Output type : GEOMETRIC osculating elements
+ Output format : 10
+ Reference frame : Ecliptic of J2000.0
+ *******************************************************************************
+ JDTDB
+ EC QR IN
+ OM W Tp
+ N MA TA
+ A AD PR
+ *******************************************************************************
+ 2451544.500000000 = A.D. 2000-Jan-01 00:00:00.0000 TDB
+ EC= 1.541574629185453E-02 QR= 6.172355997200113E-05 IN= 2.605134469886070E+01
+ OM= 8.481060425430211E+01 W = 3.423766042007272E+02 Tp= 2451544.512579276226
+ N = 1.128020331783272E+03 MA= 3.458103207628809E+02 TA= 3.453690572568199E+02
+ A = 6.268997268597137E-05 AD= 6.365638539994159E-05 PR= 3.191431837322302E-01
+ *******************************************************************************/
+ private static final double axisPhobosAU = 6.268997268597137E-05; // Semi-major axis [au]
+ private static final double eccentricityPhobos = 1.541574629185453E-02; // Eccentricity [-]
+ private static final double inclinationPhobos = 2.605134469886070E+01; // Inclination [degrees]
+ private static final double argPeriapsisPhobos = 3.423766042007272E+02; // Arg perifocus [degrees]
+ private static final double longNodePhobos = 8.481060425430211E+01; // Long asc node [degrees]
+ private static final double periapsisPassagePhobos = 2451544.512579276226; // Time of periapsis [JD]
+ private static final double meanMotionPhobos = 1.128020331783272E+03; // Mean motion [degrees/day]
+ private static final double[] PHOBOSORBITPARS = new double[]
+ {axisPhobosAU, eccentricityPhobos, inclinationPhobos, argPeriapsisPhobos,
+ longNodePhobos, periapsisPassagePhobos, meanMotionPhobos};
+
+ /*******************************************************************************
+ https://ssd.jpl.nasa.gov/horizons/app.html#/
+ Ephemeris / WWW_USER Tue Dec 28 01:42:07 2021 Pasadena, USA / Horizons
+ *******************************************************************************
+ Target body name: Deimos (402) {source: mar097}
+ Center body name: Mars Barycenter (4) {source: mar097}
+ Center-site name: BODY CENTER
+ *******************************************************************************
+ Start time : A.D. 2000-Jan-01 00:00:00.0000 TDB
+ Stop time : A.D. 2000-Jan-02 00:00:00.0000 TDB
+ Step-size : 1440 minutes
+ *******************************************************************************
+ Center geodetic : 0.00000000,0.00000000,0.0000000 {E-lon(deg),Lat(deg),Alt(km)}
+ Center cylindric: 0.00000000,0.00000000,0.0000000 {E-lon(deg),Dxy(km),Dz(km)}
+ Center radii : (undefined)
+ Keplerian GM : 9.5495487654055916E-11 au^3/d^2
+ Output units : AU-D, deg, Julian Day Number (Tp)
+ Output type : GEOMETRIC osculating elements
+ Output format : 10
+ Reference frame : Ecliptic of J2000.0
+ *******************************************************************************
+ JDTDB
+ EC QR IN
+ OM W Tp
+ N MA TA
+ A AD PR
+ *******************************************************************************
+ 2451544.500000000 = A.D. 2000-Jan-01 00:00:00.0000 TDB
+ EC= 2.419714130280285E-04 QR= 1.567750331522140E-04 IN= 2.757017394957633E+01
+ OM= 8.366378701785975E+01 W = 1.902419328274697E+02 Tp= 2451544.906276242342
+ N = 2.851287804624782E+02 MA= 2.441589505052237E+02 TA= 2.441339985159917E+02
+ A = 1.568129774099390E-04 AD= 1.568509216676640E-04 PR= 1.262587380397310E+00
+ *******************************************************************************/
+ private static final double axisDeimosAU = 1.568129774099390E-04; // Semi-major axis [au]
+ private static final double eccentricityDeimos = 2.419714130280285E-04; // Eccentricity [-]
+ private static final double inclinationDeimos = 2.757017394957633E+01; // Inclination [degrees]
+ private static final double argPeriapsisDeimos = 1.902419328274697E+02; // Arg perifocus [degrees]
+ private static final double longNodeDeimos = 8.366378701785975E+01; // Long asc node [degrees]
+ private static final double periapsisPassageDeimos = 2451544.906276242342; // Time of periapsis [JD]
+ private static final double meanMotionDeimos = 2.851287804624782E+02; // Mean motion [degrees/day]
+ private static final double[] DEIMOSORBITPARS = new double[]
+ {axisDeimosAU, eccentricityDeimos, inclinationDeimos, argPeriapsisDeimos,
+ longNodeDeimos, periapsisPassageDeimos, meanMotionDeimos};
+
// Singleton instance
private static SolarSystemParameters instance = null;
-
+
// Mass in kg for solar system bodies
- private final Map