Changeset 3974 for LMDZ6/trunk/libf/phylmd/surf_land_bucket_mod.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/surf_land_bucket_mod.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/surf_land_bucket_mod.F90

@@ -24 +24 @@
     USE cpl_mod
     USE dimphy
-    USE geometry_mod, ONLY: latitude 
+    USE geometry_mod, ONLY: longitude,latitude 
     USE mod_grid_phy_lmdz
     USE mod_phys_lmdz_para
@@ -103 +103 @@
        
 ! calculate temperature, heat capacity and conduction flux in soil
-    IF (soil_model) THEN 
-       CALL soil(dtime, is_ter, knon, snow, tsurf, tsoil, soilcap, soilflux)
+    IF (soil_model) THEN
+       CALL soil(dtime, is_ter, knon, snow, tsurf, qsol,  & 
+        & longitude(knindex(1:knon)), latitude(knindex(1:knon)), tsoil, soilcap, soilflux)
+
        DO i=1, knon
           cal(i) = RCPD / soilcap(i)