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

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

@@ -21 +21 @@
 ! ------------
 use numerics,         only: dp, qp, k4
-use layered_deposits, only: layering
+use layered_deposits, only: layering, del_layering
 
 ! DECLARATION
@@ -231 +231 @@
 implicit none
 
+! LOCAL VARIABLES
+! ---------------
+integer :: i, islope
+
 ! CODE
 ! ----
@@ -253 +257 @@
 delta_h2o_ads(:) = 0._dp
 delta_co2_ads(:) = 0._dp
+do islope = 1,nslope
+    do i = 1,ngrid
+        layerings_map(i,islope)%nb_str = 0
+        nullify(layerings_map(i,islope)%bottom,layerings_map(i,islope)%top)
+    end do
+end do
 
 END SUBROUTINE allocate_startevo_state
@@ -386 +396 @@
 !-----------------------------------------------------------------------
 
-! DECLARATION
-! -----------
-implicit none
+! DEPENDENCIES
+! ------------
+use geometry, only: ngrid, nslope
+
+! DECLARATION
+! -----------
+implicit none
+
+! LOCAL VARIABLES
+! ---------------
+integer :: i, islope
 
 ! CODE
@@ -407 +425 @@
 if (allocated(tsoil_ts)) deallocate(tsoil_ts)
 if (allocated(tsoil_ts_old)) deallocate(tsoil_ts_old)
-if (allocated(layerings_map)) deallocate(layerings_map)
+if (allocated(layerings_map)) then
+    do islope = 1,nslope
+        do i = 1,ngrid
+            call del_layering(layerings_map(i,islope))
+        end do
+    end do
+    deallocate(layerings_map)
+end if
 if (allocated(h2o_soildensity_avg)) deallocate(h2o_soildensity_avg)
 if (allocated(delta_co2_ads)) deallocate(delta_co2_ads)
@@ -424 +449 @@
 if (allocated(d_co2ice_ini)) deallocate(d_co2ice_ini)
 if (allocated(d_h2oice)) deallocate(d_h2oice)
-if (allocated(layerings_map)) deallocate(layerings_map)
 if (allocated(flux_ssice_avg)) deallocate(flux_ssice_avg)