diff --git a/examples/TidalEarth/README.rst b/examples/TidalEarth/README.rst index 494057eda..560af8945 100644 --- a/examples/TidalEarth/README.rst +++ b/examples/TidalEarth/README.rst @@ -15,8 +15,8 @@ from its host star. =================== ============ This example shows the thermal, magnetic, and orbital evolution of an -Earth-like planet using **ThermInt**, **RadHeat**, and **EqTide**. Here we use the "orbit-only" tidal model, which is highly idealized. -The interior model is 1-D and many free parameters have been +Earth-like planet using **ThermInt**, **RadHeat**, and **EqTide**. Here we use the "orbit-only" (DB15) tidal model, +which is highly idealized. The interior model is 1-D and many free parameters have been tuned to match `the modern Earth <../EarthInterior>`_. This example reproduces the results of `Driscoll & Barnes (2015) `_. The generated figures should be compared to their Figures 4 and 5. diff --git a/src/body.c b/src/body.c index 3f3adadf7..101e72ad3 100644 --- a/src/body.c +++ b/src/body.c @@ -771,15 +771,11 @@ double fdHflowSecMan(BODY *body, EVOLVE *evolve, int iBody) { double dHflowSecMan = 0; if (body[iBody].bThermint) { - dHflowSecMan += fdPowerThermint(body, iBody); - // dHflowSecMan += fdHflowSecManThermint(body,iBody); + dHflowSecMan += fdPowerThermint(body, iBody); } if (body[iBody].bEqtide) { - // dHflowSecMan -= fdTidePower(body,iBody,evolve->iEqtideModel); // formerly - // dTidalPowerMan - dHflowSecMan -= body[iBody].dTidalPowMan; - } - // Should add RadHeat here + dHflowSecMan -= fdTidePower(body,iBody,evolve->iEqtideModel); + } // XXX Should add RadHeat here return dHflowSecMan; } diff --git a/src/module.c b/src/module.c index 5eddcdc5c..70d3f5ffa 100644 --- a/src/module.c +++ b/src/module.c @@ -981,8 +981,6 @@ void VerifyModuleMultiRadheatThermint(BODY *body, UPDATE *update, OPTIONS *options, int iBody, int *iModuleProps, int *iModuleForce) { - /* This will need modification if material can move between layers */ - if (body[iBody].bThermint) { if (!body[iBody].bRadheat) { if (control->Io.iVerbose > VERBINPUT) { @@ -1007,7 +1005,11 @@ void VerifyModuleMultiEqtideThermint(BODY *body, UPDATE *update, int iBody, int *iModuleProps, int *iModuleForce) { - // Initialize mantle tidal power + if (control->Evolve.iEqtideModel == CTL) { + fprintf(stderr,"ERROR: The CTL EqTide model cannot be coupled to ThermInt.\n"); + exit(EXIT_INPUT); + } + body[iBody].dTidalPowMan = 0; if (body[iBody].bEqtide && body[iBody].bThermint) { diff --git a/tests/EqtideRadheatThermint/CPLRadheatThermint/c.in b/tests/EqtideRadheatThermint/CPLRadheatThermint/c.in new file mode 100644 index 000000000..a511a590f --- /dev/null +++ b/tests/EqtideRadheatThermint/CPLRadheatThermint/c.in @@ -0,0 +1,66 @@ +# Planet c's parameters +sName c # Body's name +saModules radheat thermint eqtide + +# Physical Properties +dMass -2.8857142857142857 +dRadius -1.367 # Radius, negative -> Earth radii +# Orbital Properties +dEcc 0.08571428571428572 +dSemi -0.03984 # Semi-major axis, negative -> AU + +# mantle parameters that we will change +#dActViscMan 500000 + + +# EQTIDE Parameters +bForceEqSpin 1 +dObliquity 0 +dRadGyra 0.5 +dTidalQ 10 # +dK2 0.299 # Love number of degree 2 (Yoder 1995) +saTidePerts sun # Names of perturbing bodies (only central body may have >1) +bFixOrbit 0 +#bTideLock 1 + +# RADHEAT Parameters +# *Num* are in numbers of atoms, negative -> Earth vals +### 40K +d40KPowerMan -1 +d40KPowerCore -1 +d40KPowerCrust -1 +### 232Th +d232ThPowerMan -1 +d232ThPowerCore -1 +d232ThPowerCrust -1 +### 235U +d235UPowerMan -1 +d235UPowerCore -1 +d235UPowerCrust -1 +### 238U +d238UPowerMan -1 +d238UPowerCore -1 +d238UPowerCrust -1 + +### THERMINT inputs. +dTMan 3500 +dTCore 5500 +dEruptEff 0.2 +dViscJumpMan 2.0 +dTrefLind 5600 +dImK2ManOrbModel 2 +dShModRef 1e5 +dStiffness 1.71e13 +dActViscMan 309285.71428571426 + +saOutputOrder -Time -TMan -TUMan -TLMan -TCMB -TCore $ + -HflowUMan -HflowMeltMan -RadPowerMan -RadPowerCore -RadPowerCrust $ + -HflowCMB -HflowSecMan -HFlowSurf $ + -TDotMan -TDotCore -TJumpLMan -TJumpUMan -RIC -RayleighMan -ViscUMan -ViscLMan $ + -MeltMassFluxMan FMeltUMan $ + -MagMom -CoreBuoyTherm -CoreBuoyCompo -CoreBuoyTotal -MagPauseRad $ + -BLUMan -BLLMan $ + K2 ImK2 SemiMajorAxis Eccentricity ShmodUMan MeanMotion ViscUManArr $ + ChiOC ChiIC MassChiOC MassChiIC MassOC MassIC $ + -RadPowerTotal SurfEnFluxEqtide -PowerEqtide TideLock TidalQ -SurfEnFluxTotal + diff --git a/tests/EqtideRadheatThermint/CPLRadheatThermint/sun.in b/tests/EqtideRadheatThermint/CPLRadheatThermint/sun.in new file mode 100644 index 000000000..532fbf19d --- /dev/null +++ b/tests/EqtideRadheatThermint/CPLRadheatThermint/sun.in @@ -0,0 +1,11 @@ +# host star parameters +sName sun +dMass 0.118 +dRadius 0.00072354 +dLuminosity 5.50944e30 +sStellarModel none +saModules stellar eqtide +saTidePerts c + +dTidalQ 1e6 +dK2 1.5 diff --git a/tests/EqtideRadheatThermint/CPLRadheatThermint/test_CPLRadheatThermint.py b/tests/EqtideRadheatThermint/CPLRadheatThermint/test_CPLRadheatThermint.py new file mode 100644 index 000000000..d5db1d185 --- /dev/null +++ b/tests/EqtideRadheatThermint/CPLRadheatThermint/test_CPLRadheatThermint.py @@ -0,0 +1,478 @@ +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": 5.522265e+40, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.system.TotEnergy": {"value": -2.036611e+40, "unit": u.Joule}, + "log.initial.system.PotEnergy": {"value": -2.036790e+40, "unit": u.Joule}, + "log.initial.system.KinEnergy": {"value": 1.817220e+36, "unit": u.Joule}, + "log.initial.system.DeltaTime": {"value": 0.000000, "unit": u.sec}, + "log.initial.sun.Mass": {"value": 2.346331e+29, "unit": u.kg}, + "log.initial.sun.Obliquity": {"value": 0.000000, "unit": u.rad}, + "log.initial.sun.PrecA": {"value": 0.000000, "unit": u.rad}, + "log.initial.sun.Xobl": {"value": 0.000000}, + "log.initial.sun.Yobl": {"value": 0.000000}, + "log.initial.sun.Zobl": {"value": 1.000000}, + "log.initial.sun.Radius": {"value": 1.082400e+08, "unit": u.m}, + "log.initial.sun.RadGyra": {"value": 0.500000}, + "log.initial.sun.RotAngMom": {"value": 4.997713e+40, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.sun.RotKinEnergy": {"value": 1.817220e+36, "unit": u.Joule}, + "log.initial.sun.RotVel": {"value": 7871.438080, "unit": u.m / u.sec}, + "log.initial.sun.BodyType": {"value": 0.000000}, + "log.initial.sun.RotRate": {"value": 7.272205e-05, "unit": 1 / u.sec}, + "log.initial.sun.RotPer": {"value": 8.640000e+04, "unit": u.sec}, + "log.initial.sun.Density": {"value": 4.417095e+04, "unit": u.kg / u.m ** 3}, + "log.initial.sun.SurfEnFluxTotal": {"value": 6.760931e-06, "unit": u.kg / u.sec ** 3}, + "log.initial.sun.TidalQ": {"value": 1.000000e+06}, + "log.initial.sun.ImK2": {"value": -1.500000e-06}, + "log.initial.sun.K2": {"value": 1.500000}, + "log.initial.sun.K2Man": {"value": 0.010000}, + "log.initial.sun.Imk2Man": {"value": 0.000000}, + "log.initial.sun.TidalQMantle": {"value": 100.000000}, + "log.initial.sun.HEcc": {"value": 0.000000}, + "log.initial.sun.HZLimitDryRunaway": {"value": 1.933828e+28, "unit": u.m}, + "log.initial.sun.HZLimRecVenus": {"value": 1.593583e+28, "unit": u.m}, + "log.initial.sun.HZLimRunaway": {"value": 2.080916e+28, "unit": u.m}, + "log.initial.sun.HZLimMoistGreenhouse": {"value": 2.108557e+28, "unit": u.m}, + "log.initial.sun.HZLimMaxGreenhouse": {"value": 3.574093e+28, "unit": u.m}, + "log.initial.sun.HZLimEarlyMars": {"value": 3.899625e+28, "unit": u.m}, + "log.initial.sun.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3}, + "log.initial.sun.KEcc": {"value": 0.000000}, + "log.initial.sun.Eccentricity": {"value": -1.000000}, + "log.initial.sun.OrbEnergy": {"value": 0.000000, "unit": u.Joule}, + "log.initial.sun.MeanMotion": {"value": -1.000000, "unit": 1 / u.sec}, + "log.initial.sun.OrbPeriod": {"value": -1.000000, "unit": u.sec}, + "log.initial.sun.SemiMajorAxis": {"value": -1.000000, "unit": u.m}, + "log.initial.sun.CriticalSemiMajorAxis": {"value": -1.000000, "unit": u.m}, + "log.initial.sun.COPP": {"value": 0.000000}, + "log.initial.sun.OrbAngMom": {"value": 0.000000, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.sun.LongP": {"value": 0.000000, "unit": u.rad}, + "log.initial.sun.LXUVTot": {"value": 7.801060e+57, "unit": u.kg / u.sec ** 3}, + "log.initial.sun.TotOrbEnergy": {"value": -2.264152e+34, "unit": u.Joule}, + "log.initial.sun.OrbPotEnergy": {"value": -1.000000, "unit": u.Joule}, + "log.initial.sun.LostEnergy": {"value": 5.562685e-309, "unit": u.Joule}, + "log.initial.sun.LostAngMom": {"value": 5.562685e-309, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.sun.LockTime": {"value": -1.000000, "unit": u.sec}, + "log.initial.sun.BodyDsemiDtEqtide": {"value": -1.000000}, + "log.initial.sun.BodyDeccDt": {"value": -1.000000}, + "log.initial.sun.DOblDtEqtide": {"value": 0.000000, "unit": u.rad / u.sec}, + "log.initial.sun.DRotPerDtEqtide": {"value": 2.696923e-20}, + "log.initial.sun.DRotRateDtEqtide": {"value": -2.269974e-29, "unit": 1 / u.sec ** 2}, + "log.initial.sun.EqRotRateDiscrete": {"value": 8.600935e-06, "unit": 1 / u.sec}, + "log.initial.sun.EqRotPerDiscrete": {"value": 7.305235e+05, "unit": u.sec}, + "log.initial.sun.EqRotRateCont": {"value": 9.201245e-06, "unit": 1 / u.sec}, + "log.initial.sun.EqRotPerCont": {"value": 6.828625e+05, "unit": u.sec}, + "log.initial.sun.EqRotPer": {"value": 7.305235e+05, "unit": u.sec}, + "log.initial.sun.EqTidePower": {"value": 0.000000, "unit": 1 / u.sec}, + "log.initial.sun.GammaRot": {"value": -1.000000, "unit": u.sec}, + "log.initial.sun.GammaOrb": {"value": -1.000000, "unit": u.sec}, + "log.initial.sun.OceanK2": {"value": 0.010000}, + "log.initial.sun.EnvTidalQ": {"value": -1.000000}, + "log.initial.sun.OceanTidalQ": {"value": -1.000000}, + "log.initial.sun.TideLock": {"value": 0.000000}, + "log.initial.sun.RotTimeEqtide": {"value": 3.203651e+24, "unit": u.sec}, + "log.initial.sun.EnvK2": {"value": 0.010000}, + "log.initial.sun.OblTimeEqtide": {"value": -1.000000}, + "log.initial.sun.PowerEqtide": {"value": 9.953877e+11, "unit": u.W}, + "log.initial.sun.SurfEnFluxEqtide": {"value": 6.760931e-06, "unit": u.kg / u.sec ** 3}, + "log.initial.sun.Luminosity": {"value": 7.801060e+60, "unit": u.W}, + "log.initial.sun.LXUVStellar": {"value": 7.801060e+57, "unit": u.W}, + "log.initial.sun.Temperature": {"value": 5778.000000, "unit": u.K}, + "log.initial.sun.LXUVFrac": {"value": 0.001000}, + "log.initial.sun.RossbyNumber": {"value": 0.078260}, + "log.initial.sun.DRotPerDtStellar": {"value": 2.889133e-12}, + "log.initial.c.Mass": {"value": 1.723402e+25, "unit": u.kg}, + "log.initial.c.Obliquity": {"value": 0.000000, "unit": u.rad}, + "log.initial.c.PrecA": {"value": 0.000000, "unit": u.rad}, + "log.initial.c.Xobl": {"value": 0.000000}, + "log.initial.c.Yobl": {"value": 0.000000}, + "log.initial.c.Zobl": {"value": 1.000000}, + "log.initial.c.Radius": {"value": 8.718863e+06, "unit": u.m}, + "log.initial.c.RadGyra": {"value": 0.500000}, + "log.initial.c.RotAngMom": {"value": 2.817034e+33, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.c.RotKinEnergy": {"value": 1.211456e+28, "unit": u.Joule}, + "log.initial.c.RotVel": {"value": 74.990373, "unit": u.m / u.sec}, + "log.initial.c.BodyType": {"value": 0.000000}, + "log.initial.c.RotRate": {"value": 8.600935e-06, "unit": 1 / u.sec}, + "log.initial.c.RotPer": {"value": 7.305235e+05, "unit": u.sec}, + "log.initial.c.Density": {"value": 6207.525742, "unit": u.kg / u.m ** 3}, + "log.initial.c.SurfEnFluxTotal": {"value": 286.139740, "unit": u.W / u.m ** 2}, + "log.initial.c.TidalQ": {"value": 27.579829}, + "log.initial.c.ImK2": {"value": -2.877378e-12}, + "log.initial.c.K2": {"value": 1.500000}, + "log.initial.c.K2Man": {"value": 1.500000}, + "log.initial.c.Imk2Man": {"value": -2.877378e-12}, + "log.initial.c.TidalQMantle": {"value": 27.579829}, + "log.initial.c.ViscUMan": {"value": 1.223055e+07, "unit": u.m ** 2 / u.sec}, + "log.initial.c.HEcc": {"value": 0.000000}, + "log.initial.c.HZLimitDryRunaway": {"value": 1.940971e+28, "unit": u.m}, + "log.initial.c.HZLimRecVenus": {"value": 1.593583e+28, "unit": u.m}, + "log.initial.c.HZLimRunaway": {"value": 2.080916e+28, "unit": u.m}, + "log.initial.c.HZLimMoistGreenhouse": {"value": 2.108557e+28, "unit": u.m}, + "log.initial.c.HZLimMaxGreenhouse": {"value": 3.574093e+28, "unit": u.m}, + "log.initial.c.HZLimEarlyMars": {"value": 3.899625e+28, "unit": u.m}, + "log.initial.c.Instellation": {"value": 1.754105e+40, "unit": u.kg / u.sec ** 3}, + "log.initial.c.KEcc": {"value": 0.085714}, + "log.initial.c.Eccentricity": {"value": 0.085714}, + "log.initial.c.OrbEnergy": {"value": -2.264152e+34, "unit": u.Joule}, + "log.initial.c.MeanMotion": {"value": 8.600935e-06, "unit": 1 / u.sec}, + "log.initial.c.OrbPeriod": {"value": 7.305235e+05, "unit": u.sec}, + "log.initial.c.SemiMajorAxis": {"value": 5.959979e+09, "unit": u.m}, + "log.initial.c.CriticalSemiMajorAxis": {"value": -1.000000, "unit": u.m}, + "log.initial.c.COPP": {"value": 0.000000}, + "log.initial.c.OrbAngMom": {"value": 5.245521e+39, "unit": (u.kg * u.m ** 2) / u.sec}, + "log.initial.c.LongP": {"value": 0.000000, "unit": u.rad}, + "log.initial.c.TotOrbEnergy": {"value": -2.264152e+34, "unit": u.Joule}, + "log.initial.c.OrbPotEnergy": {"value": -4.528304e+34, "unit": u.Joule}, + "log.initial.c.LostEnergy": {"value": 5.562685e-309, "unit": u.Joule}, + "log.initial.c.TidalRadius": {"value": 8.718863e+06, "unit": u.m}, + "log.initial.c.DsemiDtEqtide": {"value": 3.660829e-14, "unit": u.m / u.sec}, + "log.initial.c.DeccDtEqtide": {"value": 1.141343e-24, "unit": 1 / u.sec}, + "log.initial.c.DMeanMotionDtEqtide": {"value": -7.924496e-29, "unit": 1 / u.sec ** 2}, + "log.initial.c.DOrbPerDtEqtide": {"value": 6.730699e-18}, + "log.initial.c.EccTimeEqtide": {"value": 7.509951e+22, "unit": u.sec}, + "log.initial.c.SemiTimeEqtide": {"value": 1.628041e+23, "unit": u.sec}, + "log.initial.c.DHEccDtEqtide": {"value": 0.000000, "unit": 1 / u.sec}, + "log.initial.c.DKEccDtEqtide": {"value": 1.141343e-24, "unit": 1 / u.sec}, + "log.initial.c.DXoblDtEqtide": {"value": 0.000000, "unit": 1 / u.sec}, + "log.initial.c.DYoblDtEqtide": {"value": 0.000000, "unit": 1 / u.sec}, + "log.initial.c.DZoblDtEqtide": {"value": 0.000000, "unit": 1 / u.sec}, + "log.initial.c.LockTime": {"value": 0.000000, "unit": u.sec}, + "log.initial.c.BodyDsemiDtEqtide": {"value": -1.000000}, + "log.initial.c.BodyDeccDt": {"value": -1.000000}, + "log.initial.c.DOblDtEqtide": {"value": 0.000000, "unit": u.rad / u.sec}, + "log.initial.c.DRotPerDtEqtide": {"value": -4.724686e-298}, + "log.initial.c.DRotRateDtEqtide": {"value": 5.562685e-309, "unit": 1 / u.sec ** 2}, + "log.initial.c.EqRotRateDiscrete": {"value": 8.600935e-06, "unit": 1 / u.sec}, + "log.initial.c.EqRotPerDiscrete": {"value": 7.305235e+05, "unit": u.sec}, + "log.initial.c.EqRotRateCont": {"value": 9.201245e-06, "unit": 1 / u.sec}, + "log.initial.c.EqRotPerCont": {"value": 6.828625e+05, "unit": u.sec}, + "log.initial.c.EqRotPer": {"value": 7.305235e+05, "unit": u.sec}, + "log.initial.c.EqTidePower": {"value": 5.377595e+06, "unit": 1 / u.sec}, + "log.initial.c.GammaRot": {"value": -1.000000, "unit": u.sec}, + "log.initial.c.GammaOrb": {"value": -1.000000, "unit": u.sec}, + "log.initial.c.OceanK2": {"value": 0.010000}, + "log.initial.c.EnvTidalQ": {"value": -1.000000}, + "log.initial.c.OceanTidalQ": {"value": -1.000000}, + "log.initial.c.TideLock": {"value": 1.000000}, + "log.initial.c.RotTimeEqtide": {"value": 1.546184e+303, "unit": u.sec}, + "log.initial.c.EnvK2": {"value": 0.010000}, + "log.initial.c.OblTimeEqtide": {"value": -1.000000}, + "log.initial.c.PowerEqtide": {"value": 7.885913e-06, "unit": u.TW}, + "log.initial.c.SurfEnFluxEqtide": {"value": 8.255102e-09, "unit": u.kg / u.sec ** 3}, + "log.initial.c.D26AlPowerDt": {"value": -1.000000}, + "log.initial.c.D26AlNumManDt": {"value": 0.000000, "unit": 1 / u.sec}, + "log.initial.c.D40KPowerDt": {"value": -1.000000}, + "log.initial.c.D40KNumManDt": {"value": -1.694596e+26, "unit": 1 / u.sec}, + "log.initial.c.D232ThNumManDt": {"value": -9.534863e+23, "unit": 1 / u.sec}, + "log.initial.c.D238UNumManDt": {"value": -1.408812e+24, "unit": 1 / u.sec}, + "log.initial.c.D235UNumManDt": {"value": -3.089017e+24, "unit": 1 / u.sec}, + "log.initial.c.RadPowerMan": {"value": 74.591650, "unit": u.TW}, + "log.initial.c.RadPowerCore": {"value": 34.624677, "unit": u.TW}, + "log.initial.c.RadPowerCrust": {"value": 33.267780, "unit": u.TW}, + "log.initial.c.RadPowerTotal": {"value": 142.484107, "unit": u.TW}, + "log.initial.c.SurfEnFluxRadTotal": {"value": 0.149155, "unit": u.kg / u.sec ** 3}, + "log.initial.c.TMan": {"value": 3500.000000, "unit": u.K}, + "log.initial.c.TUMan": {"value": 2450.000000, "unit": u.K}, + "log.initial.c.TLMan": {"value": 4550.000000, "unit": u.K}, + "log.initial.c.TCore": {"value": 5500.000000, "unit": u.K}, + "log.initial.c.TCMB": {"value": 4400.000000, "unit": u.K}, + "log.initial.c.BLUMan": {"value": 0.020759, "unit": u.km}, + "log.initial.c.BLLMan": {"value": 0.603556, "unit": u.km}, + "log.initial.c.TJumpUMan": {"value": 2150.000000, "unit": u.K}, + "log.initial.c.TJumpLMan": {"value": 150.000000, "unit": u.K}, + "log.initial.c.SignTJumpUMan": {"value": 1.000000}, + "log.initial.c.SignTJumpLMan": {"value": -1.000000}, + "log.initial.c.ViscLMan": {"value": 2.097051e+10, "unit": u.m ** 2 / u.sec}, + "log.initial.c.ShmodUMan": {"value": 95.354299}, + "log.initial.c.FMeltUMan": {"value": 1.000000}, + "log.initial.c.FMeltLMan": {"value": 0.807856}, + "log.initial.c.MeltfactorUMan": {"value": 1.925879e+07}, + "log.initial.c.MeltfactorLMan": {"value": 2.246447e+04}, + "log.initial.c.DepthMeltMan": {"value": 2.551886e+06}, + "log.initial.c.TDepthMeltMan": {"value": 3725.932721}, + "log.initial.c.TJumpMeltMan": {"value": 2149.989620}, + "log.initial.c.MeltMassFluxMan": {"value": 9.462947e+10, "unit": u.kg / u.sec}, + "log.initial.c.ViscUManArr": {"value": 2.355457e+14, "unit": u.m ** 2 / u.sec}, + "log.initial.c.RayleighMan": {"value": 1.904937e+17, "unit": u.nd}, + "log.initial.c.ViscMMan": {"value": 1.223055e+08, "unit": u.m ** 2 / u.sec}, + "log.initial.c.TDotMan": {"value": -4.190393e-11, "unit": u.K / u.sec}, + "log.initial.c.TDotCore": {"value": 2.532624e-13, "unit": u.K / u.sec}, + "log.initial.c.HfluxUMan": {"value": 434.983048}, + "log.initial.c.HflowUMan": {"value": 2.218694e+05, "unit": u.TW}, + "log.initial.c.HfluxLMan": {"value": -2.485271}, + "log.initial.c.HflowLMan": {"value": -3.784353e+14}, + "log.initial.c.HfluxCMB": {"value": -2.485271}, + "log.initial.c.HflowCMB": {"value": -378.435255, "unit": u.TW}, + "log.initial.c.HflowLatentMan": {"value": 6.108891e+16}, + "log.initial.c.HflowMeltMan": {"value": 5.147345e+04, "unit": u.TW}, + "log.initial.c.HflowLatentIC": {"value": 0.000000}, + "log.initial.c.PowerGravIC": {"value": 0.000000}, + "log.initial.c.HflowSurf": {"value": 2.733429e+05, "unit": u.TW}, + "log.initial.c.HflowSecMan": {"value": 2.125578e+05, "unit": u.TW}, + "log.initial.c.HfluxCMBAd": {"value": 0.041915}, + "log.initial.c.HfluxCMBConv": {"value": -2.527186}, + "log.initial.c.RIC": {"value": 0.000000, "unit": u.km}, + "log.initial.c.DRICDTCMB": {"value": 0.000000}, + "log.initial.c.RICDot": {"value": 0.000000}, + "log.initial.c.ChiOC": {"value": 0.094946}, + "log.initial.c.ChiIC": {"value": 0.000000}, + "log.initial.c.ThermConductOC": {"value": 110.000000}, + "log.initial.c.MassOC": {"value": 1.941615e+24}, + "log.initial.c.MassIC": {"value": 0.000000}, + "log.initial.c.MassChiOC": {"value": 1.843482e+23}, + "log.initial.c.MassChiIC": {"value": 0.000000}, + "log.initial.c.DTChi": {"value": 0.000000}, + "log.initial.c.CoreBuoyTherm": {"value": -2.902026e-11, "unit": u.m ** 2 / u.sec ** 3}, + "log.initial.c.CoreBuoyCompo": {"value": 0.000000, "unit": u.m ** 2 / u.sec ** 3}, + "log.initial.c.CoreBuoyTotal": {"value": -2.902026e-11, "unit": u.m ** 2 / u.sec ** 3}, + "log.initial.c.GravICB": {"value": 0.000000}, + "log.initial.c.MagMom": {"value": 0.000000}, + "log.initial.c.PresSWind": {"value": 2.676100e-09}, + "log.initial.c.MagPauseRad": {"value": 0.000000}, + "log.initial.c.ViscJumpMan": {"value": 2.000000}, + "log.initial.c.EruptEff": {"value": 0.200000}, + "log.initial.c.TrefLind": {"value": 5600.000000}, + "log.initial.c.DynViscUMan": {"value": 3.057639e+08, "unit": u.Joule}, + "log.initial.c.TsolUMan": {"value": 1449.476709}, + "log.final.system.Age": {"value": 3.155760e+13, "unit": u.sec, "rtol": 1e-4}, + "log.final.system.Time": {"value": 3.155760e+13, "unit": u.sec, "rtol": 1e-4}, + "log.final.system.TotAngMom": {"value": 5.522265e+40, "unit": (u.kg * u.m ** 2) / u.sec, "rtol": 1e-4}, + "log.final.system.TotEnergy": {"value": -2.036611e+40, "unit": u.Joule, "rtol": 1e-4}, + "log.final.system.PotEnergy": {"value": -2.036790e+40, "unit": u.Joule, "rtol": 1e-4}, + "log.final.system.KinEnergy": {"value": 1.813388e+36, "unit": u.Joule, "rtol": 1e-4}, + "log.final.sun.Mass": {"value": 2.346331e+29, "unit": u.kg, "rtol": 1e-4}, + "log.final.sun.Obliquity": {"value": 0.000000, "unit": u.rad, "rtol": 1e-4}, + "log.final.sun.PrecA": {"value": 0.000000, "unit": u.rad, "rtol": 1e-4}, + "log.final.sun.Xobl": {"value": 0.000000, "rtol": 1e-4}, + "log.final.sun.Yobl": {"value": 0.000000, "rtol": 1e-4}, + "log.final.sun.Zobl": {"value": 1.000000, "rtol": 1e-4}, + "log.final.sun.Radius": {"value": 1.082400e+08, "unit": u.m, "rtol": 1e-4}, + "log.final.sun.RadGyra": {"value": 0.500000, "rtol": 1e-4}, + "log.final.sun.RotAngMom": {"value": 4.992442e+40, "unit": (u.kg * u.m ** 2) / u.sec, "rtol": 1e-4}, + "log.final.sun.RotKinEnergy": {"value": 1.813388e+36, "unit": u.Joule, "rtol": 1e-4}, + "log.final.sun.RotVel": {"value": 7863.136081, "unit": u.m / u.sec, "rtol": 1e-4}, + "log.final.sun.BodyType": {"value": 0.000000, "rtol": 1e-4}, + "log.final.sun.RotRate": {"value": 7.264535e-05, "unit": 1 / u.sec, "rtol": 1e-4}, + "log.final.sun.RotPer": {"value": 8.649122e+04, "unit": u.sec, "rtol": 1e-4}, + "log.final.sun.Density": {"value": 4.417095e+04, "unit": u.kg / u.m ** 3, "rtol": 1e-4}, + "log.final.sun.SurfEnFluxTotal": {"value": 6.752804e-06, "unit": u.kg / u.sec ** 3, "rtol": 1e-4}, + "log.final.sun.TidalQ": {"value": 1.000000e+06, "rtol": 1e-4}, + "log.final.sun.ImK2": {"value": -1.500000e-06, "rtol": 1e-4}, + "log.final.sun.K2": {"value": 1.500000, "rtol": 1e-4}, + "log.final.sun.K2Man": {"value": 0.010000, "rtol": 1e-4}, + "log.final.sun.Imk2Man": {"value": 0.000000, "rtol": 1e-4}, + "log.final.sun.TidalQMantle": {"value": 100.000000, "rtol": 1e-4}, + "log.final.sun.HEcc": {"value": 0.000000, "rtol": 1e-4}, + "log.final.sun.HZLimitDryRunaway": {"value": 1.933828e+28, "unit": u.m, "rtol": 1e-4}, + "log.final.sun.HZLimRecVenus": {"value": 1.593583e+28, "unit": u.m, "rtol": 1e-4}, + "log.final.sun.HZLimRunaway": {"value": 2.080916e+28, "unit": u.m, "rtol": 1e-4}, + "log.final.sun.HZLimMoistGreenhouse": {"value": 2.108557e+28, "unit": u.m, "rtol": 1e-4}, + "log.final.sun.HZLimMaxGreenhouse": {"value": 3.574093e+28, "unit": u.m, "rtol": 1e-4}, + "log.final.sun.HZLimEarlyMars": {"value": 3.899625e+28, "unit": u.m, "rtol": 1e-4}, + "log.final.sun.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3, "rtol": 1e-4}, + "log.final.sun.KEcc": {"value": 0.000000, "rtol": 1e-4}, + "log.final.sun.Eccentricity": {"value": -1.000000, "rtol": 1e-4}, + "log.final.sun.OrbEnergy": {"value": 0.000000, "unit": u.Joule, "rtol": 1e-4}, + "log.final.sun.MeanMotion": {"value": -1.000000, "unit": 1 / u.sec, "rtol": 1e-4}, + "log.final.sun.OrbPeriod": {"value": -1.000000, "unit": u.sec, "rtol": 1e-4}, + "log.final.sun.SemiMajorAxis": {"value": -1.000000, "unit": u.m, "rtol": 1e-4}, + "log.final.sun.CriticalSemiMajorAxis": {"value": -1.000000, "unit": u.m, "rtol": 1e-4}, + "log.final.sun.COPP": {"value": 0.000000, "rtol": 1e-4}, + "log.final.sun.OrbAngMom": {"value": 0.000000, "unit": (u.kg * u.m ** 2) / u.sec, "rtol": 1e-4}, + "log.final.sun.LongP": {"value": 0.000000, "unit": u.rad, "rtol": 1e-4}, + "log.final.sun.LXUVTot": {"value": 7.801060e+57, "unit": u.kg / u.sec ** 3, "rtol": 1e-4}, + "log.final.sun.TotOrbEnergy": {"value": -2.264152e+34, "unit": u.Joule, "rtol": 1e-4}, + "log.final.sun.OrbPotEnergy": {"value": -1.000000, "unit": u.Joule, "rtol": 1e-4}, + "log.final.sun.LostEnergy": {"value": 3.831219e+33, "unit": u.Joule, "rtol": 1e-4}, + "log.final.sun.LostAngMom": {"value": 5.271084e+37, "unit": (u.kg * u.m ** 2) / u.sec, "rtol": 1e-4}, + "log.final.sun.LockTime": {"value": -1.000000, "unit": u.sec, "rtol": 1e-4}, + "log.final.sun.BodyDsemiDtEqtide": {"value": -1.000000, "rtol": 1e-4}, + "log.final.sun.BodyDeccDt": {"value": -1.000000, "rtol": 1e-4}, + "log.final.sun.DOblDtEqtide": {"value": 0.000000, "unit": u.rad / u.sec, "rtol": 1e-4}, + "log.final.sun.DRotPerDtEqtide": {"value": 2.702621e-20, "rtol": 1e-4}, + "log.final.sun.DRotRateDtEqtide": {"value": -2.269974e-29, "unit": 1 / u.sec ** 2, "rtol": 1e-4}, + "log.final.sun.EqRotRateDiscrete": {"value": 8.600935e-06, "unit": 1 / u.sec, "rtol": 1e-4}, + "log.final.sun.EqRotPerDiscrete": {"value": 7.305235e+05, "unit": u.sec, "rtol": 1e-4}, + "log.final.sun.EqRotRateCont": {"value": 9.201245e-06, "unit": 1 / u.sec, "rtol": 1e-4}, + "log.final.sun.EqRotPerCont": {"value": 6.828625e+05, "unit": u.sec, "rtol": 1e-4}, + "log.final.sun.EqRotPer": {"value": 7.305235e+05, "unit": u.sec, "rtol": 1e-4}, + "log.final.sun.EqTidePower": {"value": 0.000000, "unit": 1 / u.sec, "rtol": 1e-4}, + "log.final.sun.GammaRot": {"value": -1.000000, "unit": u.sec, "rtol": 1e-4}, + "log.final.sun.GammaOrb": {"value": -1.000000, "unit": u.sec, "rtol": 1e-4}, + "log.final.sun.OceanK2": {"value": 0.010000, "rtol": 1e-4}, + "log.final.sun.EnvTidalQ": {"value": -1.000000, "rtol": 1e-4}, + "log.final.sun.OceanTidalQ": {"value": -1.000000, "rtol": 1e-4}, + "log.final.sun.TideLock": {"value": 0.000000, "rtol": 1e-4}, + "log.final.sun.RotTimeEqtide": {"value": 3.200273e+24, "unit": u.sec, "rtol": 1e-4}, + "log.final.sun.EnvK2": {"value": 0.010000, "rtol": 1e-4}, + "log.final.sun.OblTimeEqtide": {"value": -1.000000, "rtol": 1e-4}, + "log.final.sun.PowerEqtide": {"value": 9.941912e+11, "unit": u.W, "rtol": 1e-4}, + "log.final.sun.SurfEnFluxEqtide": {"value": 6.752804e-06, "unit": u.kg / u.sec ** 3, "rtol": 1e-4}, + "log.final.sun.Luminosity": {"value": 7.801060e+60, "unit": u.W, "rtol": 1e-4}, + "log.final.sun.LXUVStellar": {"value": 7.801060e+57, "unit": u.W, "rtol": 1e-4}, + "log.final.sun.Temperature": {"value": 5778.000000, "unit": u.K, "rtol": 1e-4}, + "log.final.sun.LXUVFrac": {"value": 0.001000, "rtol": 1e-4}, + "log.final.sun.RossbyNumber": {"value": 0.078342, "rtol": 1e-4}, + "log.final.sun.DRotPerDtStellar": {"value": 2.892183e-12, "rtol": 1e-4}, + "log.final.c.Mass": {"value": 1.723402e+25, "unit": u.kg, "rtol": 1e-4}, + "log.final.c.Obliquity": {"value": 0.000000, "unit": u.rad, "rtol": 1e-4}, + "log.final.c.PrecA": {"value": 0.000000, "unit": u.rad, "rtol": 1e-4}, + "log.final.c.Xobl": {"value": 0.000000, "rtol": 1e-4}, + "log.final.c.Yobl": {"value": 0.000000, "rtol": 1e-4}, + "log.final.c.Zobl": {"value": 1.000000, "rtol": 1e-4}, + "log.final.c.Radius": {"value": 8.718863e+06, "unit": u.m, "rtol": 1e-4}, + "log.final.c.RadGyra": {"value": 0.500000, "rtol": 1e-4}, + "log.final.c.RotAngMom": {"value": 2.817034e+33, "unit": (u.kg * u.m ** 2) / u.sec, "rtol": 1e-4}, + "log.final.c.RotKinEnergy": {"value": 1.211456e+28, "unit": u.Joule, "rtol": 1e-4}, + "log.final.c.RotVel": {"value": 74.990373, "unit": u.m / u.sec, "rtol": 1e-4}, + "log.final.c.BodyType": {"value": 0.000000, "rtol": 1e-4}, + "log.final.c.RotRate": {"value": 8.600935e-06, "unit": 1 / u.sec, "rtol": 1e-4}, + "log.final.c.RotPer": {"value": 7.305235e+05, "unit": u.sec, "rtol": 1e-4}, + "log.final.c.Density": {"value": 6207.525742, "unit": u.kg / u.m ** 3, "rtol": 1e-4}, + "log.final.c.SurfEnFluxTotal": {"value": 66.911059, "unit": u.W / u.m ** 2, "rtol": 1e-4}, + "log.final.c.TidalQ": {"value": 92.319939, "rtol": 1e-4}, + "log.final.c.ImK2": {"value": -7.853352e-12, "rtol": 1e-4}, + "log.final.c.K2": {"value": 1.500000, "rtol": 1e-4}, + "log.final.c.K2Man": {"value": 1.500000, "rtol": 1e-4}, + "log.final.c.Imk2Man": {"value": -7.853352e-12, "rtol": 1e-4}, + "log.final.c.TidalQMantle": {"value": 92.319939, "rtol": 1e-4}, + "log.final.c.ViscUMan": {"value": 3.735877e+08, "unit": u.m ** 2 / u.sec, "rtol": 1e-4}, + "log.final.c.HEcc": {"value": 0.000000, "rtol": 1e-4}, + "log.final.c.HZLimitDryRunaway": {"value": 1.940971e+28, "unit": u.m, "rtol": 1e-4}, + "log.final.c.HZLimRecVenus": {"value": 1.593583e+28, "unit": u.m, "rtol": 1e-4}, + "log.final.c.HZLimRunaway": {"value": 2.080916e+28, "unit": u.m, "rtol": 1e-4}, + "log.final.c.HZLimMoistGreenhouse": {"value": 2.108557e+28, "unit": u.m, "rtol": 1e-4}, + "log.final.c.HZLimMaxGreenhouse": {"value": 3.574093e+28, "unit": u.m, "rtol": 1e-4}, + "log.final.c.HZLimEarlyMars": {"value": 3.899625e+28, "unit": u.m, "rtol": 1e-4}, + "log.final.c.Instellation": {"value": 1.754105e+40, "unit": u.kg / u.sec ** 3, "rtol": 1e-4}, + "log.final.c.KEcc": {"value": 0.085714, "rtol": 1e-4}, + "log.final.c.Eccentricity": {"value": 0.085714, "rtol": 1e-4}, + "log.final.c.OrbEnergy": {"value": -2.264152e+34, "unit": u.Joule, "rtol": 1e-4}, + "log.final.c.MeanMotion": {"value": 8.600935e-06, "unit": 1 / u.sec, "rtol": 1e-4}, + "log.final.c.OrbPeriod": {"value": 7.305235e+05, "unit": u.sec, "rtol": 1e-4}, + "log.final.c.SemiMajorAxis": {"value": 5.959979e+09, "unit": u.m, "rtol": 1e-4}, + "log.final.c.CriticalSemiMajorAxis": {"value": -1.000000, "unit": u.m, "rtol": 1e-4}, + "log.final.c.COPP": {"value": 0.000000, "rtol": 1e-4}, + "log.final.c.OrbAngMom": {"value": 5.245521e+39, "unit": (u.kg * u.m ** 2) / u.sec, "rtol": 1e-4}, + "log.final.c.LongP": {"value": 0.000000, "unit": u.rad, "rtol": 1e-4}, + "log.final.c.TotOrbEnergy": {"value": -2.264152e+34, "unit": u.Joule, "rtol": 1e-4}, + "log.final.c.OrbPotEnergy": {"value": -4.528304e+34, "unit": u.Joule, "rtol": 1e-4}, + "log.final.c.LostEnergy": {"value": 4.656483e+20, "unit": u.Joule, "rtol": 1e-4}, + "log.final.c.TidalRadius": {"value": 8.718863e+06, "unit": u.m, "rtol": 1e-4}, + "log.final.c.DsemiDtEqtide": {"value": 3.660470e-14, "unit": u.m / u.sec, "rtol": 1e-4}, + "log.final.c.DeccDtEqtide": {"value": 1.137829e-24, "unit": 1 / u.sec, "rtol": 1e-4}, + "log.final.c.DMeanMotionDtEqtide": {"value": -7.923719e-29, "unit": 1 / u.sec ** 2, "rtol": 1e-4}, + "log.final.c.DOrbPerDtEqtide": {"value": 6.730039e-18, "rtol": 1e-4}, + "log.final.c.EccTimeEqtide": {"value": 7.533142e+22, "unit": u.sec, "rtol": 1e-4}, + "log.final.c.SemiTimeEqtide": {"value": 1.628200e+23, "unit": u.sec, "rtol": 1e-4}, + "log.final.c.DHEccDtEqtide": {"value": 0.000000, "unit": 1 / u.sec, "rtol": 1e-4}, + "log.final.c.DKEccDtEqtide": {"value": 1.137829e-24, "unit": 1 / u.sec, "rtol": 1e-4}, + "log.final.c.DXoblDtEqtide": {"value": 0.000000, "unit": 1 / u.sec, "rtol": 1e-4}, + "log.final.c.DYoblDtEqtide": {"value": 0.000000, "unit": 1 / u.sec, "rtol": 1e-4}, + "log.final.c.DZoblDtEqtide": {"value": 0.000000, "unit": 1 / u.sec, "rtol": 1e-4}, + "log.final.c.LockTime": {"value": 0.000000, "unit": u.sec, "rtol": 1e-4}, + "log.final.c.BodyDsemiDtEqtide": {"value": -1.000000, "rtol": 1e-4}, + "log.final.c.BodyDeccDt": {"value": -1.000000, "rtol": 1e-4}, + "log.final.c.DOblDtEqtide": {"value": 0.000000, "unit": u.rad / u.sec, "rtol": 1e-4}, + "log.final.c.DRotPerDtEqtide": {"value": -4.724686e-298, "rtol": 1e-4}, + "log.final.c.DRotRateDtEqtide": {"value": 5.562685e-309, "unit": 1 / u.sec ** 2, "rtol": 1e-4}, + "log.final.c.EqRotRateDiscrete": {"value": 8.600935e-06, "unit": 1 / u.sec, "rtol": 1e-4}, + "log.final.c.EqRotPerDiscrete": {"value": 7.305235e+05, "unit": u.sec, "rtol": 1e-4}, + "log.final.c.EqRotRateCont": {"value": 9.201245e-06, "unit": 1 / u.sec, "rtol": 1e-4}, + "log.final.c.EqRotPerCont": {"value": 6.828625e+05, "unit": u.sec, "rtol": 1e-4}, + "log.final.c.EqRotPer": {"value": 7.305235e+05, "unit": u.sec, "rtol": 1e-4}, + "log.final.c.EqTidePower": {"value": 1.467730e+07, "unit": 1 / u.sec, "rtol": 1e-4}, + "log.final.c.GammaRot": {"value": -1.000000, "unit": u.sec, "rtol": 1e-4}, + "log.final.c.GammaOrb": {"value": -1.000000, "unit": u.sec, "rtol": 1e-4}, + "log.final.c.OceanK2": {"value": 0.010000, "rtol": 1e-4}, + "log.final.c.EnvTidalQ": {"value": -1.000000, "rtol": 1e-4}, + "log.final.c.OceanTidalQ": {"value": -1.000000, "rtol": 1e-4}, + "log.final.c.TideLock": {"value": 1.000000, "rtol": 1e-4}, + "log.final.c.RotTimeEqtide": {"value": 1.546184e+303, "unit": u.sec, "rtol": 1e-4}, + "log.final.c.EnvK2": {"value": 0.010000, "rtol": 1e-4}, + "log.final.c.OblTimeEqtide": {"value": -1.000000, "rtol": 1e-4}, + "log.final.c.PowerEqtide": {"value": 2.152336e-05, "unit": u.TW, "rtol": 1e-4}, + "log.final.c.SurfEnFluxEqtide": {"value": 2.253101e-08, "unit": u.kg / u.sec ** 3, "rtol": 1e-4}, + "log.final.c.D26AlPowerDt": {"value": -1.000000, "rtol": 1e-4}, + "log.final.c.D26AlNumManDt": {"value": 0.000000, "unit": 1 / u.sec, "rtol": 1e-4}, + "log.final.c.D40KPowerDt": {"value": -1.000000, "rtol": 1e-4}, + "log.final.c.D40KNumManDt": {"value": -1.693664e+26, "unit": 1 / u.sec, "rtol": 1e-4}, + "log.final.c.D232ThNumManDt": {"value": -9.534391e+23, "unit": 1 / u.sec, "rtol": 1e-4}, + "log.final.c.D238UNumManDt": {"value": -1.408594e+24, "unit": 1 / u.sec, "rtol": 1e-4}, + "log.final.c.D235UNumManDt": {"value": -3.085976e+24, "unit": 1 / u.sec, "rtol": 1e-4}, + "log.final.c.RadPowerMan": {"value": 74.549698, "unit": u.TW, "rtol": 1e-4}, + "log.final.c.RadPowerCore": {"value": 34.605537, "unit": u.TW, "rtol": 1e-4}, + "log.final.c.RadPowerCrust": {"value": 33.248925, "unit": u.TW, "rtol": 1e-4}, + "log.final.c.RadPowerTotal": {"value": 142.404189, "unit": u.TW, "rtol": 1e-4}, + "log.final.c.SurfEnFluxRadTotal": {"value": 0.149071, "unit": u.kg / u.sec ** 3, "rtol": 1e-4}, + "log.final.c.TMan": {"value": 2856.678300, "unit": u.K, "rtol": 1e-4}, + "log.final.c.TUMan": {"value": 1999.674810, "unit": u.K, "rtol": 1e-4}, + "log.final.c.TLMan": {"value": 3713.681791, "unit": u.K, "rtol": 1e-4}, + "log.final.c.TCore": {"value": 5500.346073, "unit": u.K, "rtol": 1e-4}, + "log.final.c.TCMB": {"value": 4400.276858, "unit": u.K, "rtol": 1e-4}, + "log.final.c.BLUMan": {"value": 0.070181, "unit": u.km, "rtol": 1e-4}, + "log.final.c.BLLMan": {"value": 32.062434, "unit": u.km, "rtol": 1e-4}, + "log.final.c.TJumpUMan": {"value": 1699.674810, "unit": u.K, "rtol": 1e-4}, + "log.final.c.TJumpLMan": {"value": 686.595068, "unit": u.K, "rtol": 1e-4}, + "log.final.c.SignTJumpUMan": {"value": 1.000000, "rtol": 1e-4}, + "log.final.c.SignTJumpLMan": {"value": 1.000000, "rtol": 1e-4}, + "log.final.c.ViscLMan": {"value": 1.438969e+16, "unit": u.m ** 2 / u.sec, "rtol": 1e-4}, + "log.final.c.ShmodUMan": {"value": 870.127293, "rtol": 1e-4}, + "log.final.c.FMeltUMan": {"value": 1.000000, "rtol": 1e-4}, + "log.final.c.FMeltLMan": {"value": 0.000000, "rtol": 1e-4}, + "log.final.c.MeltfactorUMan": {"value": 1.925879e+07, "rtol": 1e-4}, + "log.final.c.MeltfactorLMan": {"value": 1.000000, "rtol": 1e-4}, + "log.final.c.DepthMeltMan": {"value": 1.028496e+06, "rtol": 1e-4}, + "log.final.c.TDepthMeltMan": {"value": 2513.887961, "rtol": 1e-4}, + "log.final.c.TJumpMeltMan": {"value": 1699.639720, "rtol": 1e-4}, + "log.final.c.MeltMassFluxMan": {"value": 2.799111e+10, "unit": u.kg / u.sec, "rtol": 1e-4}, + "log.final.c.ViscUManArr": {"value": 7.194847e+15, "unit": u.m ** 2 / u.sec, "rtol": 1e-4}, + "log.final.c.RayleighMan": {"value": 6.470321e+15, "unit": u.nd, "rtol": 1e-4}, + "log.final.c.ViscMMan": {"value": 3.735877e+09, "unit": u.m ** 2 / u.sec, "rtol": 1e-4}, + "log.final.c.TDotMan": {"value": -1.028371e-11, "unit": u.K / u.sec, "rtol": 1e-4}, + "log.final.c.TDotCore": {"value": 1.224863e-15, "unit": u.K / u.sec, "rtol": 1e-4}, + "log.final.c.HfluxUMan": {"value": 101.716955, "rtol": 1e-4}, + "log.final.c.HflowUMan": {"value": 5.188221e+04, "unit": u.TW, "rtol": 1e-4}, + "log.final.c.HfluxLMan": {"value": 0.214143, "rtol": 1e-4}, + "log.final.c.HflowLMan": {"value": 3.260785e+13, "rtol": 1e-4}, + "log.final.c.HfluxCMB": {"value": 0.214143, "rtol": 1e-4}, + "log.final.c.HflowCMB": {"value": 32.607845, "unit": u.TW, "rtol": 1e-4}, + "log.final.c.HflowLatentMan": {"value": 1.164734e+16, "rtol": 1e-4}, + "log.final.c.HflowMeltMan": {"value": 1.203642e+04, "unit": u.TW, "rtol": 1e-4}, + "log.final.c.HflowLatentIC": {"value": 0.000000, "rtol": 1e-4}, + "log.final.c.PowerGravIC": {"value": 0.000000, "rtol": 1e-4}, + "log.final.c.HflowSurf": {"value": 6.391863e+04, "unit": u.TW, "rtol": 1e-4}, + "log.final.c.HflowSecMan": {"value": 5.216413e+04, "unit": u.TW, "rtol": 1e-4}, + "log.final.c.HfluxCMBAd": {"value": 0.041920, "rtol": 1e-4}, + "log.final.c.HfluxCMBConv": {"value": 0.172223, "rtol": 1e-4}, + "log.final.c.RIC": {"value": 0.000000, "unit": u.km, "rtol": 1e-4}, + "log.final.c.DRICDTCMB": {"value": 0.000000, "rtol": 1e-4}, + "log.final.c.RICDot": {"value": 0.000000, "rtol": 1e-4}, + "log.final.c.ChiOC": {"value": 0.094946, "rtol": 1e-4}, + "log.final.c.ChiIC": {"value": 0.000000, "rtol": 1e-4}, + "log.final.c.ThermConductOC": {"value": 110.006921, "rtol": 1e-4}, + "log.final.c.MassOC": {"value": 1.941615e+24, "rtol": 1e-4}, + "log.final.c.MassIC": {"value": 0.000000, "rtol": 1e-4}, + "log.final.c.MassChiOC": {"value": 1.843482e+23, "rtol": 1e-4}, + "log.final.c.MassChiIC": {"value": 0.000000, "rtol": 1e-4}, + "log.final.c.DTChi": {"value": 0.000000, "rtol": 1e-4}, + "log.final.c.CoreBuoyTherm": {"value": 1.977674e-12, "unit": u.m ** 2 / u.sec ** 3, "rtol": 1e-4}, + "log.final.c.CoreBuoyCompo": {"value": 0.000000, "unit": u.m ** 2 / u.sec ** 3, "rtol": 1e-4}, + "log.final.c.CoreBuoyTotal": {"value": 1.977674e-12, "unit": u.m ** 2 / u.sec ** 3, "rtol": 1e-4}, + "log.final.c.GravICB": {"value": 0.000000, "rtol": 1e-4}, + "log.final.c.MagMom": {"value": 1.291847, "rtol": 1e-4}, + "log.final.c.PresSWind": {"value": 2.676100e-09, "rtol": 1e-4}, + "log.final.c.MagPauseRad": {"value": 1.089167, "rtol": 1e-4}, + "log.final.c.ViscJumpMan": {"value": 2.000000, "rtol": 1e-4}, + "log.final.c.EruptEff": {"value": 0.200000, "rtol": 1e-4}, + "log.final.c.TrefLind": {"value": 5600.000000, "rtol": 1e-4}, + "log.final.c.DynViscUMan": {"value": 9.339693e+09, "unit": u.Joule, "rtol": 1e-4}, + "log.final.c.TsolUMan": {"value": 1449.547670, "rtol": 1e-4}, + } +) +class Test_CPLRadheatThermint(Benchmark): + pass diff --git a/tests/EqtideRadheatThermint/CPLRadheatThermint/vpl.in b/tests/EqtideRadheatThermint/CPLRadheatThermint/vpl.in new file mode 100644 index 000000000..35e6d8b69 --- /dev/null +++ b/tests/EqtideRadheatThermint/CPLRadheatThermint/vpl.in @@ -0,0 +1,40 @@ +# Example primary input file for VPLanet +sSystemName lp890-9 # System Name +iVerbose 5 # Verbosity level +bOverwrite 1 # Allow file overwrites? + +# All space after a # is ignored, as is white space +# The first lowercase letter(s) denote the cast: b=boolean, i=int, d=double, +# s=string. An "a" indicates an array and multiple arguments are allowed/expected. + +# List of "body files" that contain body-specific parameters +saBodyFiles sun.in $ # The host star + c.in # planet + +# Array options can continue to the next line with a terminating "$". The $ can be +# at the end of the string or not. Comments are allowed afterwards. + +# Input/Output Units +sUnitMass solar # Options: gram, kg, Earth, Neptune, Jupiter, solar +sUnitLength aU # Options: cm, m, km, Earth, Jupiter, solar, AU +sUnitTime YEARS # Options: sec, day, year, Myr, Gyr +sUnitAngle d # Options: deg, rad +sUnitTemp K + +# Units specified in the primary input file are propagated into the bodies. Otherwise +# specify units on a per body basis in the body files. +# Most string arguments can be in any case and need only be unambiguous. + +# Input/Output +bDoLog 1 # Write a log file? +iDigits 6 # Maximum number of digits to right of decimal +dMinValue 1e-10 # Minimum value of eccentricity/obliquity + +# Option names must be exact in spelling and case. + +# Evolution Parameters +bDoForward 1 # Perform a forward evolution? +bVarDt 1 # Use variable timestepping? +dEta 0.001 # Coefficient for variable timestepping +dStopTime 1e6 # Stop time for evolution +dOutputTime 1e5 # Output timesteps (assuming in body files)