Changeset 2441

Timestamp:

Feb 12, 2016, 5:22:02 PM (9 years ago)

Author:

fhourdin

Message:

Introduction d'une longueur de mélange minimum dans Yamada4 (Cheruy)
Correction du Cd neutre pour chaleur et humidité k^{2/ln(z/z0m)/ln(z/z0h)}

Location:

LMDZ5/trunk/libf/phylmd

Files:

• cdrag.F90 (modified) (2 diffs)
• yamada4.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/cdrag.F90

@@ -114 +114 @@
   LOGICAL, PARAMETER    :: zxli=.FALSE. ! calcul des cdrags selon Laurent Li
   REAL, DIMENSION(klon) :: zcdn_m, zcdn_h         ! Drag coefficient in neutral conditions
+  REAL zzzcd
 !
 ! Fonctions thermodynamiques et fonctions d'instabilite
@@ -176 +177 @@
 
 
-! Coefficients CD neutres pour m et h
-     zcdn_m(i) = (CKAP/LOG(1.+zgeop1(i)/(RG*z0m(i))))**2
-     zcdn_h(i) = (CKAP/LOG(1.+zgeop1(i)/(RG*z0h(i))))**2
+! Coefficients CD neutres pour m et h : k^2/ln(z/z0) et k^2/(ln(z/z0)*ln(z/z0h))
+     zzzcd=CKAP/LOG(1.+zgeop1(i)/(RG*z0m(i)))
+     zcdn_m(i) = zzzcd*zzzcd
+     zcdn_h(i) = zzzcd*(CKAP/LOG(1.+zgeop1(i)/(RG*z0m(i))))
 
      IF (zri(i) .GT. 0.) THEN      ! situation stable

LMDZ5/trunk/libf/phylmd/yamada4.F90

@@ -6 +6 @@
   USE dimphy
   USE print_control_mod, ONLY: prt_level
+  USE ioipsl_getin_p_mod, ONLY : getin_p
+
   IMPLICIT NONE
 
@@ -75 +77 @@
   DATA first, ipas/.FALSE., 0/
   !$OMP THREADPRIVATE( first,ipas)
+  REAL,SAVE :: lmixmin=1.
+  !$OMP THREADPRIVATE(lmixmin) 
+
 
   INTEGER ig, k
@@ -107 +112 @@
   fl(zzz, zl0, zq2, zn2) = max(min(l0(ig)*kap*zlev(ig, &
     k)/(kap*zlev(ig,k)+l0(ig)),0.5*sqrt(q2(ig,k))/sqrt( &
-    max(n2(ig,k),1.E-10))), 1.)
+    max(n2(ig,k),1.E-10))), lmixmin)
 
 
@@ -116 +121 @@
     ALLOCATE (l0(klon))
     firstcall = .FALSE.
+    CALL getin_p('lmixmin',lmixmin)
   END IF
 
@@ -341 +347 @@
     DO k = 2, klev - 1
       DO ig = 1, ngrid
-        l(ig, k) = max(l(ig,k), 1.)
+        l(ig, k) = max(l(ig,k), lmixmin)
         km(ig, k) = l(ig, k)*sqrt(q2(ig,k))*sm(ig, k)
         q2(ig, k) = q2(ig, k) + dt*km(ig, k)*m2(ig, k)*(1.-rif(ig,k))