RoryBarnes · chandan-kumarsahu · Sep 21, 2022 · Sep 21, 2022
diff --git a/MagmOc_test_runs/1K_and_10K_radioactive_potassium/1_100TO_10K/LP8909_10K_100TO.png b/MagmOc_test_runs/1K_and_10K_radioactive_potassium/1_100TO_10K/LP8909_10K_100TO.png
diff --git a/MagmOc_test_runs/1K_and_10K_radioactive_potassium/1_100TO_10K/LP8909_10K_10TO.png b/MagmOc_test_runs/1K_and_10K_radioactive_potassium/1_100TO_10K/LP8909_10K_10TO.png
diff --git a/MagmOc_test_runs/1K_and_10K_radioactive_potassium/1_100TO_10K/LP8909_10K_1TO.png b/MagmOc_test_runs/1K_and_10K_radioactive_potassium/1_100TO_10K/LP8909_10K_1TO.png
diff --git a/MagmOc_test_runs/1K_and_10K_radioactive_potassium/1_100TO_10K/LP8909_10K_50TO.png b/MagmOc_test_runs/1K_and_10K_radioactive_potassium/1_100TO_10K/LP8909_10K_50TO.png
diff --git a/MagmOc_test_runs/1K_and_10K_radioactive_potassium/1_100TO_10K/LP8909_10K_5TO.png b/MagmOc_test_runs/1K_and_10K_radioactive_potassium/1_100TO_10K/LP8909_10K_5TO.png
diff --git a/MagmOc_test_runs/1K_and_10K_radioactive_potassium/1_100TO_1K/LP8909_100TO.png b/MagmOc_test_runs/1K_and_10K_radioactive_potassium/1_100TO_1K/LP8909_100TO.png
diff --git a/MagmOc_test_runs/1K_and_10K_radioactive_potassium/1_100TO_1K/LP8909_10TO.png b/MagmOc_test_runs/1K_and_10K_radioactive_potassium/1_100TO_1K/LP8909_10TO.png
diff --git a/MagmOc_test_runs/1K_and_10K_radioactive_potassium/1_100TO_1K/LP8909_15TO.png b/MagmOc_test_runs/1K_and_10K_radioactive_potassium/1_100TO_1K/LP8909_15TO.png
diff --git a/MagmOc_test_runs/1K_and_10K_radioactive_potassium/1_100TO_1K/LP8909_1TO.png b/MagmOc_test_runs/1K_and_10K_radioactive_potassium/1_100TO_1K/LP8909_1TO.png
diff --git a/MagmOc_test_runs/1K_and_10K_radioactive_potassium/1_100TO_1K/LP8909_25TO.png b/MagmOc_test_runs/1K_and_10K_radioactive_potassium/1_100TO_1K/LP8909_25TO.png
diff --git a/MagmOc_test_runs/1K_and_10K_radioactive_potassium/1_100TO_1K/LP8909_500TO.png b/MagmOc_test_runs/1K_and_10K_radioactive_potassium/1_100TO_1K/LP8909_500TO.png
diff --git a/MagmOc_test_runs/1K_and_10K_radioactive_potassium/1_100TO_1K/LP8909_50TO.png b/MagmOc_test_runs/1K_and_10K_radioactive_potassium/1_100TO_1K/LP8909_50TO.png
diff --git a/MagmOc_test_runs/1K_and_10K_radioactive_potassium/1_100TO_1K/LP8909_5TO.png b/MagmOc_test_runs/1K_and_10K_radioactive_potassium/1_100TO_1K/LP8909_5TO.png
diff --git a/MagmOc_test_runs/1M_1K_1-25TO/10TO/LP8909.LP8909.forward b/MagmOc_test_runs/1M_1K_1-25TO/10TO/LP8909.LP8909.forward
diff --git a/MagmOc_test_runs/1M_1K_1-25TO/10TO/LP8909.c.forward b/MagmOc_test_runs/1M_1K_1-25TO/10TO/LP8909.c.forward
diff --git a/MagmOc_test_runs/1M_1K_1-25TO/10TO/LP8909.in b/MagmOc_test_runs/1M_1K_1-25TO/10TO/LP8909.in
@@ -0,0 +1,21 @@
+# The host star
+sName	        LP8909                # Body's name
+saModules	stellar, eqtide          # Modules to apply, exact spelling required
+
+# Output
+saOutputOrder Age -Luminosity -LXUVStellar -Radius Temperature
+
+# Physical Parameters
+sStellarModel             baraffe     # stellar model grid
+dMass                     0.118        # stellar mass [msun]
+dAge                      5e6         # age at which stellar grids are initialized [yr] (this is a fixed parameter)
+dSatXUVFrac               0.00120226        # fsat = initial Lxuv / Lbol [dimensionless ratio]
+dSatXUVTime               3e9        # tsat = saturation time [yr]
+dXUVBeta                  1.18       # beta = exponential decay rate [dimensionless]
+
+# EQTIDE Parameters
+dTidalQ         100          # Tidal phase lag
+dK2             0.5          # Love number of degree 2
+dMaxLockDiff    0.1
+sTideModel      p2           # Tidal model, p2=CPL, t8=CTL
+saTidePerts     c         # Body name(s) of tidal perturbers
diff --git a/MagmOc_test_runs/1M_1K_1-25TO/10TO/LP8909.log b/MagmOc_test_runs/1M_1K_1-25TO/10TO/LP8909.log
diff --git a/MagmOc_test_runs/1M_1K_1-25TO/10TO/LP8909_10TO.png b/MagmOc_test_runs/1M_1K_1-25TO/10TO/LP8909_10TO.png
diff --git a/MagmOc_test_runs/1M_1K_1-25TO/10TO/README.rst b/MagmOc_test_runs/1M_1K_1-25TO/10TO/README.rst
@@ -0,0 +1,32 @@
+Magma ocean evolution of LP 890-9 and c
+============================================
+
+=========================   ========================================
+**Date**                    09/14/2022
+**Modules**                 AtmEsc, EqTide, MagmOc, RadHeat, STELLAR
+**Approx. runtime**         xxx minutes
+=========================   ========================================
+
+TODO:
+
+1. Modify sun.in to match LP 890-9's properties
+2. Modify g.in to match planet c
+3. Explore a wider range of water contents
+4. Try different cases with different tidal Q's and radiogenic inventories
+
+To run this example
+-------------------
+
+.. code-block:: bash
+
+    python makeplot.py
+
+
+Expected output
+---------------
+
+.. figure:: Trappist1g_2TO.png
+   :width: 600px
+   :align: center
+
+Caption
diff --git a/MagmOc_test_runs/1M_1K_1-25TO/10TO/b.in b/MagmOc_test_runs/1M_1K_1-25TO/10TO/b.in
@@ -0,0 +1,73 @@
+# Planet b parameters
+sName         b			                     # Body's name
+saModules     magmoc, atmesc, radheat, eqtide       # Modules to apply, exact spelling required
+sColor        13aed5
+
+# Physical Properties
+dMass         -2.3                       # Mass, negative -> Earth masses
+dRadius       -1.320                     # Radius, negative -> Earth radii
+dRotPeriod    -1                           # Rotation Period, negative -> days
+
+# MAGMOC Properties
+dMassFracFeOIni		0.0788		  # Initial mass fraction of FeO in the mantle (default: 0.0788)
+dWaterMassAtm		-2		  # Initial water mass in the system (neg: terrestrial oceans, default: -1)
+dSurfTemp		4000		  # Initial surface temperature (default: 4000K)
+dManMeltDensity		4000		  # Mantle melt density (default: 4000kg/m^3)
+bHaltMantleSolidified 	0		  # Halt when mantle solidified? (default: 0 = no, 1 = yes)
+bHaltMantleMeltFracLow  0     # Halt when melt fraction smaller than 0.4 at surface? (default: 0 = no, 1 = yes)
+bHaltAtmDesiSurfCool    0                 # Halt when atmosphere desiccated? (default: 0 = no, 1 = yes)
+bHaltEnterHabZone       0                 # Halt when planet enters habitable zone? (default: 0 = no, 1 = yes)
+bHaltAllPlanetsSolid  	0                 # Halt when all planets are soldified	(for multiple planet systems)
+bHaltAllPlanetsDesicc 	0                 # Halt when all planets are desiccated (for multiple planet systems)
+sRadioHeatModel		none              # Radiogenic heating model (default: none = RadHeat; schaefer = BSE composition)
+sMagmOcAtmModel 	grey		  # Atmospheric net flux model (default: grey; petit only for GJ1132b-H2O)
+
+# ATMESC Properties
+dXFrac                    1.0             # X-Ray/XUV absorption radius (fraction of planet radius)
+dSurfWaterMass            -1.0            # Initial surface water (Earth oceans)
+dEnvelopeMass             0               # Initial envelope mass (Earth masses)
+bHaltSurfaceDesiccated    0               # Halt when dry?
+bHaltEnvelopeGone         0               # Halt when evaporated?
+dMinSurfWaterMass         -1.e-5          # Planet is desiccated when water content drops below this (Earth oceans)
+sWaterLossModel           lbexact
+sPlanetRadiusModel        none
+bInstantO2Sink            0
+sAtmXAbsEffH2OModel       none
+dAtmXAbsEffH2O 		  0.3
+
+# RADHEAT Parameters
+# Scale with the planetary mass
+# *Num* are in numbers of atoms, negative -> Earth vals
+### 40K
+d40KPowerMan    -1
+d40KPowerCore   -0
+d40KPowerCrust  -0.0
+
+### 232Th
+d232ThPowerMan      -1
+d232ThPowerCore     -0
+d232ThPowerCrust    -0
+
+### 235U
+d235UPowerMan      -1
+d235UPowerCore     -0
+d235UPowerCrust    -0
+
+### 238U
+d238UPowerMan      -1
+d238UPowerCore     -0
+d238UPowerCrust    -0
+
+# EQTIDE Parameters
+dTidalQ         100		   # Tidal phase lag
+dK2             0.5	           # Love number of degree 2
+dMaxLockDiff    0.01
+saTidePerts     LP8909          # Body name(s) of tidal perturbers
+
+# Orbital Properties
+dSemi         -0.01875              # Semi-major axis, negative -> AU
+dEcc          0.00               # Eccentricity
+
+# Output
+saOutputOrder Time -PotTemp -SurfTemp -SolidRadius -WaterMassMOAtm -WaterMassSol -OxygenMassMOAtm -OxygenMassSol -PressWaterAtm -PressOxygenAtm -HydrogenMassSpace -OxygenMassSpace FracFe2O3Man NetFluxAtmo WaterFracMelt -RadioPower -TidalPower $
+-SemiMajorAxis -HZInnerEdge
diff --git a/MagmOc_test_runs/1M_1K_1-25TO/10TO/c.in b/MagmOc_test_runs/1M_1K_1-25TO/10TO/c.in
@@ -0,0 +1,73 @@
+# Planet c parameters
+sName         c			                     # Body's name
+saModules     magmoc, atmesc, radheat, eqtide       # Modules to apply, exact spelling required
+sColor        7a13d5
+
+# Physical Properties
+dMass         -1                       # Mass, negative -> Earth masses
+dRadius       -1.367                      # Radius, negative -> Earth radii
+dRotPeriod    -8.457                      # Rotation Period, negative -> days
+
+# MAGMOC Properties
+dMassFracFeOIni		0.0788		  # Initial mass fraction of FeO in the mantle (default: 0.0788)
+dWaterMassAtm		-10		  # Initial water mass in the system (neg: terrestrial oceans, default: -1)
+dSurfTemp		4000		  # Initial surface temperature (default: 4000K)
+dManMeltDensity		4000		  # Mantle melt density (default: 4000kg/m^3)
+bHaltMantleSolidified 	0		  # Halt when mantle solidified? (default: 0 = no, 1 = yes)
+bHaltMantleMeltFracLow  0     # Halt when melt fraction smaller than 0.4 at surface? (default: 0 = no, 1 = yes)
+bHaltAtmDesiSurfCool                     # Halt when atmosphere desiccated? (default: 0 = no, 1 = yes)
+bHaltEnterHabZone       0                 # Halt when planet enters habitable zone? (default: 0 = no, 1 = yes)
+bHaltAllPlanetsSolid  	1                 # Halt when all planets are soldified	(for multiple planet systems)
+bHaltAllPlanetsDesicc 	0                 # Halt when all planets are desiccated (for multiple planet systems)
+sRadioHeatModel		none              # Radiogenic heating model (default: none = RadHeat; schaefer = BSE composition)
+sMagmOcAtmModel 	grey		  # Atmospheric net flux model (default: grey; petit only for GJ1132b-H2O)
+
+# ATMESC Properties
+dXFrac                    1.0             # X-Ray/XUV absorption radius (fraction of planet radius)
+dSurfWaterMass            -1.0            # Initial surface water (Earth oceans)
+dEnvelopeMass             0               # Initial envelope mass (Earth masses)
+bHaltSurfaceDesiccated    0               # Halt when dry?
+bHaltEnvelopeGone         0               # Halt when evaporated?
+dMinSurfWaterMass         -1.e-5          # Planet is desiccated when water content drops below this (Earth oceans)
+sWaterLossModel           lbexact
+sPlanetRadiusModel        none
+bInstantO2Sink            0
+sAtmXAbsEffH2OModel       none
+dAtmXAbsEffH2O 		  0.3
+
+# RADHEAT Parameters
+# Scale with planetary mass
+# *Num* are in numbers of atoms, negative -> Earth vals
+### 40K
+d40KNumMan    -1.0
+d40KNumCore   -0
+d40KNumCrust  -0.0
+
+### 232Th
+d232ThNumMan      -1.0
+d232ThNumCore     -0
+d232ThNumCrust    -0
+
+### 235U
+d235UNumMan      -1.0
+d235UNumCore     -0
+d235UNumCrust    -0
+
+### 238U
+d238UNumMan      -1.0
+d238UNumCore     -0
+d238UNumCrust    -0
+
+# EQTIDE Parameters
+dTidalQ         100		   # Tidal phase lag
+dK2             0.5	           # Love number of degree 2
+dMaxLockDiff    0.01
+saTidePerts     LP8909          # Body name(s) of tidal perturbers
+
+# Orbital Properties
+dSemi         -0.03984              # Semi-major axis, negative -> AU
+dEcc          0.00                  # Eccentricity
+
+# Output
+saOutputOrder Time -PotTemp -SurfTemp -SolidRadius -WaterMassMOAtm -WaterMassSol -OxygenMassMOAtm -OxygenMassSol -PressWaterAtm -PressOxygenAtm -HydrogenMassSpace -OxygenMassSpace FracFe2O3Man NetFluxAtmo WaterFracMelt -RadioPower -TidalPower $
+-SemiMajorAxis -HZInnerEdge