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	Int
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

dLuminosityAmplitude
Description	Amplitude of luminosity oscillation for SINEWAVE stellar model
Type	Bool
Modules	STELLAR
Files	Any
Default value	0.001

dLuminosityPeriod
Description	Period of luminosity oscillation for SINEWAVE stellar model
Type	Bool
Modules	STELLAR
Files	Any
Default value	0.001

dLuminosityPhase
Description	Phase of luminosity oscillation at age=0 for SINEWAVE stellar model
Type	Bool
Modules	STELLAR
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

sFlareBandPass
Description	Currently suppressed due to execution errors.
Type	String
Modules	ALL
Files	Any
Default value	KEPLER
Allowed values	KEPLER UV GOES SXR BOLOMETRIC TESSUV

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 SINEWAVE 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, TIAN, and LS16.
Type	String
Modules	ALL
Files	Any
Default value	LBEXACT
Allowed values	LB15 LBEXACT TIAN LS16

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

CrossoverMass
Description	Crossover mass, defines oxygen is escaping if mc > mo
Custom unit	amu
Modules	AtmEsc

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

FXUVCRITDRAG
Description	Critical Drag 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

HDiffFlux
Description	The Hydrogen Diffusion Flux (m^-2 s^-1)
Modules	AtmEsc

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

HREFFLUX
Description	The Hydrogen Reference Flux (m^-2 s^-1)
Modules	AtmEsc

HRefODragMod
Description	Multiply by H ref flux to get flux with drag
Modules	AtmEsc

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

TsolUMan
Description	Upper Mantle Thermal Boundary Layer Solidus Temperature
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

WaterEscapeRegime
Description	What water escape regime the code is currently in
Modules	AtmEsc

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

XH2O
Description	Water mixing ratio in upper atmosphere
Modules	AtmEsc

XO
Description	Atomic oxygen mixing ratio in upper atmosphere
Modules	AtmEsc

XO2
Description	Molecular oxygen mixing ratio in 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