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

LMDZ6/trunk/libf/phylmd/ocean_forced_mod.F90

@@ -180 +180 @@
 !
     USE dimphy
+    USE geometry_mod, ONLY: longitude,latitude
     USE calcul_fluxs_mod
     USE surface_data,     ONLY : calice, calsno
@@ -260 +261 @@
     IF (soil_model) THEN 
 ! update tsoil and calculate soilcap and soilflux
-       CALL soil(dtime, is_sic, knon, snow, tsurf_tmp, tsoil,soilcap, soilflux)
+       CALL soil(dtime, is_sic, knon, snow, tsurf_tmp, qsol, &
+        & longitude(knindex(1:knon)), latitude(knindex(1:knon)), tsoil,soilcap, soilflux)
        cal(1:knon) = RCPD / soilcap(1:knon)
        radsol(1:knon) = radsol(1:knon)  + soilflux(1:knon)