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

Timestamp:

Apr 10, 2026, 7:17:55 PM (6 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/sorption.F90 (modified) (4 diffs)

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

@@ -18 +18 @@
 ! DEPENDENCIES
 ! ------------
-use numerics, only: dp, qp, di, k4, minieps
+use numerics, only: dp, qp, di, k4, eps
 
 ! DECLARATION
@@ -230 +230 @@
         deltam_reg_slope(ig,islope) = 0._dp
         do iloop = 1,index_breccia
-            if (TI(ig,iloop,islope) < inertie_thresold .and. abs(h2o_ice(ig,islope)) < minieps .and. abs(co2_ice(ig,islope)) < minieps) then
+            if (TI(ig,iloop,islope) < inertie_thresold .and. abs(h2o_ice(ig,islope)) < eps .and. abs(co2_ice(ig,islope)) < eps) then
                 if (iloop == 1) then
                     deltam_reg_complete(ig,iloop,islope) = (dm_h2o_regolith_slope(ig,iloop,islope) - h2o_ads_reg(ig,iloop,islope))*(layer(iloop))
@@ -347 +347 @@
     do islope = 1,nslope
         do iloop = 1,index_breccia
-            if (TI(ig,iloop,islope) < inertie_thresold .and. abs(h2o_ice(ig,islope)) < minieps .and. abs(co2_ice(ig,islope)) < minieps) then
+            if (TI(ig,iloop,islope) < inertie_thresold .and. abs(h2o_ice(ig,islope)) < eps .and. abs(co2_ice(ig,islope)) < eps) then
                 dm_co2_regolith_slope(ig,iloop,islope) = as*rho_regolith*m_theta*(1._dp - theta_h2o_ads(ig,iloop,islope))*alpha*pco2_avg(ig)/ &
                                                          (alpha*pco2_avg(ig) + sqrt(tsoil(ig,iloop,islope))*exp(beta/tsoil(ig,iloop,islope)))
             else
-                if (abs(co2_ads_reg(ig,iloop,islope)) < minieps) then !!! we are at first call
+                if (abs(co2_ads_reg(ig,iloop,islope)) < eps) then !!! we are at first call
                     dm_co2_regolith_slope(ig,iloop,islope) = as*rho_regolith*m_theta*(1._dp - theta_h2o_ads(ig,iloop,islope))*alpha*pco2_avg(ig) &
                                                              /(alpha*pco2_avg(ig)+sqrt(tsoil(ig,iloop,islope))*exp(beta/tsoil(ig,iloop,islope)))
@@ -368 +368 @@
         deltam_reg_slope(ig,islope) = 0._dp
         do iloop = 1,index_breccia
-            if (TI(ig,iloop,islope) < inertie_thresold .and. abs(h2o_ice(ig,islope)) < minieps .and. abs(co2_ice(ig,islope)) < minieps) then
+            if (TI(ig,iloop,islope) < inertie_thresold .and. abs(h2o_ice(ig,islope)) < eps .and. abs(co2_ice(ig,islope)) < eps) then
                 if (iloop == 1) then
                     deltam_reg_complete(ig,iloop,islope) = (dm_co2_regolith_slope(ig,iloop,islope) - co2_ads_reg(ig,iloop,islope))*(layer(iloop))