Merge pull request #326 from FormingWorlds/hn/tideineq

Dummy tides inequality, convergence robustness, docs
FormingWorlds · Jan 21, 2025 · c089723 · c089723
2 parents aa0f07b + 0f1ddde
commit c089723
Show file tree

Hide file tree

Showing 40 changed files with 162 additions and 344 deletions.
diff --git a/README.md b/README.md
@@ -56,20 +56,24 @@ If you make use of PROTEUS, please reference the following manuscripts, and stat
        author = {{Nicholls}, Harrison and {Pierrehumbert}, Raymond T. and {Lichtenberg}, Tim and {Soucasse}, Laurent and {Smeets}, Stef},
         title = "{Convective shutdown in the atmospheres of lava worlds}",
       journal = {Monthly Notices of the Royal Astronomical Society},
-         year = 2024,
-       month  = dec,
-         doi  = {10.1093/mnras/stae2772},
-          url = {https://arxiv.org/abs/2412.11987}
+         year = 2025,
+        month = jan,
+       volume = {536},
+       number = {3},
+        pages = {2957-2971},
+          doi = {10.1093/mnras/stae2772},
 }
 
 @ARTICLE{Nicholls2024JGRP,
-       author = {{Nicholls}, Harrison and {Lichtenberg}, Tim and {Bower}, Dan J. and and {Pierrehumbert}, Raymond T.},
-        title = "{Magma ocean evolution at arbitrary redox state}",
+       author = {{Nicholls}, Harrison and {Lichtenberg}, Tim and {Bower}, Dan J. and {Pierrehumbert}, Raymond},
+        title = "{Magma Ocean Evolution at Arbitrary Redox State}",
       journal = {Journal of Geophysical Research (Planets)},
          year = 2024,
-       month  = nov,
-         doi  = {10.1029/2024JE008576},
-          url = {https://arxiv.org/abs/2411.19137}
+        month = dec,
+       volume = {129},
+       number = {12},
+        pages = {2024JE008576},
+          doi = {10.1029/2024JE008576},
 }
 
 @ARTICLE{Lichtenberg2021JGRP,
@@ -83,10 +87,6 @@ If you make use of PROTEUS, please reference the following manuscripts, and stat
           eid = {e06711},
         pages = {e06711},
           doi = {10.1029/2020JE006711},
-archivePrefix = {arXiv},
-       eprint = {2101.10991},
- primaryClass = {astro-ph.EP},
-       adsurl = {https://ui.adsabs.harvard.edu/abs/2021JGRE..12606711L}
 }
 ```
 

diff --git a/input/default.toml b/input/default.toml
@@ -109,7 +109,7 @@ author  = "Harrison Nicholls, Tim Lichtenberg"
 
     [orbit.dummy]
         H_tide  = 1e-11  # Fixed tidal power density [W kg-1]
-        Phi_tide = 0.3    # Maximum local melt fraction to which tides can be applied.
+        Phi_tide = "<0.3"   # Tidal heating applied when inequality locally satisfied
 
 
 # Planetary structure - physics table
@@ -179,6 +179,8 @@ author  = "Harrison Nicholls, Tim Lichtenberg"
     F_initial       = 8.0E4 # Initial heat flux guess [W m-2]
     radiogenic_heat = true  # enable radiogenic heat production
     tidal_heat      = false  # enable tidal heat production
+    rheo_phi_loc    = 0.4    # Centre of rheological transition
+    rheo_phi_wid    = 0.15   # Width of rheological transition
 
     module = "spider"   # Which interior module to use
 

diff --git a/input/demos/aragog.toml b/input/demos/aragog.toml
@@ -109,7 +109,7 @@ author  = "Harrison Nicholls, Tim Lichtenberg"
 
     [orbit.dummy]
         H_tide  = 1e-7  # Fixed tidal power density [W kg-1]
-        Phi_tide = 0.3    # Maximum local melt fraction to which tides can be applied.
+        Phi_tide = "<0.3"   # Tidal heating applied when inequality locally satisfied
 
 
 # Planetary structure - physics table
@@ -179,6 +179,8 @@ author  = "Harrison Nicholls, Tim Lichtenberg"
     F_initial       = 1e6   # Initial heat flux guess [W m-2]
     radiogenic_heat = false  # enable radiogenic heat production
     tidal_heat      = true  # enable tidal heat production
+    rheo_phi_loc    = 0.4    # Centre of rheological transition
+    rheo_phi_wid    = 0.15   # Width of rheological transition
 
     module = "aragog"   # Which interior module to use
 

diff --git a/input/demos/dummy.toml b/input/demos/dummy.toml
@@ -109,7 +109,7 @@ author  = "Harrison Nicholls, Tim Lichtenberg"
 
     [orbit.dummy]
         H_tide  = 1e-9  # Fixed tidal power density [W kg-1]
-        Phi_tide = 0.3    # Maximum local melt fraction to which tides can be applied.
+        Phi_tide = "<0.3"   # Tidal heating applied when inequality locally satisfied
 
 
 # Planetary structure - physics table
@@ -179,6 +179,8 @@ author  = "Harrison Nicholls, Tim Lichtenberg"
     F_initial       = 1e6 # Initial heat flux guess [W m-2]
     radiogenic_heat = false  # enable radiogenic heat production
     tidal_heat      = true  # enable tidal heat production
+    rheo_phi_loc    = 0.4    # Centre of rheological transition
+    rheo_phi_wid    = 0.15   # Width of rheological transition
 
     module = "dummy"   # Which interior module to use
 

diff --git a/input/demos/escape.toml b/input/demos/escape.toml
@@ -109,7 +109,7 @@ author  = "Harrison Nicholls, Tim Lichtenberg"
 
     [orbit.dummy]
         H_tide  = 1e-11  # Fixed tidal power density [W kg-1]
-        Phi_tide = 0.3    # Maximum local melt fraction to which tides can be applied.
+        Phi_tide = "<0.3"   # Tidal heating applied when inequality locally satisfied
 
 
 # Planetary structure - physics table
@@ -179,6 +179,8 @@ author  = "Harrison Nicholls, Tim Lichtenberg"
     F_initial       = 8.0E4 # Initial heat flux guess [W m-2]
     radiogenic_heat = true  # enable radiogenic heat production
     tidal_heat      = false  # enable tidal heat production
+    rheo_phi_loc    = 0.4    # Centre of rheological transition
+    rheo_phi_wid    = 0.15   # Width of rheological transition
 
     module = "spider"   # Which interior module to use
 

diff --git a/input/demos/spider_tides.toml b/input/demos/spider_tides.toml
diff --git a/input/planets/hd63433d.toml b/input/planets/hd63433d.toml
@@ -109,7 +109,7 @@ author  = "Harrison Nicholls, Tim Lichtenberg"
 
     [orbit.dummy]
         H_tide  = 1e-11  # Fixed tidal power density [W kg-1]
-        Phi_tide = 0.3    # Maximum local melt fraction to which tides can be applied.
+        Phi_tide = "<0.3"   # Tidal heating applied when inequality locally satisfied
 
 
 # Planetary structure - physics table
@@ -179,6 +179,8 @@ author  = "Harrison Nicholls, Tim Lichtenberg"
     F_initial       = 1e5   # Initial heat flux guess [W m-2]
     radiogenic_heat = true  # enable radiogenic heat production
     tidal_heat      = false  # enable tidal heat production
+    rheo_phi_loc    = 0.4    # Centre of rheological transition
+    rheo_phi_wid    = 0.15   # Width of rheological transition
 
     module = "spider"   # Which interior module to use