Changeset 4881 for LMDZ6/trunk/libf/phylmd/coef_diff_turb_mod.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/coef_diff_turb_mod.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/coef_diff_turb_mod.F90

@@ -14 +14 @@
   SUBROUTINE coef_diff_turb(dtime, nsrf, knon, ni, &
        ypaprs, ypplay, yu, yv, yq, yt, yts, yqsurf, ycdragm, &
-       ycoefm, ycoefh ,yq2, ydrgpro)
+       ycoefm, ycoefh ,yq2, yeps, ydrgpro)
  
     USE dimphy
@@ -47 +47 @@
 ! Output arguments
 !****************************************************************************************
+    REAL, DIMENSION(klon,klev+1), INTENT(OUT)  :: yeps
     REAL, DIMENSION(klon,klev), INTENT(OUT)    :: ycoefh
     REAL, DIMENSION(klon,klev), INTENT(OUT)    :: ycoefm
@@ -69 +70 @@
     ykmq = 0 !ym missing init
     
-    
+    yeps(:,:) = 0.
+
 !****************************************************************************************    
 ! Calcul de coefficients de diffusion turbulent de l'atmosphere : 
@@ -182 +184 @@
           CALL yamada4(ni,nsrf,knon,dtime,RG,RD,ypaprs,yt, &
                yzlev,yzlay,yu,yv,yteta, &
-               ycdragm,yq2,ykmm,ykmn,ykmq,yustar, &
+               ycdragm,yq2,yeps,ykmm,ykmn,ykmq,yustar, &
                iflag_pbl,ydrgpro)
 !FC