Changeset 3938 for trunk/LMDZ.COMMON/libf/evolution/layering_mod.F90

Timestamp:

Oct 29, 2025, 10:02:05 AM (3 weeks ago)

Author:

jbclement

Message:

PEM:

• Correction of the process to balance the H2O flux from and into the atmosphere accross reservoirs: (i) computation of H2O amount going from/in the atmosphere is corrected (missing dt and wrong condition); (ii) computation of the scaling factor is corrected because we need to balance all the icetable/adsorbed/surface ice flux but it affects only sublimating/condensing surface ice flux; (iii) process of balancing is modified to be applied to every flux by scaling instead of applying it only to positive or negative flux by trimming; (iv) balance is now fully processed by the tendency before evolving the ice. This is done through dedicated subroutines.
• Correction of ice conversion between the layering algorithm and PEM ice variables (density factor was missing).
• Addition of H2O flux balance and related stopping criterion for the layering algorithm.
• Few updates for files in the deftank to be in line with the wiki Planeto pages.

JBC

File:

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

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

@@ -594 +594 @@
     if (h_toplag < h_patchy_dust) then ! But the dust lag is too thin to considered ice as subsurface ice
         h_toplag = 0.
-        h2o_ice = current%h_h2oice
+        h2o_ice = current%h_h2oice*rho_h2oice
     endif
 else if (current%h_co2ice > 0. .and. current%h_co2ice > current%h_h2oice) then ! No, there is more CO2 ice than H2O ice
     h_toplag = 0. ! Because it matters only for H2O ice
-    if (h_toplag < h_patchy_ice) co2_ice = current%h_co2ice ! But the dust lag is too thin to considered ice as subsurface ice
+    if (h_toplag < h_patchy_ice) co2_ice = current%h_co2ice*rho_co2ice ! But the dust lag is too thin to considered ice as subsurface ice
 endif