Changeset 3161 for trunk/LMDZ.COMMON/libf/evolution/glaciers_mod.F90

Timestamp:

Dec 19, 2023, 3:14:48 PM (11 months ago)

Author:

jbclement

Message:

PEM:

• Addition of flags defined in the "run_PEM.def" to decide to do or not CO2 & H2O ice metamorphism: 'metam_co2ice' and 'metam_h2oice' (default is false).
• The variations of infinite reservoirs ('watercap') during the PCM years are now taken into account to update H2O ice at the PEM initialization.
• 'ini_h2o_bigreservoir' is renamed into 'ini_huge_h2oice'.
• Some cleanings, in particular for the main program "pem.F90".

JBC

File:

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

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

@@ -3 +3 @@
 implicit none
 
-logical :: co2_ice_flow ! True by default, to compute co2 ice flow. Read in run_PEM.def
-logical :: h2o_ice_flow ! True by default, to compute h2o ice flow. Read in run_PEM.def
+
+! Flags for ice management
+logical :: h2oice_flow ! True by default, to compute H2O ice flow. Read in "run_PEM.def"
+logical :: co2ice_flow ! True by default, to compute CO2 ice flow. Read in "run_PEM.def"
+logical :: metam_h2oice ! False by default, to compute H2O ice metamorphism. Read in "run_PEM.def"
+logical :: metam_co2ice ! False by default, to compute CO2 ice metamorphism. Read in "run_PEM.def"
+
+! Thresholds for ice management
+real, save :: inf_h2oice_threshold   ! To consider the amount of H2O ice as an infinite reservoir
+real, save :: metam_h2oice_threshold ! To consider frost is becoming perennial H2O ice
+real, save :: metam_co2ice_threshold ! To consider frost is becoming perennial CO2 ice
 
 !=======================================================================