Changeset 4180 for trunk/LMDZ.COMMON/libf/evolution/tendencies.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/tendencies.F90 (modified) (5 diffs)

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

@@ -18 +18 @@
 ! DEPENDENCIES
 ! ------------
-use numerics, only: dp, di, k4, minieps, tol
+use numerics, only: dp, di, k4, eps
 
 ! DECLARATION
@@ -62 +62 @@
 
 ! If the difference is too small, then there is no evolution
-where (abs(d_ice) < minieps) d_ice = 0._dp
+where (abs(d_ice) < eps) d_ice = 0._dp
 
 ! If the tendency is negative but there is no ice reservoir for the PEM
-where (d_ice(:,:) < 0._dp .and. abs(perice(:,:)) < minieps) d_ice(:,:) = 0._dp
+where (d_ice(:,:) < 0._dp .and. abs(perice(:,:)) < eps) d_ice(:,:) = 0._dp
 
 END SUBROUTINE compute_tendice
@@ -170 +170 @@
 ! ARGUMENTS
 ! ---------
-real(dp), dimension(:,:),     intent(in) :: h2oice_depth_old ! Old H2O ice depth
-real(dp), dimension(:,:),     intent(in) :: h2oice_depth_new ! New H2O ice depth
-real(dp), dimension(:,:),     intent(in) :: tsurf            ! Surface temperature
-real(dp), dimension(:,:,:,:), intent(in) :: tsoil_ts_old     ! Old soil temperature time series
-real(dp), dimension(:,:,:,:), intent(in) :: tsoil_ts_new     ! New soil temperature time series
-real(dp), dimension(:,:),     intent(inout) :: flux_ssice_avg       ! Tendency of sub-surface ice
+real(dp), dimension(:,:),     intent(in)    :: h2oice_depth_old ! Old H2O ice depth
+real(dp), dimension(:,:),     intent(in)    :: h2oice_depth_new ! New H2O ice depth
+real(dp), dimension(:,:),     intent(in)    :: tsurf            ! Surface temperature
+real(dp), dimension(:,:,:,:), intent(in)    :: tsoil_ts_old     ! Old soil temperature time series
+real(dp), dimension(:,:,:,:), intent(in)    :: tsoil_ts_new     ! New soil temperature time series
+real(dp), dimension(:,:),     intent(inout) :: flux_ssice_avg   ! Tendency of sub-surface ice
 
 ! LOCAL VARIABLES
@@ -190 +190 @@
     do islope = 1,nslope
         ! Higher resistance due to growing lag layer (higher depth)
-        Rz_old = h2oice_depth_old(i,islope)*zcdv/coef_ssdif_PCM(i,islope) ! Old resistance from PCM
-        Rz_new = h2oice_depth_new(i,islope)*zcdv/coef_ssdif_PCM(i,islope) ! New resistance based on new depth
-        R_dec = Rz_old/Rz_new ! Decrease because of resistance
+        Rz_old = h2oice_depth_old(i,islope)*zcdv/max(abs(coef_ssdif_PCM(i,islope)),eps) ! Old resistance from PCM
+        Rz_new = h2oice_depth_new(i,islope)*zcdv/max(abs(coef_ssdif_PCM(i,islope)),eps) ! New resistance based on new depth
+        R_dec = Rz_old/max(Rz_new,eps) ! Decrease because of resistance
 
         ! The maxmimum of the daily averages over one year for the saturation vapor pressure at the ice table location
@@ -203 +203 @@
 
         ! Lower humidity due to growing lag layer (higher depth)
-        if (abs(psv_max_old) < tol) then
+        if (abs(psv_max_old) < eps) then
             hum_dec = 1._dp
         else