Pass tidal heating to Aragog, SPIDER, Dummy

src/proteus/config/_interior.py

@@ -32,6 +32,7 @@ class Interior:
     grain_size: float = field(validator=gt(0))
     F_initial: float = field(validator=gt(0))
     radiogenic_heat: bool = field(validator=no_radio_if_dummy)
+    tidal_heat: bool
 
     module: str = field(validator=in_(('spider', 'aragog', 'dummy')))
 

src/proteus/interior/aragog.py

@@ -122,7 +122,7 @@ def SetupAragogSolver(config:Config, hf_row:dict):
             gravitational_separation = False,
             mixing = False,
             radionuclides = config.interior.radiogenic_heat,
-            tidal = False,
+            tidal = config.interior.tidal_heat,
             )
 
     initial_condition = _InitialConditionParameters(
@@ -200,7 +200,7 @@ def _append_radnuc(_iso, _cnc):
 
     aragog_solver = Solver(param)
 
-def UpdateAragogSolver(dt:float, hf_row:dict, output_dir:str = None):
+def UpdateAragogSolver(dt:float, hf_row:dict, config:Config, output_dir:str = None):
 
     # Set solver time
     # hf_row["Time"] is in yr so do not need to scale as long as scaling time is secs_per_year
@@ -221,6 +221,14 @@ def UpdateAragogSolver(dt:float, hf_row:dict, output_dir:str = None):
     # Update boundary conditions
     aragog_solver.parameters.boundary_conditions.outer_boundary_value = hf_row["F_atm"]
 
+    # Update tidal heating
+    tidal_value= 0.0
+    if config.interior.tidal_heat:
+        if config.orbit.dummy:
+            tidal_value = config.orbit.dummy.H_tide
+    tidal_value /= aragog_solver.parameters.scalings.power_per_mass
+    aragog_solver.parameters.energy.tidal_value = tidal_value
+
     return
 
 def WriteAragogOutput(output_dir:str, time:float):

src/proteus/interior/dummy.py

@@ -94,7 +94,11 @@ def _calc_phi(tmp:float):
 
     # Subtract tidal contribution to the total heat flux.
     #    This heat energy is generated only in the mantle, not in the core.
-    output["F_tidal"] = config.orbit.dummy.H_tide * output["M_mantle"] / area
+    tidal_flux = 0.0
+    if config.interior.tidal_heat:
+        if config.orbit.dummy:
+            tidal_flux = config.orbit.dummy.H_tide * output["M_mantle"] / area
+    output["F_tidal"] = tidal_flux
 
     # Radiogenic heating not included
     output["F_radio"] = 0.0

src/proteus/interior/spider.py

@@ -270,6 +270,15 @@ def _try_spider( dirs:dict, config:Config,
     call_sequence.extend(["-MIXING    ", "1"]) # mixing (latent heat transport)
     call_sequence.extend(["-SEPARATION", "1"]) # gravitational separation of solid/melt
 
+    # Tidal heating
+    tidal_value = 0.0
+    if config.interior.tidal_heat:
+        if config.orbit.dummy:
+            tidal_value = config.orbit.dummy.H_tide
+
+        call_sequence.extend(["-HTIDAL",       "1"])
+        call_sequence.extend(["-htidal_value", "%.5e"%tidal_value])
+
     # Properties lookup data (folder relative to SPIDER src)
     folder = "lookup_data/1TPa-dK09-elec-free/"
     call_sequence.extend(["-phase_names",  "melt,solid"])