From 82e6a709d021cb1b007506c7bfe8e919d672ed4c Mon Sep 17 00:00:00 2001 From: Harrison Nicholls Date: Fri, 22 Nov 2024 14:12:23 +0000 Subject: [PATCH] Store tidal constants --- input/l9859d.toml | 6 +++--- src/proteus/utils/constants.py | 8 ++++++++ 2 files changed, 11 insertions(+), 3 deletions(-) diff --git a/input/l9859d.toml b/input/l9859d.toml index 0ae0d9b9..ffb49b6c 100644 --- a/input/l9859d.toml +++ b/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] diff --git a/src/proteus/utils/constants.py b/src/proteus/utils/constants.py index c9906e0b..d1e752b7 100644 --- a/src/proteus/utils/constants.py +++ b/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