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

Timestamp:

Feb 17, 2026, 2:45:53 PM (13 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/ice_table.F90 (modified) (4 diffs)

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

@@ -67 +67 @@
 allocate(icetable_thickness(ngrid,nslope))
 allocate(ice_porefilling(ngrid,nsoil,nslope))
+icetable_depth(:,:) = 0._dp
+icetable_thickness(:,:) = 0._dp
+ice_porefilling(:,:,:) = 0._dp
 
 END SUBROUTINE ini_ice_table
@@ -303 +306 @@
         do islope = 1,nslope
             call compute_Tice(tsoil(ig,:,islope),tsurf(ig,islope),icetable_depth(ig,islope),Tice)
-            delta_icetable(ig) = delta_icetable(ig) + porosity*rho_ice(Tice,'h2o')*(icetable_thickness(ig,islope) - icetable_thickness_old(ig,islope)) & ! convention > 0. <=> it condenses
+            delta_icetable(ig) = delta_icetable(ig) + porosity*rho_ice(Tice,'H2O')*(icetable_thickness(ig,islope) - icetable_thickness_old(ig,islope)) & ! convention > 0. <=> it condenses
                               *subslope_dist(ig,islope)/cos(def_slope_mean(islope)*pi/180._dp)
         end do
@@ -313 +316 @@
             do isoil = 1,nsoil
                 call compute_Tice(tsoil(ig,:,islope),tsurf(ig,islope),mlayer(isoil - 1),Tice)
-                delta_icetable(ig) = delta_icetable(ig) + rho_ice(Tice,'h2o')*(ice_porefilling(ig,isoil,islope) - ice_porefilling_old(ig,isoil,islope)) & ! convention > 0. <=> it condenses
+                delta_icetable(ig) = delta_icetable(ig) + rho_ice(Tice,'H2O')*(ice_porefilling(ig,isoil,islope) - ice_porefilling_old(ig,isoil,islope)) & ! convention > 0. <=> it condenses
                                   *subslope_dist(ig,islope)/cos(def_slope_mean(islope)*pi/180._dp)
             end do
@@ -586 +589 @@
 ! ----
 select case (trim(adjustl(ice)))
-    case('h2o')
+    case('H2O')
         rho = -3.5353e-4_dp*T**2 + 0.0351_dp*T + 933.5030_dp ! Rottgers 2012
-    case('co2')
+    case('CO2')
         rho = (1.72391_dp - 2.53e-4_dp*T - 2.87_dp*1.e-7_dp*T**2)*1.e3_dp ! Mangan et al. 2017
     case default