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

Timestamp:

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

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

@@ -379 +379 @@
      fracIR = 0.04*p0(k)/600.; fracDust = 0.02*p0(k)/600.
      B = Diff*bigstep*86400.*365.24/(porosity*927.)
-     !B = Diff*bigstep*86400.*365.24/(porosity*rho_ice(Tb(),'h2o'))
+     !B = Diff*bigstep*86400.*365.24/(porosity*rho_ice(Tb(),'H2O'))
 
      typeT = -9
@@ -869 +869 @@
 
   B = Diff*bigstep*86400.*365.24/(porosity*927.)
-  !B = Diff*bigstep*86400.*365.24/(porosity*rho_ice(T,'h2o'))
+  !B = Diff*bigstep*86400.*365.24/(porosity*rho_ice(T,'H2O'))
 
   ! advance ice table, avdrho>0 is retreat
@@ -889 +889 @@
         beta = (1-icefrac)/(1-porosity)/icefrac
         beta = Diff*bigstep*beta*86400*365.24/927.
-        !beta = Diff*bigstep*beta*86400*365.24/rho_ice(T,'h2o')
+        !beta = Diff*bigstep*beta*86400*365.24/rho_ice(T,'H2O')
         zdepthT = sqrt(2*beta*avdrho*18./8314. + zdepthT**2)
      endif