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

Timestamp:

Feb 16, 2026, 10:28:56 AM (2 weeks 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/frost.F90 (modified) (3 diffs)

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

@@ -164 +164 @@
 
 !=======================================================================
-SUBROUTINE compute_frost4PCM(min_h2ofrost,min_co2frost)
+SUBROUTINE compute_frost4PCM(minPCM_h2ofrost,minPCM_co2frost)
 !-----------------------------------------------------------------------
 ! NAME
@@ -190 +190 @@
 ! ARGUMENTS
 ! ---------
-real(dp), dimension(:,:), intent(in) :: min_h2ofrost, min_co2frost
+real(dp), dimension(:,:), intent(in) :: minPCM_h2ofrost, minPCM_co2frost
 
 ! CODE
@@ -198 +198 @@
 h2o_frost4PCM(:,:) = 0._dp
 co2_frost4PCM(:,:) = 0._dp
-where (h2ofrost_PCM(:,:) > 0._dp) h2o_frost4PCM(:,:) = h2ofrost_PCM(:,:) - min_h2ofrost(:,:)
-where (co2frost_PCM(:,:) > 0._dp) co2_frost4PCM(:,:) = co2frost_PCM(:,:) - min_co2frost(:,:)
+where (h2ofrost_PCM(:,:) > 0._dp) h2o_frost4PCM(:,:) = h2ofrost_PCM(:,:) - minPCM_h2ofrost(:,:)
+where (co2frost_PCM(:,:) > 0._dp) co2_frost4PCM(:,:) = co2frost_PCM(:,:) - minPCM_co2frost(:,:)
 
 END SUBROUTINE compute_frost4PCM