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

Timestamp:

May 15, 2024, 4:45:14 PM (7 months ago)

Author:

jbclement

Message:

PEM:

• Update of some default parameters;
• Correction of a bug when the ice table was too low for the PEM subsurface discretization;
• Correction of the computation for the change of sublimating ice necessary to the stopping criteria (the mistake was introduced in r3149) + simplification of the algorithm + renaming of variables with more explicit names;
• Cleaning of "soil_thermalproperties_mod.F90".

JBC

File:

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

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

@@ -12 +12 @@
 #---------- Orbital parameters ----------#
 # Do you want to follow an orbital forcing read in "obl_ecc_lsp.asc"? Default = .false.
-  evol_orbit_pem=.true.
+# evol_orbit_pem=.false.
 
 # If evol_orbit_pem=.true., number of Earth years before present to start the PEM run? Default = 0
-  year_earth_bp_ini=-1000000
+# year_earth_bp_ini=0
 
 # Do you want to vary the obliquity when following "obl_ecc_lsp.asc"? Default = .true.
-  var_obl=.true.
+# var_obl=.false.
 
 # Do you want to vary the eccentricity when following "obl_ecc_lsp.asc"? Default = .true.
-  var_ecc=.true.
+# var_ecc=.false.
 
 # Do you want to vary the ls perihelie when following "obl_ecc_lsp.asc"? Default = .true.
-  var_lsp=.true.
+# var_lsp=.false.
 
 #---------- Stopping criteria parameters ----------#
 # If evol_orbit_pem=.false., maximal number of iterations if no stopping criterion is reached? Default=100000000
-#  Max_iter_pem=100000000
+# Max_iter_pem=10
 
 # Acceptance rate of sublimating H2O ice surface change? Default = 0.2
@@ -46 +46 @@
 # soil_pem=.true.
 
-# Do you want to run with adsoprtion in the PEM? Default = .true.
+# Do you want to run with adsoprtion in the PEM? Default = .false.
 # adsorption_pem=.true.
 
@@ -90 +90 @@
 
 # Do you want the CO2 ice to flow along subslope inside a cell? Default = .true.
-# co2ice_flow=.true.
+# co2ice_flow=.false.
 
 #---------- Layering parameters ----------#