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

Timestamp:

Jul 1, 2025, 2:17:10 PM (6 months ago)

Author:

jbclement

Message:

PEM:
Addition of the flag 'impose_dust_ratio' (default: .false.) to impose a dust-to-ice ratio instead of a dust tendency for the layering algorithm. The ratio value can be set in the "run_PEM.def" by 'dust_to_ice_ratio' (default: 0.01).
JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/deftank/run_PEM.def (modified) (7 diffs)

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

@@ -14 +14 @@
 # evol_orbit_pem=.false.
 
-# If evol_orbit_pem=.true., number of Earth years before present to start the PEM run? Default = 0
+# If evol_orbit_pem=.true., number of Earth years before present to start the PEM run. Default = 0
 # year_earth_bp_ini=0
 
@@ -30 +30 @@
 
 #---------- Stopping criteria parameters ----------#
-# If evol_orbit_pem=.false., maximal number of iterations if no stopping criterion is reached? Default=100000000
+# If evol_orbit_pem=.false., maximal number of iterations if no stopping criterion is reached. Default=100000000
 # Max_iter_pem=10
 
-# 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
 
-# Acceptance rate of sublimating CO2 ice surface change? Default = 0.2
+# Acceptance rate of sublimating CO2 ice surface change. Default = 0.2
 # co2_ice_crit=0.2
 
-# Acceptance rate of pressure surface change? Default = 0.15
+# Acceptance rate of pressure surface change. Default = 0.15
 # ps_criterion = 0.15
 
-# Acceptance of change for obliquity? Default = 1.
+# Acceptance of change for obliquity. Default = 1.
 # max_change_obl=1.
 
@@ -61 +61 @@
 # reg_thprop_dependp=.false.
 
-# Value of the geothermal flux? Default = 0.
+# Value of the geothermal flux. Default = 0.
 # fluxgeo=0.
 
-# Depth at which the breccia layer begins? Default = 10 m
+# Depth at which the breccia layer begins. Default = 10 m
 # depth_breccia=10.
 
-# Depth at which the bedrock layer begins? Default = 1000 m
+# Depth at which the bedrock layer begins. Default = 1000 m
 # depth_bedrock=1000.
 
@@ -77 +77 @@
 
 #---------- Ice management parameters ----------#
-# Amount of H2O ice to initialize the huge reservoir if the variable is not present in "startfi_PEM.nc"? Default = 9200. kg.m-2 (= 10 m)
+# Amount of H2O ice to initialize the huge reservoir if the variable is not present in "startfi_PEM.nc". Default = 9200. kg.m-2 (= 10 m)
 # ini_huge_h2oice=1.e4
 
-# Threshold to consider the amount of H2O ice as an infinite reservoir? Default = 460. kg.m-2 (= 0.5 m)
+# Threshold to consider the amount of H2O ice as an infinite reservoir. Default = 460. kg.m-2 (= 0.5 m)
 # inf_h2oice_threshold=2.e3
 
@@ -86 +86 @@
 # metam_h2oice=.false.
 
-# Threshold to consider frost is becoming perennial H2O ice? Default = 460. kg.m-2 (= 0.5 m)
+# Threshold to consider frost is becoming perennial H2O ice. Default = 460. kg.m-2 (= 0.5 m)
 # metam_h2oice_threshold=5.e-2
 
@@ -95 +95 @@
 # metam_co2ice=.false.
 
-# Threshold to consider frost is becoming perennial CO2 ice? Default = 16500. kg.m-2 (= 10 m)
+# Threshold to consider frost is becoming perennial CO2 ice. Default = 16500. kg.m-2 (= 10 m)
 # metam_co2ice_threshold=16.e3
 
@@ -104 +104 @@
 # Do you want to run with the layering algorithm? Default = .false.
 # layering=.true.
+ 
+# Value of the dust tendency. Default = 1.e-3 kg.m-2.y-1
+# d_dust=1.e-3
+ 
+# Do you want to impose a dust-to-ice ratio instead of a dust tendency? Default = .false.
+# impose_dust_ratio=.true.
+ 
+# If impose_dust_ratio=.true., value of the dust-to-ice ratio. Default = 0.01
+# dust2ice_ratio=0.05
 
 # Some definitions for the physics, in file 'callphys.def'