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

Timestamp:

Mar 16, 2026, 6:09:43 PM (2 days ago)

Author:

jbclement

Message:

PEM:
Relocate Mars-specific parameters out of "planet" module and into the modules that own their physics domain.
JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/ice_table.F90 (modified) (2 diffs)

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

@@ -276 +276 @@
 ! ------------
 use geometry, only: ngrid, nslope, nsoil
-use soil,     only: mlayer
+use soil,     only: mlayer, regolith_porosity
 use slopes,   only: subslope_dist, def_slope_mean
-use planet,   only: porosity
 use maths,    only: pi
 
@@ -306 +305 @@
         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) + regolith_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