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

Timestamp:

May 21, 2025, 3:57:33 PM (7 months 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/stopping_crit_mod.F90 (modified) (2 diffs)

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

@@ -69 +69 @@
 endif
 
-!if (abs(h2oice_ini_surf) < 1.e-5) stopPEM = 0
-
 END SUBROUTINE stopping_crit_h2o_ice
 
@@ -133 +131 @@
 endif
 
-!if (abs(co2ice_sublim_surf_ini) < 1.e-5) stopPEM = 0
-
 if (ps_avg_global < ps_avg_global_ini*(1. - ps_criterion)) then
     stopPEM = 4