Changeset 4074 for trunk/LMDZ.COMMON/libf/evolution/pem.F90

Timestamp:

Feb 17, 2026, 2:45:53 PM (9 days ago)

Author:

jbclement

Message:

PEM:

• Correct management of H2O ice tendency in 1D when there is not enough ice anymore.
• Clean initialization of allocatable module arrays (especially needed when no slope)
• One more renaming for consistency + few small updates thoughout the code.

JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/pem.F90 (modified) (1 diff)

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

@@ -246 +246 @@
         if (h2o_ice(i,islope) > 0._dp) then
             is_h2oice_ini(i,islope) = .true.
-            if (d_co2ice(i,islope) < 0._dp) h2oice_sublim_coverage_ini = h2oice_sublim_coverage_ini + cell_area(i)*subslope_dist(i,islope)/cos(pi*def_slope_mean(islope)/180._dp)
+            if (d_h2oice(i,islope) < 0._dp) h2oice_sublim_coverage_ini = h2oice_sublim_coverage_ini + cell_area(i)*subslope_dist(i,islope)/cos(pi*def_slope_mean(islope)/180._dp)
         end if
     end do