Changeset 4135 for trunk/LMDZ.COMMON/libf/evolution/soil_therm_inertia.F90

Timestamp:

Mar 16, 2026, 6:09:43 PM (4 weeks ago)

Author:

jbclement

Message:

PEM:
Relocate Mars-specific parameters out of "planet" module and into the modules that own their physics domain.
JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/soil_therm_inertia.F90 (modified) (4 diffs)

trunk/LMDZ.COMMON/libf/evolution/soil_therm_inertia.F90

@@ -26 +26 @@
 ! PARAMETERS
 ! ----------
+real(dp), parameter          :: TI_regolith_avg = 250._dp      ! Average of observed thermal inertia for regolith [SI], Mellon et al. 2000
+real(dp), parameter          :: TI_breccia = 750._dp           ! Thermal inertia for breccia [SI], Mellon et al. 2000; Wood et al. 2008
+real(dp), parameter          :: TI_bedrock = 2300._dp          ! Thermal inertia for bedrock [SI], Mellon et al. 2000; Wood et al. 2008
+real(dp), parameter, private :: P610 = 610._dp                 ! Reference mean surface pressure [Pa]
 real(dp), parameter, private :: reg_inertie_thresold = 370._dp ! Above this thermal inertia, the regolith has too much cementation to be dependant on the pressure [J/m^2/K/s^1/2]
 real(dp), parameter, private :: reg_inertie_minvalue = 50._dp  ! Minimum value of the Thermal Inertia at low pressure (Piqueux & Christensen 2009) [J/m^2/K/s^1/2]
@@ -56 +60 @@
 ! DEPENDENCIES
 ! ------------
-use planet, only: porosity
+use soil, only: regolith_porosity
 
 ! DECLARATION
@@ -78 +82 @@
     ice_thermalinertia = inertie_purewaterice
 else
-    ice_thermalinertia = sqrt(surf_thermalinertia**2 + porosity*pore_filling*inertie_purewaterice**2) ! Siegler et al., 2012
+    ice_thermalinertia = sqrt(surf_thermalinertia**2 + regolith_porosity*pore_filling*inertie_purewaterice**2) ! Siegler et al., 2012
 end if
 
@@ -106 +110 @@
 use geometry, only: ngrid, nslope, nsoil
 use soil,     only: volcapa, layer, inertiedat, depth_breccia, depth_bedrock, index_breccia, index_bedrock, reg_thprop_dependp
-use planet,   only: TI_breccia, TI_bedrock, TI_regolith_avg, P610
 use display,  only: print_msg, LVL_NFO