Changeset 4071 for trunk/LMDZ.COMMON/libf/evolution/glaciers.F90

Timestamp:

Feb 16, 2026, 10:28:56 AM (12 days ago)

Author:

jbclement

Message:

PEM:

• Making the computation of ice tendencies more reliable by doing it after 'read_startpem' to know exactly where perennial ice is (no matter if there is a "startpem.nc" or not). Moreover, when there is no "startpem.nc", location of perennial ice depends now on 'watercaptag' and on huge amount of frost. This prevents negative ice tendency while there is no ice which can happen with weird PCM inputs (i.e. 'watercaptag' = F & 'watercap' /= 0 & no "stratpem.nc").
• Few small safeguards throughout the code.

JBC

File:

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

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

@@ -18 +18 @@
 ! DEPENDENCIES
 ! ------------
-use numerics, only: dp, di, k4
+use numerics, only: dp, di, k4, minieps
 
 ! DECLARATION
@@ -108 +108 @@
 ! ARGUMENTS
 ! ---------
-real(dp),    dimension(:,:),            intent(in)    :: vmr_co2_PEM     ! Grid points x Time field : VMR of CO2 in the first layer [mol/mol]
-real(dp),    dimension(:,:),            intent(in)    :: ps_PCM          ! Grid points x Time field: surface pressure given by the PCM [Pa]
-real(dp),                               intent(in)    :: ps_avg_glob_ini ! Global averaged surface pressure at the beginning [Pa]
-real(dp),                               intent(in)    :: ps_avg_global   ! Global averaged surface pressure during the PEM iteration [Pa]
-real(dp),    dimension(:,:),            intent(inout) :: co2ice          ! Grid points x Slope field: CO2 ice on the subgrid slopes [kg/m^2]
-logical(k4), dimension(:,:),            intent(out)   :: is_co2ice_flow  ! Flag to see if there is flow on the subgrid slopes
+real(dp),    dimension(:,:), intent(in)    :: vmr_co2_PEM     ! Grid points x Time field : VMR of CO2 in the first layer [mol/mol]
+real(dp),    dimension(:,:), intent(in)    :: ps_PCM          ! Grid points x Time field: surface pressure given by the PCM [Pa]
+real(dp),                    intent(in)    :: ps_avg_glob_ini ! Global averaged surface pressure at the beginning [Pa]
+real(dp),                    intent(in)    :: ps_avg_global   ! Global averaged surface pressure during the PEM iteration [Pa]
+real(dp),    dimension(:,:), intent(inout) :: co2ice          ! Grid points x Slope field: CO2 ice on the subgrid slopes [kg/m^2]
+logical(k4), dimension(:,:), intent(out)   :: is_co2ice_flow  ! Flag to see if there is flow on the subgrid slopes
 
 ! LOCAL VARIABLES
@@ -301 +301 @@
                     iaval = islope + 1
                 end if
-                do while (iaval /= iflat .and. subslope_dist(ig,iaval) == 0)
+                do while (iaval /= iflat .and. abs(subslope_dist(ig,iaval)) < minieps)
                     if (iaval > iflat) then
                         iaval = iaval - 1