Changeset 3161 for trunk/LMDZ.COMMON/libf/evolution/deftank

Timestamp:

Dec 19, 2023, 3:14:48 PM (13 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/deftank/run_PEM.def (modified) (3 diffs)

trunk/LMDZ.COMMON/libf/evolution/deftank/run_PEM.def

@@ -28 +28 @@
 #  Max_iter_pem=100000000
 
-# Acceptance rate of sublimating H2o ice surface change? Default = 0.2
+# Acceptance rate of sublimating H2O ice surface change? Default = 0.2
 # h2o_ice_crit=0.2
 
@@ -47 +47 @@
 # adsorption_pem=.true.
 
-# Do you want the glaciers to flow along subslope inside a mesh? Default = .true.
-# co2glaciersflow=.true.
-
 # Do you want to modify the soil thermal properties with the pressure? Default = .false.
 # reg_thprop_dependp=.false.
 
 #---------- Layering parameters ----------#
-# Value of the geothermal flux? Default= 0.
+# Value of the geothermal flux? Default = 0.
 # fluxgeo=0.
 
@@ -70 +67 @@
 
 #---------- Ice management parameters ----------#
-# Amount of h2o ice to initialize the big reservoir if the variable is not present in "startfi_PEM.nc"? Default = 1.e4
-# ini_h2o_bigreservoir=1.e4
+# Amount of H2O ice to initialize the huge reservoir if the variable is not present in "startfi_PEM.nc"? Default = 1.e4
+# ini_huge_h2oice=1.e4
 
 # Threshold to consider the amount of H2O ice as an infinite reservoir? Default = 2.e3
 # inf_h2oice_threshold=2.e3
 
+# Do you want H2O frost to transform into perennial H2O ice? Default = .false.
+# metam_h2oice=.false.
+
 # Threshold to consider frost is becoming perennial H2O ice? Default = 5.e-2
 # metam_h2oice_threshold=5.e-2
+
+# Do you want the H2O ice to flow along subslope inside a cell? Default = .true.
+# h2oice_flow=.true.
+
+# Do you want CO2 frost to transform into perennial CO2 ice? Default = .false.
+# metam_co2ice=.false.
 
 # Threshold to consider frost is becoming perennial CO2 ice? Default = 16.e3
 # metam_co2ice_threshold=16.e3
 
+# Do you want the CO2 ice to flow along subslope inside a cell? Default = .true.
+# co2ice_flow=.true.
+
 # Some definitions for the physics, in file 'callphys.def'
 INCLUDEDEF=callphys.def