From 3b719b3176c98c3daedeef6cd0dc6559c87c41ef Mon Sep 17 00:00:00 2001 From: Rory Barnes Date: Fri, 9 Feb 2024 11:08:10 -0800 Subject: [PATCH] Added tests for DistOrb-RD4 and DistOrb-LL2. --- src/distorb.c | 5 +- tests/Distorb/LL2/earth.in | 19 + tests/Distorb/LL2/george.in | 16 + tests/Distorb/LL2/jupiter.in | 17 + tests/Distorb/LL2/mars.in | 17 + tests/Distorb/LL2/mercury.in | 22 + tests/Distorb/LL2/neptune.in | 16 + tests/Distorb/LL2/saturn.in | 15 + tests/Distorb/LL2/sun.in | 6 + tests/Distorb/LL2/test_LL2.py | 825 ++++++++++++++++++++++++++++++++++ tests/Distorb/LL2/venus.in | 17 + tests/Distorb/LL2/vpl.in | 19 + tests/Distorb/RD4/earth.in | 19 + tests/Distorb/RD4/george.in | 16 + tests/Distorb/RD4/jupiter.in | 17 + tests/Distorb/RD4/mars.in | 17 + tests/Distorb/RD4/mercury.in | 22 + tests/Distorb/RD4/neptune.in | 16 + tests/Distorb/RD4/saturn.in | 15 + tests/Distorb/RD4/sun.in | 6 + tests/Distorb/RD4/test_RD4.py | 825 ++++++++++++++++++++++++++++++++++ tests/Distorb/RD4/venus.in | 17 + tests/Distorb/RD4/vpl.in | 19 + 23 files changed, 1979 insertions(+), 4 deletions(-) create mode 100644 tests/Distorb/LL2/earth.in create mode 100644 tests/Distorb/LL2/george.in create mode 100644 tests/Distorb/LL2/jupiter.in create mode 100644 tests/Distorb/LL2/mars.in create mode 100644 tests/Distorb/LL2/mercury.in create mode 100644 tests/Distorb/LL2/neptune.in create mode 100644 tests/Distorb/LL2/saturn.in create mode 100644 tests/Distorb/LL2/sun.in create mode 100644 tests/Distorb/LL2/test_LL2.py create mode 100644 tests/Distorb/LL2/venus.in create mode 100644 tests/Distorb/LL2/vpl.in create mode 100644 tests/Distorb/RD4/earth.in create mode 100644 tests/Distorb/RD4/george.in create mode 100644 tests/Distorb/RD4/jupiter.in create mode 100644 tests/Distorb/RD4/mars.in create mode 100644 tests/Distorb/RD4/mercury.in create mode 100644 tests/Distorb/RD4/neptune.in create mode 100644 tests/Distorb/RD4/saturn.in create mode 100644 tests/Distorb/RD4/sun.in create mode 100644 tests/Distorb/RD4/test_RD4.py create mode 100644 tests/Distorb/RD4/venus.in create mode 100644 tests/Distorb/RD4/vpl.in diff --git a/src/distorb.c b/src/distorb.c index 6bb32f1ff..0a0875df8 100644 --- a/src/distorb.c +++ b/src/distorb.c @@ -2559,10 +2559,7 @@ void HessEigen(double **amat, int origsize, double real[], double imag[]) { q /= lrcorner; r /= lrcorner; } - } else { - // fprintf(stderr,"ERROR: k=m in distorb.c:HessEigen."); - // exit(EXIT_INT); - } + } value = sqrt(p * p + q * q + r * r); s = (double)fiSign(p) * value; if (s != 0.0) { diff --git a/tests/Distorb/LL2/earth.in b/tests/Distorb/LL2/earth.in new file mode 100644 index 000000000..8a35e4ed6 --- /dev/null +++ b/tests/Distorb/LL2/earth.in @@ -0,0 +1,19 @@ +sName Earth +saModules distorb distrot +dMass 3.0018090452e-06 +dRadius -1. +dRotPeriod -1 +dObliquity 23.44 #no Moons yet in vplanet.... +dSemi 1.00000321 +dEcc 0.01671327 +dInc 5e-5 +dLongP 102.948601 +dLongA 348.73936 +bGRCorr 1 +saOutputOrder Time Ecce Inc ArgP LongA Obliq PrecA -TotEn -TotAngMom +dDynEllip 3.28e-3 +dPrecA 0.0 +bForcePrecRate 1 +bCalcDynEllip 0 +#bOverrideMaxEcc 1 +#dHaltMaxEcc 0.4 diff --git a/tests/Distorb/LL2/george.in b/tests/Distorb/LL2/george.in new file mode 100644 index 000000000..6469da757 --- /dev/null +++ b/tests/Distorb/LL2/george.in @@ -0,0 +1,16 @@ +sName George +saModules distorb +dMass 4.3634170854e-05 +dRadius -1.50000 +dRotPeriod -1.00000 +dObliquity 0.00000 +dSemi 19.19203990 +dEcc 0.04713385 +dInc 0.76986 +dLongP 170.994339 +dLongA 74.22988 +#dDfcrit 0.1 +bGRCorr 1 +saOutputOrder Time Ecce Inc ArgP LongA +#bOverrideMaxEcc 1 +#dHaltMaxEcc 0.4 diff --git a/tests/Distorb/LL2/jupiter.in b/tests/Distorb/LL2/jupiter.in new file mode 100644 index 000000000..32d145dc8 --- /dev/null +++ b/tests/Distorb/LL2/jupiter.in @@ -0,0 +1,17 @@ +sName Jupiter +saModules distorb +dMass 9.5407035176e-04 +dRadius -1.50000 +dRotPeriod -1.00000 +dObliquity 0.00000 +dSemi 5.20880408 +dEcc 0.04932699 +dInc 1.30530 +dLongP 15.1623 +dLongA 100.55615 +#dDfcrit 0.1 +bGRCorr 1 +#sOrbitModel ll2 +saOutputOrder Time Ecce Inc ArgP LongA +#bOverrideMaxEcc 1 +#dHaltMaxEcc 0.4 diff --git a/tests/Distorb/LL2/mars.in b/tests/Distorb/LL2/mars.in new file mode 100644 index 000000000..dbb36dd5e --- /dev/null +++ b/tests/Distorb/LL2/mars.in @@ -0,0 +1,17 @@ +sName Mars +saModules distorb distrot +dMass 3.2253768844e-07 +dRadius -0.53202 +dRotPeriod -1.0259233 +dObliquity 25.189417 #obliquity, dynell from Armstrong 2004, Bouquillon & Souchay 1999 +dSemi 1.52366290 +dEcc 0.09341266 +dInc 1.85061 +dLongP 336.040919 +dLongA 49.57854 +bGRCorr 1 +saOutputOrder Time Ecce Inc ArgP LongA Obliq PrecA PrecFNat +dDynEllip 5.363e-3 +dPrecA 0.0 +#bOverrideMaxEcc 1 +#dHaltMaxEcc 0.4 diff --git a/tests/Distorb/LL2/mercury.in b/tests/Distorb/LL2/mercury.in new file mode 100644 index 000000000..2c52f03e1 --- /dev/null +++ b/tests/Distorb/LL2/mercury.in @@ -0,0 +1,22 @@ +sName Mercury +saModules distorb +dMass 1.6592964824e-07 +dRadius -1.00 +dRotPeriod -1.00000 +dObliquity 23.5 +dSemi 0.38709897 +dEcc 0.20563056 +dInc 7.00487 +dLongP 77.456453 +dLongA 48.33167 +bGRCorr 1 +#sOrbitModel rd4 +sOrbitModel ll2 +saOutputOrder Time Ecce Inc ArgP LongA +#bOverrideMaxEcc 1 +#dHaltMaxEcc 0.4 +#bOutputLapl 1 +bHaltHillStab 0 +bHaltCloseEnc 1 +bInvPlane 0 +bOutputEigen 0 diff --git a/tests/Distorb/LL2/neptune.in b/tests/Distorb/LL2/neptune.in new file mode 100644 index 000000000..fecc292a5 --- /dev/null +++ b/tests/Distorb/LL2/neptune.in @@ -0,0 +1,16 @@ +sName Neptune +saModules distorb +dMass 5.1472361809e-05 +dRadius -1.50000 +dRotPeriod -1.00000 +dObliquity 0.00000 +dSemi 30.07050641 +dEcc 0.00857659 +dInc 1.76917 +dLongP 44.633872 +dLongA 131.72169 +#dDfcrit 0.1 +bGRCorr 1 +saOutputOrder Time Ecce Inc ArgP LongA +#bOverrideMaxEcc 1 +#dHaltMaxEcc 0.4 diff --git a/tests/Distorb/LL2/saturn.in b/tests/Distorb/LL2/saturn.in new file mode 100644 index 000000000..9ca4cd5b3 --- /dev/null +++ b/tests/Distorb/LL2/saturn.in @@ -0,0 +1,15 @@ +sName Saturn +saModules distorb +dMass 2.8565829146e-04 +dRadius -1.00 #radius, rotper, and obl don't matter without distrot and/or eqtide +dRotPeriod -1.00000 +dObliquity 23.5 +dSemi 9.53999265 +dEcc 0.05434133 +dInc 2.48446 +dLongP 92.193712 +dLongA 113.71504 +bGRCorr 1 +saOutputOrder Time Ecce Inc ArgP LongA +#bOverrideMaxEcc 1 +#dHaltMaxEcc 0.4 diff --git a/tests/Distorb/LL2/sun.in b/tests/Distorb/LL2/sun.in new file mode 100644 index 000000000..62b1f88fc --- /dev/null +++ b/tests/Distorb/LL2/sun.in @@ -0,0 +1,6 @@ +# sun parameters +sName sun +dMass 1 +dSemi 0 +dEcc 0 +dRadius 0.00135 diff --git a/tests/Distorb/LL2/test_LL2.py b/tests/Distorb/LL2/test_LL2.py new file mode 100644 index 000000000..2c2f5ca7e --- /dev/null +++ b/tests/Distorb/LL2/test_LL2.py @@ -0,0 +1,825 @@ +from benchmark import Benchmark, benchmark +import astropy.units as u + +@benchmark( + { + "log.initial.system.Age": {"value": 0.000000, "unit": u.sec}, + "log.initial.system.Time": {"value": 0.000000, "unit": u.sec}, + "log.initial.system.TotAngMom": {"value": 3.277987e+43, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.system.TotEnergy": {"value": -7.839546e+41, "unit": u.Joule}, + "log.initial.system.PotEnergy": {"value": -7.839908e+41, "unit": u.Joule}, + "log.initial.system.KinEnergy": {"value": 5.361272e+37, "unit": u.Joule}, + "log.initial.system.DeltaTime": {"value": 0.000000, "unit": u.sec}, + "log.initial.sun.Mass": {"value": 1.988416e+30, "unit": u.kg}, + "log.initial.sun.Radius": {"value": 2.019571e+08, "unit": u.m}, + "log.initial.sun.RadGyra": {"value": 0.500000}, + "log.initial.sun.BodyType": {"value": 0.000000}, + "log.initial.sun.Density": {"value": 5.762900e+04, "unit": u.kg / u.m ** 3}, + "log.initial.sun.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.initial.sun.HZLimRecVenus": {"value": -1.000000}, + "log.initial.sun.HZLimRunaway": {"value": -1.000000}, + "log.initial.sun.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.initial.sun.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.initial.sun.HZLimEarlyMars": {"value": -1.000000}, + "log.initial.sun.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.initial.Mercury.Mass": {"value": 3.299372e+23, "unit": u.kg}, + "log.initial.Mercury.Radius": {"value": 6.378100e+06, "unit": u.m}, + "log.initial.Mercury.RadGyra": {"value": 0.500000}, + "log.initial.Mercury.RotKinEnergy": {"value": 8.872707e+27, "unit": u.Joule}, + "log.initial.Mercury.RotVel": {"value": 463.828521, "unit": u.m / u.sec}, + "log.initial.Mercury.BodyType": {"value": 0.000000}, + "log.initial.Mercury.Density": {"value": 303.576857, "unit": u.kg / u.m ** 3}, + "log.initial.Mercury.HEcc": {"value": 0.200722}, + "log.initial.Mercury.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.initial.Mercury.HZLimRecVenus": {"value": -1.000000}, + "log.initial.Mercury.HZLimRunaway": {"value": -1.000000}, + "log.initial.Mercury.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.initial.Mercury.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.initial.Mercury.HZLimEarlyMars": {"value": -1.000000}, + "log.initial.Mercury.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.initial.Mercury.KEcc": {"value": 0.044659}, + "log.initial.Mercury.Eccentricity": {"value": 0.205631}, + "log.initial.Mercury.OrbEnergy": {"value": -3.780641e+32, "unit": u.Joule}, + "log.initial.Mercury.MeanMotion": {"value": 8.266749e-07, "unit": 1 / u.sec}, + "log.initial.Mercury.OrbPeriod": {"value": 7.600552e+06, "unit": u.sec}, + "log.initial.Mercury.SemiMajorAxis": {"value": 5.790918e+10, "unit": u.m}, + "log.initial.Mercury.OrbAngMom": {"value": 8.951156e+38, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.Mercury.ArgP": {"value": 0.508323, "unit": u.rad}, + "log.initial.Mercury.Inc": {"value": 0.122258, "unit": u.rad}, + "log.initial.Mercury.LongA": {"value": 0.843547, "unit": u.rad}, + "log.initial.Mercury.MeanLongitude": {"value": 1.351870, "unit": u.rad}, + "log.initial.Mercury.LongP": {"value": 1.351870, "unit": u.rad}, + "log.initial.Mercury.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.initial.Mercury.OrbPotEnergy": {"value": -7.561283e+32, "unit": u.Joule}, + "log.initial.Mercury.Sinc": {"value": 0.061091}, + "log.initial.Mercury.Pinc": {"value": 0.045635}, + "log.initial.Mercury.Qinc": {"value": 0.040614}, + "log.initial.Mercury.DIncDtDistOrb": {"value": 0.122411, "unit": u.rad / u.sec}, + "log.initial.Mercury.DSincDtDistOrb": {"value": 0.061091, "unit": 1 / u.sec}, + "log.initial.Mercury.DLongADtDistOrb": {"value": -4.943962e-17, "unit": u.rad / u.sec}, + "log.initial.Mercury.DLongPDtDistOrb": {"value": -5.137340e-17, "unit": u.rad / u.sec}, + "log.initial.Mercury.DHeccDtDistOrb": {"value": 0.200722, "unit": 1 / u.sec}, + "log.initial.Mercury.DKeccDtDistOrb": {"value": 0.044659, "unit": 1 / u.sec}, + "log.initial.Mercury.DPincDtDistOrb": {"value": 0.045635, "unit": 1 / u.sec}, + "log.initial.Mercury.DQincDtDistOrb": {"value": 0.040614, "unit": 1 / u.sec}, + "log.initial.Venus.Mass": {"value": 4.864625e+24, "unit": u.kg}, + "log.initial.Venus.Obliquity": {"value": 3.095516, "unit": u.rad}, + "log.initial.Venus.PrecA": {"value": 0.000000, "unit": u.rad}, + "log.initial.Venus.Xobl": {"value": 0.046060}, + "log.initial.Venus.Yobl": {"value": 0.000000}, + "log.initial.Venus.Zobl": {"value": -0.998939}, + "log.initial.Venus.Radius": {"value": 6.058557e+06, "unit": u.m}, + "log.initial.Venus.RadGyra": {"value": 0.500000}, + "log.initial.Venus.RotAngMom": {"value": 1.332159e+31, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.Venus.RotKinEnergy": {"value": 1.987717e+24, "unit": u.Joule}, + "log.initial.Venus.RotVel": {"value": 1.807997, "unit": u.m / u.sec}, + "log.initial.Venus.BodyType": {"value": 0.000000}, + "log.initial.Venus.RotRate": {"value": 2.984203e-07, "unit": 1 / u.sec}, + "log.initial.Venus.RotPer": {"value": 2.105482e+07, "unit": u.sec}, + "log.initial.Venus.Density": {"value": 5222.195213, "unit": u.kg / u.m ** 3}, + "log.initial.Venus.HEcc": {"value": 0.005070}, + "log.initial.Venus.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.initial.Venus.HZLimRecVenus": {"value": -1.000000}, + "log.initial.Venus.HZLimRunaway": {"value": -1.000000}, + "log.initial.Venus.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.initial.Venus.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.initial.Venus.HZLimEarlyMars": {"value": -1.000000}, + "log.initial.Venus.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.initial.Venus.KEcc": {"value": -0.004493}, + "log.initial.Venus.Eccentricity": {"value": 0.006775}, + "log.initial.Venus.OrbEnergy": {"value": -2.983094e+33, "unit": u.Joule}, + "log.initial.Venus.MeanMotion": {"value": 3.236385e-07, "unit": 1 / u.sec}, + "log.initial.Venus.OrbPeriod": {"value": 1.941421e+07, "unit": u.sec}, + "log.initial.Venus.SemiMajorAxis": {"value": 1.082092e+11, "unit": u.m}, + "log.initial.Venus.COPP": {"value": 0.000234}, + "log.initial.Venus.OrbAngMom": {"value": 1.843430e+40, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.Venus.ArgP": {"value": 0.957634, "unit": u.rad}, + "log.initial.Venus.Inc": {"value": 0.059249, "unit": u.rad}, + "log.initial.Venus.LongA": {"value": 1.338331, "unit": u.rad}, + "log.initial.Venus.MeanLongitude": {"value": 2.295964, "unit": u.rad}, + "log.initial.Venus.LongP": {"value": 2.295964, "unit": u.rad}, + "log.initial.Venus.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.initial.Venus.OrbPotEnergy": {"value": -5.966187e+33, "unit": u.Joule}, + "log.initial.Venus.Sinc": {"value": 0.029620}, + "log.initial.Venus.Pinc": {"value": 0.028823}, + "log.initial.Venus.Qinc": {"value": 0.006824}, + "log.initial.Venus.DIncDtDistOrb": {"value": 0.059266, "unit": u.rad / u.sec}, + "log.initial.Venus.DSincDtDistOrb": {"value": 0.029620, "unit": 1 / u.sec}, + "log.initial.Venus.DLongADtDistOrb": {"value": 2.775558e-17, "unit": u.rad / u.sec}, + "log.initial.Venus.DLongPDtDistOrb": {"value": 2.428396e-15, "unit": u.rad / u.sec}, + "log.initial.Venus.DHeccDtDistOrb": {"value": 0.005070, "unit": 1 / u.sec}, + "log.initial.Venus.DKeccDtDistOrb": {"value": -0.004493, "unit": 1 / u.sec}, + "log.initial.Venus.DPincDtDistOrb": {"value": 0.028823, "unit": 1 / u.sec}, + "log.initial.Venus.DQincDtDistOrb": {"value": 0.006824, "unit": 1 / u.sec}, + "log.initial.Venus.CassiniOne": {"value": 0.819026}, + "log.initial.Venus.CassiniTwo": {"value": 0.573757}, + "log.initial.Venus.DOblDtDistRot": {"value": -8.754474e-14, "unit": u.rad}, + "log.initial.Venus.DPrecADtDistRot": {"value": -6.331626e-12, "unit": u.rad / u.sec}, + "log.initial.Venus.DXoblDtDistRot": {"value": 8.782448e-14, "unit": 1 / u.sec}, + "log.initial.Venus.DYoblDtDistRot": {"value": -2.916372e-13, "unit": 1 / u.sec}, + "log.initial.Venus.DZoblDtDistRot": {"value": -4.049528e-15, "unit": 1 / u.sec}, + "log.initial.Venus.OblTimeDistRot": {"value": 3.588557e+13, "unit": u.rad}, + "log.initial.Venus.PrecATimeDistRot": {"value": 9.923494e+11, "unit": u.sec}, + "log.initial.Venus.XoblTimeDistRot": {"value": 1.138635e+13, "unit": 1 / u.sec}, + "log.initial.Venus.YoblTimeDistRot": {"value": 3.428918e+12, "unit": 1 / u.sec}, + "log.initial.Venus.ZoblTimeDistRot": {"value": 2.469424e+14, "unit": 1 / u.sec}, + "log.initial.Venus.DynEllip": {"value": 1.312500e-05}, + "log.initial.Venus.PrecFNat": {"value": -6.903194e-12, "unit": u.rad / u.sec}, + "log.initial.Earth.Mass": {"value": 5.968845e+24, "unit": u.kg}, + "log.initial.Earth.Obliquity": {"value": 0.409105, "unit": u.rad}, + "log.initial.Earth.PrecA": {"value": 0.000000, "unit": u.rad}, + "log.initial.Earth.Xobl": {"value": 0.397789}, + "log.initial.Earth.Yobl": {"value": 0.000000}, + "log.initial.Earth.Zobl": {"value": 0.917477}, + "log.initial.Earth.Radius": {"value": 6.378100e+06, "unit": u.m}, + "log.initial.Earth.RadGyra": {"value": 0.500000}, + "log.initial.Earth.RotAngMom": {"value": 4.414475e+33, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.Earth.RotKinEnergy": {"value": 1.605149e+29, "unit": u.Joule}, + "log.initial.Earth.RotVel": {"value": 463.828521, "unit": u.m / u.sec}, + "log.initial.Earth.BodyType": {"value": 0.000000}, + "log.initial.Earth.RotRate": {"value": 7.272205e-05, "unit": 1 / u.sec}, + "log.initial.Earth.RotPer": {"value": 8.640000e+04, "unit": u.sec}, + "log.initial.Earth.Density": {"value": 5491.964602, "unit": u.kg / u.m ** 3}, + "log.initial.Earth.HEcc": {"value": 0.016288}, + "log.initial.Earth.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.initial.Earth.HZLimRecVenus": {"value": -1.000000}, + "log.initial.Earth.HZLimRunaway": {"value": -1.000000}, + "log.initial.Earth.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.initial.Earth.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.initial.Earth.HZLimEarlyMars": {"value": -1.000000}, + "log.initial.Earth.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.initial.Earth.KEcc": {"value": -0.003745}, + "log.initial.Earth.Eccentricity": {"value": 0.016713}, + "log.initial.Earth.OrbEnergy": {"value": -2.647556e+33, "unit": u.Joule}, + "log.initial.Earth.MeanMotion": {"value": 1.990977e-07, "unit": 1 / u.sec}, + "log.initial.Earth.OrbPeriod": {"value": 3.155830e+07, "unit": u.sec}, + "log.initial.Earth.SemiMajorAxis": {"value": 1.495984e+11, "unit": u.m}, + "log.initial.Earth.COPP": {"value": 0.006479}, + "log.initial.Earth.OrbAngMom": {"value": 2.659183e+40, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.Earth.ArgP": {"value": 1.993327, "unit": u.rad}, + "log.initial.Earth.Inc": {"value": 8.726646e-07, "unit": u.rad}, + "log.initial.Earth.LongA": {"value": 6.086650, "unit": u.rad}, + "log.initial.Earth.MeanLongitude": {"value": 1.796792, "unit": u.rad}, + "log.initial.Earth.LongP": {"value": 1.796792, "unit": u.rad}, + "log.initial.Earth.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.initial.Earth.OrbPotEnergy": {"value": -5.295112e+33, "unit": u.Joule}, + "log.initial.Earth.Sinc": {"value": 4.363323e-07}, + "log.initial.Earth.Pinc": {"value": -8.520368e-08}, + "log.initial.Earth.Qinc": {"value": 4.279325e-07}, + "log.initial.Earth.DIncDtDistOrb": {"value": 8.726646e-07, "unit": u.rad / u.sec}, + "log.initial.Earth.DSincDtDistOrb": {"value": 4.363323e-07, "unit": 1 / u.sec}, + "log.initial.Earth.DLongADtDistOrb": {"value": -2.387424e-12, "unit": u.rad / u.sec}, + "log.initial.Earth.DLongPDtDistOrb": {"value": 9.649399e-17, "unit": u.rad / u.sec}, + "log.initial.Earth.DHeccDtDistOrb": {"value": 0.016288, "unit": 1 / u.sec}, + "log.initial.Earth.DKeccDtDistOrb": {"value": -0.003745, "unit": 1 / u.sec}, + "log.initial.Earth.DPincDtDistOrb": {"value": -8.520368e-08, "unit": 1 / u.sec}, + "log.initial.Earth.DQincDtDistOrb": {"value": 4.279325e-07, "unit": 1 / u.sec}, + "log.initial.Earth.CassiniOne": {"value": 0.953009}, + "log.initial.Earth.CassiniTwo": {"value": 0.302943}, + "log.initial.Earth.DOblDtDistRot": {"value": -7.271557e-14, "unit": u.rad}, + "log.initial.Earth.DPrecADtDistRot": {"value": 2.446527e-12, "unit": u.rad / u.sec}, + "log.initial.Earth.DXoblDtDistRot": {"value": -6.671488e-14, "unit": 1 / u.sec}, + "log.initial.Earth.DYoblDtDistRot": {"value": 9.732005e-13, "unit": 1 / u.sec}, + "log.initial.Earth.DZoblDtDistRot": {"value": 2.892542e-14, "unit": 1 / u.sec}, + "log.initial.Earth.OblTimeDistRot": {"value": 4.320385e+13, "unit": u.rad}, + "log.initial.Earth.PrecATimeDistRot": {"value": 2.568205e+12, "unit": u.sec}, + "log.initial.Earth.XoblTimeDistRot": {"value": 1.498916e+13, "unit": 1 / u.sec}, + "log.initial.Earth.YoblTimeDistRot": {"value": 1.027537e+12, "unit": 1 / u.sec}, + "log.initial.Earth.ZoblTimeDistRot": {"value": 3.457167e+13, "unit": 1 / u.sec}, + "log.initial.Earth.DynEllip": {"value": 0.003280}, + "log.initial.Earth.PrecFNat": {"value": 2.461543e-12, "unit": u.rad / u.sec}, + "log.initial.Mars.Mass": {"value": 6.413391e+23, "unit": u.kg}, + "log.initial.Mars.Obliquity": {"value": 0.439638, "unit": u.rad}, + "log.initial.Mars.PrecA": {"value": 0.000000, "unit": u.rad}, + "log.initial.Mars.Xobl": {"value": 0.425612}, + "log.initial.Mars.Yobl": {"value": 0.000000}, + "log.initial.Mars.Zobl": {"value": 0.904906}, + "log.initial.Mars.Radius": {"value": 3.393277e+06, "unit": u.m}, + "log.initial.Mars.RadGyra": {"value": 0.500000}, + "log.initial.Mars.RotAngMom": {"value": 1.308632e+32, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.Mars.RotKinEnergy": {"value": 4.638086e+27, "unit": u.Joule}, + "log.initial.Mars.RotVel": {"value": 240.530700, "unit": u.m / u.sec}, + "log.initial.Mars.BodyType": {"value": 0.000000}, + "log.initial.Mars.RotRate": {"value": 7.088449e-05, "unit": 1 / u.sec}, + "log.initial.Mars.RotPer": {"value": 8.863977e+04, "unit": u.sec}, + "log.initial.Mars.Density": {"value": 3918.693171, "unit": u.kg / u.m ** 3}, + "log.initial.Mars.HEcc": {"value": -0.037933}, + "log.initial.Mars.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.initial.Mars.HZLimRecVenus": {"value": -1.000000}, + "log.initial.Mars.HZLimRunaway": {"value": -1.000000}, + "log.initial.Mars.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.initial.Mars.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.initial.Mars.HZLimEarlyMars": {"value": -1.000000}, + "log.initial.Mars.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.initial.Mars.KEcc": {"value": 0.085364}, + "log.initial.Mars.Eccentricity": {"value": 0.093413}, + "log.initial.Mars.OrbEnergy": {"value": -1.867046e+32, "unit": u.Joule}, + "log.initial.Mars.MeanMotion": {"value": 1.058606e-07, "unit": 1 / u.sec}, + "log.initial.Mars.OrbPeriod": {"value": 5.935338e+07, "unit": u.sec}, + "log.initial.Mars.SemiMajorAxis": {"value": 2.279367e+11, "unit": u.m}, + "log.initial.Mars.COPP": {"value": -0.016145}, + "log.initial.Mars.OrbAngMom": {"value": 3.511943e+39, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.Mars.ArgP": {"value": 4.999712, "unit": u.rad}, + "log.initial.Mars.Inc": {"value": 0.032299, "unit": u.rad}, + "log.initial.Mars.LongA": {"value": 0.865309, "unit": u.rad}, + "log.initial.Mars.MeanLongitude": {"value": -0.418165, "unit": u.rad}, + "log.initial.Mars.LongP": {"value": 5.865020, "unit": u.rad}, + "log.initial.Mars.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.initial.Mars.OrbPotEnergy": {"value": -3.734093e+32, "unit": u.Joule}, + "log.initial.Mars.Sinc": {"value": 0.016149}, + "log.initial.Mars.Pinc": {"value": 0.012294}, + "log.initial.Mars.Qinc": {"value": 0.010471}, + "log.initial.Mars.DIncDtDistOrb": {"value": 0.032302, "unit": u.rad / u.sec}, + "log.initial.Mars.DSincDtDistOrb": {"value": 0.016149, "unit": 1 / u.sec}, + "log.initial.Mars.DLongADtDistOrb": {"value": 1.040834e-17, "unit": u.rad / u.sec}, + "log.initial.Mars.DLongPDtDistOrb": {"value": 2.705084e-17, "unit": u.rad / u.sec}, + "log.initial.Mars.DHeccDtDistOrb": {"value": -0.037933, "unit": 1 / u.sec}, + "log.initial.Mars.DKeccDtDistOrb": {"value": 0.085364, "unit": 1 / u.sec}, + "log.initial.Mars.DPincDtDistOrb": {"value": 0.012294, "unit": 1 / u.sec}, + "log.initial.Mars.DQincDtDistOrb": {"value": 0.010471, "unit": 1 / u.sec}, + "log.initial.Mars.CassiniOne": {"value": 0.056899}, + "log.initial.Mars.CassiniTwo": {"value": 0.998380}, + "log.initial.Mars.DOblDtDistRot": {"value": 8.402460e-15, "unit": u.rad}, + "log.initial.Mars.DPrecADtDistRot": {"value": 1.308952e-12, "unit": u.rad / u.sec}, + "log.initial.Mars.DXoblDtDistRot": {"value": 7.603434e-15, "unit": 1 / u.sec}, + "log.initial.Mars.DYoblDtDistRot": {"value": 5.571058e-13, "unit": 1 / u.sec}, + "log.initial.Mars.DZoblDtDistRot": {"value": -3.576189e-15, "unit": 1 / u.sec}, + "log.initial.Mars.OblTimeDistRot": {"value": 3.738896e+14, "unit": u.rad}, + "log.initial.Mars.PrecATimeDistRot": {"value": 4.800166e+12, "unit": u.sec}, + "log.initial.Mars.XoblTimeDistRot": {"value": 1.315195e+14, "unit": 1 / u.sec}, + "log.initial.Mars.YoblTimeDistRot": {"value": 1.794991e+12, "unit": 1 / u.sec}, + "log.initial.Mars.ZoblTimeDistRot": {"value": 2.796273e+14, "unit": 1 / u.sec}, + "log.initial.Mars.DynEllip": {"value": 0.005363}, + "log.initial.Mars.PrecFNat": {"value": 1.166082e-12, "unit": u.rad / u.sec}, + "log.initial.Jupiter.Mass": {"value": 1.897089e+27, "unit": u.kg}, + "log.initial.Jupiter.Radius": {"value": 9.567150e+06, "unit": u.m}, + "log.initial.Jupiter.RadGyra": {"value": 0.500000}, + "log.initial.Jupiter.RotKinEnergy": {"value": 1.147876e+32, "unit": u.Joule}, + "log.initial.Jupiter.RotVel": {"value": 695.742781, "unit": u.m / u.sec}, + "log.initial.Jupiter.BodyType": {"value": 0.000000}, + "log.initial.Jupiter.Density": {"value": 5.171914e+05, "unit": u.kg / u.m ** 3}, + "log.initial.Jupiter.HEcc": {"value": 0.012902}, + "log.initial.Jupiter.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.initial.Jupiter.HZLimRecVenus": {"value": -1.000000}, + "log.initial.Jupiter.HZLimRunaway": {"value": -1.000000}, + "log.initial.Jupiter.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.initial.Jupiter.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.initial.Jupiter.HZLimEarlyMars": {"value": -1.000000}, + "log.initial.Jupiter.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.initial.Jupiter.KEcc": {"value": 0.047610}, + "log.initial.Jupiter.Eccentricity": {"value": 0.049327}, + "log.initial.Jupiter.OrbEnergy": {"value": -1.615496e+35, "unit": u.Joule}, + "log.initial.Jupiter.MeanMotion": {"value": 1.675590e-08, "unit": 1 / u.sec}, + "log.initial.Jupiter.OrbPeriod": {"value": 3.749835e+08, "unit": u.sec}, + "log.initial.Jupiter.SemiMajorAxis": {"value": 7.792260e+11, "unit": u.m}, + "log.initial.Jupiter.OrbAngMom": {"value": 1.925924e+43, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.Jupiter.ArgP": {"value": 4.792781, "unit": u.rad}, + "log.initial.Jupiter.Inc": {"value": 0.022782, "unit": u.rad}, + "log.initial.Jupiter.LongA": {"value": 1.755036, "unit": u.rad}, + "log.initial.Jupiter.MeanLongitude": {"value": 0.264632, "unit": u.rad}, + "log.initial.Jupiter.LongP": {"value": 0.264632, "unit": u.rad}, + "log.initial.Jupiter.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.initial.Jupiter.OrbPotEnergy": {"value": -3.230992e+35, "unit": u.Joule}, + "log.initial.Jupiter.Sinc": {"value": 0.011391}, + "log.initial.Jupiter.Pinc": {"value": 0.011198}, + "log.initial.Jupiter.Qinc": {"value": -0.002087}, + "log.initial.Jupiter.DIncDtDistOrb": {"value": 0.022783, "unit": u.rad / u.sec}, + "log.initial.Jupiter.DSincDtDistOrb": {"value": 0.011391, "unit": 1 / u.sec}, + "log.initial.Jupiter.DLongADtDistOrb": {"value": -1.734723e-16, "unit": u.rad / u.sec}, + "log.initial.Jupiter.DLongPDtDistOrb": {"value": 6.613633e-18, "unit": u.rad / u.sec}, + "log.initial.Jupiter.DHeccDtDistOrb": {"value": 0.012902, "unit": 1 / u.sec}, + "log.initial.Jupiter.DKeccDtDistOrb": {"value": 0.047610, "unit": 1 / u.sec}, + "log.initial.Jupiter.DPincDtDistOrb": {"value": 0.011198, "unit": 1 / u.sec}, + "log.initial.Jupiter.DQincDtDistOrb": {"value": -0.002087, "unit": 1 / u.sec}, + "log.initial.Saturn.Mass": {"value": 5.680075e+26, "unit": u.kg}, + "log.initial.Saturn.Radius": {"value": 6.378100e+06, "unit": u.m}, + "log.initial.Saturn.RadGyra": {"value": 0.500000}, + "log.initial.Saturn.RotKinEnergy": {"value": 1.527492e+31, "unit": u.Joule}, + "log.initial.Saturn.RotVel": {"value": 463.828521, "unit": u.m / u.sec}, + "log.initial.Saturn.BodyType": {"value": 0.000000}, + "log.initial.Saturn.Density": {"value": 5.226266e+05, "unit": u.kg / u.m ** 3}, + "log.initial.Saturn.HEcc": {"value": 0.054302}, + "log.initial.Saturn.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.initial.Saturn.HZLimRecVenus": {"value": -1.000000}, + "log.initial.Saturn.HZLimRunaway": {"value": -1.000000}, + "log.initial.Saturn.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.initial.Saturn.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.initial.Saturn.HZLimEarlyMars": {"value": -1.000000}, + "log.initial.Saturn.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.initial.Saturn.KEcc": {"value": -0.002080}, + "log.initial.Saturn.Eccentricity": {"value": 0.054341}, + "log.initial.Saturn.OrbEnergy": {"value": -2.640963e+34, "unit": u.Joule}, + "log.initial.Saturn.MeanMotion": {"value": 6.757836e-09, "unit": 1 / u.sec}, + "log.initial.Saturn.OrbPeriod": {"value": 9.297629e+08, "unit": u.sec}, + "log.initial.Saturn.SemiMajorAxis": {"value": 1.427163e+12, "unit": u.m}, + "log.initial.Saturn.OrbAngMom": {"value": 7.804453e+42, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.Saturn.ArgP": {"value": 5.907567, "unit": u.rad}, + "log.initial.Saturn.Inc": {"value": 0.043362, "unit": u.rad}, + "log.initial.Saturn.LongA": {"value": 1.984702, "unit": u.rad}, + "log.initial.Saturn.MeanLongitude": {"value": 1.609084, "unit": u.rad}, + "log.initial.Saturn.LongP": {"value": 1.609084, "unit": u.rad}, + "log.initial.Saturn.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.initial.Saturn.OrbPotEnergy": {"value": -5.281926e+34, "unit": u.Joule}, + "log.initial.Saturn.Sinc": {"value": 0.021679}, + "log.initial.Saturn.Pinc": {"value": 0.019849}, + "log.initial.Saturn.Qinc": {"value": -0.008719}, + "log.initial.Saturn.DIncDtDistOrb": {"value": 0.043369, "unit": u.rad / u.sec}, + "log.initial.Saturn.DSincDtDistOrb": {"value": 0.021679, "unit": 1 / u.sec}, + "log.initial.Saturn.DLongADtDistOrb": {"value": 2.775558e-17, "unit": u.rad / u.sec}, + "log.initial.Saturn.DLongPDtDistOrb": {"value": -1.194791e-16, "unit": u.rad / u.sec}, + "log.initial.Saturn.DHeccDtDistOrb": {"value": 0.054302, "unit": 1 / u.sec}, + "log.initial.Saturn.DKeccDtDistOrb": {"value": -0.002080, "unit": 1 / u.sec}, + "log.initial.Saturn.DPincDtDistOrb": {"value": 0.019849, "unit": 1 / u.sec}, + "log.initial.Saturn.DQincDtDistOrb": {"value": -0.008719, "unit": 1 / u.sec}, + "log.initial.George.Mass": {"value": 8.676288e+25, "unit": u.kg}, + "log.initial.George.Radius": {"value": 9.567150e+06, "unit": u.m}, + "log.initial.George.RadGyra": {"value": 0.500000}, + "log.initial.George.RotKinEnergy": {"value": 5.249784e+30, "unit": u.Joule}, + "log.initial.George.RotVel": {"value": 695.742781, "unit": u.m / u.sec}, + "log.initial.George.BodyType": {"value": 0.000000}, + "log.initial.George.Density": {"value": 2.365362e+04, "unit": u.kg / u.m ** 3}, + "log.initial.George.HEcc": {"value": 0.007378}, + "log.initial.George.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.initial.George.HZLimRecVenus": {"value": -1.000000}, + "log.initial.George.HZLimRunaway": {"value": -1.000000}, + "log.initial.George.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.initial.George.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.initial.George.HZLimEarlyMars": {"value": -1.000000}, + "log.initial.George.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.initial.George.KEcc": {"value": -0.046553}, + "log.initial.George.Eccentricity": {"value": 0.047134}, + "log.initial.George.OrbEnergy": {"value": -2.005253e+33, "unit": u.Joule}, + "log.initial.George.MeanMotion": {"value": 2.368075e-09, "unit": 1 / u.sec}, + "log.initial.George.OrbPeriod": {"value": 2.653288e+09, "unit": u.sec}, + "log.initial.George.SemiMajorAxis": {"value": 2.871088e+12, "unit": u.m}, + "log.initial.George.OrbAngMom": {"value": 1.691690e+42, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.George.ArgP": {"value": 1.688858, "unit": u.rad}, + "log.initial.George.Inc": {"value": 0.013437, "unit": u.rad}, + "log.initial.George.LongA": {"value": 1.295556, "unit": u.rad}, + "log.initial.George.MeanLongitude": {"value": 2.984414, "unit": u.rad}, + "log.initial.George.LongP": {"value": 2.984414, "unit": u.rad}, + "log.initial.George.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.initial.George.OrbPotEnergy": {"value": -4.010505e+33, "unit": u.Joule}, + "log.initial.George.Sinc": {"value": 0.006718}, + "log.initial.George.Pinc": {"value": 0.006465}, + "log.initial.George.Qinc": {"value": 0.001826}, + "log.initial.George.DIncDtDistOrb": {"value": 0.013437, "unit": u.rad / u.sec}, + "log.initial.George.DSincDtDistOrb": {"value": 0.006718, "unit": 1 / u.sec}, + "log.initial.George.DLongADtDistOrb": {"value": -5.551115e-17, "unit": u.rad / u.sec}, + "log.initial.George.DLongPDtDistOrb": {"value": -2.081668e-17, "unit": u.rad / u.sec}, + "log.initial.George.DHeccDtDistOrb": {"value": 0.007378, "unit": 1 / u.sec}, + "log.initial.George.DKeccDtDistOrb": {"value": -0.046553, "unit": 1 / u.sec}, + "log.initial.George.DPincDtDistOrb": {"value": 0.006465, "unit": 1 / u.sec}, + "log.initial.George.DQincDtDistOrb": {"value": 0.001826, "unit": 1 / u.sec}, + "log.initial.Neptune.Mass": {"value": 1.023485e+26, "unit": u.kg}, + "log.initial.Neptune.Radius": {"value": 9.567150e+06, "unit": u.m}, + "log.initial.Neptune.RadGyra": {"value": 0.500000}, + "log.initial.Neptune.RotKinEnergy": {"value": 6.192825e+30, "unit": u.Joule}, + "log.initial.Neptune.RotVel": {"value": 695.742781, "unit": u.m / u.sec}, + "log.initial.Neptune.BodyType": {"value": 0.000000}, + "log.initial.Neptune.Density": {"value": 2.790262e+04, "unit": u.kg / u.m ** 3}, + "log.initial.Neptune.HEcc": {"value": 0.006026}, + "log.initial.Neptune.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.initial.Neptune.HZLimRecVenus": {"value": -1.000000}, + "log.initial.Neptune.HZLimRunaway": {"value": -1.000000}, + "log.initial.Neptune.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.initial.Neptune.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.initial.Neptune.HZLimEarlyMars": {"value": -1.000000}, + "log.initial.Neptune.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.initial.Neptune.KEcc": {"value": 0.006103}, + "log.initial.Neptune.Eccentricity": {"value": 0.008577}, + "log.initial.Neptune.OrbEnergy": {"value": -1.509722e+33, "unit": u.Joule}, + "log.initial.Neptune.MeanMotion": {"value": 1.207446e-09, "unit": 1 / u.sec}, + "log.initial.Neptune.OrbPeriod": {"value": 5.203698e+09, "unit": u.sec}, + "log.initial.Neptune.SemiMajorAxis": {"value": 4.498484e+12, "unit": u.m}, + "log.initial.Neptune.OrbAngMom": {"value": 2.500594e+42, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.Neptune.ArgP": {"value": 4.763216, "unit": u.rad}, + "log.initial.Neptune.Inc": {"value": 0.030878, "unit": u.rad}, + "log.initial.Neptune.LongA": {"value": 2.298977, "unit": u.rad}, + "log.initial.Neptune.MeanLongitude": {"value": 0.779008, "unit": u.rad}, + "log.initial.Neptune.LongP": {"value": 0.779008, "unit": u.rad}, + "log.initial.Neptune.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.initial.Neptune.OrbPotEnergy": {"value": -3.019443e+33, "unit": u.Joule}, + "log.initial.Neptune.Sinc": {"value": 0.015438}, + "log.initial.Neptune.Pinc": {"value": 0.011523}, + "log.initial.Neptune.Qinc": {"value": -0.010274}, + "log.initial.Neptune.DIncDtDistOrb": {"value": 0.030880, "unit": u.rad / u.sec}, + "log.initial.Neptune.DSincDtDistOrb": {"value": 0.015438, "unit": 1 / u.sec}, + "log.initial.Neptune.DLongADtDistOrb": {"value": 5.551115e-17, "unit": u.rad / u.sec}, + "log.initial.Neptune.DLongPDtDistOrb": {"value": -5.551115e-17, "unit": u.rad / u.sec}, + "log.initial.Neptune.DHeccDtDistOrb": {"value": 0.006026, "unit": 1 / u.sec}, + "log.initial.Neptune.DKeccDtDistOrb": {"value": 0.006103, "unit": 1 / u.sec}, + "log.initial.Neptune.DPincDtDistOrb": {"value": 0.011523, "unit": 1 / u.sec}, + "log.initial.Neptune.DQincDtDistOrb": {"value": -0.010274, "unit": 1 / u.sec}, + "log.final.system.Age": {"value": 1.577880e+11, "unit": u.sec}, + "log.final.system.Time": {"value": 1.577880e+11, "unit": u.sec}, + "log.final.system.TotAngMom": {"value": 3.277987e+43, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.final.system.TotEnergy": {"value": -7.839546e+41, "unit": u.Joule}, + "log.final.system.PotEnergy": {"value": -7.839908e+41, "unit": u.Joule}, + "log.final.system.KinEnergy": {"value": 5.361272e+37, "unit": u.Joule}, + "log.final.system.DeltaTime": {"value": 1.168800e+09, "unit": u.sec}, + "log.final.sun.Mass": {"value": 1.988416e+30, "unit": u.kg}, + "log.final.sun.Radius": {"value": 2.019571e+08, "unit": u.m}, + "log.final.sun.RadGyra": {"value": 0.500000}, + "log.final.sun.BodyType": {"value": 0.000000}, + "log.final.sun.Density": {"value": 5.762900e+04, "unit": u.kg / u.m ** 3}, + "log.final.sun.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.final.sun.HZLimRecVenus": {"value": -1.000000}, + "log.final.sun.HZLimRunaway": {"value": -1.000000}, + "log.final.sun.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.final.sun.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.final.sun.HZLimEarlyMars": {"value": -1.000000}, + "log.final.sun.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.final.Mercury.Mass": {"value": 3.299372e+23, "unit": u.kg}, + "log.final.Mercury.Radius": {"value": 6.378100e+06, "unit": u.m}, + "log.final.Mercury.RadGyra": {"value": 0.500000}, + "log.final.Mercury.RotKinEnergy": {"value": 8.872707e+27, "unit": u.Joule}, + "log.final.Mercury.RotVel": {"value": 463.828521, "unit": u.m / u.sec}, + "log.final.Mercury.BodyType": {"value": 0.000000}, + "log.final.Mercury.Density": {"value": 303.576857, "unit": u.kg / u.m ** 3}, + "log.final.Mercury.HEcc": {"value": 0.205108}, + "log.final.Mercury.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.final.Mercury.HZLimRecVenus": {"value": -1.000000}, + "log.final.Mercury.HZLimRunaway": {"value": -1.000000}, + "log.final.Mercury.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.final.Mercury.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.final.Mercury.HZLimEarlyMars": {"value": -1.000000}, + "log.final.Mercury.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.final.Mercury.KEcc": {"value": 0.015919}, + "log.final.Mercury.Eccentricity": {"value": 0.205724}, + "log.final.Mercury.OrbEnergy": {"value": -3.780641e+32, "unit": u.Joule}, + "log.final.Mercury.MeanMotion": {"value": 8.266749e-07, "unit": 1 / u.sec}, + "log.final.Mercury.OrbPeriod": {"value": 7.600552e+06, "unit": u.sec}, + "log.final.Mercury.SemiMajorAxis": {"value": 5.790918e+10, "unit": u.m}, + "log.final.Mercury.OrbAngMom": {"value": 8.950976e+38, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.final.Mercury.ArgP": {"value": 0.748997, "unit": u.rad}, + "log.final.Mercury.Inc": {"value": 0.119685, "unit": u.rad}, + "log.final.Mercury.LongA": {"value": 0.744341, "unit": u.rad}, + "log.final.Mercury.MeanLongitude": {"value": 1.493338, "unit": u.rad}, + "log.final.Mercury.LongP": {"value": 1.493338, "unit": u.rad}, + "log.final.Mercury.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.final.Mercury.OrbPotEnergy": {"value": -7.561283e+32, "unit": u.Joule}, + "log.final.Mercury.Sinc": {"value": 0.059807}, + "log.final.Mercury.Pinc": {"value": 0.040518}, + "log.final.Mercury.Qinc": {"value": 0.043990}, + "log.final.Mercury.DIncDtDistOrb": {"value": 0.119816, "unit": u.rad / u.sec}, + "log.final.Mercury.DSincDtDistOrb": {"value": 0.059801, "unit": 1 / u.sec}, + "log.final.Mercury.DLongADtDistOrb": {"value": -0.000466, "unit": u.rad / u.sec}, + "log.final.Mercury.DLongPDtDistOrb": {"value": 0.000660, "unit": u.rad / u.sec}, + "log.final.Mercury.DHeccDtDistOrb": {"value": 0.205118, "unit": 1 / u.sec}, + "log.final.Mercury.DKeccDtDistOrb": {"value": 0.015784, "unit": 1 / u.sec}, + "log.final.Mercury.DPincDtDistOrb": {"value": 0.040494, "unit": 1 / u.sec}, + "log.final.Mercury.DQincDtDistOrb": {"value": 0.044005, "unit": 1 / u.sec}, + "log.final.Venus.Mass": {"value": 4.864625e+24, "unit": u.kg}, + "log.final.Venus.Obliquity": {"value": 3.087670, "unit": u.rad}, + "log.final.Venus.PrecA": {"value": 5.411501, "unit": u.rad}, + "log.final.Venus.Xobl": {"value": 0.034712}, + "log.final.Venus.Yobl": {"value": -0.041285}, + "log.final.Venus.Zobl": {"value": -0.999333}, + "log.final.Venus.Radius": {"value": 6.058557e+06, "unit": u.m}, + "log.final.Venus.RadGyra": {"value": 0.500000}, + "log.final.Venus.RotAngMom": {"value": 1.332159e+31, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.final.Venus.RotKinEnergy": {"value": 1.987717e+24, "unit": u.Joule}, + "log.final.Venus.RotVel": {"value": 1.807997, "unit": u.m / u.sec}, + "log.final.Venus.BodyType": {"value": 0.000000}, + "log.final.Venus.RotRate": {"value": 2.984203e-07, "unit": 1 / u.sec}, + "log.final.Venus.RotPer": {"value": 2.105482e+07, "unit": u.sec}, + "log.final.Venus.Density": {"value": 5222.195213, "unit": u.kg / u.m ** 3}, + "log.final.Venus.HEcc": {"value": 0.003782}, + "log.final.Venus.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.final.Venus.HZLimRecVenus": {"value": -1.000000}, + "log.final.Venus.HZLimRunaway": {"value": -1.000000}, + "log.final.Venus.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.final.Venus.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.final.Venus.HZLimEarlyMars": {"value": -1.000000}, + "log.final.Venus.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.final.Venus.KEcc": {"value": -0.002863}, + "log.final.Venus.Eccentricity": {"value": 0.004744}, + "log.final.Venus.OrbEnergy": {"value": -2.983094e+33, "unit": u.Joule}, + "log.final.Venus.MeanMotion": {"value": 3.236385e-07, "unit": 1 / u.sec}, + "log.final.Venus.OrbPeriod": {"value": 1.941421e+07, "unit": u.sec}, + "log.final.Venus.SemiMajorAxis": {"value": 1.082092e+11, "unit": u.m}, + "log.final.Venus.COPP": {"value": 0.000204}, + "log.final.Venus.OrbAngMom": {"value": 1.843452e+40, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.final.Venus.ArgP": {"value": 1.130238, "unit": u.rad}, + "log.final.Venus.Inc": {"value": 0.057031, "unit": u.rad}, + "log.final.Venus.LongA": {"value": 1.088424, "unit": u.rad}, + "log.final.Venus.MeanLongitude": {"value": 2.218661, "unit": u.rad}, + "log.final.Venus.LongP": {"value": 2.218661, "unit": u.rad}, + "log.final.Venus.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.final.Venus.OrbPotEnergy": {"value": -5.966187e+33, "unit": u.Joule}, + "log.final.Venus.Sinc": {"value": 0.028512}, + "log.final.Venus.Pinc": {"value": 0.025258}, + "log.final.Venus.Qinc": {"value": 0.013226}, + "log.final.Venus.DIncDtDistOrb": {"value": 0.057030, "unit": u.rad / u.sec}, + "log.final.Venus.DSincDtDistOrb": {"value": 0.028503, "unit": 1 / u.sec}, + "log.final.Venus.DLongADtDistOrb": {"value": -0.001200, "unit": u.rad / u.sec}, + "log.final.Venus.DLongPDtDistOrb": {"value": -0.000785, "unit": u.rad / u.sec}, + "log.final.Venus.DHeccDtDistOrb": {"value": 0.003779, "unit": 1 / u.sec}, + "log.final.Venus.DKeccDtDistOrb": {"value": -0.002855, "unit": 1 / u.sec}, + "log.final.Venus.DPincDtDistOrb": {"value": 0.025235, "unit": 1 / u.sec}, + "log.final.Venus.DQincDtDistOrb": {"value": 0.013252, "unit": 1 / u.sec}, + "log.final.Venus.CassiniOne": {"value": 0.259910}, + "log.final.Venus.CassiniTwo": {"value": 0.965633}, + "log.final.Venus.DOblDtDistRot": {"value": 2.269808e-15, "unit": u.rad}, + "log.final.Venus.DPrecADtDistRot": {"value": -5.122749e-12, "unit": u.rad / u.sec}, + "log.final.Venus.DXoblDtDistRot": {"value": -2.129629e-13, "unit": 1 / u.sec}, + "log.final.Venus.DYoblDtDistRot": {"value": -1.760723e-13, "unit": 1 / u.sec}, + "log.final.Venus.DZoblDtDistRot": {"value": 1.232433e-16, "unit": 1 / u.sec}, + "log.final.Venus.OblTimeDistRot": {"value": 1.384078e+15, "unit": u.rad}, + "log.final.Venus.PrecATimeDistRot": {"value": 1.226526e+12, "unit": u.sec}, + "log.final.Venus.XoblTimeDistRot": {"value": 4.695654e+12, "unit": 1 / u.sec}, + "log.final.Venus.YoblTimeDistRot": {"value": 5.679485e+12, "unit": 1 / u.sec}, + "log.final.Venus.ZoblTimeDistRot": {"value": 8.114030e+15, "unit": 1 / u.sec}, + "log.final.Venus.DynEllip": {"value": 1.312500e-05}, + "log.final.Venus.PrecFNat": {"value": -6.905677e-12, "unit": u.rad / u.sec}, + "log.final.Earth.Mass": {"value": 5.968845e+24, "unit": u.kg}, + "log.final.Earth.Obliquity": {"value": 0.398052, "unit": u.rad}, + "log.final.Earth.PrecA": {"value": 0.386029, "unit": u.rad}, + "log.final.Earth.Xobl": {"value": 0.359099}, + "log.final.Earth.Yobl": {"value": 0.145945}, + "log.final.Earth.Zobl": {"value": 0.921818}, + "log.final.Earth.Radius": {"value": 6.378100e+06, "unit": u.m}, + "log.final.Earth.RadGyra": {"value": 0.500000}, + "log.final.Earth.RotAngMom": {"value": 4.414475e+33, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.final.Earth.RotKinEnergy": {"value": 1.605149e+29, "unit": u.Joule}, + "log.final.Earth.RotVel": {"value": 463.828521, "unit": u.m / u.sec}, + "log.final.Earth.BodyType": {"value": 0.000000}, + "log.final.Earth.RotRate": {"value": 7.272205e-05, "unit": 1 / u.sec}, + "log.final.Earth.RotPer": {"value": 8.640000e+04, "unit": u.sec}, + "log.final.Earth.Density": {"value": 5491.964602, "unit": u.kg / u.m ** 3}, + "log.final.Earth.HEcc": {"value": 0.012529}, + "log.final.Earth.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.final.Earth.HZLimRecVenus": {"value": -1.000000}, + "log.final.Earth.HZLimRunaway": {"value": -1.000000}, + "log.final.Earth.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.final.Earth.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.final.Earth.HZLimEarlyMars": {"value": -1.000000}, + "log.final.Earth.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.final.Earth.KEcc": {"value": -0.007039}, + "log.final.Earth.Eccentricity": {"value": 0.014370}, + "log.final.Earth.OrbEnergy": {"value": -2.647556e+33, "unit": u.Joule}, + "log.final.Earth.MeanMotion": {"value": 1.990977e-07, "unit": 1 / u.sec}, + "log.final.Earth.OrbPeriod": {"value": 3.155830e+07, "unit": u.sec}, + "log.final.Earth.SemiMajorAxis": {"value": 1.495984e+11, "unit": u.m}, + "log.final.Earth.COPP": {"value": 0.004856}, + "log.final.Earth.OrbAngMom": {"value": 2.659280e+40, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.final.Earth.ArgP": {"value": 5.524282, "unit": u.rad}, + "log.final.Earth.Inc": {"value": 0.010982, "unit": u.rad}, + "log.final.Earth.LongA": {"value": 2.841562, "unit": u.rad}, + "log.final.Earth.MeanLongitude": {"value": 2.082658, "unit": u.rad}, + "log.final.Earth.LongP": {"value": 2.082658, "unit": u.rad}, + "log.final.Earth.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.final.Earth.OrbPotEnergy": {"value": -5.295112e+33, "unit": u.Joule}, + "log.final.Earth.Sinc": {"value": 0.005491}, + "log.final.Earth.Pinc": {"value": 0.001623}, + "log.final.Earth.Qinc": {"value": -0.005245}, + "log.final.Earth.DIncDtDistOrb": {"value": 0.011030, "unit": u.rad / u.sec}, + "log.final.Earth.DSincDtDistOrb": {"value": 0.005515, "unit": 1 / u.sec}, + "log.final.Earth.DLongADtDistOrb": {"value": -0.000982, "unit": u.rad / u.sec}, + "log.final.Earth.DLongPDtDistOrb": {"value": 0.001383, "unit": u.rad / u.sec}, + "log.final.Earth.DHeccDtDistOrb": {"value": 0.012508, "unit": 1 / u.sec}, + "log.final.Earth.DKeccDtDistOrb": {"value": -0.007050, "unit": 1 / u.sec}, + "log.final.Earth.DPincDtDistOrb": {"value": 0.001635, "unit": 1 / u.sec}, + "log.final.Earth.DQincDtDistOrb": {"value": -0.005267, "unit": 1 / u.sec}, + "log.final.Earth.CassiniOne": {"value": 0.957336}, + "log.final.Earth.CassiniTwo": {"value": 0.288978}, + "log.final.Earth.DOblDtDistRot": {"value": -6.736846e-14, "unit": u.rad}, + "log.final.Earth.DPrecADtDistRot": {"value": 2.451482e-12, "unit": u.rad / u.sec}, + "log.final.Earth.DXoblDtDistRot": {"value": -4.153135e-13, "unit": 1 / u.sec}, + "log.final.Earth.DYoblDtDistRot": {"value": 8.569418e-13, "unit": 1 / u.sec}, + "log.final.Earth.DZoblDtDistRot": {"value": 2.611358e-14, "unit": 1 / u.sec}, + "log.final.Earth.OblTimeDistRot": {"value": 4.663299e+13, "unit": u.rad}, + "log.final.Earth.PrecATimeDistRot": {"value": 2.563015e+12, "unit": u.sec}, + "log.final.Earth.XoblTimeDistRot": {"value": 2.407820e+12, "unit": 1 / u.sec}, + "log.final.Earth.YoblTimeDistRot": {"value": 1.166940e+12, "unit": 1 / u.sec}, + "log.final.Earth.ZoblTimeDistRot": {"value": 3.829425e+13, "unit": 1 / u.sec}, + "log.final.Earth.DynEllip": {"value": 0.003280}, + "log.final.Earth.PrecFNat": {"value": 2.472919e-12, "unit": u.rad / u.sec}, + "log.final.Mars.Mass": {"value": 6.413391e+23, "unit": u.kg}, + "log.final.Mars.Obliquity": {"value": 0.440047, "unit": u.rad}, + "log.final.Mars.PrecA": {"value": 0.208098, "unit": u.rad}, + "log.final.Mars.Xobl": {"value": 0.416791}, + "log.final.Mars.Yobl": {"value": 0.088007}, + "log.final.Mars.Zobl": {"value": 0.904732}, + "log.final.Mars.Radius": {"value": 3.393277e+06, "unit": u.m}, + "log.final.Mars.RadGyra": {"value": 0.500000}, + "log.final.Mars.RotAngMom": {"value": 1.308632e+32, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.final.Mars.RotKinEnergy": {"value": 4.638086e+27, "unit": u.Joule}, + "log.final.Mars.RotVel": {"value": 240.530700, "unit": u.m / u.sec}, + "log.final.Mars.BodyType": {"value": 0.000000}, + "log.final.Mars.RotRate": {"value": 7.088449e-05, "unit": 1 / u.sec}, + "log.final.Mars.RotPer": {"value": 8.863977e+04, "unit": u.sec}, + "log.final.Mars.Density": {"value": 3918.693171, "unit": u.kg / u.m ** 3}, + "log.final.Mars.HEcc": {"value": -0.004457}, + "log.final.Mars.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.final.Mars.HZLimRecVenus": {"value": -1.000000}, + "log.final.Mars.HZLimRunaway": {"value": -1.000000}, + "log.final.Mars.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.final.Mars.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.final.Mars.HZLimEarlyMars": {"value": -1.000000}, + "log.final.Mars.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.final.Mars.KEcc": {"value": 0.097550}, + "log.final.Mars.Eccentricity": {"value": 0.097652}, + "log.final.Mars.OrbEnergy": {"value": -1.867046e+32, "unit": u.Joule}, + "log.final.Mars.MeanMotion": {"value": 1.058606e-07, "unit": 1 / u.sec}, + "log.final.Mars.OrbPeriod": {"value": 5.935338e+07, "unit": u.sec}, + "log.final.Mars.SemiMajorAxis": {"value": 2.279367e+11, "unit": u.m}, + "log.final.Mars.COPP": {"value": -0.001899}, + "log.final.Mars.OrbAngMom": {"value": 3.510508e+39, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.final.Mars.ArgP": {"value": 5.651330, "unit": u.rad}, + "log.final.Mars.Inc": {"value": 0.024189, "unit": u.rad}, + "log.final.Mars.LongA": {"value": 0.586198, "unit": u.rad}, + "log.final.Mars.MeanLongitude": {"value": -0.045657, "unit": u.rad}, + "log.final.Mars.LongP": {"value": 6.237528, "unit": u.rad}, + "log.final.Mars.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.final.Mars.OrbPotEnergy": {"value": -3.734093e+32, "unit": u.Joule}, + "log.final.Mars.Sinc": {"value": 0.012094}, + "log.final.Mars.Pinc": {"value": 0.006690}, + "log.final.Mars.Qinc": {"value": 0.010075}, + "log.final.Mars.DIncDtDistOrb": {"value": 0.024148, "unit": u.rad / u.sec}, + "log.final.Mars.DSincDtDistOrb": {"value": 0.012073, "unit": 1 / u.sec}, + "log.final.Mars.DLongADtDistOrb": {"value": -0.001532, "unit": u.rad / u.sec}, + "log.final.Mars.DLongPDtDistOrb": {"value": 0.001708, "unit": u.rad / u.sec}, + "log.final.Mars.DHeccDtDistOrb": {"value": -0.004291, "unit": 1 / u.sec}, + "log.final.Mars.DKeccDtDistOrb": {"value": 0.097576, "unit": 1 / u.sec}, + "log.final.Mars.DPincDtDistOrb": {"value": 0.006663, "unit": 1 / u.sec}, + "log.final.Mars.DQincDtDistOrb": {"value": 0.010068, "unit": 1 / u.sec}, + "log.final.Mars.CassiniOne": {"value": 0.214771}, + "log.final.Mars.CassiniTwo": {"value": 0.976664}, + "log.final.Mars.DOblDtDistRot": {"value": -4.353171e-15, "unit": u.rad}, + "log.final.Mars.DPrecADtDistRot": {"value": 1.330659e-12, "unit": u.rad / u.sec}, + "log.final.Mars.DXoblDtDistRot": {"value": -1.209612e-13, "unit": 1 / u.sec}, + "log.final.Mars.DYoblDtDistRot": {"value": 5.537936e-13, "unit": 1 / u.sec}, + "log.final.Mars.DZoblDtDistRot": {"value": 1.854371e-15, "unit": 1 / u.sec}, + "log.final.Mars.OblTimeDistRot": {"value": 7.216792e+14, "unit": u.rad}, + "log.final.Mars.PrecATimeDistRot": {"value": 4.721858e+12, "unit": u.sec}, + "log.final.Mars.XoblTimeDistRot": {"value": 8.267111e+12, "unit": 1 / u.sec}, + "log.final.Mars.YoblTimeDistRot": {"value": 1.805727e+12, "unit": 1 / u.sec}, + "log.final.Mars.ZoblTimeDistRot": {"value": 5.392666e+14, "unit": 1 / u.sec}, + "log.final.Mars.DynEllip": {"value": 0.005363}, + "log.final.Mars.PrecFNat": {"value": 1.167289e-12, "unit": u.rad / u.sec}, + "log.final.Jupiter.Mass": {"value": 1.897089e+27, "unit": u.kg}, + "log.final.Jupiter.Radius": {"value": 9.567150e+06, "unit": u.m}, + "log.final.Jupiter.RadGyra": {"value": 0.500000}, + "log.final.Jupiter.RotKinEnergy": {"value": 1.147876e+32, "unit": u.Joule}, + "log.final.Jupiter.RotVel": {"value": 695.742781, "unit": u.m / u.sec}, + "log.final.Jupiter.BodyType": {"value": 0.000000}, + "log.final.Jupiter.Density": {"value": 5.171914e+05, "unit": u.kg / u.m ** 3}, + "log.final.Jupiter.HEcc": {"value": 0.023092}, + "log.final.Jupiter.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.final.Jupiter.HZLimRecVenus": {"value": -1.000000}, + "log.final.Jupiter.HZLimRunaway": {"value": -1.000000}, + "log.final.Jupiter.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.final.Jupiter.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.final.Jupiter.HZLimEarlyMars": {"value": -1.000000}, + "log.final.Jupiter.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.final.Jupiter.KEcc": {"value": 0.049706}, + "log.final.Jupiter.Eccentricity": {"value": 0.054808}, + "log.final.Jupiter.OrbEnergy": {"value": -1.615496e+35, "unit": u.Joule}, + "log.final.Jupiter.MeanMotion": {"value": 1.675590e-08, "unit": 1 / u.sec}, + "log.final.Jupiter.OrbPeriod": {"value": 3.749835e+08, "unit": u.sec}, + "log.final.Jupiter.SemiMajorAxis": {"value": 7.792260e+11, "unit": u.m}, + "log.final.Jupiter.OrbAngMom": {"value": 1.925373e+43, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.final.Jupiter.ArgP": {"value": 4.794915, "unit": u.rad}, + "log.final.Jupiter.Inc": {"value": 0.022522, "unit": u.rad}, + "log.final.Jupiter.LongA": {"value": 1.923184, "unit": u.rad}, + "log.final.Jupiter.MeanLongitude": {"value": 0.434914, "unit": u.rad}, + "log.final.Jupiter.LongP": {"value": 0.434914, "unit": u.rad}, + "log.final.Jupiter.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.final.Jupiter.OrbPotEnergy": {"value": -3.230992e+35, "unit": u.Joule}, + "log.final.Jupiter.Sinc": {"value": 0.011261}, + "log.final.Jupiter.Pinc": {"value": 0.010569}, + "log.final.Jupiter.Qinc": {"value": -0.003886}, + "log.final.Jupiter.DIncDtDistOrb": {"value": 0.022529, "unit": u.rad / u.sec}, + "log.final.Jupiter.DSincDtDistOrb": {"value": 0.011264, "unit": 1 / u.sec}, + "log.final.Jupiter.DLongADtDistOrb": {"value": 0.000759, "unit": u.rad / u.sec}, + "log.final.Jupiter.DLongPDtDistOrb": {"value": 0.000858, "unit": u.rad / u.sec}, + "log.final.Jupiter.DHeccDtDistOrb": {"value": 0.023144, "unit": 1 / u.sec}, + "log.final.Jupiter.DKeccDtDistOrb": {"value": 0.049706, "unit": 1 / u.sec}, + "log.final.Jupiter.DPincDtDistOrb": {"value": 0.010568, "unit": 1 / u.sec}, + "log.final.Jupiter.DQincDtDistOrb": {"value": -0.003896, "unit": 1 / u.sec}, + "log.final.Saturn.Mass": {"value": 5.680075e+26, "unit": u.kg}, + "log.final.Saturn.Radius": {"value": 6.378100e+06, "unit": u.m}, + "log.final.Saturn.RadGyra": {"value": 0.500000}, + "log.final.Saturn.RotKinEnergy": {"value": 1.527492e+31, "unit": u.Joule}, + "log.final.Saturn.RotVel": {"value": 463.828521, "unit": u.m / u.sec}, + "log.final.Saturn.BodyType": {"value": 0.000000}, + "log.final.Saturn.Density": {"value": 5.226266e+05, "unit": u.kg / u.m ** 3}, + "log.final.Saturn.HEcc": {"value": 0.035518}, + "log.final.Saturn.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.final.Saturn.HZLimRecVenus": {"value": -1.000000}, + "log.final.Saturn.HZLimRunaway": {"value": -1.000000}, + "log.final.Saturn.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.final.Saturn.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.final.Saturn.HZLimEarlyMars": {"value": -1.000000}, + "log.final.Saturn.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.final.Saturn.KEcc": {"value": -0.017420}, + "log.final.Saturn.Eccentricity": {"value": 0.039560}, + "log.final.Saturn.OrbEnergy": {"value": -2.640963e+34, "unit": u.Joule}, + "log.final.Saturn.MeanMotion": {"value": 6.757836e-09, "unit": 1 / u.sec}, + "log.final.Saturn.OrbPeriod": {"value": 9.297629e+08, "unit": u.sec}, + "log.final.Saturn.SemiMajorAxis": {"value": 1.427163e+12, "unit": u.m}, + "log.final.Saturn.OrbAngMom": {"value": 7.809883e+42, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.final.Saturn.ArgP": {"value": 0.265956, "unit": u.rad}, + "log.final.Saturn.Inc": {"value": 0.043735, "unit": u.rad}, + "log.final.Saturn.LongA": {"value": 1.760834, "unit": u.rad}, + "log.final.Saturn.MeanLongitude": {"value": 2.026790, "unit": u.rad}, + "log.final.Saturn.LongP": {"value": 2.026790, "unit": u.rad}, + "log.final.Saturn.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.final.Saturn.OrbPotEnergy": {"value": -5.281926e+34, "unit": u.Joule}, + "log.final.Saturn.Sinc": {"value": 0.021866}, + "log.final.Saturn.Pinc": {"value": 0.021472}, + "log.final.Saturn.Qinc": {"value": -0.004130}, + "log.final.Saturn.DIncDtDistOrb": {"value": 0.043735, "unit": u.rad / u.sec}, + "log.final.Saturn.DSincDtDistOrb": {"value": 0.021862, "unit": 1 / u.sec}, + "log.final.Saturn.DLongADtDistOrb": {"value": -0.001040, "unit": u.rad / u.sec}, + "log.final.Saturn.DLongPDtDistOrb": {"value": 0.002124, "unit": u.rad / u.sec}, + "log.final.Saturn.DHeccDtDistOrb": {"value": 0.035415, "unit": 1 / u.sec}, + "log.final.Saturn.DKeccDtDistOrb": {"value": -0.017464, "unit": 1 / u.sec}, + "log.final.Saturn.DPincDtDistOrb": {"value": 0.021473, "unit": 1 / u.sec}, + "log.final.Saturn.DQincDtDistOrb": {"value": -0.004107, "unit": 1 / u.sec}, + "log.final.George.Mass": {"value": 8.676288e+25, "unit": u.kg}, + "log.final.George.Radius": {"value": 9.567150e+06, "unit": u.m}, + "log.final.George.RadGyra": {"value": 0.500000}, + "log.final.George.RotKinEnergy": {"value": 5.249784e+30, "unit": u.Joule}, + "log.final.George.RotVel": {"value": 695.742781, "unit": u.m / u.sec}, + "log.final.George.BodyType": {"value": 0.000000}, + "log.final.George.Density": {"value": 2.365362e+04, "unit": u.kg / u.m ** 3}, + "log.final.George.HEcc": {"value": 0.004113}, + "log.final.George.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.final.George.HZLimRecVenus": {"value": -1.000000}, + "log.final.George.HZLimRunaway": {"value": -1.000000}, + "log.final.George.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.final.George.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.final.George.HZLimEarlyMars": {"value": -1.000000}, + "log.final.George.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.final.George.KEcc": {"value": -0.045830}, + "log.final.George.Eccentricity": {"value": 0.046014}, + "log.final.George.OrbEnergy": {"value": -2.005253e+33, "unit": u.Joule}, + "log.final.George.MeanMotion": {"value": 2.368075e-09, "unit": 1 / u.sec}, + "log.final.George.OrbPeriod": {"value": 2.653288e+09, "unit": u.sec}, + "log.final.George.SemiMajorAxis": {"value": 2.871088e+12, "unit": u.m}, + "log.final.George.OrbAngMom": {"value": 1.691778e+42, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.final.George.ArgP": {"value": 1.675289, "unit": u.rad}, + "log.final.George.Inc": {"value": 0.012182, "unit": u.rad}, + "log.final.George.LongA": {"value": 1.376807, "unit": u.rad}, + "log.final.George.MeanLongitude": {"value": 3.052096, "unit": u.rad}, + "log.final.George.LongP": {"value": 3.052096, "unit": u.rad}, + "log.final.George.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.final.George.OrbPotEnergy": {"value": -4.010505e+33, "unit": u.Joule}, + "log.final.George.Sinc": {"value": 0.006091}, + "log.final.George.Pinc": {"value": 0.005976}, + "log.final.George.Qinc": {"value": 0.001174}, + "log.final.George.DIncDtDistOrb": {"value": 0.012177, "unit": u.rad / u.sec}, + "log.final.George.DSincDtDistOrb": {"value": 0.006088, "unit": 1 / u.sec}, + "log.final.George.DLongADtDistOrb": {"value": 0.000452, "unit": u.rad / u.sec}, + "log.final.George.DLongPDtDistOrb": {"value": 0.000308, "unit": u.rad / u.sec}, + "log.final.George.DHeccDtDistOrb": {"value": 0.004098, "unit": 1 / u.sec}, + "log.final.George.DKeccDtDistOrb": {"value": -0.045827, "unit": 1 / u.sec}, + "log.final.George.DPincDtDistOrb": {"value": 0.005975, "unit": 1 / u.sec}, + "log.final.George.DQincDtDistOrb": {"value": 0.001171, "unit": 1 / u.sec}, + "log.final.Neptune.Mass": {"value": 1.023485e+26, "unit": u.kg}, + "log.final.Neptune.Radius": {"value": 9.567150e+06, "unit": u.m}, + "log.final.Neptune.RadGyra": {"value": 0.500000}, + "log.final.Neptune.RotKinEnergy": {"value": 6.192825e+30, "unit": u.Joule}, + "log.final.Neptune.RotVel": {"value": 695.742781, "unit": u.m / u.sec}, + "log.final.Neptune.BodyType": {"value": 0.000000}, + "log.final.Neptune.Density": {"value": 2.790262e+04, "unit": u.kg / u.m ** 3}, + "log.final.Neptune.HEcc": {"value": 0.006336}, + "log.final.Neptune.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.final.Neptune.HZLimRecVenus": {"value": -1.000000}, + "log.final.Neptune.HZLimRunaway": {"value": -1.000000}, + "log.final.Neptune.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.final.Neptune.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.final.Neptune.HZLimEarlyMars": {"value": -1.000000}, + "log.final.Neptune.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.final.Neptune.KEcc": {"value": 0.006138}, + "log.final.Neptune.Eccentricity": {"value": 0.008822}, + "log.final.Neptune.OrbEnergy": {"value": -1.509722e+33, "unit": u.Joule}, + "log.final.Neptune.MeanMotion": {"value": 1.207446e-09, "unit": 1 / u.sec}, + "log.final.Neptune.OrbPeriod": {"value": 5.203698e+09, "unit": u.sec}, + "log.final.Neptune.SemiMajorAxis": {"value": 4.498484e+12, "unit": u.m}, + "log.final.Neptune.OrbAngMom": {"value": 2.500588e+42, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.final.Neptune.ArgP": {"value": 4.790721, "unit": u.rad}, + "log.final.Neptune.Inc": {"value": 0.031082, "unit": u.rad}, + "log.final.Neptune.LongA": {"value": 2.293691, "unit": u.rad}, + "log.final.Neptune.MeanLongitude": {"value": 0.801226, "unit": u.rad}, + "log.final.Neptune.LongP": {"value": 0.801226, "unit": u.rad}, + "log.final.Neptune.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.final.Neptune.OrbPotEnergy": {"value": -3.019443e+33, "unit": u.Joule}, + "log.final.Neptune.Sinc": {"value": 0.015540}, + "log.final.Neptune.Pinc": {"value": 0.011654}, + "log.final.Neptune.Qinc": {"value": -0.010281}, + "log.final.Neptune.DIncDtDistOrb": {"value": 0.031085, "unit": u.rad / u.sec}, + "log.final.Neptune.DSincDtDistOrb": {"value": 0.015541, "unit": 1 / u.sec}, + "log.final.Neptune.DLongADtDistOrb": {"value": -2.576494e-05, "unit": u.rad / u.sec}, + "log.final.Neptune.DLongPDtDistOrb": {"value": 0.000114, "unit": u.rad / u.sec}, + "log.final.Neptune.DHeccDtDistOrb": {"value": 0.006337, "unit": 1 / u.sec}, + "log.final.Neptune.DKeccDtDistOrb": {"value": 0.006138, "unit": 1 / u.sec}, + "log.final.Neptune.DPincDtDistOrb": {"value": 0.011654, "unit": 1 / u.sec}, + "log.final.Neptune.DQincDtDistOrb": {"value": -0.010281, "unit": 1 / u.sec}, + } +) +class Test_LL2(Benchmark): + pass diff --git a/tests/Distorb/LL2/venus.in b/tests/Distorb/LL2/venus.in new file mode 100644 index 000000000..b3996f95d --- /dev/null +++ b/tests/Distorb/LL2/venus.in @@ -0,0 +1,17 @@ +sName Venus +saModules distorb distrot +dMass 2.4464824121e-06 +dRadius -0.9499 +dRotPeriod -243.69 +dObliquity 177.36 +dSemi 0.72333377 +dEcc 0.00677478 +dInc 3.39471 +dLongP 131.54907 +dLongA 76.68069 +bGRCorr 1 +saOutputOrder Time Ecce Inc ArgP LongA Obliq PrecA +dDynEllip 1.3125e-5 #based on J2 from Yoder 1995 +dPrecA 0.0 +#bOverrideMaxEcc 1 +#dHaltMaxEcc 0.4 diff --git a/tests/Distorb/LL2/vpl.in b/tests/Distorb/LL2/vpl.in new file mode 100644 index 000000000..d732ee11d --- /dev/null +++ b/tests/Distorb/LL2/vpl.in @@ -0,0 +1,19 @@ +sSystemName solarsys +iVerbose 5 +iDigits 6 +bOverwrite 1 +sUnitMass solar +sUnitLength au +sUnitTime y +sUnitAngle d +bDoLog 1 +saBodyFiles sun.in mercury.in venus.in earth.in mars.in jupiter.in $ + saturn.in george.in neptune.in +#saBodyFiles sun.in venus.in earth.in mars.in jupiter.in saturn.in george.in neptune.in +#saBodyFiles sun.in mercury.in jupiter.in saturn.in + +bDoForward 1 +bVarDt 1 +dEta 0.01 +dStopTime 5000 +dOutputTime 1000 diff --git a/tests/Distorb/RD4/earth.in b/tests/Distorb/RD4/earth.in new file mode 100644 index 000000000..486903832 --- /dev/null +++ b/tests/Distorb/RD4/earth.in @@ -0,0 +1,19 @@ +sName Earth +saModules distorb distrot +dMass 3.0018090452e-06 +dRadius -1. +dRotPeriod -1 +dObliquity 23.44 #no Moons yet in vplanet.... +dSemi 1.00000321 +dEcc 0.01671327 +dInc 5e-5 +dLongP 102.948601 +dLongA 348.73936 +# bGRCorr 0 +saOutputOrder Time Ecce Inc ArgP LongA Obliq PrecA -TotEn -TotAngMom +dDynEllip 3.28e-3 +dPrecA 0.0 +bForcePrecRate 1 +bCalcDynEllip 0 +#bOverrideMaxEcc 1 +#dHaltMaxEcc 0.4 diff --git a/tests/Distorb/RD4/george.in b/tests/Distorb/RD4/george.in new file mode 100644 index 000000000..956f86877 --- /dev/null +++ b/tests/Distorb/RD4/george.in @@ -0,0 +1,16 @@ +sName George +saModules distorb +dMass 4.3634170854e-05 +dRadius -1.50000 +dRotPeriod -1.00000 +dObliquity 0.00000 +dSemi 19.19203990 +dEcc 0.04713385 +dInc 0.76986 +dLongP 170.994339 +dLongA 74.22988 +#dDfcrit 0.1 +#bGRCorr 0 +saOutputOrder Time Ecce Inc ArgP LongA +#bOverrideMaxEcc 1 +#dHaltMaxEcc 0.4 diff --git a/tests/Distorb/RD4/jupiter.in b/tests/Distorb/RD4/jupiter.in new file mode 100644 index 000000000..158243ee3 --- /dev/null +++ b/tests/Distorb/RD4/jupiter.in @@ -0,0 +1,17 @@ +sName Jupiter +saModules distorb +dMass 9.5407035176e-04 +dRadius -1.50000 +dRotPeriod -1.00000 +dObliquity 0.00000 +dSemi 5.20880408 +dEcc 0.04932699 +dInc 1.30530 +dLongP 15.1623 +dLongA 100.55615 +#dDfcrit 0.1 +#bGRCorr 0 +#sOrbitModel ll2 +saOutputOrder Time Ecce Inc ArgP LongA +#bOverrideMaxEcc 1 +#dHaltMaxEcc 0.4 diff --git a/tests/Distorb/RD4/mars.in b/tests/Distorb/RD4/mars.in new file mode 100644 index 000000000..81f2802e5 --- /dev/null +++ b/tests/Distorb/RD4/mars.in @@ -0,0 +1,17 @@ +sName Mars +saModules distorb distrot +dMass 3.2253768844e-07 +dRadius -0.53202 +dRotPeriod -1.0259233 +dObliquity 25.189417 #obliquity, dynell from Armstrong 2004, Bouquillon & Souchay 1999 +dSemi 1.52366290 +dEcc 0.09341266 +dInc 1.85061 +dLongP 336.040919 +dLongA 49.57854 +# bGRCorr 0 +saOutputOrder Time Ecce Inc ArgP LongA Obliq PrecA PrecFNat +dDynEllip 5.363e-3 +dPrecA 0.0 +#bOverrideMaxEcc 1 +#dHaltMaxEcc 0.4 diff --git a/tests/Distorb/RD4/mercury.in b/tests/Distorb/RD4/mercury.in new file mode 100644 index 000000000..e8bad6f62 --- /dev/null +++ b/tests/Distorb/RD4/mercury.in @@ -0,0 +1,22 @@ +sName Mercury +saModules distorb +dMass 1.6592964824e-07 +dRadius -1.00 +dRotPeriod -1.00000 +dObliquity 23.5 +dSemi 0.38709897 +dEcc 0.20563056 +dInc 7.00487 +dLongP 77.456453 +dLongA 48.33167 +bGRCorr 1 +sOrbitModel rd4 +#sOrbitModel ll2 +saOutputOrder Time Ecce Inc ArgP LongA +#bOverrideMaxEcc 1 +#dHaltMaxEcc 0.4 +#bOutputLapl 1 +bHaltHillStab 0 +bHaltCloseEnc 1 +bInvPlane 0 +bOutputEigen 0 diff --git a/tests/Distorb/RD4/neptune.in b/tests/Distorb/RD4/neptune.in new file mode 100644 index 000000000..b59e59a24 --- /dev/null +++ b/tests/Distorb/RD4/neptune.in @@ -0,0 +1,16 @@ +sName Neptune +saModules distorb +dMass 5.1472361809e-05 +dRadius -1.50000 +dRotPeriod -1.00000 +dObliquity 0.00000 +dSemi 30.07050641 +dEcc 0.00857659 +dInc 1.76917 +dLongP 44.633872 +dLongA 131.72169 +#dDfcrit 0.1 +#bGRCorr 0 +saOutputOrder Time Ecce Inc ArgP LongA +#bOverrideMaxEcc 1 +#dHaltMaxEcc 0.4 diff --git a/tests/Distorb/RD4/saturn.in b/tests/Distorb/RD4/saturn.in new file mode 100644 index 000000000..41bc9365b --- /dev/null +++ b/tests/Distorb/RD4/saturn.in @@ -0,0 +1,15 @@ +sName Saturn +saModules distorb +dMass 2.8565829146e-04 +dRadius -1.00 #radius, rotper, and obl don't matter without distrot and/or eqtide +dRotPeriod -1.00000 +dObliquity 23.5 +dSemi 9.53999265 +dEcc 0.05434133 +dInc 2.48446 +dLongP 92.193712 +dLongA 113.71504 +# bGRCorr 0 +saOutputOrder Time Ecce Inc ArgP LongA +#bOverrideMaxEcc 1 +#dHaltMaxEcc 0.4 diff --git a/tests/Distorb/RD4/sun.in b/tests/Distorb/RD4/sun.in new file mode 100644 index 000000000..62b1f88fc --- /dev/null +++ b/tests/Distorb/RD4/sun.in @@ -0,0 +1,6 @@ +# sun parameters +sName sun +dMass 1 +dSemi 0 +dEcc 0 +dRadius 0.00135 diff --git a/tests/Distorb/RD4/test_RD4.py b/tests/Distorb/RD4/test_RD4.py new file mode 100644 index 000000000..7fe51fe38 --- /dev/null +++ b/tests/Distorb/RD4/test_RD4.py @@ -0,0 +1,825 @@ +from benchmark import Benchmark, benchmark +import astropy.units as u + +@benchmark( + { + "log.initial.system.Age": {"value": 0.000000, "unit": u.sec}, + "log.initial.system.Time": {"value": 0.000000, "unit": u.sec}, + "log.initial.system.TotAngMom": {"value": 3.277987e+43, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.system.TotEnergy": {"value": -7.839546e+41, "unit": u.Joule}, + "log.initial.system.PotEnergy": {"value": -7.839908e+41, "unit": u.Joule}, + "log.initial.system.KinEnergy": {"value": 5.361272e+37, "unit": u.Joule}, + "log.initial.system.DeltaTime": {"value": 0.000000, "unit": u.sec}, + "log.initial.sun.Mass": {"value": 1.988416e+30, "unit": u.kg}, + "log.initial.sun.Radius": {"value": 2.019571e+08, "unit": u.m}, + "log.initial.sun.RadGyra": {"value": 0.500000}, + "log.initial.sun.BodyType": {"value": 0.000000}, + "log.initial.sun.Density": {"value": 5.762900e+04, "unit": u.kg / u.m ** 3}, + "log.initial.sun.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.initial.sun.HZLimRecVenus": {"value": -1.000000}, + "log.initial.sun.HZLimRunaway": {"value": -1.000000}, + "log.initial.sun.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.initial.sun.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.initial.sun.HZLimEarlyMars": {"value": -1.000000}, + "log.initial.sun.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.initial.Mercury.Mass": {"value": 3.299372e+23, "unit": u.kg}, + "log.initial.Mercury.Radius": {"value": 6.378100e+06, "unit": u.m}, + "log.initial.Mercury.RadGyra": {"value": 0.500000}, + "log.initial.Mercury.RotKinEnergy": {"value": 8.872707e+27, "unit": u.Joule}, + "log.initial.Mercury.RotVel": {"value": 463.828521, "unit": u.m / u.sec}, + "log.initial.Mercury.BodyType": {"value": 0.000000}, + "log.initial.Mercury.Density": {"value": 303.576857, "unit": u.kg / u.m ** 3}, + "log.initial.Mercury.HEcc": {"value": 0.200722}, + "log.initial.Mercury.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.initial.Mercury.HZLimRecVenus": {"value": -1.000000}, + "log.initial.Mercury.HZLimRunaway": {"value": -1.000000}, + "log.initial.Mercury.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.initial.Mercury.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.initial.Mercury.HZLimEarlyMars": {"value": -1.000000}, + "log.initial.Mercury.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.initial.Mercury.KEcc": {"value": 0.044659}, + "log.initial.Mercury.Eccentricity": {"value": 0.205631}, + "log.initial.Mercury.OrbEnergy": {"value": -3.780641e+32, "unit": u.Joule}, + "log.initial.Mercury.MeanMotion": {"value": 8.266749e-07, "unit": 1 / u.sec}, + "log.initial.Mercury.OrbPeriod": {"value": 7.600552e+06, "unit": u.sec}, + "log.initial.Mercury.SemiMajorAxis": {"value": 5.790918e+10, "unit": u.m}, + "log.initial.Mercury.OrbAngMom": {"value": 8.951156e+38, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.Mercury.ArgP": {"value": 0.508323, "unit": u.rad}, + "log.initial.Mercury.Inc": {"value": 0.122258, "unit": u.rad}, + "log.initial.Mercury.LongA": {"value": 0.843547, "unit": u.rad}, + "log.initial.Mercury.MeanLongitude": {"value": 1.351870, "unit": u.rad}, + "log.initial.Mercury.LongP": {"value": 1.351870, "unit": u.rad}, + "log.initial.Mercury.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.initial.Mercury.OrbPotEnergy": {"value": -7.561283e+32, "unit": u.Joule}, + "log.initial.Mercury.Sinc": {"value": 0.061091}, + "log.initial.Mercury.Pinc": {"value": 0.045635}, + "log.initial.Mercury.Qinc": {"value": 0.040614}, + "log.initial.Mercury.DIncDtDistOrb": {"value": -3.236872e-14, "unit": u.rad / u.sec}, + "log.initial.Mercury.DSincDtDistOrb": {"value": -1.615413e-14, "unit": 1 / u.sec}, + "log.initial.Mercury.DLongADtDistOrb": {"value": -6.892307e-13, "unit": u.rad / u.sec}, + "log.initial.Mercury.DLongPDtDistOrb": {"value": 8.053323e-13, "unit": u.rad / u.sec}, + "log.initial.Mercury.DHeccDtDistOrb": {"value": 4.211759e-14, "unit": 1 / u.sec}, + "log.initial.Mercury.DKeccDtDistOrb": {"value": -1.602794e-13, "unit": 1 / u.sec}, + "log.initial.Mercury.DPincDtDistOrb": {"value": -4.005988e-14, "unit": 1 / u.sec}, + "log.initial.Mercury.DQincDtDistOrb": {"value": 2.071370e-14, "unit": 1 / u.sec}, + "log.initial.Venus.Mass": {"value": 4.864625e+24, "unit": u.kg}, + "log.initial.Venus.Obliquity": {"value": 3.095516, "unit": u.rad}, + "log.initial.Venus.PrecA": {"value": 0.000000, "unit": u.rad}, + "log.initial.Venus.Xobl": {"value": 0.046060}, + "log.initial.Venus.Yobl": {"value": 0.000000}, + "log.initial.Venus.Zobl": {"value": -0.998939}, + "log.initial.Venus.Radius": {"value": 6.058557e+06, "unit": u.m}, + "log.initial.Venus.RadGyra": {"value": 0.500000}, + "log.initial.Venus.RotAngMom": {"value": 1.332159e+31, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.Venus.RotKinEnergy": {"value": 1.987717e+24, "unit": u.Joule}, + "log.initial.Venus.RotVel": {"value": 1.807997, "unit": u.m / u.sec}, + "log.initial.Venus.BodyType": {"value": 0.000000}, + "log.initial.Venus.RotRate": {"value": 2.984203e-07, "unit": 1 / u.sec}, + "log.initial.Venus.RotPer": {"value": 2.105482e+07, "unit": u.sec}, + "log.initial.Venus.Density": {"value": 5222.195213, "unit": u.kg / u.m ** 3}, + "log.initial.Venus.HEcc": {"value": 0.005070}, + "log.initial.Venus.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.initial.Venus.HZLimRecVenus": {"value": -1.000000}, + "log.initial.Venus.HZLimRunaway": {"value": -1.000000}, + "log.initial.Venus.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.initial.Venus.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.initial.Venus.HZLimEarlyMars": {"value": -1.000000}, + "log.initial.Venus.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.initial.Venus.KEcc": {"value": -0.004493}, + "log.initial.Venus.Eccentricity": {"value": 0.006775}, + "log.initial.Venus.OrbEnergy": {"value": -2.983094e+33, "unit": u.Joule}, + "log.initial.Venus.MeanMotion": {"value": 3.236385e-07, "unit": 1 / u.sec}, + "log.initial.Venus.OrbPeriod": {"value": 1.941421e+07, "unit": u.sec}, + "log.initial.Venus.SemiMajorAxis": {"value": 1.082092e+11, "unit": u.m}, + "log.initial.Venus.COPP": {"value": 0.000234}, + "log.initial.Venus.OrbAngMom": {"value": 1.843430e+40, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.Venus.ArgP": {"value": 0.957634, "unit": u.rad}, + "log.initial.Venus.Inc": {"value": 0.059249, "unit": u.rad}, + "log.initial.Venus.LongA": {"value": 1.338331, "unit": u.rad}, + "log.initial.Venus.MeanLongitude": {"value": 2.295964, "unit": u.rad}, + "log.initial.Venus.LongP": {"value": 2.295964, "unit": u.rad}, + "log.initial.Venus.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.initial.Venus.OrbPotEnergy": {"value": -5.966187e+33, "unit": u.Joule}, + "log.initial.Venus.Sinc": {"value": 0.029620}, + "log.initial.Venus.Pinc": {"value": 0.028823}, + "log.initial.Venus.Qinc": {"value": 0.006824}, + "log.initial.Venus.DIncDtDistOrb": {"value": -4.771288e-15, "unit": u.rad / u.sec}, + "log.initial.Venus.DSincDtDistOrb": {"value": -2.384597e-15, "unit": 1 / u.sec}, + "log.initial.Venus.DLongADtDistOrb": {"value": -1.520890e-12, "unit": u.rad / u.sec}, + "log.initial.Venus.DLongPDtDistOrb": {"value": 2.616045e-14, "unit": u.rad / u.sec}, + "log.initial.Venus.DHeccDtDistOrb": {"value": -1.139880e-14, "unit": 1 / u.sec}, + "log.initial.Venus.DKeccDtDistOrb": {"value": 9.865404e-15, "unit": 1 / u.sec}, + "log.initial.Venus.DPincDtDistOrb": {"value": -1.269872e-14, "unit": 1 / u.sec}, + "log.initial.Venus.DQincDtDistOrb": {"value": 4.328781e-14, "unit": 1 / u.sec}, + "log.initial.Venus.CassiniOne": {"value": 0.819026}, + "log.initial.Venus.CassiniTwo": {"value": 0.573757}, + "log.initial.Venus.DOblDtDistRot": {"value": -8.616948e-14, "unit": u.rad}, + "log.initial.Venus.DPrecADtDistRot": {"value": -6.349869e-12, "unit": u.rad / u.sec}, + "log.initial.Venus.DXoblDtDistRot": {"value": 8.644482e-14, "unit": 1 / u.sec}, + "log.initial.Venus.DYoblDtDistRot": {"value": -2.924774e-13, "unit": 1 / u.sec}, + "log.initial.Venus.DZoblDtDistRot": {"value": -3.985913e-15, "unit": 1 / u.sec}, + "log.initial.Venus.OblTimeDistRot": {"value": 3.645830e+13, "unit": u.rad}, + "log.initial.Venus.PrecATimeDistRot": {"value": 9.894984e+11, "unit": u.sec}, + "log.initial.Venus.XoblTimeDistRot": {"value": 1.156807e+13, "unit": 1 / u.sec}, + "log.initial.Venus.YoblTimeDistRot": {"value": 3.419067e+12, "unit": 1 / u.sec}, + "log.initial.Venus.ZoblTimeDistRot": {"value": 2.508836e+14, "unit": 1 / u.sec}, + "log.initial.Venus.DynEllip": {"value": 1.312500e-05}, + "log.initial.Venus.PrecFNat": {"value": -6.903194e-12, "unit": u.rad / u.sec}, + "log.initial.Earth.Mass": {"value": 5.968845e+24, "unit": u.kg}, + "log.initial.Earth.Obliquity": {"value": 0.409105, "unit": u.rad}, + "log.initial.Earth.PrecA": {"value": 0.000000, "unit": u.rad}, + "log.initial.Earth.Xobl": {"value": 0.397789}, + "log.initial.Earth.Yobl": {"value": 0.000000}, + "log.initial.Earth.Zobl": {"value": 0.917477}, + "log.initial.Earth.Radius": {"value": 6.378100e+06, "unit": u.m}, + "log.initial.Earth.RadGyra": {"value": 0.500000}, + "log.initial.Earth.RotAngMom": {"value": 4.414475e+33, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.Earth.RotKinEnergy": {"value": 1.605149e+29, "unit": u.Joule}, + "log.initial.Earth.RotVel": {"value": 463.828521, "unit": u.m / u.sec}, + "log.initial.Earth.BodyType": {"value": 0.000000}, + "log.initial.Earth.RotRate": {"value": 7.272205e-05, "unit": 1 / u.sec}, + "log.initial.Earth.RotPer": {"value": 8.640000e+04, "unit": u.sec}, + "log.initial.Earth.Density": {"value": 5491.964602, "unit": u.kg / u.m ** 3}, + "log.initial.Earth.HEcc": {"value": 0.016288}, + "log.initial.Earth.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.initial.Earth.HZLimRecVenus": {"value": -1.000000}, + "log.initial.Earth.HZLimRunaway": {"value": -1.000000}, + "log.initial.Earth.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.initial.Earth.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.initial.Earth.HZLimEarlyMars": {"value": -1.000000}, + "log.initial.Earth.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.initial.Earth.KEcc": {"value": -0.003745}, + "log.initial.Earth.Eccentricity": {"value": 0.016713}, + "log.initial.Earth.OrbEnergy": {"value": -2.647556e+33, "unit": u.Joule}, + "log.initial.Earth.MeanMotion": {"value": 1.990977e-07, "unit": 1 / u.sec}, + "log.initial.Earth.OrbPeriod": {"value": 3.155830e+07, "unit": u.sec}, + "log.initial.Earth.SemiMajorAxis": {"value": 1.495984e+11, "unit": u.m}, + "log.initial.Earth.COPP": {"value": 0.006479}, + "log.initial.Earth.OrbAngMom": {"value": 2.659183e+40, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.Earth.ArgP": {"value": 1.993327, "unit": u.rad}, + "log.initial.Earth.Inc": {"value": 8.726646e-07, "unit": u.rad}, + "log.initial.Earth.LongA": {"value": 6.086650, "unit": u.rad}, + "log.initial.Earth.MeanLongitude": {"value": 1.796792, "unit": u.rad}, + "log.initial.Earth.LongP": {"value": 1.796792, "unit": u.rad}, + "log.initial.Earth.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.initial.Earth.OrbPotEnergy": {"value": -5.295112e+33, "unit": u.Joule}, + "log.initial.Earth.Sinc": {"value": 4.363323e-07}, + "log.initial.Earth.Pinc": {"value": -8.520368e-08}, + "log.initial.Earth.Qinc": {"value": 4.279325e-07}, + "log.initial.Earth.DIncDtDistOrb": {"value": -7.164122e-14, "unit": u.rad / u.sec}, + "log.initial.Earth.DSincDtDistOrb": {"value": -3.582061e-14, "unit": 1 / u.sec}, + "log.initial.Earth.DLongADtDistOrb": {"value": -8.868033e-09, "unit": u.rad / u.sec}, + "log.initial.Earth.DLongPDtDistOrb": {"value": 1.744957e-12, "unit": u.rad / u.sec}, + "log.initial.Earth.DHeccDtDistOrb": {"value": -1.966253e-14, "unit": 1 / u.sec}, + "log.initial.Earth.DKeccDtDistOrb": {"value": -2.540401e-14, "unit": 1 / u.sec}, + "log.initial.Earth.DPincDtDistOrb": {"value": 3.199859e-15, "unit": 1 / u.sec}, + "log.initial.Earth.DQincDtDistOrb": {"value": -3.588662e-14, "unit": 1 / u.sec}, + "log.initial.Earth.CassiniOne": {"value": 0.953009}, + "log.initial.Earth.CassiniTwo": {"value": 0.302943}, + "log.initial.Earth.DOblDtDistRot": {"value": -7.177324e-14, "unit": u.rad}, + "log.initial.Earth.DPrecADtDistRot": {"value": 7.711339e-12, "unit": u.rad / u.sec}, + "log.initial.Earth.DXoblDtDistRot": {"value": -6.585030e-14, "unit": 1 / u.sec}, + "log.initial.Earth.DYoblDtDistRot": {"value": 3.067482e-12, "unit": 1 / u.sec}, + "log.initial.Earth.DZoblDtDistRot": {"value": 2.855057e-14, "unit": 1 / u.sec}, + "log.initial.Earth.OblTimeDistRot": {"value": 4.377109e+13, "unit": u.rad}, + "log.initial.Earth.PrecATimeDistRot": {"value": 8.147982e+11, "unit": u.sec}, + "log.initial.Earth.XoblTimeDistRot": {"value": 1.518596e+13, "unit": 1 / u.sec}, + "log.initial.Earth.YoblTimeDistRot": {"value": 3.260003e+11, "unit": 1 / u.sec}, + "log.initial.Earth.ZoblTimeDistRot": {"value": 3.502557e+13, "unit": 1 / u.sec}, + "log.initial.Earth.DynEllip": {"value": 0.003280}, + "log.initial.Earth.PrecFNat": {"value": 2.461543e-12, "unit": u.rad / u.sec}, + "log.initial.Mars.Mass": {"value": 6.413391e+23, "unit": u.kg}, + "log.initial.Mars.Obliquity": {"value": 0.439638, "unit": u.rad}, + "log.initial.Mars.PrecA": {"value": 0.000000, "unit": u.rad}, + "log.initial.Mars.Xobl": {"value": 0.425612}, + "log.initial.Mars.Yobl": {"value": 0.000000}, + "log.initial.Mars.Zobl": {"value": 0.904906}, + "log.initial.Mars.Radius": {"value": 3.393277e+06, "unit": u.m}, + "log.initial.Mars.RadGyra": {"value": 0.500000}, + "log.initial.Mars.RotAngMom": {"value": 1.308632e+32, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.Mars.RotKinEnergy": {"value": 4.638086e+27, "unit": u.Joule}, + "log.initial.Mars.RotVel": {"value": 240.530700, "unit": u.m / u.sec}, + "log.initial.Mars.BodyType": {"value": 0.000000}, + "log.initial.Mars.RotRate": {"value": 7.088449e-05, "unit": 1 / u.sec}, + "log.initial.Mars.RotPer": {"value": 8.863977e+04, "unit": u.sec}, + "log.initial.Mars.Density": {"value": 3918.693171, "unit": u.kg / u.m ** 3}, + "log.initial.Mars.HEcc": {"value": -0.037933}, + "log.initial.Mars.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.initial.Mars.HZLimRecVenus": {"value": -1.000000}, + "log.initial.Mars.HZLimRunaway": {"value": -1.000000}, + "log.initial.Mars.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.initial.Mars.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.initial.Mars.HZLimEarlyMars": {"value": -1.000000}, + "log.initial.Mars.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.initial.Mars.KEcc": {"value": 0.085364}, + "log.initial.Mars.Eccentricity": {"value": 0.093413}, + "log.initial.Mars.OrbEnergy": {"value": -1.867046e+32, "unit": u.Joule}, + "log.initial.Mars.MeanMotion": {"value": 1.058606e-07, "unit": 1 / u.sec}, + "log.initial.Mars.OrbPeriod": {"value": 5.935338e+07, "unit": u.sec}, + "log.initial.Mars.SemiMajorAxis": {"value": 2.279367e+11, "unit": u.m}, + "log.initial.Mars.COPP": {"value": -0.016145}, + "log.initial.Mars.OrbAngMom": {"value": 3.511943e+39, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.Mars.ArgP": {"value": 4.999712, "unit": u.rad}, + "log.initial.Mars.Inc": {"value": 0.032299, "unit": u.rad}, + "log.initial.Mars.LongA": {"value": 0.865309, "unit": u.rad}, + "log.initial.Mars.MeanLongitude": {"value": -0.418165, "unit": u.rad}, + "log.initial.Mars.LongP": {"value": 5.865020, "unit": u.rad}, + "log.initial.Mars.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.initial.Mars.OrbPotEnergy": {"value": -3.734093e+32, "unit": u.Joule}, + "log.initial.Mars.Sinc": {"value": 0.016149}, + "log.initial.Mars.Pinc": {"value": 0.012294}, + "log.initial.Mars.Qinc": {"value": 0.010471}, + "log.initial.Mars.DIncDtDistOrb": {"value": -4.499455e-14, "unit": u.rad / u.sec}, + "log.initial.Mars.DSincDtDistOrb": {"value": -2.249434e-14, "unit": 1 / u.sec}, + "log.initial.Mars.DLongADtDistOrb": {"value": -1.618093e-12, "unit": u.rad / u.sec}, + "log.initial.Mars.DLongPDtDistOrb": {"value": 2.428805e-12, "unit": u.rad / u.sec}, + "log.initial.Mars.DHeccDtDistOrb": {"value": 1.954142e-13, "unit": 1 / u.sec}, + "log.initial.Mars.DKeccDtDistOrb": {"value": 1.189524e-13, "unit": 1 / u.sec}, + "log.initial.Mars.DPincDtDistOrb": {"value": -3.406796e-14, "unit": 1 / u.sec}, + "log.initial.Mars.DQincDtDistOrb": {"value": 5.307548e-15, "unit": 1 / u.sec}, + "log.initial.Mars.CassiniOne": {"value": 0.056899}, + "log.initial.Mars.CassiniTwo": {"value": 0.998380}, + "log.initial.Mars.DOblDtDistRot": {"value": 1.060610e-14, "unit": u.rad}, + "log.initial.Mars.DPrecADtDistRot": {"value": 1.311792e-12, "unit": u.rad / u.sec}, + "log.initial.Mars.DXoblDtDistRot": {"value": 9.597522e-15, "unit": 1 / u.sec}, + "log.initial.Mars.DYoblDtDistRot": {"value": 5.583147e-13, "unit": 1 / u.sec}, + "log.initial.Mars.DZoblDtDistRot": {"value": -4.514086e-15, "unit": 1 / u.sec}, + "log.initial.Mars.OblTimeDistRot": {"value": 2.962061e+14, "unit": u.rad}, + "log.initial.Mars.PrecATimeDistRot": {"value": 4.789772e+12, "unit": u.sec}, + "log.initial.Mars.XoblTimeDistRot": {"value": 1.041936e+14, "unit": 1 / u.sec}, + "log.initial.Mars.YoblTimeDistRot": {"value": 1.791105e+12, "unit": 1 / u.sec}, + "log.initial.Mars.ZoblTimeDistRot": {"value": 2.215288e+14, "unit": 1 / u.sec}, + "log.initial.Mars.DynEllip": {"value": 0.005363}, + "log.initial.Mars.PrecFNat": {"value": 1.166082e-12, "unit": u.rad / u.sec}, + "log.initial.Jupiter.Mass": {"value": 1.897089e+27, "unit": u.kg}, + "log.initial.Jupiter.Radius": {"value": 9.567150e+06, "unit": u.m}, + "log.initial.Jupiter.RadGyra": {"value": 0.500000}, + "log.initial.Jupiter.RotKinEnergy": {"value": 1.147876e+32, "unit": u.Joule}, + "log.initial.Jupiter.RotVel": {"value": 695.742781, "unit": u.m / u.sec}, + "log.initial.Jupiter.BodyType": {"value": 0.000000}, + "log.initial.Jupiter.Density": {"value": 5.171914e+05, "unit": u.kg / u.m ** 3}, + "log.initial.Jupiter.HEcc": {"value": 0.012902}, + "log.initial.Jupiter.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.initial.Jupiter.HZLimRecVenus": {"value": -1.000000}, + "log.initial.Jupiter.HZLimRunaway": {"value": -1.000000}, + "log.initial.Jupiter.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.initial.Jupiter.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.initial.Jupiter.HZLimEarlyMars": {"value": -1.000000}, + "log.initial.Jupiter.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.initial.Jupiter.KEcc": {"value": 0.047610}, + "log.initial.Jupiter.Eccentricity": {"value": 0.049327}, + "log.initial.Jupiter.OrbEnergy": {"value": -1.615496e+35, "unit": u.Joule}, + "log.initial.Jupiter.MeanMotion": {"value": 1.675590e-08, "unit": 1 / u.sec}, + "log.initial.Jupiter.OrbPeriod": {"value": 3.749835e+08, "unit": u.sec}, + "log.initial.Jupiter.SemiMajorAxis": {"value": 7.792260e+11, "unit": u.m}, + "log.initial.Jupiter.OrbAngMom": {"value": 1.925924e+43, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.Jupiter.ArgP": {"value": 4.792781, "unit": u.rad}, + "log.initial.Jupiter.Inc": {"value": 0.022782, "unit": u.rad}, + "log.initial.Jupiter.LongA": {"value": 1.755036, "unit": u.rad}, + "log.initial.Jupiter.MeanLongitude": {"value": 0.264632, "unit": u.rad}, + "log.initial.Jupiter.LongP": {"value": 0.264632, "unit": u.rad}, + "log.initial.Jupiter.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.initial.Jupiter.OrbPotEnergy": {"value": -3.230992e+35, "unit": u.Joule}, + "log.initial.Jupiter.Sinc": {"value": 0.011391}, + "log.initial.Jupiter.Pinc": {"value": 0.011198}, + "log.initial.Jupiter.Qinc": {"value": -0.002087}, + "log.initial.Jupiter.DIncDtDistOrb": {"value": -1.115896e-14, "unit": u.rad / u.sec}, + "log.initial.Jupiter.DSincDtDistOrb": {"value": -5.579119e-15, "unit": 1 / u.sec}, + "log.initial.Jupiter.DLongADtDistOrb": {"value": 9.914453e-13, "unit": u.rad / u.sec}, + "log.initial.Jupiter.DLongPDtDistOrb": {"value": 9.811600e-13, "unit": u.rad / u.sec}, + "log.initial.Jupiter.DHeccDtDistOrb": {"value": 5.722680e-14, "unit": 1 / u.sec}, + "log.initial.Jupiter.DKeccDtDistOrb": {"value": 2.613989e-14, "unit": 1 / u.sec}, + "log.initial.Jupiter.DPincDtDistOrb": {"value": -7.553601e-15, "unit": 1 / u.sec}, + "log.initial.Jupiter.DQincDtDistOrb": {"value": -1.007998e-14, "unit": 1 / u.sec}, + "log.initial.Saturn.Mass": {"value": 5.680075e+26, "unit": u.kg}, + "log.initial.Saturn.Radius": {"value": 6.378100e+06, "unit": u.m}, + "log.initial.Saturn.RadGyra": {"value": 0.500000}, + "log.initial.Saturn.RotKinEnergy": {"value": 1.527492e+31, "unit": u.Joule}, + "log.initial.Saturn.RotVel": {"value": 463.828521, "unit": u.m / u.sec}, + "log.initial.Saturn.BodyType": {"value": 0.000000}, + "log.initial.Saturn.Density": {"value": 5.226266e+05, "unit": u.kg / u.m ** 3}, + "log.initial.Saturn.HEcc": {"value": 0.054302}, + "log.initial.Saturn.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.initial.Saturn.HZLimRecVenus": {"value": -1.000000}, + "log.initial.Saturn.HZLimRunaway": {"value": -1.000000}, + "log.initial.Saturn.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.initial.Saturn.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.initial.Saturn.HZLimEarlyMars": {"value": -1.000000}, + "log.initial.Saturn.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.initial.Saturn.KEcc": {"value": -0.002080}, + "log.initial.Saturn.Eccentricity": {"value": 0.054341}, + "log.initial.Saturn.OrbEnergy": {"value": -2.640963e+34, "unit": u.Joule}, + "log.initial.Saturn.MeanMotion": {"value": 6.757836e-09, "unit": 1 / u.sec}, + "log.initial.Saturn.OrbPeriod": {"value": 9.297629e+08, "unit": u.sec}, + "log.initial.Saturn.SemiMajorAxis": {"value": 1.427163e+12, "unit": u.m}, + "log.initial.Saturn.OrbAngMom": {"value": 7.804453e+42, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.Saturn.ArgP": {"value": 5.907567, "unit": u.rad}, + "log.initial.Saturn.Inc": {"value": 0.043362, "unit": u.rad}, + "log.initial.Saturn.LongA": {"value": 1.984702, "unit": u.rad}, + "log.initial.Saturn.MeanLongitude": {"value": 1.609084, "unit": u.rad}, + "log.initial.Saturn.LongP": {"value": 1.609084, "unit": u.rad}, + "log.initial.Saturn.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.initial.Saturn.OrbPotEnergy": {"value": -5.281926e+34, "unit": u.Joule}, + "log.initial.Saturn.Sinc": {"value": 0.021679}, + "log.initial.Saturn.Pinc": {"value": 0.019849}, + "log.initial.Saturn.Qinc": {"value": -0.008719}, + "log.initial.Saturn.DIncDtDistOrb": {"value": 1.440129e-14, "unit": u.rad / u.sec}, + "log.initial.Saturn.DSincDtDistOrb": {"value": 7.198953e-15, "unit": 1 / u.sec}, + "log.initial.Saturn.DLongADtDistOrb": {"value": -1.445905e-12, "unit": u.rad / u.sec}, + "log.initial.Saturn.DLongPDtDistOrb": {"value": 2.529529e-12, "unit": u.rad / u.sec}, + "log.initial.Saturn.DHeccDtDistOrb": {"value": -9.349643e-14, "unit": 1 / u.sec}, + "log.initial.Saturn.DKeccDtDistOrb": {"value": -1.339773e-13, "unit": 1 / u.sec}, + "log.initial.Saturn.DPincDtDistOrb": {"value": 1.919813e-14, "unit": 1 / u.sec}, + "log.initial.Saturn.DQincDtDistOrb": {"value": 2.580392e-14, "unit": 1 / u.sec}, + "log.initial.George.Mass": {"value": 8.676288e+25, "unit": u.kg}, + "log.initial.George.Radius": {"value": 9.567150e+06, "unit": u.m}, + "log.initial.George.RadGyra": {"value": 0.500000}, + "log.initial.George.RotKinEnergy": {"value": 5.249784e+30, "unit": u.Joule}, + "log.initial.George.RotVel": {"value": 695.742781, "unit": u.m / u.sec}, + "log.initial.George.BodyType": {"value": 0.000000}, + "log.initial.George.Density": {"value": 2.365362e+04, "unit": u.kg / u.m ** 3}, + "log.initial.George.HEcc": {"value": 0.007378}, + "log.initial.George.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.initial.George.HZLimRecVenus": {"value": -1.000000}, + "log.initial.George.HZLimRunaway": {"value": -1.000000}, + "log.initial.George.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.initial.George.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.initial.George.HZLimEarlyMars": {"value": -1.000000}, + "log.initial.George.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.initial.George.KEcc": {"value": -0.046553}, + "log.initial.George.Eccentricity": {"value": 0.047134}, + "log.initial.George.OrbEnergy": {"value": -2.005253e+33, "unit": u.Joule}, + "log.initial.George.MeanMotion": {"value": 2.368075e-09, "unit": 1 / u.sec}, + "log.initial.George.OrbPeriod": {"value": 2.653288e+09, "unit": u.sec}, + "log.initial.George.SemiMajorAxis": {"value": 2.871088e+12, "unit": u.m}, + "log.initial.George.OrbAngMom": {"value": 1.691690e+42, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.George.ArgP": {"value": 1.688858, "unit": u.rad}, + "log.initial.George.Inc": {"value": 0.013437, "unit": u.rad}, + "log.initial.George.LongA": {"value": 1.295556, "unit": u.rad}, + "log.initial.George.MeanLongitude": {"value": 2.984414, "unit": u.rad}, + "log.initial.George.LongP": {"value": 2.984414, "unit": u.rad}, + "log.initial.George.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.initial.George.OrbPotEnergy": {"value": -4.010505e+33, "unit": u.Joule}, + "log.initial.George.Sinc": {"value": 0.006718}, + "log.initial.George.Pinc": {"value": 0.006465}, + "log.initial.George.Qinc": {"value": 0.001826}, + "log.initial.George.DIncDtDistOrb": {"value": -9.237547e-15, "unit": u.rad / u.sec}, + "log.initial.George.DSincDtDistOrb": {"value": -4.618669e-15, "unit": 1 / u.sec}, + "log.initial.George.DLongADtDistOrb": {"value": 4.153309e-13, "unit": u.rad / u.sec}, + "log.initial.George.DLongPDtDistOrb": {"value": 4.490508e-13, "unit": u.rad / u.sec}, + "log.initial.George.DHeccDtDistOrb": {"value": -2.218761e-14, "unit": 1 / u.sec}, + "log.initial.George.DKeccDtDistOrb": {"value": 4.782430e-15, "unit": 1 / u.sec}, + "log.initial.George.DPincDtDistOrb": {"value": -3.686480e-15, "unit": 1 / u.sec}, + "log.initial.George.DQincDtDistOrb": {"value": -3.940522e-15, "unit": 1 / u.sec}, + "log.initial.Neptune.Mass": {"value": 1.023485e+26, "unit": u.kg}, + "log.initial.Neptune.Radius": {"value": 9.567150e+06, "unit": u.m}, + "log.initial.Neptune.RadGyra": {"value": 0.500000}, + "log.initial.Neptune.RotKinEnergy": {"value": 6.192825e+30, "unit": u.Joule}, + "log.initial.Neptune.RotVel": {"value": 695.742781, "unit": u.m / u.sec}, + "log.initial.Neptune.BodyType": {"value": 0.000000}, + "log.initial.Neptune.Density": {"value": 2.790262e+04, "unit": u.kg / u.m ** 3}, + "log.initial.Neptune.HEcc": {"value": 0.006026}, + "log.initial.Neptune.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.initial.Neptune.HZLimRecVenus": {"value": -1.000000}, + "log.initial.Neptune.HZLimRunaway": {"value": -1.000000}, + "log.initial.Neptune.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.initial.Neptune.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.initial.Neptune.HZLimEarlyMars": {"value": -1.000000}, + "log.initial.Neptune.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.initial.Neptune.KEcc": {"value": 0.006103}, + "log.initial.Neptune.Eccentricity": {"value": 0.008577}, + "log.initial.Neptune.OrbEnergy": {"value": -1.509722e+33, "unit": u.Joule}, + "log.initial.Neptune.MeanMotion": {"value": 1.207446e-09, "unit": 1 / u.sec}, + "log.initial.Neptune.OrbPeriod": {"value": 5.203698e+09, "unit": u.sec}, + "log.initial.Neptune.SemiMajorAxis": {"value": 4.498484e+12, "unit": u.m}, + "log.initial.Neptune.OrbAngMom": {"value": 2.500594e+42, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.Neptune.ArgP": {"value": 4.763216, "unit": u.rad}, + "log.initial.Neptune.Inc": {"value": 0.030878, "unit": u.rad}, + "log.initial.Neptune.LongA": {"value": 2.298977, "unit": u.rad}, + "log.initial.Neptune.MeanLongitude": {"value": 0.779008, "unit": u.rad}, + "log.initial.Neptune.LongP": {"value": 0.779008, "unit": u.rad}, + "log.initial.Neptune.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.initial.Neptune.OrbPotEnergy": {"value": -3.019443e+33, "unit": u.Joule}, + "log.initial.Neptune.Sinc": {"value": 0.015438}, + "log.initial.Neptune.Pinc": {"value": 0.011523}, + "log.initial.Neptune.Qinc": {"value": -0.010274}, + "log.initial.Neptune.DIncDtDistOrb": {"value": 1.231584e-15, "unit": u.rad / u.sec}, + "log.initial.Neptune.DSincDtDistOrb": {"value": 6.157186e-16, "unit": 1 / u.sec}, + "log.initial.Neptune.DLongADtDistOrb": {"value": -3.314407e-14, "unit": u.rad / u.sec}, + "log.initial.Neptune.DLongPDtDistOrb": {"value": 1.268089e-13, "unit": u.rad / u.sec}, + "log.initial.Neptune.DHeccDtDistOrb": {"value": 1.900488e-15, "unit": 1 / u.sec}, + "log.initial.Neptune.DKeccDtDistOrb": {"value": 3.769281e-16, "unit": 1 / u.sec}, + "log.initial.Neptune.DPincDtDistOrb": {"value": 8.000991e-16, "unit": 1 / u.sec}, + "log.initial.Neptune.DQincDtDistOrb": {"value": -2.785158e-17, "unit": 1 / u.sec}, + "log.final.system.Age": {"value": 1.577880e+11, "unit": u.sec}, + "log.final.system.Time": {"value": 1.577880e+11, "unit": u.sec}, + "log.final.system.TotAngMom": {"value": 3.277986e+43, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.final.system.TotEnergy": {"value": -7.839546e+41, "unit": u.Joule}, + "log.final.system.PotEnergy": {"value": -7.839908e+41, "unit": u.Joule}, + "log.final.system.KinEnergy": {"value": 5.361272e+37, "unit": u.Joule}, + "log.final.system.DeltaTime": {"value": 1.434436e+09, "unit": u.sec}, + "log.final.sun.Mass": {"value": 1.988416e+30, "unit": u.kg}, + "log.final.sun.Radius": {"value": 2.019571e+08, "unit": u.m}, + "log.final.sun.RadGyra": {"value": 0.500000}, + "log.final.sun.BodyType": {"value": 0.000000}, + "log.final.sun.Density": {"value": 5.762900e+04, "unit": u.kg / u.m ** 3}, + "log.final.sun.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.final.sun.HZLimRecVenus": {"value": -1.000000}, + "log.final.sun.HZLimRunaway": {"value": -1.000000}, + "log.final.sun.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.final.sun.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.final.sun.HZLimEarlyMars": {"value": -1.000000}, + "log.final.sun.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.final.Mercury.Mass": {"value": 3.299372e+23, "unit": u.kg}, + "log.final.Mercury.Radius": {"value": 6.378100e+06, "unit": u.m}, + "log.final.Mercury.RadGyra": {"value": 0.500000}, + "log.final.Mercury.RotKinEnergy": {"value": 8.872707e+27, "unit": u.Joule}, + "log.final.Mercury.RotVel": {"value": 463.828521, "unit": u.m / u.sec}, + "log.final.Mercury.BodyType": {"value": 0.000000}, + "log.final.Mercury.Density": {"value": 303.576857, "unit": u.kg / u.m ** 3}, + "log.final.Mercury.HEcc": {"value": 0.205843}, + "log.final.Mercury.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.final.Mercury.HZLimRecVenus": {"value": -1.000000}, + "log.final.Mercury.HZLimRunaway": {"value": -1.000000}, + "log.final.Mercury.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.final.Mercury.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.final.Mercury.HZLimEarlyMars": {"value": -1.000000}, + "log.final.Mercury.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.final.Mercury.KEcc": {"value": 0.016917}, + "log.final.Mercury.Eccentricity": {"value": 0.206537}, + "log.final.Mercury.OrbEnergy": {"value": -3.780641e+32, "unit": u.Joule}, + "log.final.Mercury.MeanMotion": {"value": 8.266749e-07, "unit": 1 / u.sec}, + "log.final.Mercury.OrbPeriod": {"value": 7.600552e+06, "unit": u.sec}, + "log.final.Mercury.SemiMajorAxis": {"value": 5.790918e+10, "unit": u.m}, + "log.final.Mercury.OrbAngMom": {"value": 8.949411e+38, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.final.Mercury.ArgP": {"value": 0.758125, "unit": u.rad}, + "log.final.Mercury.Inc": {"value": 0.117269, "unit": u.rad}, + "log.final.Mercury.LongA": {"value": 0.730670, "unit": u.rad}, + "log.final.Mercury.MeanLongitude": {"value": 1.488795, "unit": u.rad}, + "log.final.Mercury.LongP": {"value": 1.488795, "unit": u.rad}, + "log.final.Mercury.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.final.Mercury.OrbPotEnergy": {"value": -7.561283e+32, "unit": u.Joule}, + "log.final.Mercury.Sinc": {"value": 0.058601}, + "log.final.Mercury.Pinc": {"value": 0.039109}, + "log.final.Mercury.Qinc": {"value": 0.043642}, + "log.final.Mercury.DIncDtDistOrb": {"value": -3.036467e-14, "unit": u.rad / u.sec}, + "log.final.Mercury.DSincDtDistOrb": {"value": -1.515625e-14, "unit": 1 / u.sec}, + "log.final.Mercury.DLongADtDistOrb": {"value": -7.437719e-13, "unit": u.rad / u.sec}, + "log.final.Mercury.DLongPDtDistOrb": {"value": 7.980477e-13, "unit": u.rad / u.sec}, + "log.final.Mercury.DHeccDtDistOrb": {"value": 1.853390e-14, "unit": 1 / u.sec}, + "log.final.Mercury.DKeccDtDistOrb": {"value": -1.638587e-13, "unit": 1 / u.sec}, + "log.final.Mercury.DPincDtDistOrb": {"value": -4.257435e-14, "unit": 1 / u.sec}, + "log.final.Mercury.DQincDtDistOrb": {"value": 1.780054e-14, "unit": 1 / u.sec}, + "log.final.Venus.Mass": {"value": 4.864625e+24, "unit": u.kg}, + "log.final.Venus.Obliquity": {"value": 3.087758, "unit": u.rad}, + "log.final.Venus.PrecA": {"value": 5.408448, "unit": u.rad}, + "log.final.Venus.Xobl": {"value": 0.034529}, + "log.final.Venus.Yobl": {"value": -0.041324}, + "log.final.Venus.Zobl": {"value": -0.999328}, + "log.final.Venus.Radius": {"value": 6.058557e+06, "unit": u.m}, + "log.final.Venus.RadGyra": {"value": 0.500000}, + "log.final.Venus.RotAngMom": {"value": 1.332159e+31, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.final.Venus.RotKinEnergy": {"value": 1.987717e+24, "unit": u.Joule}, + "log.final.Venus.RotVel": {"value": 1.807997, "unit": u.m / u.sec}, + "log.final.Venus.BodyType": {"value": 0.000000}, + "log.final.Venus.RotRate": {"value": 2.984203e-07, "unit": 1 / u.sec}, + "log.final.Venus.RotPer": {"value": 2.105482e+07, "unit": u.sec}, + "log.final.Venus.Density": {"value": 5222.195213, "unit": u.kg / u.m ** 3}, + "log.final.Venus.HEcc": {"value": 0.003720}, + "log.final.Venus.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.final.Venus.HZLimRecVenus": {"value": -1.000000}, + "log.final.Venus.HZLimRunaway": {"value": -1.000000}, + "log.final.Venus.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.final.Venus.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.final.Venus.HZLimEarlyMars": {"value": -1.000000}, + "log.final.Venus.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.final.Venus.KEcc": {"value": -0.002882}, + "log.final.Venus.Eccentricity": {"value": 0.004706}, + "log.final.Venus.OrbEnergy": {"value": -2.983094e+33, "unit": u.Joule}, + "log.final.Venus.MeanMotion": {"value": 3.236385e-07, "unit": 1 / u.sec}, + "log.final.Venus.OrbPeriod": {"value": 1.941421e+07, "unit": u.sec}, + "log.final.Venus.SemiMajorAxis": {"value": 1.082092e+11, "unit": u.m}, + "log.final.Venus.COPP": {"value": 0.000200}, + "log.final.Venus.OrbAngMom": {"value": 1.843452e+40, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.final.Venus.ArgP": {"value": 1.138318, "unit": u.rad}, + "log.final.Venus.Inc": {"value": 0.057128, "unit": u.rad}, + "log.final.Venus.LongA": {"value": 1.091682, "unit": u.rad}, + "log.final.Venus.MeanLongitude": {"value": 2.230000, "unit": u.rad}, + "log.final.Venus.LongP": {"value": 2.230000, "unit": u.rad}, + "log.final.Venus.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.final.Venus.OrbPotEnergy": {"value": -5.966187e+33, "unit": u.Joule}, + "log.final.Venus.Sinc": {"value": 0.028560}, + "log.final.Venus.Pinc": {"value": 0.025344}, + "log.final.Venus.Qinc": {"value": 0.013166}, + "log.final.Venus.DIncDtDistOrb": {"value": -2.188034e-14, "unit": u.rad / u.sec}, + "log.final.Venus.DSincDtDistOrb": {"value": -1.093571e-14, "unit": 1 / u.sec}, + "log.final.Venus.DLongADtDistOrb": {"value": -1.609241e-12, "unit": u.rad / u.sec}, + "log.final.Venus.DLongPDtDistOrb": {"value": -9.443380e-13, "unit": u.rad / u.sec}, + "log.final.Venus.DHeccDtDistOrb": {"value": -5.773409e-15, "unit": 1 / u.sec}, + "log.final.Venus.DKeccDtDistOrb": {"value": 1.009562e-14, "unit": 1 / u.sec}, + "log.final.Venus.DPincDtDistOrb": {"value": -3.089155e-14, "unit": 1 / u.sec}, + "log.final.Venus.DQincDtDistOrb": {"value": 3.574359e-14, "unit": 1 / u.sec}, + "log.final.Venus.CassiniOne": {"value": 0.258011}, + "log.final.Venus.CassiniTwo": {"value": 0.966142}, + "log.final.Venus.DOblDtDistRot": {"value": 1.579691e-15, "unit": u.rad}, + "log.final.Venus.DPrecADtDistRot": {"value": -5.151856e-12, "unit": u.rad / u.sec}, + "log.final.Venus.DXoblDtDistRot": {"value": -2.139114e-13, "unit": 1 / u.sec}, + "log.final.Venus.DYoblDtDistRot": {"value": -1.766696e-13, "unit": 1 / u.sec}, + "log.final.Venus.DZoblDtDistRot": {"value": 8.562977e-17, "unit": 1 / u.sec}, + "log.final.Venus.OblTimeDistRot": {"value": 1.988739e+15, "unit": u.rad}, + "log.final.Venus.PrecATimeDistRot": {"value": 1.219596e+12, "unit": u.sec}, + "log.final.Venus.XoblTimeDistRot": {"value": 4.674832e+12, "unit": 1 / u.sec}, + "log.final.Venus.YoblTimeDistRot": {"value": 5.660283e+12, "unit": 1 / u.sec}, + "log.final.Venus.ZoblTimeDistRot": {"value": 1.167818e+16, "unit": 1 / u.sec}, + "log.final.Venus.DynEllip": {"value": 1.312500e-05}, + "log.final.Venus.PrecFNat": {"value": -6.905640e-12, "unit": u.rad / u.sec}, + "log.final.Earth.Mass": {"value": 5.968845e+24, "unit": u.kg}, + "log.final.Earth.Obliquity": {"value": 0.398931, "unit": u.rad}, + "log.final.Earth.PrecA": {"value": 1.226713, "unit": u.rad}, + "log.final.Earth.Xobl": {"value": 0.131032}, + "log.final.Earth.Yobl": {"value": 0.365665}, + "log.final.Earth.Zobl": {"value": 0.921477}, + "log.final.Earth.Radius": {"value": 6.378100e+06, "unit": u.m}, + "log.final.Earth.RadGyra": {"value": 0.500000}, + "log.final.Earth.RotAngMom": {"value": 4.414475e+33, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.final.Earth.RotKinEnergy": {"value": 1.605149e+29, "unit": u.Joule}, + "log.final.Earth.RotVel": {"value": 463.828521, "unit": u.m / u.sec}, + "log.final.Earth.BodyType": {"value": 0.000000}, + "log.final.Earth.RotRate": {"value": 7.272205e-05, "unit": 1 / u.sec}, + "log.final.Earth.RotPer": {"value": 8.640000e+04, "unit": u.sec}, + "log.final.Earth.Density": {"value": 5491.964602, "unit": u.kg / u.m ** 3}, + "log.final.Earth.HEcc": {"value": 0.012499}, + "log.final.Earth.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.final.Earth.HZLimRecVenus": {"value": -1.000000}, + "log.final.Earth.HZLimRunaway": {"value": -1.000000}, + "log.final.Earth.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.final.Earth.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.final.Earth.HZLimEarlyMars": {"value": -1.000000}, + "log.final.Earth.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.final.Earth.KEcc": {"value": -0.006991}, + "log.final.Earth.Eccentricity": {"value": 0.014322}, + "log.final.Earth.OrbEnergy": {"value": -2.647556e+33, "unit": u.Joule}, + "log.final.Earth.MeanMotion": {"value": 1.990977e-07, "unit": 1 / u.sec}, + "log.final.Earth.OrbPeriod": {"value": 3.155830e+07, "unit": u.sec}, + "log.final.Earth.SemiMajorAxis": {"value": 1.495984e+11, "unit": u.m}, + "log.final.Earth.COPP": {"value": 0.004855}, + "log.final.Earth.OrbAngMom": {"value": 2.659282e+40, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.final.Earth.ArgP": {"value": 5.518696, "unit": u.rad}, + "log.final.Earth.Inc": {"value": 0.010890, "unit": u.rad}, + "log.final.Earth.LongA": {"value": 2.845289, "unit": u.rad}, + "log.final.Earth.MeanLongitude": {"value": 2.080800, "unit": u.rad}, + "log.final.Earth.LongP": {"value": 2.080800, "unit": u.rad}, + "log.final.Earth.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.final.Earth.OrbPotEnergy": {"value": -5.295112e+33, "unit": u.Joule}, + "log.final.Earth.Sinc": {"value": 0.005445}, + "log.final.Earth.Pinc": {"value": 0.001590}, + "log.final.Earth.Qinc": {"value": -0.005208}, + "log.final.Earth.DIncDtDistOrb": {"value": 6.554951e-14, "unit": u.rad / u.sec}, + "log.final.Earth.DSincDtDistOrb": {"value": 3.277427e-14, "unit": 1 / u.sec}, + "log.final.Earth.DLongADtDistOrb": {"value": -1.305094e-12, "unit": u.rad / u.sec}, + "log.final.Earth.DLongPDtDistOrb": {"value": 1.863720e-12, "unit": u.rad / u.sec}, + "log.final.Earth.DHeccDtDistOrb": {"value": -2.764468e-14, "unit": 1 / u.sec}, + "log.final.Earth.DKeccDtDistOrb": {"value": -1.511986e-14, "unit": 1 / u.sec}, + "log.final.Earth.DPincDtDistOrb": {"value": 1.636611e-14, "unit": 1 / u.sec}, + "log.final.Earth.DQincDtDistOrb": {"value": -2.927114e-14, "unit": 1 / u.sec}, + "log.final.Earth.CassiniOne": {"value": 0.420045}, + "log.final.Earth.CassiniTwo": {"value": 0.907503}, + "log.final.Earth.DOblDtDistRot": {"value": -5.056231e-14, "unit": u.rad}, + "log.final.Earth.DPrecADtDistRot": {"value": 7.830724e-12, "unit": u.rad / u.sec}, + "log.final.Earth.DXoblDtDistRot": {"value": -2.879142e-12, "unit": 1 / u.sec}, + "log.final.Earth.DYoblDtDistRot": {"value": 9.822112e-13, "unit": 1 / u.sec}, + "log.final.Earth.DZoblDtDistRot": {"value": 1.964009e-14, "unit": 1 / u.sec}, + "log.final.Earth.OblTimeDistRot": {"value": 6.213309e+13, "unit": u.rad}, + "log.final.Earth.PrecATimeDistRot": {"value": 8.023760e+11, "unit": u.sec}, + "log.final.Earth.XoblTimeDistRot": {"value": 3.473256e+11, "unit": 1 / u.sec}, + "log.final.Earth.YoblTimeDistRot": {"value": 1.018111e+12, "unit": 1 / u.sec}, + "log.final.Earth.ZoblTimeDistRot": {"value": 5.091626e+13, "unit": 1 / u.sec}, + "log.final.Earth.DynEllip": {"value": 0.003280}, + "log.final.Earth.PrecFNat": {"value": 2.471999e-12, "unit": u.rad / u.sec}, + "log.final.Mars.Mass": {"value": 6.413391e+23, "unit": u.kg}, + "log.final.Mars.Obliquity": {"value": 0.440202, "unit": u.rad}, + "log.final.Mars.PrecA": {"value": 0.208647, "unit": u.rad}, + "log.final.Mars.Xobl": {"value": 0.416880}, + "log.final.Mars.Yobl": {"value": 0.088265}, + "log.final.Mars.Zobl": {"value": 0.904666}, + "log.final.Mars.Radius": {"value": 3.393277e+06, "unit": u.m}, + "log.final.Mars.RadGyra": {"value": 0.500000}, + "log.final.Mars.RotAngMom": {"value": 1.308632e+32, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.final.Mars.RotKinEnergy": {"value": 4.638086e+27, "unit": u.Joule}, + "log.final.Mars.RotVel": {"value": 240.530700, "unit": u.m / u.sec}, + "log.final.Mars.BodyType": {"value": 0.000000}, + "log.final.Mars.RotRate": {"value": 7.088449e-05, "unit": 1 / u.sec}, + "log.final.Mars.RotPer": {"value": 8.863977e+04, "unit": u.sec}, + "log.final.Mars.Density": {"value": 3918.693171, "unit": u.kg / u.m ** 3}, + "log.final.Mars.HEcc": {"value": -0.004053}, + "log.final.Mars.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.final.Mars.HZLimRecVenus": {"value": -1.000000}, + "log.final.Mars.HZLimRunaway": {"value": -1.000000}, + "log.final.Mars.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.final.Mars.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.final.Mars.HZLimEarlyMars": {"value": -1.000000}, + "log.final.Mars.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.final.Mars.KEcc": {"value": 0.097597}, + "log.final.Mars.Eccentricity": {"value": 0.097681}, + "log.final.Mars.OrbEnergy": {"value": -1.867046e+32, "unit": u.Joule}, + "log.final.Mars.MeanMotion": {"value": 1.058606e-07, "unit": 1 / u.sec}, + "log.final.Mars.OrbPeriod": {"value": 5.935338e+07, "unit": u.sec}, + "log.final.Mars.SemiMajorAxis": {"value": 2.279367e+11, "unit": u.m}, + "log.final.Mars.COPP": {"value": -0.001727}, + "log.final.Mars.OrbAngMom": {"value": 3.510498e+39, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.final.Mars.ArgP": {"value": 5.667007, "unit": u.rad}, + "log.final.Mars.Inc": {"value": 0.024169, "unit": u.rad}, + "log.final.Mars.LongA": {"value": 0.574677, "unit": u.rad}, + "log.final.Mars.MeanLongitude": {"value": -0.041502, "unit": u.rad}, + "log.final.Mars.LongP": {"value": 6.241683, "unit": u.rad}, + "log.final.Mars.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.final.Mars.OrbPotEnergy": {"value": -3.734093e+32, "unit": u.Joule}, + "log.final.Mars.Sinc": {"value": 0.012084}, + "log.final.Mars.Pinc": {"value": 0.006569}, + "log.final.Mars.Qinc": {"value": 0.010143}, + "log.final.Mars.DIncDtDistOrb": {"value": -5.827632e-14, "unit": u.rad / u.sec}, + "log.final.Mars.DSincDtDistOrb": {"value": -2.913603e-14, "unit": 1 / u.sec}, + "log.final.Mars.DLongADtDistOrb": {"value": -2.143373e-12, "unit": u.rad / u.sec}, + "log.final.Mars.DLongPDtDistOrb": {"value": 2.349706e-12, "unit": u.rad / u.sec}, + "log.final.Mars.DHeccDtDistOrb": {"value": 2.283055e-13, "unit": 1 / u.sec}, + "log.final.Mars.DKeccDtDistOrb": {"value": 3.405735e-14, "unit": 1 / u.sec}, + "log.final.Mars.DPincDtDistOrb": {"value": -3.757798e-14, "unit": 1 / u.sec}, + "log.final.Mars.DQincDtDistOrb": {"value": -1.037687e-14, "unit": 1 / u.sec}, + "log.final.Mars.CassiniOne": {"value": 0.219290}, + "log.final.Mars.CassiniTwo": {"value": 0.975660}, + "log.final.Mars.DOblDtDistRot": {"value": -4.741784e-15, "unit": u.rad}, + "log.final.Mars.DPrecADtDistRot": {"value": 1.333064e-12, "unit": u.rad / u.sec}, + "log.final.Mars.DXoblDtDistRot": {"value": -1.218600e-13, "unit": 1 / u.sec}, + "log.final.Mars.DYoblDtDistRot": {"value": 5.548395e-13, "unit": 1 / u.sec}, + "log.final.Mars.DZoblDtDistRot": {"value": 2.020578e-15, "unit": 1 / u.sec}, + "log.final.Mars.OblTimeDistRot": {"value": 6.625339e+14, "unit": u.rad}, + "log.final.Mars.PrecATimeDistRot": {"value": 4.713341e+12, "unit": u.sec}, + "log.final.Mars.XoblTimeDistRot": {"value": 8.206136e+12, "unit": 1 / u.sec}, + "log.final.Mars.YoblTimeDistRot": {"value": 1.802323e+12, "unit": 1 / u.sec}, + "log.final.Mars.ZoblTimeDistRot": {"value": 4.949078e+14, "unit": 1 / u.sec}, + "log.final.Mars.DynEllip": {"value": 0.005363}, + "log.final.Mars.PrecFNat": {"value": 1.167214e-12, "unit": u.rad / u.sec}, + "log.final.Jupiter.Mass": {"value": 1.897089e+27, "unit": u.kg}, + "log.final.Jupiter.Radius": {"value": 9.567150e+06, "unit": u.m}, + "log.final.Jupiter.RadGyra": {"value": 0.500000}, + "log.final.Jupiter.RotKinEnergy": {"value": 1.147876e+32, "unit": u.Joule}, + "log.final.Jupiter.RotVel": {"value": 695.742781, "unit": u.m / u.sec}, + "log.final.Jupiter.BodyType": {"value": 0.000000}, + "log.final.Jupiter.Density": {"value": 5.171914e+05, "unit": u.kg / u.m ** 3}, + "log.final.Jupiter.HEcc": {"value": 0.023210}, + "log.final.Jupiter.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.final.Jupiter.HZLimRecVenus": {"value": -1.000000}, + "log.final.Jupiter.HZLimRunaway": {"value": -1.000000}, + "log.final.Jupiter.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.final.Jupiter.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.final.Jupiter.HZLimEarlyMars": {"value": -1.000000}, + "log.final.Jupiter.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.final.Jupiter.KEcc": {"value": 0.049745}, + "log.final.Jupiter.Eccentricity": {"value": 0.054893}, + "log.final.Jupiter.OrbEnergy": {"value": -1.615496e+35, "unit": u.Joule}, + "log.final.Jupiter.MeanMotion": {"value": 1.675590e-08, "unit": 1 / u.sec}, + "log.final.Jupiter.OrbPeriod": {"value": 3.749835e+08, "unit": u.sec}, + "log.final.Jupiter.SemiMajorAxis": {"value": 7.792260e+11, "unit": u.m}, + "log.final.Jupiter.OrbAngMom": {"value": 1.925364e+43, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.final.Jupiter.ArgP": {"value": 4.790758, "unit": u.rad}, + "log.final.Jupiter.Inc": {"value": 0.022601, "unit": u.rad}, + "log.final.Jupiter.LongA": {"value": 1.928985, "unit": u.rad}, + "log.final.Jupiter.MeanLongitude": {"value": 0.436557, "unit": u.rad}, + "log.final.Jupiter.LongP": {"value": 0.436557, "unit": u.rad}, + "log.final.Jupiter.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.final.Jupiter.OrbPotEnergy": {"value": -3.230992e+35, "unit": u.Joule}, + "log.final.Jupiter.Sinc": {"value": 0.011300}, + "log.final.Jupiter.Pinc": {"value": 0.010583}, + "log.final.Jupiter.Qinc": {"value": -0.003962}, + "log.final.Jupiter.DIncDtDistOrb": {"value": 9.025020e-15, "unit": u.rad / u.sec}, + "log.final.Jupiter.DSincDtDistOrb": {"value": 4.512222e-15, "unit": 1 / u.sec}, + "log.final.Jupiter.DLongADtDistOrb": {"value": 1.061829e-12, "unit": u.rad / u.sec}, + "log.final.Jupiter.DLongPDtDistOrb": {"value": 1.178614e-12, "unit": u.rad / u.sec}, + "log.final.Jupiter.DHeccDtDistOrb": {"value": 7.115393e-14, "unit": 1 / u.sec}, + "log.final.Jupiter.DKeccDtDistOrb": {"value": -5.144369e-16, "unit": 1 / u.sec}, + "log.final.Jupiter.DPincDtDistOrb": {"value": 1.920879e-17, "unit": 1 / u.sec}, + "log.final.Jupiter.DQincDtDistOrb": {"value": -1.281949e-14, "unit": 1 / u.sec}, + "log.final.Saturn.Mass": {"value": 5.680075e+26, "unit": u.kg}, + "log.final.Saturn.Radius": {"value": 6.378100e+06, "unit": u.m}, + "log.final.Saturn.RadGyra": {"value": 0.500000}, + "log.final.Saturn.RotKinEnergy": {"value": 1.527492e+31, "unit": u.Joule}, + "log.final.Saturn.RotVel": {"value": 463.828521, "unit": u.m / u.sec}, + "log.final.Saturn.BodyType": {"value": 0.000000}, + "log.final.Saturn.Density": {"value": 5.226266e+05, "unit": u.kg / u.m ** 3}, + "log.final.Saturn.HEcc": {"value": 0.035123}, + "log.final.Saturn.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.final.Saturn.HZLimRecVenus": {"value": -1.000000}, + "log.final.Saturn.HZLimRunaway": {"value": -1.000000}, + "log.final.Saturn.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.final.Saturn.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.final.Saturn.HZLimEarlyMars": {"value": -1.000000}, + "log.final.Saturn.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.final.Saturn.KEcc": {"value": -0.017651}, + "log.final.Saturn.Eccentricity": {"value": 0.039309}, + "log.final.Saturn.OrbEnergy": {"value": -2.640963e+34, "unit": u.Joule}, + "log.final.Saturn.MeanMotion": {"value": 6.757836e-09, "unit": 1 / u.sec}, + "log.final.Saturn.OrbPeriod": {"value": 9.297629e+08, "unit": u.sec}, + "log.final.Saturn.SemiMajorAxis": {"value": 1.427163e+12, "unit": u.m}, + "log.final.Saturn.OrbAngMom": {"value": 7.809960e+42, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.final.Saturn.ArgP": {"value": 0.283838, "unit": u.rad}, + "log.final.Saturn.Inc": {"value": 0.043583, "unit": u.rad}, + "log.final.Saturn.LongA": {"value": 1.752644, "unit": u.rad}, + "log.final.Saturn.MeanLongitude": {"value": 2.036482, "unit": u.rad}, + "log.final.Saturn.LongP": {"value": 2.036482, "unit": u.rad}, + "log.final.Saturn.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.final.Saturn.OrbPotEnergy": {"value": -5.281926e+34, "unit": u.Joule}, + "log.final.Saturn.Sinc": {"value": 0.021790}, + "log.final.Saturn.Pinc": {"value": 0.021431}, + "log.final.Saturn.Qinc": {"value": -0.003941}, + "log.final.Saturn.DIncDtDistOrb": {"value": -1.134630e-14, "unit": u.rad / u.sec}, + "log.final.Saturn.DSincDtDistOrb": {"value": -5.671805e-15, "unit": 1 / u.sec}, + "log.final.Saturn.DLongADtDistOrb": {"value": -1.473580e-12, "unit": u.rad / u.sec}, + "log.final.Saturn.DLongPDtDistOrb": {"value": 2.966731e-12, "unit": u.rad / u.sec}, + "log.final.Saturn.DHeccDtDistOrb": {"value": -1.425473e-13, "unit": 1 / u.sec}, + "log.final.Saturn.DKeccDtDistOrb": {"value": -5.888032e-14, "unit": 1 / u.sec}, + "log.final.Saturn.DPincDtDistOrb": {"value": 2.285947e-16, "unit": 1 / u.sec}, + "log.final.Saturn.DQincDtDistOrb": {"value": 3.260560e-14, "unit": 1 / u.sec}, + "log.final.George.Mass": {"value": 8.676288e+25, "unit": u.kg}, + "log.final.George.Radius": {"value": 9.567150e+06, "unit": u.m}, + "log.final.George.RadGyra": {"value": 0.500000}, + "log.final.George.RotKinEnergy": {"value": 5.249784e+30, "unit": u.Joule}, + "log.final.George.RotVel": {"value": 695.742781, "unit": u.m / u.sec}, + "log.final.George.BodyType": {"value": 0.000000}, + "log.final.George.Density": {"value": 2.365362e+04, "unit": u.kg / u.m ** 3}, + "log.final.George.HEcc": {"value": 0.004073}, + "log.final.George.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.final.George.HZLimRecVenus": {"value": -1.000000}, + "log.final.George.HZLimRunaway": {"value": -1.000000}, + "log.final.George.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.final.George.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.final.George.HZLimEarlyMars": {"value": -1.000000}, + "log.final.George.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.final.George.KEcc": {"value": -0.045828}, + "log.final.George.Eccentricity": {"value": 0.046009}, + "log.final.George.OrbEnergy": {"value": -2.005253e+33, "unit": u.Joule}, + "log.final.George.MeanMotion": {"value": 2.368075e-09, "unit": 1 / u.sec}, + "log.final.George.OrbPeriod": {"value": 2.653288e+09, "unit": u.sec}, + "log.final.George.SemiMajorAxis": {"value": 2.871088e+12, "unit": u.m}, + "log.final.George.OrbAngMom": {"value": 1.691778e+42, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.final.George.ArgP": {"value": 1.674515, "unit": u.rad}, + "log.final.George.Inc": {"value": 0.012161, "unit": u.rad}, + "log.final.George.LongA": {"value": 1.378432, "unit": u.rad}, + "log.final.George.MeanLongitude": {"value": 3.052947, "unit": u.rad}, + "log.final.George.LongP": {"value": 3.052947, "unit": u.rad}, + "log.final.George.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.final.George.OrbPotEnergy": {"value": -4.010505e+33, "unit": u.Joule}, + "log.final.George.Sinc": {"value": 0.006081}, + "log.final.George.Pinc": {"value": 0.005968}, + "log.final.George.Qinc": {"value": 0.001162}, + "log.final.George.DIncDtDistOrb": {"value": -6.761146e-15, "unit": u.rad / u.sec}, + "log.final.George.DSincDtDistOrb": {"value": -3.380510e-15, "unit": 1 / u.sec}, + "log.final.George.DLongADtDistOrb": {"value": 6.267136e-13, "unit": u.rad / u.sec}, + "log.final.George.DLongPDtDistOrb": {"value": 4.235494e-13, "unit": u.rad / u.sec}, + "log.final.George.DHeccDtDistOrb": {"value": -1.993962e-14, "unit": 1 / u.sec}, + "log.final.George.DKeccDtDistOrb": {"value": 4.227316e-15, "unit": 1 / u.sec}, + "log.final.George.DPincDtDistOrb": {"value": -2.589623e-15, "unit": 1 / u.sec}, + "log.final.George.DQincDtDistOrb": {"value": -4.386728e-15, "unit": 1 / u.sec}, + "log.final.Neptune.Mass": {"value": 1.023485e+26, "unit": u.kg}, + "log.final.Neptune.Radius": {"value": 9.567150e+06, "unit": u.m}, + "log.final.Neptune.RadGyra": {"value": 0.500000}, + "log.final.Neptune.RotKinEnergy": {"value": 6.192825e+30, "unit": u.Joule}, + "log.final.Neptune.RotVel": {"value": 695.742781, "unit": u.m / u.sec}, + "log.final.Neptune.BodyType": {"value": 0.000000}, + "log.final.Neptune.Density": {"value": 2.790262e+04, "unit": u.kg / u.m ** 3}, + "log.final.Neptune.HEcc": {"value": 0.006340}, + "log.final.Neptune.HZLimitDryRunaway": {"value": -1.000000, "unit": u.m}, + "log.final.Neptune.HZLimRecVenus": {"value": -1.000000}, + "log.final.Neptune.HZLimRunaway": {"value": -1.000000}, + "log.final.Neptune.HZLimMoistGreenhouse": {"value": -1.000000}, + "log.final.Neptune.HZLimMaxGreenhouse": {"value": -1.000000}, + "log.final.Neptune.HZLimEarlyMars": {"value": -1.000000}, + "log.final.Neptune.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.final.Neptune.KEcc": {"value": 0.006139}, + "log.final.Neptune.Eccentricity": {"value": 0.008825}, + "log.final.Neptune.OrbEnergy": {"value": -1.509722e+33, "unit": u.Joule}, + "log.final.Neptune.MeanMotion": {"value": 1.207446e-09, "unit": 1 / u.sec}, + "log.final.Neptune.OrbPeriod": {"value": 5.203698e+09, "unit": u.sec}, + "log.final.Neptune.SemiMajorAxis": {"value": 4.498484e+12, "unit": u.m}, + "log.final.Neptune.OrbAngMom": {"value": 2.500588e+42, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.final.Neptune.ArgP": {"value": 4.791147, "unit": u.rad}, + "log.final.Neptune.Inc": {"value": 0.031082, "unit": u.rad}, + "log.final.Neptune.LongA": {"value": 2.293563, "unit": u.rad}, + "log.final.Neptune.MeanLongitude": {"value": 0.801525, "unit": u.rad}, + "log.final.Neptune.LongP": {"value": 0.801525, "unit": u.rad}, + "log.final.Neptune.TotOrbEnergy": {"value": -1.976696e+35, "unit": u.Joule}, + "log.final.Neptune.OrbPotEnergy": {"value": -3.019443e+33, "unit": u.Joule}, + "log.final.Neptune.Sinc": {"value": 0.015540}, + "log.final.Neptune.Pinc": {"value": 0.011655}, + "log.final.Neptune.Qinc": {"value": -0.010279}, + "log.final.Neptune.DIncDtDistOrb": {"value": 1.344875e-15, "unit": u.rad / u.sec}, + "log.final.Neptune.DSincDtDistOrb": {"value": 6.723561e-16, "unit": 1 / u.sec}, + "log.final.Neptune.DLongADtDistOrb": {"value": -3.600781e-14, "unit": u.rad / u.sec}, + "log.final.Neptune.DLongPDtDistOrb": {"value": 1.576014e-13, "unit": u.rad / u.sec}, + "log.final.Neptune.DHeccDtDistOrb": {"value": 2.059479e-15, "unit": 1 / u.sec}, + "log.final.Neptune.DKeccDtDistOrb": {"value": 5.822243e-17, "unit": 1 / u.sec}, + "log.final.Neptune.DPincDtDistOrb": {"value": 8.743871e-16, "unit": 1 / u.sec}, + "log.final.Neptune.DQincDtDistOrb": {"value": -2.507269e-17, "unit": 1 / u.sec}, + } +) +class Test_RD4(Benchmark): + pass diff --git a/tests/Distorb/RD4/venus.in b/tests/Distorb/RD4/venus.in new file mode 100644 index 000000000..847e10b32 --- /dev/null +++ b/tests/Distorb/RD4/venus.in @@ -0,0 +1,17 @@ +sName Venus +saModules distorb distrot +dMass 2.4464824121e-06 +dRadius -0.9499 +dRotPeriod -243.69 +dObliquity 177.36 +dSemi 0.72333377 +dEcc 0.00677478 +dInc 3.39471 +dLongP 131.54907 +dLongA 76.68069 +#bGRCorr 0 +saOutputOrder Time Ecce Inc ArgP LongA Obliq PrecA +dDynEllip 1.3125e-5 #based on J2 from Yoder 1995 +dPrecA 0.0 +#bOverrideMaxEcc 1 +#dHaltMaxEcc 0.4 diff --git a/tests/Distorb/RD4/vpl.in b/tests/Distorb/RD4/vpl.in new file mode 100644 index 000000000..d732ee11d --- /dev/null +++ b/tests/Distorb/RD4/vpl.in @@ -0,0 +1,19 @@ +sSystemName solarsys +iVerbose 5 +iDigits 6 +bOverwrite 1 +sUnitMass solar +sUnitLength au +sUnitTime y +sUnitAngle d +bDoLog 1 +saBodyFiles sun.in mercury.in venus.in earth.in mars.in jupiter.in $ + saturn.in george.in neptune.in +#saBodyFiles sun.in venus.in earth.in mars.in jupiter.in saturn.in george.in neptune.in +#saBodyFiles sun.in mercury.in jupiter.in saturn.in + +bDoForward 1 +bVarDt 1 +dEta 0.01 +dStopTime 5000 +dOutputTime 1000