Options and Outputs

Module Lookup

Module name

Description

atmesc

Atmospheric escape of terrestrial or gaseous planets.

binary

Orbital evolution of circumbinary planets.

distorb

Orbital evolution of multi-planet systems via 2nd or 4th order secular theory.

distrot

Evolution of a planet’s rotational axis.

eqtide

Tidal evolution of a two-body systems.

flare

XUV emission of stellar flares.

galhabit

Evolution of the galactic habitat of a system.

magmoc

Thermal and geochemical evolution of a magma ocean

poise

Evolution of planetary climate using an energy balance model.

radheat

Radiogenic heating of planetary interiors.

spinbody

Orbital evolution using an N-body integrator.

stellar

Evolution of stars before and during the main sequence.

thermint

Post-magma ocean thermal evolution of terrestrial planet interiors.

Options Lookup

The contents of this page can be generated by typing vplanet -H in a terminal.

VPLanet is a general purpose planetary evolution integrator. From the command line, enter one optional command line option and one file name, e.g.:

> vplanet vpl.in

where vpl.in is the “primary” input file consisting of options and a list of files that contain the initial conditions for all bodies in a system. The code then simulates planetary systen evolution forward or backward in time. This onboard help provides information regarding the input and output files, command line options, and each option and output parameter. For more more information, see https://virtualplanetarylaboratory.github.io/vplanet, or consult the examples directory.

Support for VPLanet has been provided through grants by NASA, the NSF.

Command Line Options

Description

-v, -verbose

– Maximum verbosity, i.e. display all warnings and updates.

-q, -quiet

– No verbosity, i.e. nothing printed to device.

-h, -help

– Display short help.

-H, -Help

– Display extended help.

Input File Structure

  • Options must be the first string on the line and must be written exactly as shown below; the options are case-sensitive.

  • All characters to the right of a # sign are treated as a comment and are ignored. Blank lines are also ignored.

  • Options that take an array of arguments may span multiple lines if a $ is placed at the end of the line.

  • Output parameters (the arguments to saOutputOrder) are not case sensitive, and need only enough characters to be unambiguous.

Output File Structure

Output files consist columns of data separated by white space in the same order as listed in saOutputOrder.

Input Options

bAccuracyMode

Description

Re-invert matrix every EBM time step?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bAlbedoZA

Description

Use albedo based on zenith angle

Type

Bool

Modules

ALL

Files

Any

Default value

0

bAtmEscAuto

Description

Let AtmEsc determine the proper atmospheric escape regime for
an eroding H envelope. If the planetary radius exceeds the
Roche lobe, the unbound material is assumed to be Bondi-limited
and mass loss proceeds following Equation 4 from Owen & Wu
(2016). For planetary radius less than the Roche limit, if the
incident XUV flux exceeds the critical flux (Equation A25 from
Luger et. al. 2015), the mass loss is
radiation/recombination-limited and scales as XUV^0.5.
Otherwise, the loss is energy-limited and is linearly
proportional to the incident XUV flux.

Type

Bool

Modules

ALL

Files

Any

Default value

0

bCalcAB

Description

Calculate A and B in OLR function, from (T & pCO2)

Type

Bool

Modules

ALL

Files

Any

Default value

0

bColdStart

Description

Start from snowball Earth conditions

Type

Bool

Modules

ALL

Files

Any

Default value

0

bDiffRot

Description

Adjust heat diffusion for rotation rate?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bDiscreteRot

Description

In the CPL model of EqTide, this option toggles between two
ways in which the rotation rate behaves once it has damped to
equilibrium. The rigorous model only has two states: 1:1 and
3:2 frequency ratios, but this is due to truncation in the
solution; the equilibrium value is a continuous function of
eccentricity. If this option is set to 1, then only the 1:1 and
3:2 states are availeble, with a boundary at e = sqrt(1/19); if
set to 0, the rotation rate is (1+9.5e^2) times the mean
motion.

Type

Bool

Modules

ALL

Files

Any

Default value

1

bDoBackward

Description

Do Backward Integration?

Type

Bool

Modules

ALL

Files

Any

Default value

No

bDoForward

Description

Do Forward Integration?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bDoLog

Description

Write Log File?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bElevFB

Description

Use elevation feedback for ice sheet ablation?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bEnvTides

Description

Include effects of gaseous envelope tides?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bEvolveRG

Description

Set this flag to 0 to ignore the role of mass concentration in
stellar evolution. Only useful for testing purposes.

Type

Bool

Modules

STELLAR

Files

Any

Default value

1

bFixOrbit

Description

In EqTide, setting this value to 1 holds the eccentricity and
semi-major axis fixed during the evolution. This option can be
helpful for testing or for faking planet-planet perturbations
that maintain an eccentricity. (Although a bFixEcc option would
probably be better!)

Type

Bool

Modules

ALL

Files

Any

Default value

0

bForceEcc

Description

Force Eccentricity to evolve sinusoidally?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bForceEqSpin

Description

Set this option to 1 to force the rotation rate to the
equilibrium tidal value in EqTide. If set to 0, the rotation
rate may “bounce” back and forth across the equilibrium value
due to numerical issues. Setting to 0 may result in a
simulation that requires orders of magnitude more computational
time. The threshold to force the equilibrium value is set by
dMaxLockDiff, but the rotation rate will only be fixed if the
derivative on the other side of the equilibrium value send the
spin rate toward equilibrium. This latter functionality enables
pre-MS stars to not tidally lock if the torque due to
contraction exceeds the tidal torque.

Type

Bool

Modules

ALL

Files

Any

Default value

0

bForceObliq

Description

Force obliquity to evolve sinusoidally?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bForcePrecRate

Description

In DisRot, set the axial precession rate to a fixed value. This
option can mimic the forcing of a natural satellite, or be used
for testing.

Type

Bool

Modules

ALL

Files

Any

Default value

0

bGalacTides

Description

Include galactic tides

Type

Bool

Modules

ALL

Files

Any

Default value

1

bGRCorr

Description

Use general relativity correction

Type

Bool

Modules

DistOrb DistRot SpiNBody

Files

Body Only

Default value

0

bHadley

Description

Enable Hadley circulation

Type

Bool

Modules

ALL

Files

Any

Default value

0

bHaltAllPlanetsDesicc

Description

Halt when all planets desiccated?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bHaltAllPlanetsSolid

Description

Halt when all planets solidified?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bHaltAtmDesiSurfCool

Description

Halt when atmosphere desiccated and surface below 1000K?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bHaltCloseEnc

Description

Halt the code if the apocenter of an interior planets is less
than 4 mutual Hill radii from the pericenter of an outer
planet.

Type

Bool

Modules

ALL

Files

Any

Default value

0

bHaltDblSync

Description

Halt at Double Synchronous State?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bHaltEndBaraffeGrid

Description

The BARRAFFE stellar model will only compute parameters until
the end of the main sequence. Setting this flag to 1 will halt
the code if the end of the model grid is reached.

Type

Bool

Modules

ALL

Files

Any

Default value

1

bHaltEnterHabZone

Description

Halt when planet enters habitable zone?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bHaltEnvelopeGone

Description

Halt When Envelope Evaporates?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bHaltHillStab

Description

If two planets come within the Hill stability criterion, the
code willhalt if this parameter is set

Type

Bool

Modules

ALL

Files

Any

Default value

0

bHaltHolmanUnstable

Description

Halt when CBP is Holman-Wiegert Unstable?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bHaltMantleMeltFracLow

Description

Halt when mantle mostly solidified?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bHaltMantleSolidified

Description

Halt when mantle solidified?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bHaltMerge

Description

Halt at Merge

Type

Bool

Modules

DistOrb EqTide SpiNBody STELLAR

Files

Any

Default value

If eqtide or distorb called 1, else 0

bHaltPosDeDt

Description

Halt if de/dt > 0?

Type

Bool

Modules

DistOrb EqTide SpiNBody

Files

Any

Default value

0

bHaltRocheLobe

Description

Halt when Roche lobe crossing occurs?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bHaltSurfaceDesiccated

Description

Halt at Desiccation?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bHaltSyncRot

Description

Halt if the rotation becomes syncrhonous?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bHaltTideLock

Description

Halt if Tide-Locked?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bHostBinary

Description

Include 3 body interactions with binary star

Type

Bool

Modules

ALL

Files

Any

Default value

0

bIceSheets

Description

Include ice sheets?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bInstantO2Sink

Description

If set to 1, then all oxygen released by photolysis is
immediately removed from the atmosphere. This mimics rapid
surface oxidation.

Type

Bool

Modules

ALL

Files

Any

Default value

0

bInvPlane

Description

Convert input coordinates to invariable plane coordinates

Type

Bool

Modules

ALL

Files

Any

Default value

0

bMantleTides

Description

Include effects of mantle tides?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bMEPDiff

Description

Calculate diffusion from max entropy production (D=B/4)?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bOceanTides

Description

Include tidal dissapation due to oceans?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bOptManQuasiSol

Description

Solidify when melt frac = 0.4?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bOutputEigen

Description

Output Eigenvalues?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bOutputEnc

Description

Output stellar encounter information?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bOutputLapl

Description

Write special files that contain the eigenvalues and
eigenvectors of the system with DistOrb. In the LL2 solution,
these are already computed. In the RD4 solution, they are
computed at time of output. These can be useful for
interpretating results.

Type

Bool

Modules

ALL

Files

Any

Default value

0

bOverrideMaxEcc

Description

Override default maximum eccentricity (MAXORBDISTORB) in
DistOrb?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bOverwrite

Description

Permit file overwrite?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bRadialMigr

Description

Use radial migration?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bReadOrbitData

Description

Rather than calculate orbital evolution with DistOrb or
SpiNBody, users may read in a previously run simulation. See
sFileOrbitData for more information.

Type

Bool

Modules

ALL

Files

Any

Default value

0

bReadOrbitOblData

Description

Read in orbital and obliquity data and use with poise

Type

Bool

Modules

ALL

Files

Any

Default value

0

bRossbyCut

Description

Van Saders, J. et al. (2019, ApJ, 872, 128) find that when the
stellar Rossby number exceeds 2.08, then the magnetic braking
is quenched. This flag enforces that behavior.

Type

Bool

Modules

STELLAR

Files

Any

Default value

0

bSeaIceModel

Description

model sea ice dynamics and heat flow?

Type

Bool

Modules

ALL

Files

Any

Default value

1

bSkipSeasEnabled

Description

Run annual before seasonal and allow skip seas?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bStellarEnc

Description

Model stellar encounters?

Type

Bool

Modules

ALL

Files

Any

Default value

1

bStopWaterLossInHZ

Description

If set to 1, then all water photolysis and hydrogen escape will
not occur for a planet in the habitable zone, defined to be
when the instellation .is less than the runaway greenhouse
threshold.

Type

Bool

Modules

ALL

Files

Any

Default value

1

bTimeEvolVelDisp

Description

Scale velocity dispersion of stars with sqrt(t)?

Type

Bool

Modules

ALL

Files

Any

Default value

1

bUseBondiLimited

Description

Force the atmospheric erosion of a H envelope to be
Bondi-limited where the mass loss is regulated by the sound
speed at the sonic point following Equation 4 from Owen & Wu
(2016). Note we compute the sound speed by assuming an
isothermal atmosphere composed of diatomic hydrogen and that
the local temperature is set by the local blackbody equlibrium
temperature (Owen & Wu 2016, Equation 2).

Type

Bool

Modules

ALL

Files

Any

Default value

0

bUseEnergyLimited

Description

Force the atmospheric erosion of a H envelope to be
energy-limited, i.e. the mass loss is directly proportional to
the incident XUV flux (Equation 5 from Luger et al. (2015)).

Type

Bool

Modules

ALL

Files

Any

Default value

0

bUseOrbParams

Description

Flag to use orbital parameters as inputs

Type

Bool

