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

Timestamp:

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

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

@@ -124 +124 @@
 call print_msg("> Flow of CO2 glaciers")
 call computeTcondCO2(vmr_co2_PEM,ps_PCM,ps_avg_glob_ini,ps_avg_global,Tcond)
-call compute_hmaxglaciers(Tcond,"co2",hmax)
+call compute_hmaxglaciers(Tcond,"CO2",hmax)
 call transfer_ice_duringflow(hmax,Tcond,co2ice,is_co2ice_flow)
 
@@ -169 +169 @@
 ! ----
 call print_msg("> Flow of H2O glaciers")
-call compute_hmaxglaciers(Tice,"h2o",hmax)
+call compute_hmaxglaciers(Tice,"H2O",hmax)
 call transfer_ice_duringflow(hmax,Tice,h2oice,is_h2oice_flow)
 
@@ -222 +222 @@
 ! ----
 select case (trim(adjustl(name_ice)))
-    case('h2o')
+    case('H2O')
         tau_d = 1.e5_dp
-    case('co2')
+    case('CO2')
         tau_d = 5.e3_dp
     case default
@@ -294 +294 @@
         if (islope /= iflat) then ! ice can be infinite on flat ground
 ! First: check that CO2 ice must flow (excess of ice on the slope), ice can accumulate infinitely  on flat ground
-            if (qice(ig,islope) >= rho_ice(Tice(ig,islope),'h2o')*hmax(ig,islope)*cos(pi*def_slope_mean(islope)/180._dp)) then
+            if (qice(ig,islope) >= rho_ice(Tice(ig,islope),'H2O')*hmax(ig,islope)*cos(pi*def_slope_mean(islope)/180._dp)) then
 ! Second: determine the flatest slopes possible:
                 if (islope > iflat) then
@@ -308 +308 @@
                     end if
                 end do
-                qice(ig,iaval) = qice(ig,iaval) + (qice(ig,islope) - rho_ice(Tice(ig,islope),'h2o')*hmax(ig,islope)*cos(pi*def_slope_mean(islope)/180._dp)) &
+                qice(ig,iaval) = qice(ig,iaval) + (qice(ig,islope) - rho_ice(Tice(ig,islope),'H2O')*hmax(ig,islope)*cos(pi*def_slope_mean(islope)/180._dp)) &
                                *subslope_dist(ig,islope)/subslope_dist(ig,iaval)*cos(pi*def_slope_mean(iaval)/180._dp)/cos(pi*def_slope_mean(islope)/180._dp)
 
-                qice(ig,islope) = rho_ice(Tice(ig,islope),'h2o')*hmax(ig,islope)*cos(pi*def_slope_mean(islope)/180._dp)
+                qice(ig,islope) = rho_ice(Tice(ig,islope),'H2O')*hmax(ig,islope)*cos(pi*def_slope_mean(islope)/180._dp)
                 flag_flow(ig,islope) = .true.
             end if ! co2ice > hmax