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

Timestamp:

May 21, 2025, 3:57:33 PM (5 weeks ago)

Author:

jbclement

Message:

PEM:
Revision of the layering structure and algorithm:

• 'stratum' components are now expressed in height
• deletion of the (redundant) thickness feature of 'stratum'
• 'stratif' in the PEM main program is renamed 'deposits'
• addition of several procedures to get useful information about a stratum (major component or thickness)
• all subsurface layers are now represented in the layering structure by strata with negative top elevation
• simplification of the different situations arising from the input tendencies
• porosity is determined for the entire stratum (and not anymore for each component) based on dominant component
• sublimation of CO2 and H2O ice is now handled simultaneously (more realistic) in a stratum
• linking the layering algorithm with the PEM initilization/finalization regarding PCM data and with the PEM stopping criteria
• making separate cases for glaciers vs layering management
• H2O sublimation flux correction is now handled with the layering when a dust lag layer layer is growing
• update of 'h2o_ice_depth' and 'zdqsdif' accordingly at the PEM end for the PCM

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 = icetable_depth_old*zcdv/coef_diff ! Old resistance from PCM
 Rz_new = icetable_depth_new*zcdv/coef_diff ! New resistance based on new depth
-R_dec = Rz_new/Rz_old ! Decrease because of resistance
+if (abs(Rz_old) < 1.e-10) then
+    R_dec = 1.
+else
+    R_dec = Rz_new/Rz_old ! Decrease because of resistance
+endif
 
 ! The maxmimum of the daily averages over one year for the saturation vapor pressure at the ice table location
@@ -103 +107 @@
 
 ! Lower humidity due to growing lag layer (higher depth)
+print*, psv_max_old,psv_max_new
 hum_dec = psv_max_old/psv_max_new ! Decrease because of lower water vapor pressure at the new depth
 
 ! Flux correction (decrease)
+print*, 'aaaaah', R_dec,hum_dec
 d_h2oice = d_h2oice/R_dec/hum_dec