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

Timestamp:

Apr 10, 2026, 7:17:55 PM (9 hours ago)

Author:

jbclement

Message:

PEM:

• Rework layering-related logic, especially clarify interactions between surface and subsurface water tendencies and disable CO2 lag resistance (inconsistent without updating pressure and mass balance + PCM).
• Prevent simultaneous activation of layering and ice flows.
• Add warning when flux_geo /= 0 while soil is disabled.
• Add new utility function "is_lvl_enabled" for displaying.
• Replace deprecated 'minieps' with 'eps'/'tol'.

JBC

File:

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

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

@@ -18 +18 @@
 ! DEPENDENCIES
 ! ------------
-use numerics, only: dp, di, k4, minieps
+use numerics, only: dp, di, k4, eps
 
 ! DECLARATION
@@ -143 +143 @@
     regolith_inertia(:,islope) = inertiedat(:,1)
     do ig = 1,ngrid
-        if (abs(h2o_ice(ig,islope)) < minieps) regolith_inertia(ig,islope) = TI_regolith_avg
+        if (abs(h2o_ice(ig,islope)) < eps) regolith_inertia(ig,islope) = TI_regolith_avg
         if (reg_thprop_dependp) then
             if (TI(ig,1,islope) < reg_inertie_thresold) then
@@ -191 +191 @@
         if (icetable_depth(ig,islope) > -1.e-6_dp) then
         ! 3.0 Case where it is perennial ice
-            if (icetable_depth(ig,islope) < minieps) then
+            if (icetable_depth(ig,islope) < eps) then
                 call get_ice_TI(.true.,1._dp,0._dp,ice_inertia)
                 do isoil = 1,nsoil