Changeset 3780 for LMDZ6/trunk/libf/phylmd/ocean_slab_mod.F90

Timestamp:

Oct 22, 2020, 2:50:18 PM (4 years ago)

Author:

evignon

Message:

Premiere comission Etienne: changements pour le 1D (forcage en Ts au dessus des continents) et inclusion drag arbres dans yamada4_num=6

File:

• LMDZ6/trunk/libf/phylmd/ocean_slab_mod.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/ocean_slab_mod.F90

@@ -421 +421 @@
 !    
 !****************************************************************************************
-    cal(:)      = 0. ! infinite thermal inertia
-    beta(:)     = 1. ! wet surface
-    dif_grnd(:) = 0. ! no diffusion into ground
+    !cal(:)      = 0. ! infinite thermal inertia
+    !beta(:)     = 1. ! wet surface
+    !dif_grnd(:) = 0. ! no diffusion into ground
+    ! EV: use calbeta
+    CALL calbeta(dtime, is_oce, knon, snow,qsurf, beta, cal, dif_grnd)
+
+
     
 ! Suppose zero surface speed
@@ -742 +746 @@
 ! set beta, cal, compute conduction fluxes inside ice/snow 
     slab_bilg(:)=0.
-    dif_grnd(:)=0.
-    beta(:) = 1.
+    !dif_grnd(:)=0.
+    !beta(:) = 1.
+    ! EV: use calbeta to calculate beta and then recalculate properly cal
+    CALL calbeta(dtime, is_sic, knon, snow, qsol, beta, cal, dif_grnd)
+
+
     DO i=1,knon
     ki=knindex(i)