Changeset 4881 for LMDZ6/trunk/libf/phylmd/yamada4.F90

Timestamp:

Mar 28, 2024, 6:06:40 PM (7 weeks ago)

Author:

evignon

Message:

extraction plus propre de la dissipation de TKE

File:

• LMDZ6/trunk/libf/phylmd/yamada4.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/yamada4.F90

@@ -2 +2 @@
 
 SUBROUTINE yamada4(ni, nsrf, ngrid, dt, g, rconst, plev, temp, zlev, zlay, u, v, teta, &
-    cd, tke, km, kn, kq, ustar, iflag_pbl, drgpro)
+    cd, tke, eps, km, kn, kq, ustar, iflag_pbl, drgpro)
 
   USE dimphy, only : klev,klon
-  USE phys_local_var_mod, only: tke_dissip,wprime
+  USE phys_local_var_mod, only: wprime
   USE yamada_ini_mod, only : new_yamada4,yamada4_num,hboville
   USE yamada_ini_mod, only : prt_level, lunout,pbl_lmixmin_alpha,b1,kap,viscom,viscoh
@@ -69 +69 @@
   ! Outputs:
   !==========
+  ! eps: tke dissipation rate 
   ! km : diffusivite turbulente de quantite de mouvement (au bas de chaque
   ! couche)
@@ -95 +96 @@
   REAL cd(klon)
   REAL tke(klon, klev+1)
+  REAL eps(klon,klev+1)
   REAL unsdz(klon, klev)
   REAL unsdzdec(klon, klev+1)
@@ -686 +688 @@
 
 ! Diagnostics
- tke_dissip(1:ngrid,:,nsrf)=0.
+ 
+ eps(:,:)=dissip(:,:)
  wprime(1:ngrid,:,nsrf)=0.
  DO k=2,klev
@@ -692 +695 @@
        jg=ni(ig)
        wprime(jg,k,nsrf)=sqrt(MAX(1./3*q2(ig,k),0.))
-       tke_dissip(jg,k,nsrf)=dissip(ig,k)
     ENDDO
  ENDDO
+
  
 !=============================================================================