Store tidal constants

input/l9859d.toml

@@ -25,7 +25,7 @@ author  = "Harrison Nicholls, Tim Lichtenberg"
 [params]
     # output files
     [params.out]
-        path        = "l9859d"
+        path        = "l9859d_notides"
         logging     = "INFO"
         plot_mod    = 5      # Plotting frequency, 0: wait until completion | n: every n iterations
         plot_fmt    = "png"  # Plotting image file format, "png" or "pdf" recommended
@@ -64,7 +64,7 @@ author  = "Harrison Nicholls, Tim Lichtenberg"
 
         # solidification
         [params.stop.solid]
-            enabled  = false
+            enabled  = true
             phi_crit = 0.005  # non-dim., model will terminate when global melt fraction < phi_crit
 
         # radiative equilibrium
@@ -113,7 +113,7 @@ author  = "Harrison Nicholls, Tim Lichtenberg"
     module  = "dummy"
 
     [orbit.dummy]
-        H_tide  = 1e-9  # Fixed tidal power density [W kg-1]
+        H_tide  = 0.0  # Fixed tidal power density [W kg-1]
 
 # Planetary structure - physics table
 [struct]

src/proteus/utils/constants.py

@@ -19,6 +19,14 @@
 mol             = 6.02214076e+23        # mol definition
 ocean_moles     = 7.68894973907177e+22  # moles of H2 (or H2O) in one present-day Earth ocean
 
+# Earth heat flux, globally averaged [W m-2]
+# https://se.copernicus.org/articles/1/5/2010/
+F_earth         = 47.2e12 / (4 * 3.141 * R_earth * R_earth)
+
+# Io radius [m] and tidal heat flux [W m-2]
+R_io            = 1.82149e6     # https://ssd.jpl.nasa.gov/sats/phys_par/#refs
+F_io            = 2.24          # https://www.nature.com/articles/nature08108
+
 # https://nssdc.gsfc.nasa.gov/planetary/factsheet/
 M_jupiter       = 1.898e27              # kg
 R_jupiter       = 1.42984e8 / 2         # m