Changeset 3974 for LMDZ6/trunk/libf/phylmd/conf_phys_m.F90

Timestamp:

Jul 30, 2021, 4:27:30 PM (3 years ago)

Author:

asima

Message:

In the "bucket" soil model, soil thermal inertia (formerly a constant value) can also be now a function of soil moisture "qsol"
(F Cheruy's idea; results from "Stage L3" Eve Rebouillat, with E Vignon, A Sima, F Cheruy).

The function depends on iflag_inertie, read from physiq.def (via conf_phys_m.F90), as follows :
(only in physiq.def_NPv6.1 for the time being)

---

# Controle de la definition de l'inertie thermique du sol dans le modele "bucket":
# 0 (default) : constant
# 1 : fonction lineaire de qsol
# 2 : 2 fonctions lineaires de qsol :
# une pour "sable" (boite lat/lon pour Sahara), une pour le reste (limons et argile)
# 3 : fonction lineaire de qsol, valeurs entre 900 pour qsol=0 et 2000 pour qsol=150

iflag_inertie=0

---

Relevant changes are in soil.F90 and physiq.def_NPv6.1 ; the other files are adjusted accordingly.

File:

• LMDZ6/trunk/libf/phylmd/conf_phys_m.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/conf_phys_m.F90

@@ -183 +183 @@
     REAL,SAVE :: exposant_glace_omp
     REAL,SAVE :: rei_min_omp, rei_max_omp
-    INTEGER,SAVE :: iflag_sic_omp
+    INTEGER,SAVE :: iflag_sic_omp, iflag_inertie_omp
     REAL,SAVE :: inertie_sol_omp,inertie_sno_omp,inertie_sic_omp
     REAL,SAVE :: inertie_lic_omp
@@ -1325 +1325 @@
     CALL getin('iflag_sic',iflag_sic_omp)
     !
+    !Config Key  = iflag_inertie
+    !Config Desc =
+    !Config Def  = 0
+    !Config Help =
+    !
+    iflag_inertie_omp = 0
+    CALL getin('iflag_inertie',iflag_inertie_omp)
+    !
     !Config Key  = inertie_sic
     !Config Desc =  
@@ -2395 +2403 @@
     albsno0 = albsno0_omp
     iflag_sic = iflag_sic_omp
+    iflag_inertie = iflag_inertie_omp
     inertie_sol = inertie_sol_omp
     inertie_sic = inertie_sic_omp
@@ -2881 +2890 @@
     WRITE(lunout,*) ' albsno0 = ', albsno0
     WRITE(lunout,*) ' iflag_sic = ', iflag_sic
+    WRITE(lunout,*) ' iflag_inertie = ', iflag_inertie
     WRITE(lunout,*) ' inertie_sol = ', inertie_sol
     WRITE(lunout,*) ' inertie_sic = ', inertie_sic