Modules

ALL

Files

Any

Default value

0

bUseRRLimited

Description

Force the atmospheric erosion of a H envelope to be
radiation/recombination-limited, i.e. the mass loss is directly
proportional to the sqrt of the incident XUV flux (Equation 13
from Luger et al. (2015) and Murray-Clay et al. (2009).

Type

Bool

Modules

ALL

Files

Any

Default value

0

bUseTidalRadius

Description

Fix radius used for CPL tidal equations?

Type

Bool

Modules

ALL

Files

Any

Default value

0

bVarDt

Description

Use variable timestep?

Type

Bool

Modules

ALL

Files

Any

Default value

None

d232ThMassCore

Description

Initial Core Mass of 232Th

Type

Double

Dimension(s)

mass

Modules

ALL

Files

Any

Default value

Primordial Earth

d232ThMassCrust

Description

Initial Crust Mass of 232Th

Type

Double

Dimension(s)

mass

Modules

ALL

Files

Any

Default value

Primordial Earth

d232ThMassMan

Description

Initial Mass of 232Th

Type

Double

Dimension(s)

mass

Modules

ALL

Files

Any

Default value

Primordial Earth

d232ThNumCore

Description

Initial Core Number of 232Th Atoms

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

Primordial Earth

d232ThNumCrust

Description

Initial Crust Number of 232Th Atoms

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

Primordial Earth

d232ThNumMan

Description

Initial Number of 232Th Atoms

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

Primordial Earth

d232ThPowerCore

Description

Initial Core Power Production from 232Th Atoms

Type

Double

Dimension(s)

energy/time

Modules

ALL

Files

Any

Default value

Primordial Earth

d232ThPowerCrust

Description

Initial Crust Power Production from 232Th Atoms

Type

Double

Dimension(s)

energy/time

Modules

ALL

Files

Any

Default value

Primordial Earth

d232ThPowerMan

Description

Initial Power Production from 232Th Atoms

Type

Double

Dimension(s)

energy/time

Modules

ALL

Files

Any

Default value

Primordial Earth

d235UMassCore

Description

Initial Core Mass of 235U

Type

Double

Dimension(s)

mass

Modules

ALL

Files

Any

Default value

Primordial Earth

d235UMassCrust

Description

Initial Crust Mass of 235U

Type

Double

Dimension(s)

mass

Modules

ALL

Files

Any

Default value

Primordial Earth

d235UMassMan

Description

Initial Mass of 235U

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

Primordial Earth

d235UNumCore

Description

Initial Core Number of 235U Atoms

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

Primordial Earth

d235UNumCrust

Description

Initial Crust Number of 235U Atoms

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

Primordial Earth

d235UNumMan

Description

Initial Number of 235U Atoms

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

Primordial Earth

d235UPowerCore

Description

Initial Core Power Production from 235U Atoms

Type

Double

Dimension(s)

energy/time

Modules

ALL

Files

Any

Default value

Primordial Earth

d235UPowerCrust

Description

Initial Crust Power Production from 235U Atoms

Type

Double

Dimension(s)

energy/time

Modules

ALL

Files

Any

Default value

Primordial Earth

d235UPowerMan

Description

Initial Power Production from 235U Atoms

Type

Double

Dimension(s)

energy/time

Modules

ALL

Files

Any

Default value

Primordial Earth

d238UMassCore

Description

Initial Core Mass of 238U

Type

Double

Dimension(s)

mass

Modules

ALL

Files

Any

Default value

Primordial Earth

d238UMassCrust

Description

Initial Crust Mass of 238U

Type

Double

Dimension(s)

mass

Modules

ALL

Files

Any

Default value

Primordial Earth

d238UMassMan

Description

Initial Mass of 238U

Type

Double

Dimension(s)

mass

Modules

ALL

Files

Any

Default value

Primordial Earth

d238UNumCore

Description

Initial Core Number of 238U Atoms

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

Primordial Earth

d238UNumCrust

Description

Initial Crust Number of 238U Atoms

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

Primordial Earth

d238UNumMan

Description

Initial Number of 238U Atoms

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

Primordial Earth

d238UPowerCore

Description

Initial Core Power Production from 238U Atoms

Type

Double

Dimension(s)

energy/time

Modules

ALL

Files

Any

Default value

Primordial Earth

d238UPowerCrust

Description

Initial Crust Power Production from 238U Atoms

Type

Double

Dimension(s)

energy/time

Modules

ALL

Files

Any

Default value

Primordial Earth

d238UPowerMan

Description

Initial Power Production from 238U Atoms

Type

Double

Dimension(s)

energy/time

Modules

ALL

Files

Any

Default value

Primordial Earth

d26AlMassCore

Description

Initial Core Mass of 26Al

Type

Double

Dimension(s)

mass

Modules

ALL

Files

Any

Default value

0

d26AlMassMan

Description

Initial Mantle Mass of 26Al

Type

Double

Dimension(s)

mass

Modules

ALL

Files

Any

Default value

0

d26AlNumCore

Description

Initial Core Number of 26Al Atoms

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

0

d26AlNumMan

Description

Initial Mantle Number of 26Al Atoms

Type

Double

Dimension(s)

mass

Modules

ALL

Files

Any

Default value

0

d26AlPowerCore

Description

Initial Core Power Production from 26Al Atoms

Type

Double

Dimension(s)

energy/time

Modules

ALL

Files

Any

Default value

0

d26AlPowerMan

Description

Initial Mantle Power Production from 26Al Atoms

Type

Double

Dimension(s)

energy/time

Modules

ALL

Files

Any

Default value

0

d40KMassCore

Description

Initial Core Mass of 40K

Type

Double

Dimension(s)

mass

Modules

ALL

Files

Any

Default value

Primordial Earth Units

d40KMassCrust

Description

Initial Crust Mass of 40K

Type

Double

Dimension(s)

mass

Modules

ALL

Files

Any

Default value

Primordial Earth

d40KMassMan

Description

Initial Mantle Mass of 40K

Type

Double

Dimension(s)

mass

Modules

ALL

Files

Any

Default value

Primordial Earth

d40KNumCore

Description

Initial Core Number of 40K Atoms

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

Primordial Earth

d40KNumCrust

Description

Initial Crust Number of 40K Atoms

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

Primordial Earth

d40KNumMan

Description

Initial Mantle Number of 40K Atoms

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

Primordial Earth

d40KPowerCore

Description

Initial Core Power Production from 40K Atoms

Type

Double

Dimension(s)

energy/time

Modules

ALL

Files

Any

Default value

Primordial Earth

d40KPowerCrust

Description

Initial Crust Power Production from 40K Atoms

Type

Double

Dimension(s)

energy/time

Modules

ALL

Files

Any

Default value

Primordial Earth

d40KPowerMan

Description

Initial Mantle Power Production from 40K Atoms

Type

Double

Dimension(s)

energy/time

Modules

ALL

Files

Any

Default value

Primordial Earth

dAblateFF

Description

Ice ablation fudge factor

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

2.3

dActViscMan

Description

Mantle viscosity activation energy, E_nu in eq (159) of Barnes
et al (2020). Physically it is a measure of how sensitive
viscosity is to temperature.

Type

Double

Dimension(s)

pressure

Modules

ALL

Files

Any

Default value

Default is ACTVISCMAN

dAdJumpC2CMB

Description

Adiabatic temp jump from ave core to CMB

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

Default is ADJUMPC2CMB

dAdJumpM2LM

Description

Constant used to relate TLMan to TMan (TLMan=AdJumpM2LM*TMan),
analogous to AdJumpM2UM.

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

Default is ADJUMPM2LM

dAdJumpM2UM

Description

Constant used to relate TUMan to TMan (TUMan=AdJumpM2UM*TMan),
eta_UM below eq (8) in Driscoll & Bercovici (2014).

Type

Double

Dimension(s)

temperature

Modules

ALL

Files

Any

Default value

Default is ADJUMPM2UM

dAge

Description

System Age

Type

Double

Custom unit

Gyr

Dimension(s)

time

Modules

ALL

Files

Any

Default value

0

dAlbedoGlobal

Description

Globally averaged albedo

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Body Only

Default value

0.3

dAlbedoLand

Description

albedo of dry land

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

0.363

dAlbedoWater

Description

albedo of open water

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

0.263

dAtmGasConst

Description

Atmospheric Gas Constant

Type

Double

Dimension(s)

energy/temperature/mass

Modules

ALL

Files

Any

Default value

4124

dAtmXAbsEffH

Description

XUV absoprtion efficiency parameter, epsilon_{XUV}, in Eq. (A1)
in Barnes et al. (2019). Must lie in the range [0,1].

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

0.15

dAtmXAbsEffH2O

Description

XUV absoprtion efficiency parameter for water vapor as defined
in Luger & Barnes (2015, AsBio, 15, 57). Must lie in range
[0,1].

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

0.30

dCBPM0

Description

Circumbinary planet initial mean anomaly for use in the BINARY
module

Type

Double

Dimension(s)

angle

Modules

ALL

Files

Any

Default value

0.0 degrees

dCBPPsi

Description

Circumbinary planet initial R, phi oscillation phase angle

Type

Double

Dimension(s)

angle

Modules

ALL

Files

Any

Default value

0.0 degrees

dCBPZeta

Description

Circumbinary planet initial z oscillation phase angle

Type

Double

Dimension(s)

angle

Modules

ALL

Files

Any

Default value

0.0 degrees

dCO2MassMOAtm

Description

Initial CO2 mass in the system

Type

Double

Dimension(s)

mass

Modules

ALL

Files

Any

Default value

0 TO

dCosObl

Description

Planet formation simulations predict that an isotropic
distribution of rotational angular momentum vectors is a
typical outcome. This result is identical to a uniform
distribution in cosine obliquity. Use this option to sample a
realistic distribution of initial obliquities.

Type

Double

Dimension(s)

nd

Modules

BINARY DistRot EqTide POISE SpiNBody

Files

Body Only

Default value

0.5

dDAdCore

Description

Length scale for core iron adiabatic temperature profile, D_N
in eq (174) of Barnes et al (2020). It determines the curvature
of the core iron adiabat as a function of radius so that
smaller values cause the core adiabat to drop faster with
radius.

Type

Double

Dimension(s)

length

Modules

ALL

Files

Any

Default value

Default is DADCORE

dDepthMO

Description

Initial depth of the magma ocean

Type

Double

Dimension(s)

length

Modules

ALL

Files

Any

Default value

core radius

dDfcrit

Description

When running DistOrb with other modules that modify the
semi-major axis, set his argument to be the maximum relative
change in the Laplace coefficients before recalculating the
values. Setting this value to be low can cause the simulation
to run very slowly, with negligible gain in accuracy.

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

0.1

dDiffusion

Description

Heat diffusion coefficient

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

0.44

dDMDensity

Description

Local dark matter density

Type

Double

Dimension(s)

mass/length^3

Modules

ALL

Files

Any

Default value

0.01 Msun pc^3

dDTChiRef

Description

Core liquidus (melting temperature) depression due to the
presence of light elements in the core. See eq (175) in Barnes
et al (2020).

Type

Double

Dimension(s)

temperature

Modules

ALL

Files

Any

Default value

Value in thermint.h

dEcc

Description

Orbital Eccentricity

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Body Only

Default value

0

dEccAmp

Description

Amplitude of forced eccentricity oscill

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

0.1

dEccPer

Description

Period of forced eccentricity oscill

Type

Double

Dimension(s)

time

Modules

ALL

Files

Any

Default value

50000

dElecCondCore

Description

Electrical conductivity of core

Type

Double

Dimension(s)

time^3*ampere^2/mass/length

Modules

ALL

Files

Any

Default value

Default is ELECCONDCORE

dEncounterRad

Description

Radius at which stellar encounters occur

Type

Double

Dimension(s)

length

Modules

ALL

Files

Any

Default value

206265 AU

dEnergyBin

Description

Number of energies consider between the minimum and maximum
energies to calculate the luminosity by flares

Type

Double

Dimension(s)


Modules

ALL

Files

Any

Default value

100 energies between dFlareMinEnergy and dFlareMaxEnergy

dEnvelopeMass

Description

Initial Envelope Mass

Type

Double

Dimension(s)

mass

Modules

ALL

Files

Any

Default value

0

dEruptEff

Description

Fraction of heat in mantle melt that escapes to the surface.
See eq (187) in Barnes et al (2020).

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

ERUPTEFF

dEta

Description

The timestep will be set to dEta times the smallest
instantaneous timescale, i.e. min(x/(dx/dt) where x represents
the primary variables.

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

1

dFixIceLat

Description

Force ice cap latitude to this value

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

None

dFixMeltfactorUMan

Description

Boolean: 1 tells the code to used a fixed upper mantle
viscosity melt factor, 0 tells the code to compute it from the
local temperature (default).

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

Default is FIXMELTFACTORUMAN

dFlareMaxEnergy

Description

Maximum Flare Energy to consider

Type

Double

Dimension(s)


Modules

ALL

Files

Any

Default value

10^29 J

dFlareMinEnergy

Description

Minimum Flare Energy to consider

Type

Double

Dimension(s)


Modules

ALL

Files

Any

Default value

10^26 J

dFlareSlope

Description

Slope for flare frequency distribution. The user needs to input
the module of the value for this parameter. The negative signal
can be use only for input the value in flares/day 1/log10(erg).

Type

Double

Dimension(s)


Modules

ALL

Files

Any

Default value

-0.68 (Proxima Centauri)

dFlowTemp

Description

Temperature of the hydrodynamic flow

Type

Double

Dimension(s)

temperature

Modules

ALL

Files

Any

Default value

400

dFreeEcc

Description

The free eccentricity of a circumbinary planet with the BINARY
module. Must be in the range [0,1).

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

0.0

dFreeInc

Description

The free inclination of a circumbinary planet with the BINARY
module. Must be in the range [0,pi).

Type

Double

Dimension(s)

angle

Modules

ALL

Files

Any

Default value

0.0 degrees

dFrzTSeaIce

Description

Temp of sea ice formation in seasonal model

Type

Double

Dimension(s)

temperature

Modules

ALL

Files

Any

Default value

-2 deg C

dGalacDensity

Description

Density of galactic environment

Type

Double

Dimension(s)

mass/length^3

Modules

ALL

Files

Any

Default value

0.102

dGasDensity

Description

Local ISM density

Type

Double

Dimension(s)

mass/length^3

Modules

ALL

Files

Any

Default value

0.05 Msun pc^3

dHalt232ThPower

Description

Minimum 232Th Power

Type

Double

Dimension(s)

energy/time

Modules

ALL

Files

Any

Default value

0

dHalt235UPower

Description

Minimum 235U Power

Type

Double

Dimension(s)

energy/time

Modules

ALL

Files

Any

Default value

0

dHalt238UPower

Description

Minimum 238U Power

Type

Double

Dimension(s)

energy/time

Modules

ALL

Files

Any

Default value

0

dHalt40KPower

Description

Minimum 40K Power

Type

Double

Dimension(s)

energy/time

Modules

ALL

Files

Any

Default value

0

dHaltMaxEcc

Description

Maximum eccentricity value that halts ntegration

Type

Double

Dimension(s)

nd

Modules

DistOrb EqTide

Files

Any

Default value

1

dHaltMaxMutualInc

Description

The execution halts when dHaltMaxMutualInc is reached. The
mutual (or relative) inclination is the value of the angle
between the orbital angular momentum vector of two bodies. If
set to 0, then the mutual inclination will not be checked every
timestep for halts, i.e. the code should run faster. For
DistOrb calculations with the RD4 null model, mutual
inclinations above 35.00 degrees, especially if the
eccentricities are significant, should be interpreted
cautiously. For the DistOrb-null model, values above 10.00
degrees are suspect. SpiNBody is accurate for any value.

Type

Double

Dimension(s)

angle

Modules

DistOrb SpiNBody

Files

Any

Default value

0 [not checked]

dHaltMinEcc

Description

Minimum Eccentricity Value that Halts Integration

Type

Double

Dimension(s)

nd

Modules

DistOrb EqTide SpiNBody

Files

Any

Default value

-1

dHaltMinObl

Description

Minimum Obliquity Value that Halts Integration

Type

Double

Dimension(s)

angle

Modules

DistOrb EqTide SpiNBody

Files

Any

Default value

-1 degrees

dHaltMinSemi

Description

Minimum Semi-Major Axis Value that Halts Integration

Type

Double

Custom unit

au

Dimension(s)

length

Modules

EqTide SpiNBody

Files

Any

Default value

0

dHaltMinTCore

Description

Halt at Minimum Core Temperature

Type

Double

Dimension(s)

temperature

Modules

ALL

Files

Any

Default value

0 K

dHaltMinTMan

Description

Halt at Minimum Mantle Temperature

Type

Double

Dimension(s)

temperature

Modules

ALL

Files

Any

Default value

0 K

dHaltRadPower

Description

Minimum Total Radgiogenic Power

Type

Double

Dimension(s)

energy/time

Modules

ALL

Files

Any

Default value

0

dHeatCapAnn

Description

Surface heat capacity in annual model

Type

Double

Dimension(s)

energy/temperature

Modules

ALL

Files

Any

Default value

0.2

dHeatCapLand

Description

Land heat capacity in seasonal model

Type

Double

Dimension(s)

energy/temperature

Modules

ALL

Files

Any

Default value

1.42e7

dHeatCapWater

Description

Water heat capacity per meter in seasonal model

Type

Double

Dimension(s)

energy/temperature

Modules

ALL

Files

Any

Default value

4.2e6

dHostBinArgP

Description

Arg periapse of host binary

Type

Double

Dimension(s)

angle

Modules

ALL

Files

Any

Default value

0.0

dHostBinEcc

Description

eccentricity of host binary

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

0.51

dHostBinInc

Description

inclination of host binary

Type

Double

Dimension(s)

angle

Modules

ALL

Files

Any

Default value

60.0 deg

dHostBinLongA

Description

Long of ascending node of host binary

Type

Double

Dimension(s)

angle

Modules

ALL

Files

Any

Default value

0.0

dHostBinMass1

Description

mass of larger host binary star

Type

Double

Dimension(s)

mass

Modules

ALL

Files

Any

Default value

1.1 Msun

dHostBinSemi

Description

Semi-major of host binary

Type

Double

Dimension(s)

length

Modules

ALL

Files

Any

Default value

17.57 AU

dIceAlbedo

Description

Albedo of ice

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

0.6

dIceDepRate

Description

Deposition rate of ice/snow to form ice sheets

Type

Double

Dimension(s)

length/time

Modules

ALL

Files

Any

Default value

2.9e-5

dImK2ManOrbModel

Description

ImK2Man model to use in orbital equations (1=ImK2(T), 2=k2/Q(T)

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

Default is IMK2MANORBMODEL

dInitIceHeight

Description

Sets initial ice sheet height

Type

Double

Dimension(s)

length

Modules

ALL

Files

Any

Default value

50

dInitIceLat

Description

Sets initial ice sheet latitude

Type

Double

Dimension(s)

angle

Modules

ALL

Files

Any

Default value

90

dJeansTime

Description

Time at which flow transitions to Jeans escape

Type

Double

Dimension(s)

time

Modules

ALL

Files

Any

Default value

1 Gyr

dK2

Description

Love Number of Degree 2

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

1

dLandFrac

Description

Fraction of land on the planetary surface

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

0.34

dLapseR

Description

Dry adiabatic lapse rate (for elev feedback)

Type

Double

Dimension(s)

temperature/length

Modules

ALL

Files

Any

Default value

9.8e-3 C/m

dLL13K0

Description

Lee+Leung 2013 Radial Epicyclic Frequency

Type

Double

Dimension(s)

time^-1

Modules

ALL

Files

Any

Default value

1 /yr

dLL13N0

Description

Lee+Leung 2013 Mean Motion

Type

Double

Dimension(s)

time^-1

Modules

ALL

Files

Any

Default value

1 /yr

dLL13PhiAB

Description

Binary Initial Mean Anomaly for use in the BINARY module

Type

Double

Dimension(s)

angle

Modules

ALL

Files

Any

Default value

0.0 degrees

dLL13V0

Description

Lee+Leung 2013 Radial Epicyclic Frequency

Type

Double

Dimension(s)

time^-1

Modules

ALL

Files

Any

Default value

1 /yr

dLuminosity

Description

Initial Luminosity

Type

Double

Dimension(s)

energy/time

Modules

ALL

Files

Any

Default value

0

dLXUV

Description

Total XUV Luminosity – Unsupported!

Type

Double

Dimension(s)

energy/time

Modules

STELLAR

Files

Body Only

Default value

-1

dLXUVFlareConst

Description

XUV luminosity of flares

Type

Double

Dimension(s)

energy/time

Modules

ALL

Files

Any

Default value

10^22 Watts or 10^29 erg/s

dMagMomCoef

Description

Coefficient in front of magnetic moment scaling law, gamma_d in
eq (192) in Barnes et al (2020).

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

Default is MAGMOMCOEF

dManHFlowPref

Description

Coefficient in front of HflowUMan, value can be anything. If
StagLid>0 then ManHFlowPref=HFLOWREDUCTSTAG. Otherwise default
behavior is 1.

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

Default is MANHFLOWPREF

dManMeltDensity

Description

Density of the molten mantle

Type

Double

Dimension(s)

mass/length^3

Modules

ALL

Files

Any

Default value

4000 kg/m^3

dMass

Description

Mass

Type

Double

Custom unit

Mearth

Dimension(s)

mass

Modules

ALL

Files

Body Only

Default value

1 Earth Mass

dMassFracFeOIni

Description

Initial Mass Fraction of FeO in the mantle

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

BSE Earth: 0.0788

dMaxLockDiff

Description

In EqTide, this is the maximum relative difference between the
actual spin rate and the equilibrium spin rate without tidal
locking. See bForceEqSpin for more information. bForceEqSpin
must be set to 1 to enable this option.

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

0

dMeanA

Description

Mean anomaly

Type

Double

Dimension(s)

angle

Modules

ALL

Files

Any

Default value

0.0

dMeanMotion

Description

Orbital Mean Motion

Type

Double

Custom unit

/Year

Dimension(s)

time^-1

Modules

ALL

Files

Body Only

Default value

1 /yr

dMeltfactorLMan

Description

Use to fix the lower mantle viscosity melt factor, defined as
epsilon_phase in eq (159,160) of Barnes et al (2020). Note that
this is an option. Default behavior is to compute this factor
from the local temperature.

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

1.000000

dMeltfactorUMan

Description

Use to fix the upper mantle viscosity melt factor, defined as
epsilon_phase in eq (159,160) of Barnes et al (2020). Note that
this is an option. Default behavior is to compute this factor
from the local temperature.

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

1.000000

dMinEnvelopeMass

Description

Minimum Envelope Mass

Type

Double

Dimension(s)

mass

Modules

ALL

Files

Any

Default value

1.e-8 Earth

dMinIceSheetHeight

Description

The minimum thickness of permanent ice in a latitude bin for it
to be labeled ice-covered. In some cases, such as rapid
thawing, a latituden can have a very low value of ice height,
primarily for numerical reasons. Parameter forces unphysically
small values of the ice height to be ignored.

Type

Double

Dimension(s)

length

Modules

ALL

Files

Any

Default value

0.001

dMinKTide

Description

Minimum value for stellar gravitaitonal enhancement of mass
loss

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

0.1

dMinStellarApproach

Description

Minimum close approach distance to primary

Type

Double

Dimension(s)

length

Modules

ALL

Files

Any

Default value

1 AU

dMinSurfWaterMass

Description

Minimum Surface Water Mass

Type

Double

Dimension(s)

mass

Modules

ALL

Files

Any

Default value

1.e-5 TO

dMinValue

Description

Minimum Non-Zero Value of Eccentricity and Obliquities

Type

Double

Dimension(s)

nd

Modules

DistOrb DistRot EqTide POISE SpiNBody

Files

Any

Default value

0

dMixingDepth

Description

Mixing depth of ocean in seasonal model

Type

Double

Dimension(s)

length

Modules

ALL

Files

Any

Default value

70

dNuLandWater

Description

Coefficient of land-ocean heat flux

Type

Double

Dimension(s)

energy/length^2

Modules

ALL

Files

Any

Default value

0.81

dObliqAmp

Description

Amplitude of forced obliquity oscill

Type

Double

Dimension(s)

angle

Modules

ALL

Files

Any

Default value

50 deg

dObliqPer

Description

Period of forced obliquity oscill

Type

Double

Dimension(s)

time

Modules

ALL

Files

Any

Default value

50000

dObliquity

Description

Obliquity

Type

Double

Dimension(s)

angle

Modules

BINARY DistRot EqTide POISE SpiNBody

Files

Body Only

Default value

0

dOrbPeriod

Description

Orbital Period

Type

Double

Custom unit

Days

Dimension(s)

time

Modules

ALL

Files

Body Only

Default value

1 year

dOutputTime

Description

Output Interval

Type

Double

Custom unit

Years

Dimension(s)

time

Modules

ALL

Files

Any

Default value

1 year

dOxygenMantleMass

Description

Initial Oxygen Mass in the Mantle

Type

Double

Dimension(s)

mass

Modules

ALL

Files

Any

Default value

0

dOxygenMass

Description

The initial oxygen mass in the atmosphere.

Type

Double

Dimension(s)

mass

Modules

ALL

Files

Any

Default value

0

dpCO2

Description

Partial pressure of CO2 in atmosphere

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

3.3e-4

dPositionXSpiNBody

Description

X position of the body

Type

Double

Dimension(s)

length

Modules

ALL

Files

Any

Default value

0

dPositionYSpiNBody

Description

Y position of the body

Type

Double

Dimension(s)

length

Modules

ALL

Files

Any

Default value

0

dPositionZSpiNBody

Description

Z position of the body

Type

Double

Dimension(s)

length

Modules

ALL

Files

Any

Default value

0

dPrecRate

Description

Fixed rate of axial precession (angle/s)

Type

Double

Dimension(s)

angle/time

Modules

ALL

Files

Any

Default value

7.7261e-12

dPresSWind

Description

Stellar Wind Pressure at body

Type

Double

Dimension(s)

pressure

Modules

ALL

Files

Any

Default value

Default is EPRESSWIND

dPresXUV

Description

Pressure at base of thermosphere

Type

Double

Dimension(s)

pressure

Modules

ALL

Files

Any

Default value

5 Pa

dRadGyra

Description

Radius of Gyration; moment of inertia constant

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Body Only

Default value

0.5

dRadius

Description

Radius

Type

Double

Custom unit

Rearth

Dimension(s)

length

Modules

ALL

Files

Body Only

Default value

1 Earth Radius

dRefHeight

Description

Reference height of atmos temp (for elev feedback)

Type

Double

Dimension(s)

length

Modules

ALL

Files

Any

Default value

1000 m

dRForm

Description

Galactic formation radius

Type

Double

Dimension(s)

length

Modules

ALL

Files

Any

Default value

4.5 kpc

dRotPeriod

Description

Rotation Period

Type

Double

Custom unit

Days

Dimension(s)

time

Modules

ALL

Files

Body Only

Default value

1 Day

dRotRate

Description

Rotational Angular Frequency

Type

Double

Custom unit

/Day

Dimension(s)

time^-1

Modules

ALL

Files

Body Only

Default value

2*pi/day

dRotVel

Description

Rotational Velocity

Type

Double

Custom unit

km/s

Dimension(s)

mass/time

Modules

ALL

Files

Body Only

Default value

0

dSatXUVFrac

Description

After formation stars emit a nearly constant amount of XUV
radiation for a time called the “saturated” phase. This
parameter sets that value relative to the total (bolometric)
luminosity. Must lie in range [0,1].

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

1e-3

dSatXUVTime

Description

The time a star will remain in its “saturated” phase.

Type

Bool

Modules

ALL

Files

Any

Default value

0.1 Gyr

dSeaIceConduct

Description

Heat conductivity of sea ice

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

2

dSeasOutputTime

Description

Output interval for seasonal parameters

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

0

dSemi

Description

Semi-Major Axis

Type

Double

Custom unit

AU

Dimension(s)

length

Modules

ALL

Files

Body Only

Default value

1 AU

dShModRef

Description

Reference kinematic mantle shear modulus coefficient, mu_ref in
eq (160) of Barnes et al (2020). Analogous to reference
viscosity ViscRef.

Type

Double

Dimension(s)

pressure

Modules

ALL

Files

Any

Default value

Default is SHMODREF

dSpecMomInertia

Description

Specific moment of inertia of polar axis

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

0.33

dSpinUpTol

Description

Tolerance for spin up phase

Type

Double

Dimension(s)

temperature

Modules

ALL

Files

Any

Default value

0.1 deg C

dStagLid

Description

Stagnant Lid Switch

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

Default is STAGLID

dStarScaleL

Description

Stellar radial scale length in MW

Type

Double

Dimension(s)

length

Modules

ALL

Files

Any

Default value

2.4 kpc

dStiffness

Description

Effective stiffness of mantle

Type

Double

Dimension(s)

pressure

Modules

ALL

Files

Any

Default value

Default is STIFFNESS

dStopTime

Description

Integration Stop Time

Type

Double

Custom unit

Years

Dimension(s)

time

Modules

ALL

Files

Any

Default value

10 years

dSurfAlbedo

Description

Albedo of (ice-free) surface

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

0.3

dSurfTemp

Description

Initial surface temp

Type

Double

Dimension(s)

temperature

Modules

ALL

Files

Any

Default value

4000 K

dSurfWaterMass

Description

Initial Surface Water Mass

Type

Double

Dimension(s)

mass

Modules

ALL

Files

Any

Default value

0

dSyncEcc

Description

Minimum Eccentricity for Non-Synchronous Rotation

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

0

dTCore

Description

Initial Core Temperature

Type

Double

Dimension(s)

temperature

Modules

ALL

Files

Any

Default value

6000 K

dTemperature

Description

Initial effective temperature

Type

Double

Dimension(s)

temperature

Modules

ALL

Files

Any

Default value

TSUN

dTGlobalInit

Description

Estimate of initial global temperature

Type

Double

Dimension(s)

temperature

Modules

ALL

Files

Any

Default value

14.85

dThermTemp

Description

The thermal temperature of a planet heated by radiation from
the primary. The user may set a value which will then remain
constant for the simulation. If the user does not specify a
value, then it will be calculated automatically from the
formula T = (F(1-A)/sigma)^0.25, where F is the incident
radiation, A is albedo and sigma is the Steffan-Boltzman
constant.

Type

Double

Dimension(s)

temperature

Modules

ALL

Files

Any

Default value

400

dTidalQ

Description

Tidal Quality Factor

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

1e6

dTidalQEnv

Description

Envelope Tidal Quality Factor

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

1.0e4

dTidalQMantle

Description

Tidal Q of Mantle

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

100

dTidalQOcean

Description

Ocean Tidal Quality Factor

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

12

dTidalRadius

Description

Eqtide Tidal Radius

Type

Double

Dimension(s)

length

Modules

ALL

Files

Any

Default value

1 Earth Radius

dTidalTau

Description

Tidal Time Lag

Type

Double

Dimension(s)

time

Modules

ALL

Files

Any

Default value

1 Second

dTimeStep

Description

Integration Timestep

Type

Double

Dimension(s)

time

Modules

ALL

Files

Any

Default value

1 year

dTMan

Description

Initial Mantle Temperature

Type

Double

Dimension(s)

temperature

Modules

ALL

Files

Any

Default value

3000 K

dTMigration

Description

Time of radial migration

Type

Double

Dimension(s)

time

Modules

ALL

Files

Any

Default value

3 Gy

dTSurf

Description

Initial Surface Temperature

Type

Double

Dimension(s)

temperature

Modules

ALL

Files

Any

Default value

300 K

dVelXSpiNBody

Description

X velocity of the body

Type

Double

Dimension(s)

length/time

Modules

ALL

Files

Any

Default value

0

dVelYSpiNBody

Description

Y velocity of the body

Type

Double

Dimension(s)

length/time

Modules

ALL

Files

Any

Default value

0

dVelZSpiNBody

Description

Z velocity of the body

Type

Double

Dimension(s)

length/time

Modules

ALL

Files

Any

Default value

0

dViscJumpMan

Description

Numerical factor that is used to get the lower mantle thermal
boundary viscosity visc_UMan from the upper mantle viscosity
visc_LM by visc_LM=ViscJumpMan*visc_UM. See eq (163) in Barnes
et al (2020).

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

VISCJUMPMAN

dViscMeltB

Description

Constant used to compute the reduction in mantle viscosity due
to the presence of melt. See eq (161) in Barnes et al (2020)
and Costa et al (2009) for values for different materials.

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

Default is VISCMELTB

dViscMeltDelta

Description

Constant used to compute the reduction in mantle viscosity due
to the presence of melt. See eq (161) in Barnes et al (2020)
and Costa et al (2009) for values for different materials.

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

Default is VISCMELTDELTA

dViscMeltGamma

Description

Constant used to compute the reduction in mantle viscosity due
to the presence of melt. See eq (162) in Barnes et al (2020)
and Costa et al (2009) for values for different materials.

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

Default is VISCMELTGAMMA

dViscMeltPhis

Description

Viscosity-melt reduction coefficient, “phi*” in Eq. 8 of
Driscoll Bercovici (2015).

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

Default is VISCMELTPHIS

dViscMeltXi

Description

Constant used to compute the reduction in mantle viscosity due
to the presence of melt. See eq (162) in Barnes et al (2020)
and Costa et al (2009) for values for different materials.

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

Default is VISCMELTXI

dViscRef

Description

Coefficient in front of viscosity equation see eq (159) in
Barnes et al. (2020).

Type

Double

Dimension(s)

length^2/time

Modules

ALL

Files

Any

Default value

1

dViscUMan

Description

Upper mantle viscosity

Type

Double

Dimension(s)

length^2/time

Modules

DistRot ThermInt

Files

Body Only

Default value

0

dWaterMassAtm

Description

Initial Water Mass in the atmosphere

Type

Double

Dimension(s)

mass

Modules

ALL

Files

Any

Default value

1 Terrestrial Ocean

dWaterPartCoeff

Description

Water partition coefficient between melt and solid

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

0.01

dXUVBeta

Description

After the “saturation” phase, the ratio of the XUV to total
luminosity will follow a power law followinfg this exponent.
Units are gigayears.

Type

Double

Dimension(s)

nd

Modules

ALL

Files

Any

Default value

1.23

iBodyType

Description

BodyType

Type

Int

Modules

ALL

Files

Body Only

Default value

0 Planet

iDigits

Description

For all floating point output, print this many number of digits
after the decimal point

Type

Int

Modules

ALL

Files

Any

Default value

4

iIceDt

Description

Minimum ice sheet timestep (unit orbital period)

Type

Int

Modules

ALL

Files

Any

Default value

5

iLatCellNum

Description

Number of latitude cells used in climate model

Type

Int

Modules

ALL

Files

Any

Default value

50

iNStepInYear

Description

Number of time-steps/year in seasonal model

Type

Int

Modules

ALL

Files

Any

Default value

60

iNumYears

Description

Number of years to run seasonal model

Type

Int

Modules

ALL

Files

Any

Default value

10

iOLRModel

Description

Outgoing longwave rad model

Type

Int

Modules

ALL

Files

Any

Default value

sms09

iRandSeed

Description

Seed for random number generator (stellar encounters)

Type

Int

Modules

ALL

Files

Any

Default value

42

iReRunSeas

Description

how often to rerun seasonal in ice sheet model, in number of
orbital periods

Type

Int

Modules

ALL

Files

Any

Default value

500

iSciNot

Description

Logarithm to Change from Standard to Scientific Notation

Type

Int

Modules

ALL

Files

Any

Default value

4

iVerbose

Description

Set how much text is written to the screen. 0 = no output, 1 =
only errors, 2 = progress updates, 3 = statements about input
choices, 4 = information about unit choices, 5 = all possible
output. Note that levels 0 and 5 can be set at execution with
the -q and -v options, respectively.

Type

Int

Modules

ALL

Files

Any

Default value

3

saBodyFiles

Description

Input files for each body

Type

String-Array

Modules

ALL

Files

Primary Only

Default value

None

saGridOutput

Description

Gridded Output Parameter(s)

Type

String-Array

Modules

POISE

Files

Body Only

Default value

None

saModules

Description

List of names of modules to be applied to the body. Spelling
must be exact, but any capitalization works

Type

String-Array

Modules

ALL

Files

Primary Only

Default value

none

saOutputOrder

Description

Output Parameter(s)

Type

String-Array

Modules

ALL

Files

Body Only

Default value

None

saTidePerts

Description

Names of bodies engaged in tidal evolution

Type

String-Array

Modules

ALL

Files

Any

Default value

none

sAtmXAbsEffH2OModel

Description

If BOLMONT16 is selected, then the value of dAtmXAbsEffH2O will
follow the model of Bolmont et al. (2017, MNRAS, 464, 3728).
NONE will not change the input value for dAtmXAbsEffH2O.

Type

String

Modules

ALL

Files

Any

Default value

NONE

Allowed values

BOLMONT16 NONE

sClimateModel

Description

Use annual or seasonal model

Type

String

Modules

ALL

Files

Any

Default value

ann

sColor

Description

Hexadecimal color code for the body to be used in vplot

Type

String

Modules

ALL

Files

Body Only

Default value

000000

sFileOrbitData

Description

File containing pre-computed orbital data. The file must have
the following format: Time SemiMajorAxis Eccentricity
Inclination ArgPericenter LongAscNode MeanAnomaly. The units of
those parameters is assumed to be the same as body being
simulated. When using this option, the integration must used a
fixed timestep (bVarDt = 0), and the timestep (dTimeStep) must
equal the cadence in the file, with time units in the
sFileOrbitData file assumed to be the same as for the body
file. See bReadOrbitData for more information.

Type

String

Modules

ALL

Files

Any

Default value

orbit.txt

sFileOrbitOblData

Description

Name of file containing orbit and obliquity time series

Type

String

Modules

ALL

Files

Any

Default value

Obl_data.txt

sGeography

Description

Type of land distribution

Type

String

Modules

ALL

Files

Any

Default value

uni3

sHZModel

Description

If KOPPARAPU13 is selected then the Recent Venus, Runaway
Greenhouse, Maximum Greenhouse, and Early Mars habitable zone
limits will be calculated from Kopparapu, R. et al. (2013, ApJ,
765, 131).

Type

String

Modules

ALL

Files

Any

Default value

Kopparapu13

sIntegrationMethod

Description

Integration Method: Euler, Runge-Kutta4 (Default =
Runge-Kutta4)

Type

String

Modules

ALL

Files

Any

Default value

Runge-Kutta4

sLogFile

Description

Log File Name

Type

String

Modules

ALL

Files

Any

Default value

null

sMagBrakingModel

Description

Magnetic braking model.

Type

String

Modules

ALL

Files

Any

Default value

REINERS

Allowed values

REINERS, SKUMANICH, MATT, NONE

sMagmOcAtmModel

Description

Atmospheric flux model: Grey or Petit

Type

String

Modules

ALL

Files

Any

Default value

GREY

sMassRad

Description

Mass-Radius Relationship for Central Body: GS99 RH00 BO06
Sotin07

Type

String

Modules

BINARY DistRot EqTide POISE SpiNBody STELLAR

Files

Body Only

Default value

None

sOrbitModel

Description

If RD4 is selected for sOrbitModel, the code will halt if an
eccentricity reaches ~0.69 unless this flag is set to true.

Type

String

Modules

ALL

Files

Any

Default value

rd4

sOutFile

Description

Name of Output File

Type

String

Modules

ALL

Files

Any

Default value

cSystemName.backward

sPlanetRadiusModel

Description

If LOPEZ12 is selected, the planet radius will follow the model
in Lopez et al. (2012, ApJ, 761, 59). PROXCENB will use the
model for Proxima b in Barnes et al. (2016, arXiv:1608.06919).
LEHMER17 is the Lehmer & Catling (2017, ApJ, 845, 130). NONE
will cause the radius to remain constant.

Type

String

Modules

ALL

Files

Any

Default value

NONE

Allowed values

LOPEZ12 PROXCENB LEHMER17 NONE.

sRadioHeatModel

Description

Radiogenic heating model

Type

String

Modules

ALL

Files

Any

Default value

NONE

sStellarModel

Description

If REINERS is selected, the stellar magnetic braking model of
Reiners & Mohanty (2012, ApJ, 746, 43) is used to modify the
rotation rate. SKUMANICH uses the model from Skumanich, A.
(1972, ApJ, 171, 565). MATT uses the model from Matt, S. et al.
(2015, ApJ, 799, 23). NONE applies no magnetic torque.

Type

String

Modules

ALL

Files

Any

Default value

BARAFFE

Allowed values

BARAFFE PROXIMA NONE

sSystemName

Description

System Name

Type

String

Modules

ALL

Files

Primary Only

Default value

None - must be supplied

sTideModel

Description

Tidal Model: p2 [constant-phase-lag, 2nd order] t8
[constant-time-lag, 8th order]

Type

String

Modules

ALL

Files

Any

Default value

p2

sUnitAngle

Description

Angle Units: Degrees Radians

Type

String

Modules

ALL

Files

Any

Default value

Radians

sUnitLength

Description

Length Units: cm m km Solar Earth Jupiter AU

Type

String

Modules

ALL

Files

Any

Default value

cm

sUnitMass

Description

Mass Units: Grams Kilograms Solar Earth Jupiter Saturn

Type

String

Modules

ALL

Files

Any

Default value

grams

sUnitTemp

Description

Temperature Units: Kelvin Celsius Farenheit

Type

String

Modules

ALL

Files

Any

Default value

Kelvin

sUnitTime

Description

Time Units: Seconds, Days Years Myr Gyr

Type

String

Modules

ALL

Files

Any

Default value

Seconds

sWaterLossModel

Description

The water loss rate will be determined by the selected model.
The options are LB15, LBEXACT, and TIAN.

Type

String

Modules

ALL

Files

Any

Default value

LBEXACT

Allowed values

LB15 LBEXACT TIAN

sWindModel

Description

If REINERS is selected, the stellar wind model of Reiners and
Mohanty (2012, ApJ, 746, 43) is used to modify the rotation
rate.

Type

String

Modules

ALL

Files

Any

Default value

REINERS

sXUVModel

Description

This parameter sets the XUV evolution model used in STELLAR.
Setting this to RIBAS (default) will evolve the XUV luminosity
according to the saturated power law of Ribas et al (2005, ApJ,
611, 680), while setting it to REINERS will use the empirical
relations of Reiners, Schussler and Passegger (2014, ApJ, 794,
144). Please note that the latter model has not been fully
vetted. Users may also set this parameter to NONE, in which
case the XUV luminosity will remain constant.

Type

String

Modules

STELLAR

Files

Any

Default value

RIBAS

Allowed values

RIBAS REINERS NONE

Output Parameters

These options follow the argument saOutputOrder.

232ThEnFlux

Description

Surface Energy Flux from 232Th

Custom unit

W/m^2

Modules

RadHeat

232ThMassCore

Description

Total Mass of 232Th

Custom unit

Mearth

Modules

RadHeat

232ThMassCrust

Description

Total Mass of 232Th

Custom unit

Mearth

Modules

RadHeat

232ThMassMan

Description

Total Mass of 232Th

Custom unit

Mearth

Modules

RadHeat

232ThNumCore

Description

Total Number of 232Th Atoms

Custom unit

Initial Primordial Earth Number

Modules

RadHeat

232ThNumCrust

Description

Total Number of 232Th Atoms

Custom unit

Initial Primordial Earth Number

Modules

RadHeat

232ThNumMan

Description

Total Number of 232Th Atoms

Custom unit

Initial Primordial Earth Number

Modules

RadHeat

232ThPowerCore

Description

Total Power Generated by 232Th

Custom unit

TW

Modules

RadHeat

232ThPowerCrust

Description

Total Power Generated by 232Th

Custom unit

TW

Modules

RadHeat

232ThPowerMan

Description

Total Power Generated by 232Th

Custom unit

TW

Modules

RadHeat

232ThPowerTotal

Description

Total Power from Decay of 232Th Atoms

Custom unit

TW

Modules

RadHeat

232ThTimescale

Description

Timescale for 232Th Power Generation

Custom unit

Gyr

Modules

RadHeat

235UEnFlux

Description

Surface Energy Flux due to 235U

Custom unit

W/m^2

Modules

RadHeat

235UMassCore

Description

Total Core Mass of 235U

Custom unit

Mearth

Modules

RadHeat

235UMassCrust

Description

Total Crust Mass of 235U

Custom unit

Mearth

Modules

RadHeat

235UMassMan

Description

Total Mass of 235U

Custom unit

Mearth

Modules

RadHeat

235UNumCore

Description

Total Core Number of 235U Atoms

Custom unit

Initial Primordial Earth Number

Modules

RadHeat

235UNumCrust

Description

Total Crust Number of 235U Atoms

Custom unit

Initial Primordial Earth Number

Modules

RadHeat

235UNumMan

Description

Total Number of 235U Atoms

Custom unit

Initial Primordial Earth Number

Modules

RadHeat

235UPowerCore

Description

Total Core Power Generated by 235U

Custom unit

TW

Modules

RadHeat

235UPowerCrust

Description

Total Crust Power Generated by 235U

Custom unit

TW

Modules

RadHeat

235UPowerMan

Description

Total Power Generated by 235U

Custom unit

TW

Modules

RadHeat

235UPowerTotal

Description

Total Power from Decay of 235U Atoms

Custom unit

TW

Modules

RadHeat

235UTimescale

Description

Timescale for 235U Power Generation

Custom unit

Gyr

Modules

RadHeat

238UEnFlux

Description

Surface Energy Flux due to 238U

Custom unit

W/m^2

Modules

RadHeat

238UMassCore

Description

Total Core Mass of 238U

Custom unit

Mearth

Modules

RadHeat

238UMassCrust

Description

Total Crust Mass of 238U

Custom unit

Mearth

Modules

RadHeat

238UMassMan

Description

Total Mass of 238U

Custom unit

Mearth

Modules

RadHeat

238UNumCore

Description

Total Core Number of 238U Atoms

Custom unit

Initial Primordial Earth Number

Modules

RadHeat

238UNumCrust

Description

Total Crust Number of 238U Atoms

Custom unit

Initial Primordial Earth Number

Modules

RadHeat

238UNumMan

Description

Total Number of 238U Atoms

Custom unit

Initial Primordial Earth Number

Modules

RadHeat

238UPowerCore

Description

Total Core Power Generated by 238U

Custom unit

TW

Modules

RadHeat

238UPowerCrust

Description

Total Crust Power Generated by 238U

Custom unit

TW

Modules

RadHeat

238UPowerMan

Description

Total Power Generated by 238U

Custom unit

TW

Modules

RadHeat

238UPowerTotal

Description

Total Power from Decay of 238U Atoms

Custom unit

TW

Modules

RadHeat

238UTimescale

Description

Timescale for 238U Power Generation

Custom unit

Gyr

Modules

RadHeat

26AlEnFlux

Description

Surface Energy Flux from 26Al

Custom unit

W/m^2

Modules

RadHeat

26AlMassCore

Description

Mass of core in 26Al

Custom unit

Mearth

Modules

RadHeat

26AlMassMan

Description

Mass of Mantle in 26Al

Custom unit

Mearth

Modules

RadHeat

26AlNumCore

Description

Total Core Number of 26Al Atoms

Custom unit

Initial Primordial Earth Number

Modules

RadHeat

26AlNumMan

Description

Total Number of 26Al Atoms

Custom unit

Initial Primordial Earth Number

Modules

RadHeat

26AlPowerCore

Description

Core Power Generated by 26Al

Custom unit

TW

Modules

RadHeat

26AlPowerMan

Description

Total Power Generated by 26Al

Custom unit

TW

Modules

RadHeat

26AlPowerTotal

Description

Total Power from Decay of 26Al Atoms

Custom unit

TW

Modules

RadHeat

26AlTimescale

Description

Timescale for 26Al Power Generation

Custom unit

Gyr

Modules

RadHeat

40KEnFlux

Description

Surface Energy Flux from 40K

Custom unit

W/m^2

Modules

RadHeat

40KMassCore

Description

Total Core Mass of 40K

Custom unit

Mearth

Modules

RadHeat

40KMassCrust

Description

Total Crust Mass of 40K

Custom unit

Mearth

Modules

RadHeat

40KMassMan

Description

Total Mass of 40K

Custom unit

Mearth

Modules

RadHeat

40KNumCore

Description

Total Core Number of 40K Atoms

Custom unit

Initial Primordial Earth Number

Modules

RadHeat

40KNumCrust

Description

Total Crust Number of 40K Atoms

Custom unit

Initial Primordial Earth Number

Modules

RadHeat

40KNumMan

Description

Total Number of 40K Atoms

Custom unit

Initial Primordial Earth Number

Modules

RadHeat

40KPowerCore

Description

Core Power Generated by 40K

Custom unit

TW

Modules

RadHeat

40KPowerCrust

Description

Crust Power Generated by 40K

Custom unit

TW

Modules

RadHeat

40KPowerMan

Description

Total Power Generated by 40K

Custom unit

TW

Modules

RadHeat

40KPowerTotal

Description

Total Power from Decay of 40K Atoms

Custom unit

TW

Modules

RadHeat

40KTimescale

Description

Timescale for 40K Power Generation

Custom unit

Gyr

Modules

RadHeat

Age

Description

System Age

Custom unit

Gyr

Modules


AlbedoGlobal

Description

Global mean bond albedo from POISE

Modules

POISE

AlbedoLat

Description

Surface albedo by latitude.

Modules

POISE

AngMX

Description

X component of angular mom vector

Modules

GalHabit

AngMY

Description

Y component of angular mom vector

Modules

GalHabit

AngMZ

Description

Z component of angular mom vector

Modules

GalHabit

AnnInsol

Description

Annual insolation by latitude.

Custom unit

W/m^2

Modules

POISE

AreaIceCov

Description

Fractional area ice covered

Custom unit


Modules

POISE

AtmGasConst

Description

Atmospheric gas constant

Custom unit

J / K kg

Modules

AtmEsc

AtmXAbsEffH2O

Description

XUV atmospheric escape efficiency for H2O

Modules

AtmEsc

BedrockH

Description

height/depth of bedrock

Custom unit

m

Modules

POISE

BinPriPhi

Description

Binary primary star azimuthal angle in BINARY

Custom unit

Deg

Modules

BINARY

BinPriR

Description

Radial position of binary primary star for module BINARY

Custom unit

AU

Modules

BINARY

BinSecPhi

Description

Binary secondary star azimuthal angle in BINARY

Custom unit

Deg

Modules

BINARY

BinSecR

Description

Radial position of binary secondary star for module BINARY

Custom unit

AU

Modules

BINARY

BLLMan

Description

Boundary Layer Thickness Lower Mantle

Custom unit

km

Modules

ThermInt

BLUMan

Description

Boundary Layer Thickness Upper Mantle

Custom unit

km

Modules

ThermInt

BodyDeccDt

Description

In EqTide, the total change in eccentricity is due to towrues
on both bodies. This output parameter returns the change due to
current body.

Custom unit

/Gyr

Modules

EqTide

BodyDsemiDtEqtide

Description

In EqTide, the total change in semi-major axis is due to
dissipation in both bodies. This output parameter returns the
change due to current body.

Custom unit

AU/Gyr

Modules

EqTide

BodyType

Description

Type of Body (0 == planet)

Modules


BondiRadius

Description

Bondi Radius

Custom unit

Rearth

Modules

AtmEsc

CassiniOne

Description

First Cassini parameter (misalignment of Cassini state vectors)

Modules

DistRot

CassiniTwo

Description

Second Cassini parameter (alignment of Cassini state vectors

Modules

DistRot

CBPPhi

Description

CBP Orbital Azimuthal Angle in BINARY

Custom unit

Deg

Modules

BINARY

ChiIC

Description

Light Element Concentration in Inner Core

Custom unit

nd

Modules

ThermInt

ChiOC

Description

Light Element Concentration in Outer Core

Custom unit

nd

Modules

ThermInt

CO2FracMelt

Description

CO2 mass fraction in magma ocean

Custom unit

null

Modules

MagmOc

CO2MassMOAtm

Description

CO2 mass in magma ocean and atmosphere

Custom unit

kg

Modules

MagmOc

CO2MassSol

Description

CO2 mass in solidified mantle

Custom unit

kg

Modules

MagmOc

COPP

Description

eccentriciy * sin(longitude of pericenter + precession angle) *
sin(obliquity)

Modules

BINARY DistRot EqTide POISE SpiNBody

CoreBuoyCompo

Description

Core Compositional Buoyancy Flux

Custom unit

m^2/s^3

Modules

ThermInt

CoreBuoyTherm

Description

Core Thermal Buoyancy Flux

Custom unit

m^2/s^3

Modules

ThermInt

CoreBuoyTotal

Description

Core Total Buoyancy Flux

Custom unit

m^2/s^3

Modules

ThermInt

CriticalSemiMajorAxis

Description

For a circumbinary planet, semi-major axes below a critical
value result in unstable orbits. This output parameter prints
the instantaneous value of that critical distance.

Custom unit

AU

Modules

BINARY EqTide STELLAR

D232ThNumManDt

Description

Time Rate of Change of the number of 232Th nuclei

Custom unit

/Gyr

Modules

RadHeat

D235UNumManDt

Description

Time Rate of Change of the Number of 235U Nuclei

Custom unit

/Gyr

Modules

RadHeat

D238UNumManDt

Description

Time Rate of Change of the Number of 238U Nuclei

Custom unit

/Gyr

Modules

RadHeat

D26AlNumManDt

Description

Time Rate of Change of 26Al Heat Generation

Custom unit

/Gyr

Modules

RadHeat

D26AlPowerDt

Description

Time Rate of Change of 26Al Power Generation

Custom unit

TW/Gyr

Modules

RadHeat

D40KNumManDt

Description

Time Rate of Change of 40K Heat Generation

Custom unit

/Gyr

Modules

RadHeat

D40KPowerDt

Description

Time Rate of Change of 40K Power Generation

Custom unit

TW/Gyr

Modules

RadHeat

DeccDtEqtide

Description

Total de/dt from EqTide

Custom unit

/Gyr

Modules

EqTide

DeltaTime

Description

Average Timestep Over Last Output Interval

Custom unit

years

Modules


Density

Description

Average Density

Custom unit

EarthDensity

Modules


DEnvMassDt

Description

Envelope Mass Loss Rate

Custom unit

Mearth/Myr

Modules

AtmEsc

DepthMeltMan

Description

Depth to base of Upper Mantle Melt Region

Custom unit

m

Modules

ThermInt

DIceMassDt

Description

derivative of mass of ice sheets/area by latitude

Custom unit

kg/m^2/s

Modules

POISE

DivFlux

Description

Divergence of flux (flow into adjacent cells) by latitude

Custom unit

W/m^2

Modules

POISE

DMeanMotionDtEqtide

Description

Total dMeanMotion/dt from EqTide

Modules

EqTide

DOblDtEqtide

Description

Time Rate of Change of Obliquity in EqTide

Custom unit

deg/Gyr

Modules

EqTide

DOrbPerDtEqtide

Description

Total dOrbPer/dt from EqTide

Custom unit

days/Gyr

Modules

EqTide

DRICDTCMB

Description

d(R_ic)/d(T_cmb)

Custom unit

m/K

Modules

ThermInt

DRotPerDtEqtide

Description

Time Rate of Change of Rotation Period in EqTide

Custom unit

days/Myr

Modules

EqTide

DRotPerDtStellar

Description

Time Rate of Change of Rotation Period in STELLAR

Custom unit

days/Myr

Modules

STELLAR

DRotRateDtEqtide

Description

Time Rate of Change of Rotational Frequency in EqTide

Modules

EqTide

DsemiDtEqtide

Description

Total da/dt from EqTide

Custom unit

AU/Gyr

Modules

EqTide

DTChi

Description

Core Liquidus Depression

Custom unit

K

Modules

ThermInt

DynEllip

Description

dynamical ellipticity of planet

Modules

DistRot

Eccentricity

Description

Orbital Eccentricity

Modules

BINARY DistOrb EqTide GalHabit POISE SpiNBody

EccTimeEqtide

Description

Timescale for Eccentricity Evolution (e/[de/dt]) in EqTide

Custom unit

years

Modules

EqTide

EccX

Description

X component of ecc vector

Modules

GalHabit

EccY

Description

Y component of ecc vector

Modules

GalHabit

EccZ

Description

Z component of ecc vector

Modules

GalHabit

EnergyResL

Description

Energy residual on land

Custom unit

W/m^2

Modules

POISE

EnergyResW

Description

Energy residual over water

Custom unit

W/m^2

Modules

POISE

EnvelopeMass

Description

Envelope mass

Custom unit

Mearth

Modules

AtmEsc

EnvTidalQ

Description

Envelope Tidal Q

Modules

EqTide

EqRotPer

Description

Equilibrium Rotation Period

Custom unit

days

Modules

EqTide

EqRotPerCont

Description

CPL2 Continuous Equilibrium Rotation Period

Custom unit

days

Modules

EqTide

EqRotPerDiscrete

Description

CPL2 Discrete Equilibrium Spin Period

Custom unit

days

Modules

EqTide

EqRotRateCont

Description

CPL2 Continuous Equilibrium Spin Rate

Modules

EqTide

EqRotRateDiscrete

Description

CPL2 Discrete Equilibrium Spin Rate

Custom unit

/day

Modules

EqTide

EqTidePower

Description

Equilibrium Power from Tides

Custom unit

/day

Modules

EqTide

EruptEff

Description

Mantle Melt Eruption Efficiency

Custom unit

nd

Modules

ThermInt

EtaO

Description

Oxygen eta parameter (Luger and Barnes 2015)

Modules

AtmEsc

FlareEnergy1

Description

First value of flare energy range

Custom unit

ergs

Modules

FLARE

FlareEnergy2

Description

Second value of flare energy range

Custom unit

ergs

Modules

FLARE

FlareEnergy3

Description

Third value of flare energy range

Custom unit

ergs

Modules

FLARE

FlareEnergy4

Description

Fourth value of flare energy range

Custom unit

ergs

Modules

FLARE

FlareEnergyMax

Description

Maximum flare energy

Custom unit

ergs

Modules

FLARE

FlareEnergyMid

Description

Middle flare energy in the range of energy values

Custom unit

ergs

Modules

FLARE

FlareEnergyMin

Description

Minimum flare energy

Custom unit

ergs

Modules

FLARE

FlareFreq1

Description

First value of flare frequency range

Custom unit

/day

Modules

FLARE

FlareFreq2

Description

Second value of flare frequency range

Custom unit

/day

Modules

FLARE

FlareFreq3

Description

Third value of flare frequency range

Custom unit

/day

Modules

FLARE

FlareFreq4

Description

Fourth value of flare frequency range

Custom unit

/day

Modules

FLARE

FlareFreqMax

Description

Frequency of the flares with the highest energy

Custom unit

/day

Modules

FLARE

FlareFreqMid

Description

Frequency of the flares with the middle energy in the energy
range

Custom unit

/day

Modules

FLARE

FlareFreqMin

Description

Frequency of the flares with the lowest energy

Custom unit

/day

Modules

FLARE

FluxIn

Description

Incoming flux by latitude

Custom unit

W/m^2

Modules

POISE

FluxInGlobal

Description

Global mean flux in (insol*(1-albedo)) from POISE

Custom unit

W/m^2

Modules

POISE

FluxMerid

Description

Total meridional (northward) flux by latitude

Custom unit

PW

Modules

POISE

FluxOut

Description

Outgoing (spaceward) flux by latitude

Custom unit

W/m^2

Modules

POISE

FluxOutGlobal

Description

Global mean flux out from POISE

Custom unit

W/m^2

Modules

POISE

FMeltLMan

Description

Melt Fraction Lower Mantle

Custom unit

nd

Modules

ThermInt

FMeltUMan

Description

Melt Fraction Upper Mantle

Custom unit

nd

Modules

ThermInt

FracFe2O3Man

Description

Fe2O3 mass fraction in magma ocean

Custom unit

null

Modules

MagmOc

FVelDisp

Description

Scaling factor for velocity dispersion

Modules

GalHabit

FXUV

Description

XUV flux

Custom unit

W/m^2

Modules

AtmEsc

GammaOrb

Description

Gamma_Orbital

Modules

EqTide

GammaRot

Description

Gamma_Rotation

Modules

EqTide

GravICB

Description

ICB Gravity

Custom unit

m/s^2

Modules

ThermInt

HEscapeRegime

Description

Integer flag for H envelope escape regime

Modules

AtmEsc

HflowCMB

Description

Heat Flow Core-Mantle Boundary

Custom unit

TW

Modules

ThermInt

HflowLatentIC

Description

Latent Heat Release at ICB

Custom unit

TW

Modules

ThermInt

HflowLatentMan

Description

Latent Heat Flow Mantle

Custom unit

TW

Modules

ThermInt

HflowLMan

Description

Heat Flow Lower Mantle

Custom unit

TW

Modules

ThermInt

HflowMeltMan

Description

Melt Heat Flow Mantle

Custom unit

TW

Modules

ThermInt

HflowSecMan

Description

Mantle Secular Heat Flow

Custom unit

TW

Modules

ThermInt

HflowSurf

Description

Heat Flow Surface Total

Custom unit

TW

Modules

ThermInt

HflowUMan

Description

Heat Flow Upper Mantle

Custom unit

TW

Modules

ThermInt

HfluxCMB

Description

Heat Flux Core-Mantle Boundary

Custom unit

W/m^2

Modules

ThermInt

HfluxCMBAd

Description

Adiabatic Heat Flux Core-Mantle Boundary

Custom unit

W/m^2

Modules

ThermInt

HfluxCMBConv

Description

Super-Adiabatic (convective) Heat Flux Core-Mantle Boundary

Custom unit

W/m^2

Modules

ThermInt

HfluxLMan

Description

Heat Flux Lower Mantle

Custom unit

W/m^2

Modules

ThermInt

HfluxUMan

Description

Heat Flux Upper Mantle

Custom unit

W/m^2

Modules

ThermInt

HydrogenMassSpace

Description

H mass lost to space

Custom unit

kg

Modules

MagmOc

HZInnerEdge

Description

Inner Edge of the Habitable Zone (Runaway Greenhouse)

Custom unit

AU

Modules

MagmOc

HZLimEarlyMars

Description

Early Mars habitable zone limit from Kopparapu et al. (2013).
The value is determined by the total luminosity and average
effective temperature of all interior bodies.

Custom unit

AU

Modules


HZLimitDryRunaway

Description

Minimum distance where a “dry” planet can be habitable,
following Abe et al. (2011), or ~415 W/m^2 net instellation.

Custom unit

AU

Modules


HZLimMaxGreenhouse

Description

Maximum greenhouse habitable zone limit from Kopparapu et al.
(2013). The value is determined by the total luminosity and
average effective temperature of all interior bodies.

Custom unit

AU

Modules


HZLimMoistGreenhouse

Description

Moist greenhouse habitable zone limit from Kopparapu et al.
(2013). The value is determined by the total luminosity and
average effective temperature of all interior bodies.

Custom unit

AU

Modules


HZLimRecVenus

Description

Recent Venus habitable zone limit from Kopparapu et al. (2013).
The value is determined by the total luminosity and average
effective temperature of all interior bodies.

Custom unit

AU

Modules


HZLimRunaway

Description

Runaway greenhouse habitable zone limit from Kopparapu et al.
(2013). The value is determined by the total luminosity and
average effective temperature of all interior bodies.

Custom unit

AU

Modules


IceAblate

Description

ice decay per orbit (ablation only)

Custom unit

m/orbit

Modules

POISE

IceAccum

Description

Ice growth per orbit (accumulation only)

Custom unit

m/orbit

Modules

POISE

IceBeltLand

Description

Does the planet have a land ice belt?

Modules

POISE

IceBeltNorthLatLand

Description

If a land ice belt is present, return the latitude of its
northern edge. If not present, return 0. Note that some ice
belts may in fact have a northern edge at the equator.

Modules

POISE

IceBeltNorthLatSea

Description

If a sea ice belt is present, return the latitude of its
northern edge. If not present, return 0. Note that some ice
belts may in fact have a northern edge at the equator.

Modules

POISE

IceBeltSea

Description

Does the planet have a sea ice belt?

Modules

POISE

IceBeltSouthLatLand

Description

If a land ice belt is present, return the latitude of its
southern edge. If not present, return 0. Note that some ice
belts may in fact have a southern edge at the equator.

Modules

POISE

IceBeltSouthLatSea

Description

If a sea ice belt is present, return the latitude of its
southern edge. If not present, return 0. Note that some ice
belts may in fact have a southern edge at the equator.

Modules

POISE

IceCapNorthLand

Description

Does the planet have a northern polar ice cap on land?

Modules

POISE

IceCapNorthLatLand

Description

If a northern land ice cap is present, return the latitude of
its southern edge. If not present, return +100 degrees.

Modules

POISE

IceCapNorthLatSea

Description

If a northern sea ice cap is present, return the latitude of
its southern edge. If not present, return +100 degrees.

Modules

POISE

IceCapNorthSea

Description

Does the planet have a northern polar sea ice cap

Modules

POISE

IceCapSouthLand

Description

Does the planet have a southern polar ice cap on land?

Modules

POISE

IceCapSouthLatLand

Description

If a southern land ice cap is present, return the latitude of
its northern edge. If not present, return -100 degrees.

Modules

POISE

IceCapSouthLatSea

Description

If a southern sea ice cap is present, return the latitude of
its northern edge. If not present, return -100 degrees.

Modules

POISE

IceCapSouthSea

Description

Does the planet have a southern polar sea ice cap?

Modules

POISE

IceFlow

Description

flow of ice sheets/area by latitude

Custom unit

m/s

Modules

POISE

IceFree

Description

Is the planet free of sea and land ice?

Modules

POISE

IceHeight

Description

Height of ice sheets

Custom unit

m

Modules

POISE

IceMass

Description

Mass of ice sheets/area by latitude

Custom unit

kg/m^2

Modules

POISE

ImK2

Description

Im(k_2)

Modules

EqTide ThermInt

Imk2Man

Description

Imaginary Love Number k2 Mantle

Custom unit

nd

Modules

EqTide ThermInt

Instellation

Description

Total instellation on a planet from all luminous, interior
bodies.

Custom unit

W/m^2

Modules


K2

Description

k_2

Modules

EqTide ThermInt

K2Man

Description

Real Love Number k2 Mantle

Custom unit

nd

Modules

EqTide ThermInt

KTide

Description

Mass loss enhancement due to stellar gravity

Modules

AtmEsc

Latitude

Description

Latitude.

Custom unit

deg

Modules

POISE

LL13K0

Description

Leung+Lee 2013 Radial epicyclic frequency

Custom unit

1/year

Modules

BINARY

LL13N0

Description

CBP Mean Motion in BINARY

Custom unit

1/year

Modules

BINARY

LL13V0

Description

CBP vertical epicyclic frequency in BINARY

Custom unit

1/year

Modules

BINARY

LockTime

Description

Time when body tidally locked.

Custom unit

Myr

Modules

EqTide

LostAngMom

Description

Lost Angular Momentum due to Magnetic Braking

Custom unit

kg*m^2/s

Modules

STELLAR

Luminosity

Description

Luminosity

Custom unit

LSUN

Modules

STELLAR

LXUVFlare

Description

XUV Luminosity from flares

Custom unit

LSUN

Modules

FLARE

LXUVFrac

Description

Fraction of luminosity in XUV

Modules

STELLAR

LXUVStellar

Description

Base X-ray/XUV Luminosity

Custom unit

LSUN

Modules

STELLAR

LXUVTot

Description

Total XUV Luminosity

Custom unit

LSUN

Modules

AtmEsc FLARE STELLAR

MagMom

Description

Core Magnetic Moment

Custom unit

EMAGMOM

Modules

ThermInt

MagPauseRad

Description

Magnetopause Radius

Custom unit

EMAGPAUSERAD

Modules

ThermInt

Mass

Description

Mass

Custom unit

Mearth

Modules


MassChiIC

Description

Mass of Chi in Inner Core

Custom unit

EMASSIC_CHI

Modules

ThermInt

MassChiOC

Description

Mass of Chi in Outer Core

Custom unit

EMASSOC_CHI

Modules

ThermInt

MassIC

Description

Mass of Inner Core

Custom unit

EMASSIC

Modules

ThermInt

MassOC

Description

Mass of Outer Core

Custom unit

EMASSOC

Modules

ThermInt

MeanAnomaly

Description

Orbital Mean Anomaly

Custom unit

Deg

Modules

SpiNBody

MeanLongitude

Description

Orbital Mean Longitude

Custom unit

Deg

Modules

DistOrb SpiNBody

MeanMotion

Description

Orbital Mean Motion

Custom unit

/day

Modules

AtmEsc BINARY DistOrb EqTide POISE SpiNBody

MeltfactorLMan

Description

Melt Factor Lower Mantle

Custom unit

nd

Modules

ThermInt

MeltfactorUMan

Description

Melt Factor Upper Mantle

Custom unit

nd

Modules

ThermInt

MeltFraction

Description

Melt fraction of the magma ocean

Custom unit

null

Modules

MagmOc

MeltMassFluxMan

Description

Mantle Melt Mass Flux

Custom unit

kg/s

Modules

ThermInt

NBadImpulse

Description

cumulative number of questionable encounters

Modules

GalHabit

NEncounters

Description

cumulative number of stellar encounters

Modules

GalHabit

NetFluxAtmo

Description

Atmospheric Net Flux (OLR-ASR)

Custom unit

W/m^2

Modules

MagmOc

Obliquity

Description

Obliquity

Custom unit

deg

Modules

DistRot EqTide POISE

OblTimeDistRot

Description

Obliquity timescale in DistRot

Custom unit

year

Modules

DistRot

OblTimeEqtide

Description

Timescale for Obliquity Evolution in EqTide

Custom unit

years

Modules

EqTide

OceanTidalQ

Description

Ocean Tidal Q

Modules

EqTide

OrbAngMom

Description

Orbital Angular Momentum

Custom unit

kg*m^2/s

Modules

BINARY DistOrb EqTide SpiNBody

OrbPeriod

Description

Orbital Period

Custom unit

days

Modules

AtmEsc BINARY DistOrb EqTide POISE SpiNBody

OxygenMantleMass

Description

Mass of oxygen in mantle

Custom unit

bars

Modules

AtmEsc

OxygenMass

Description

Oxygen mass in the atmosphere

Custom unit

bars

Modules

AtmEsc

OxygenMassMOAtm

Description

Oxygen mass in magma ocean and atmosphere

Custom unit

kg

Modules

MagmOc

OxygenMassSol

Description

Oxygen Mass in solidified mantle

Custom unit

kg

Modules

MagmOc

OxygenMassSpace

Description

O atoms lost to space

Custom unit

kg

Modules

MagmOc

PeakInsol

Description

Peak insolation by latitude.

Custom unit

W/m^2

Modules

POISE

PeriQ

Description

Pericenter distance

Custom unit

AU

Modules

GalHabit

PlanckBAvg

Description

Annually averaged Planck B coeff

Modules

POISE

PlanetRadius

Description

Planet radius

Custom unit

Rearth

Modules

AtmEsc

PotTemp

Description

Potential temperature magma ocean

Custom unit

K

Modules

MagmOc

PowerEqtide

Description

Internal Power due to Tides in EqTide

Custom unit

TW

Modules

EqTide

PowerGravIC

Description

Gravitational Power Release at ICB

Custom unit

TW

Modules

ThermInt

PrecA

Description

The precession angle is orthogonal to the obliquity and is
measured from the vernal equinox. This angle is a “dog-leg”
angle as shown in Fig. 30 of Barnes et al. (2020).

Custom unit

Deg

Modules

DistRot EqTide POISE

PrecATimeDistRot

Description

Axial precession timescale in DistRot

Custom unit

year

Modules

DistRot

PrecFNat

Description

natural precession freq of planet

Modules

DistRot

PressCO2Atm

Description

CO2 pressure in atmosphere

Custom unit

bar

Modules

MagmOc

PressOxygenAtm

Description

Oxygen pressure in atmosphere

Custom unit

bar

Modules

MagmOc

PresSurf

Description

Surface pressure due to atmosphere

Custom unit

GPa

Modules

AtmEsc

PressWaterAtm

Description

Water pressure in atmosphere

Custom unit

bar

Modules

MagmOc

PresSWind

Description

Stellar wind pressure at body

Custom unit

EPRESSWIND

Modules

ThermInt

PresXUV

Description

Pressure at base of thermosphere

Custom unit

Pa

Modules

AtmEsc

RadGyra

Description

Radius of Gyration/Moment of Inertia Constant

Modules


RadioPower

Description

Power from radiogenic heating in the mantle

Custom unit

TW

Modules

MagmOc

Radius

Description

Radius

Custom unit

Rearth

Modules


RadPowerCore

Description

Total Power Generated by Radiogenic Nuclides in Core

Custom unit

TW

Modules

RadHeat

RadPowerCrust

Description

Total Power Generated by Radiogenic Nuclides in the crust

Custom unit

TW

Modules

RadHeat

RadPowerMan

Description

Total Power Generated by Radiogenic Nuclides in the mantle

Custom unit

TW

Modules

RadHeat

RadPowerTotal

Description

Total Power Generated by Radiogenic Nuclides in Total (M+C)

Custom unit

TW

Modules

RadHeat

RadSolid

Description

Radius to the solid surface

Custom unit

Rearth

Modules

AtmEsc

RadXUV

Description

XUV radius separating hydro. dyn. escape and equilibrium

Custom unit

Rearth

Modules

AtmEsc

RayleighMan

Description

Mantle Rayleigh Number

Custom unit

nd

Modules

ThermInt

RGLimit

Description

Runaway greenhouse semi-major axis

Custom unit

AU

Modules

AtmEsc

RIC

Description

Inner Core Radius

Custom unit

km

Modules

ThermInt

RICDot

Description

Inner Core growth rate

Custom unit

m/s

Modules

ThermInt

RocheRadius

Description

Roche Lobe Radius

Custom unit

Rearth

Modules

AtmEsc

RossbyNumber

Description

Rossby Number

Modules

STELLAR

RotAngMom

Description

Rotational Angular Momentum

Custom unit

kg*m^2/s

Modules

BINARY DistRot EqTide STELLAR

RotPer

Description

Rotational Period

Custom unit

days

Modules

DistRot EqTide STELLAR

RotRate

Description

Rotational Frequency

Custom unit

/day

Modules

DistRot EqTide STELLAR

RotTimeEqtide

Description

Timescale for Rotational Evolution in EqTide

Custom unit

years

Modules

EqTide

RotVel

Description

Rotational Velocity

Custom unit

km/s

Modules

DistOrb EqTide STELLAR

RRCriticalFlux

Description

Critical XUV Flux that separates RR and energy-limited escape

Custom unit

W/m^2

Modules

AtmEsc

ScaleHeight

Description

Scale height in Lehmer & Catling (2016) model

Custom unit

km

Modules

AtmEsc

SemiMajorAxis

Description

Semi-major Axis

Custom unit

AU

Modules

AtmEsc BINARY DistOrb EqTide GalHabit POISE SpiNBody

SemiTimeEqtide

Description

Timescale for Semi-major Axis Evolution (a/[da/dt]) in EqTide

Custom unit

years

Modules

EqTide

ShmodUMan

Description

Shear Modulus Upper Mantle

Custom unit

Pa

Modules

ThermInt

SignTJumpLMan

Description

Sign of Lower Mantle Temperature Jump

Custom unit

K

Modules

ThermInt

SignTJumpUMan

Description

Sign of Upper Mantle Temperature Jump

Custom unit

K

Modules

ThermInt

SkipSeas

Description

Is Seasonal model skipped due to RGH or snowball?

Modules

POISE

Snowball

Description

Is the planet in a snowball state?

Modules

POISE

SnowballLand

Description

Is all land covered in ice?

Modules

POISE

SnowballSea

Description

Is all sea covered in ice?

Modules

POISE

SolidRadius

Description

Solidification radius of the mantle

Custom unit

Rearth

Modules

MagmOc

SurfEnFluxEqtide

Description

Surface Energy Flux due to Tides in EqTide

Custom unit

W/m^2

Modules

EqTide

SurfEnFluxRadTotal

Description

Total Surface Heat Flux Generated by Radiogenic Nuclides

Custom unit

W/m^2

Modules

RadHeat

SurfEnFluxTotal

Description

Total Surface Energy Flux

Custom unit

W/m^2

Modules

EqTide RadHeat ThermInt

SurfTemp

Description

Surface temperature magma ocean

Custom unit

K

Modules

MagmOc

SurfWaterMass

Description

Surface water mass

Custom unit

TO

Modules

AtmEsc

TCMB

Description

CMB Temperature

Custom unit

K

Modules

ThermInt

TCore

Description

Core Temperature

Custom unit

K

Modules

ThermInt

TDepthMeltMan

Description

Temp at base of Upper Mantle Melt Region

Custom unit

K

Modules

ThermInt

TDotCore

Description

Change in Core Temperature

Custom unit

K/s

Modules

ThermInt

TDotMan

Description

Change in Mantle Temperature

Custom unit

K/s

Modules

ThermInt

Temperature

Description

Effective Temperature

Modules

STELLAR

TempLat

Description

Surface temperature by latitude.

Custom unit

Celsius

Modules

POISE

TempMaxLand

Description

Maximum surface temperature on land

Custom unit

Celsius

Modules

POISE

TempMaxLat

Description

Maximum surface temperature over a year by latitude.

Custom unit

Celsius

Modules

POISE

TempMaxWater

Description

Maximum surface temperature on water

Custom unit

Celsius

Modules

POISE

TempMinLat

Description

Minimum surface temperature over a year by latitude.

Custom unit

Celsius

Modules

POISE

TGlobal

Description

Global mean temperature from POISE

Custom unit

Celsius

Modules

POISE

ThermConductOC

Description

Thermal Conductivity of the Outer Core

Custom unit

W/m/K

Modules

ThermInt

ThermTemp

Description

Isothermal atmospheric temperature

Custom unit

K

Modules

AtmEsc

TidalPower

Description

Power from tidal heating in the mantle

Custom unit

TW

Modules

MagmOc

TidalQ

Description

Tidal Q

Modules

EqTide ThermInt

TidalRadius

Description

In EqTide, the value of the radius that enters the tidal
evolution equations. This parameter is distinct from others,
such as the radius of XUV absorption or a pressure level in an
atmosphere.

Custom unit

Rearth

Modules

EqTide

TidalTau

Description

Tidal Time Lag

Custom unit

sec

Modules

EqTide

TideLock

Description

Tidally Locked?

Modules

EqTide

Time

Description

Simulation Time

Custom unit

Gyr

Modules


TJumpLMan

Description

Lower Mantle Temperature Jump

Custom unit

K

Modules

ThermInt

TJumpMeltMan

Description

Temp Jump across Upper Mantle Melt Region

Custom unit

K

Modules

ThermInt

TJumpUMan

Description

Upper Mantle Temperature Jump

Custom unit

K

Modules

ThermInt

TLMan

Description

Lower Mantle Temperature

Custom unit

K

Modules

ThermInt

TMan

Description

Mantle Temperature

Custom unit

K

Modules

ThermInt

TotAngMom

Description

Total Angular Momentum

Custom unit

kg*m^2/s

Modules

BINARY DistOrb DistRot EqTide SpiNBody STELLAR

TotEnergy

Description

Total System Energy

Custom unit

ergs

Modules

DistOrb EqTide SpiNBody STELLAR

TotIceBalance

Description

Global total ice balance in ice sheets (this time step)

Custom unit

kg

Modules

POISE

TotIceFlow

Description

Global total ice flow in ice sheets (should = 0)

Custom unit

kg

Modules

POISE

TotIceMass

Description

Global total ice mass in ice sheets

Custom unit

kg

Modules

POISE

TrefLind

Description

Core Liquidus Lindeman Reference

Custom unit

K

Modules

ThermInt

TUMan

Description

Upper Mantle Temperature

Custom unit

K

Modules

ThermInt

ViscJumpMan

Description

Mantle Viscosity Ratio

Custom unit

nd

Modules

ThermInt

ViscLMan

Description

Lower Mantle Viscosity

Custom unit

m^2/s

Modules

ThermInt

ViscMMan

Description

Average (mid) Mantle Viscosity

Custom unit

m^2/s

Modules

ThermInt

ViscUMan

Description

Upper Mantle Viscosity

Custom unit

m^2/s

Modules

ThermInt

ViscUManArr

Description

Upper Mantle Arrhenius Viscosity

Custom unit

m^2/s

Modules

ThermInt

WaterFracMelt

Description

water mass fraction in magma ocean

Custom unit

null

Modules

MagmOc

WaterMassMOAtm

Description

Water mass in magma ocean and atmosphere

Custom unit

TO

Modules

MagmOc

WaterMassSol

Description

Water mass in solidified mantle

Custom unit

TO

Modules

MagmOc

XO

Description

Atomic oxygen mixing ratio in upper atmosphere

Modules

AtmEsc

XoblTimeDistRot

Description

Xobl timescale in DistRot

Custom unit

year

Modules

DistRot

YoblTimeDistRot

Description

Yobl timescale in DistRot

Custom unit

year

Modules

DistRot

ZoblTimeDistRot

Description

Zobl timescale in DistRot

Custom unit

year

Modules

DistRot