Changeset 3256 for trunk/LMDZ.COMMON/libf

Timestamp:

Mar 7, 2024, 12:01:29 PM (10 months ago)

Author:

jbclement

Message:

PEM:
Adaptation of threshold values for ice management (in particular 'inf_h2oice_threshold') to more realistic values.
JBC

Location:

trunk/LMDZ.COMMON/libf/evolution

Files:

• changelog.txt (modified) (1 diff)
• conf_pem.F90 (modified) (3 diffs)
• deftank/run_PEM.def (modified) (3 diffs)
• glaciers_mod.F90 (modified) (1 diff)

trunk/LMDZ.COMMON/libf/evolution/changelog.txt

@@ -244 +244 @@
 == 13/02/2024 == JBC
 Small syntax correction of the if condition in "launch_pem.sh" following r3214.
+
+== 07/03/2024 == JBC
+Adaptation of threshold values for ice management (in particular 'inf_h2oice_threshold') to more realistic values.

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

@@ -150 +150 @@
 
 !#---------- Ice management parameters ----------#
-ini_huge_h2oice = 1.e4
+ini_huge_h2oice = 9200. ! kg.m-2 (= 10 m)
 call getin('ini_huge_h2oice',ini_huge_h2oice)
 
-inf_h2oice_threshold = 2.e3
+inf_h2oice_threshold = 460. ! kg.m-2 (= 0.5 m)
 call getin('inf_h2oice_threshold',inf_h2oice_threshold)
 
@@ -159 +159 @@
 call getin('metam_h2oice',metam_h2oice)
 
-metam_h2oice_threshold = 5.e-2
+metam_h2oice_threshold = 460. ! kg.m-2 (= 0.5 m)
 call getin('metam_h2oice_threshold',metam_h2oice_threshold)
 
@@ -168 +168 @@
 call getin('metam_co2ice',metam_co2ice)
 
-metam_co2ice_threshold = 16.e3
+metam_co2ice_threshold = 16500. ! kg.m-2 (= 10 m)
 call getin('metam_co2ice_threshold',metam_co2ice_threshold)
 

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

@@ -67 +67 @@
 
 #---------- Ice management parameters ----------#
-# Amount of H2O ice to initialize the huge reservoir if the variable is not present in "startfi_PEM.nc"? Default = 1.e4
+# 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 = 2.e3
+# 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
 
@@ -76 +76 @@
 # metam_h2oice=.false.
 
-# Threshold to consider frost is becoming perennial H2O ice? Default = 5.e-2
+# Threshold to consider frost is becoming perennial H2O ice? Default = 460. kg.m-2 (= 0.5 m)
 # metam_h2oice_threshold=5.e-2
 
@@ -85 +85 @@
 # metam_co2ice=.false.
 
-# Threshold to consider frost is becoming perennial CO2 ice? Default = 16.e3
+# Threshold to consider frost is becoming perennial CO2 ice? Default = 16500. kg.m-2 (= 10 m)
 # metam_co2ice_threshold=16.e3
 

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

@@ -11 +11 @@
 
 ! 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
+real, save :: inf_h2oice_threshold   ! To consider the amount of H2O ice as an infinite reservoir [kg.m-2]
+real, save :: metam_h2oice_threshold ! To consider frost is becoming perennial H2O ice [kg.m-2]
+real, save :: metam_co2ice_threshold ! To consider frost is becoming perennial CO2 ice [kg.m-2]
 
 !=======================================================================