Changeset 4071 for trunk/LMDZ.COMMON/libf/evolution/surf_ice.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/surf_ice.F90 (modified) (3 diffs)

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

@@ -200 +200 @@
 
 !=======================================================================
-SUBROUTINE build4PCM_perice(h2o_ice,co2_ice,is_h2o_perice,h2oice_PCM,co2ice_PCM)
+SUBROUTINE build4PCM_perice(h2o_ice,co2_ice,is_h2o_perice,h2o_ice4PCM,co2_ice4PCM)
 !-----------------------------------------------------------------------
 ! NAME
@@ -230 +230 @@
 ! ---------
 real(dp),    dimension(:,:), intent(inout) :: h2o_ice, co2_ice
-logical(k4), dimension(:),   intent(out)   :: is_h2o_perice          ! H2O perennial ice flag
-real(dp),    dimension(:,:), intent(out)   :: h2oice_PCM, co2ice_PCM ! Ice for PCM
+logical(k4), dimension(:),   intent(out)   :: is_h2o_perice            ! H2O perennial ice flag
+real(dp),    dimension(:,:), intent(out)   :: h2o_ice4PCM, co2_ice4PCM ! Ice for PCM
 
 ! LOCAL VARIABLES
@@ -240 +240 @@
 ! ----
 call print_msg('> Building surface ice/frost for the PCM')
-co2ice_PCM(:,:) = co2_ice(:,:)
-h2oice_PCM(:,:) = 0._dp ! Because in the Mars PCM, only the variation of perennial H2O ice is monitored, not the absolute quantity
+co2_ice4PCM(:,:) = co2_ice(:,:)
+h2o_ice4PCM(:,:) = 0._dp ! Because in the Mars PCM, only the variation of perennial H2O ice is monitored, not the absolute quantity
 do i = 1,ngrid
     ! Is H2O ice still considered as an infinite reservoir for the PCM?