Changeset 3809 for trunk/LMDZ.COMMON/libf/evolution/recomp_tend_mod.F90

Timestamp:

Jun 16, 2025, 5:11:50 PM (3 weeks ago)

Author:

jbclement

Message:

PEM:

• Correction to detect whether the stratum below the dust lag is made of ice or not.
• Correction for the initialization of surface ice and 'h2o_ice_depth' according to the layerings map.
• Correction to handle flow of glaciers with the layering algorithm.
• Simplification of "recomp_tend_h2o", making it more robust.
• Making the removing of a stratum more robust.
• Plotting the orbital variations over the simulated time in "visu_evol_layering.py".
• Lowering the value of dust tendency.
• Addition of "_co2" in the variables name for the CO2 flux correction.

JBC

File:

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

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

@@ -92 +92 @@
 Rz_old = h2oice_depth_old*zcdv/coef_diff ! Old resistance from PCM
 Rz_new = h2oice_depth_new*zcdv/coef_diff ! New resistance based on new depth
-R_dec = 1.
-if (Rz_new >= Rz_old) R_dec = Rz_new/Rz_old ! Decrease because of resistance
+R_dec = Rz_old/Rz_new ! Decrease because of resistance
 
 ! The maxmimum of the daily averages over one year for the saturation vapor pressure at the ice table location
@@ -104 +103 @@
 
 ! Lower humidity due to growing lag layer (higher depth)
-hum_dec = psv_max_old/psv_max_new ! Decrease because of lower water vapor pressure at the new depth
+if (abs(psv_max_old) < 1.e2*epsilon(1.)) then
+    hum_dec = 1.
+else
+    hum_dec = psv_max_new/psv_max_old ! Decrease because of lower water vapor pressure at the new depth
+endif
 
 ! Flux correction (decrease)
-d_h2oice = d_h2oice/R_dec/hum_dec
+d_h2oice = d_h2oice*R_dec*hum_dec
 
 END SUBROUTINE recomp_tend_h2o