Changeset 4881 for LMDZ6/trunk/libf/phylmdiso/pbl_surface_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/phylmdiso/pbl_surface_mod.F90 (modified) (15 diffs)

LMDZ6/trunk/libf/phylmdiso/pbl_surface_mod.F90

@@ -24 +24 @@
   USE climb_wind_mod,      ONLY : climb_wind_down, climb_wind_up
   USE coef_diff_turb_mod,  ONLY : coef_diff_turb
-  USE lmdz_call_atke,      ONLY :  call_atke
+  USE lmdz_call_atke,      ONLY : call_atke
   USE ioipsl_getin_p_mod,  ONLY : getin_p
   USE cdrag_mod
@@ -308 +308 @@
 !jyg<
 !!       zxfluxt,   zxfluxq,   q2m,      flux_q, tke,   &
-       zxfluxt,   zxfluxq, zxfluxqbs,   q2m, flux_q, flux_qbs, tke_x,  &
+       zxfluxt,   zxfluxq, zxfluxqbs,   q2m, flux_q, flux_qbs, tke_x, eps_x, &
 !>jyg
 !!! nrlmd+jyg le 02/05/2011 et le 20/02/2012
@@ -506 +506 @@
 ! Output variables
 !****************************************************************************************
+    REAL, DIMENSION(klon,klev+1,nbsrf+1), INTENT(OUT)   :: eps_x      ! TKE dissipation rate
+
     REAL, DIMENSION(klon),        INTENT(OUT)       :: cdragh     ! drag coefficient for T and Q
     REAL, DIMENSION(klon),        INTENT(OUT)       :: cdragm     ! drag coefficient for wind
@@ -765 +767 @@
     REAL, DIMENSION(klon,klev)         :: delp
     REAL, DIMENSION(klon,klev+1)       :: ypaprs
-    REAL, DIMENSION(klon,klev+1)       :: ytke
+    REAL, DIMENSION(klon,klev+1)       :: ytke, yeps
     REAL, DIMENSION(klon,nsoilmx)      :: ytsoil
 !FC 
@@ -785 +787 @@
     REAL, DIMENSION(klon)              :: ywake_s, ywake_cstar, ywake_dens
 !!!
-    REAL, DIMENSION(klon,klev+1)       :: ytke_x, ytke_w
+    REAL, DIMENSION(klon,klev+1)       :: ytke_x, ytke_w, yeps_x, yeps_w
     REAL, DIMENSION(klon,klev+1)       :: ywake_dltke
     REAL, DIMENSION(klon,klev)         :: yu_x, yv_x, yu_w, yv_w
@@ -1261 +1263 @@
 !!    tke(:,:,is_ave)=0.
     tke_x(:,:,is_ave)=0.
+    eps_x(:,:,is_ave)=0.
 
     wake_dltke(:,:,is_ave)=0.
@@ -1295 +1298 @@
 
     ytke=0.
+    yeps=0.
     yri0(:)=0.
 !FC 
@@ -1308 +1312 @@
 !!! nrlmd+jyg le 02/05/2011 et le 20/02/2012
     ytke_x=0.     ; ytke_w=0.        ; ywake_dltke=0.
+    yeps_x=0.     ; yeps_w=0.
     y_d_t_x=0.    ; y_d_t_w=0.       ; y_d_q_x=0.      ; y_d_q_w=0.
 !!    d_t_w=0.      ; d_q_w=0.         
@@ -1979 +1984 @@
         CALL call_atke(dtime,knon,klev,ycdragm(1:knon), ycdragh(1:knon),yus0(1:knon),yvs0(1:knon),yts(1:knon), &
                   yu(1:knon,:),yv(1:knon,:),yt(1:knon,:),yq(1:knon,:),ypplay(1:knon,:),ypaprs(1:knon,:),       &
-                  ytke(1:knon,:),ycoefm(1:knon,:), ycoefh(1:knon,:))
+                  ytke(1:knon,:),yeps(1:knon,:), ycoefm(1:knon,:), ycoefh(1:knon,:))
 
         ELSE
@@ -1985 +1990 @@
         CALL coef_diff_turb(dtime, nsrf, knon, ni,  &
             ypaprs, ypplay, yu, yv, yq, yt, yts, yqsurf, ycdragm, &
-            ycoefm, ycoefh, ytke, y_treedrg)
+            ycoefm, ycoefh, ytke, yeps, y_treedrg)
 !            ycoefm, ycoefh, ytke)
 !FC y_treedrg ajoute
@@ -2023 +2028 @@
         CALL call_atke(dtime,knon,klev,ycdragm_x(1:knon),ycdragh_x(1:knon),yus0(1:knon),yvs0(1:knon),yts_x(1:knon),    &
                        yu_x(1:knon,:),yv_x(1:knon,:),yt_x(1:knon,:),yq_x(1:knon,:),ypplay(1:knon,:),ypaprs(1:knon,:),  &
-                       ytke_x(1:knon,:),ycoefm_x(1:knon,:), ycoefh_x(1:knon,:))
+                       ytke_x(1:knon,:),yeps_x(1:knon,:),ycoefm_x(1:knon,:), ycoefh_x(1:knon,:))
 
         ELSE
@@ -2029 +2034 @@
         CALL coef_diff_turb(dtime, nsrf, knon, ni,  &
             ypaprs, ypplay, yu_x, yv_x, yq_x, yt_x, yts_x, yqsurf_x, ycdragm_x, &
-            ycoefm_x, ycoefh_x, ytke_x,y_treedrg)
+            ycoefm_x, ycoefh_x, ytke_x,yeps_x,y_treedrg)
 !            ycoefm_x, ycoefh_x, ytke_x)
 !FC doit on le mettre ( on ne l utilise pas si il y a du spliting)
@@ -2063 +2068 @@
         CALL call_atke(dtime,knon,klev,ycdragm_w(1:knon),ycdragh_w(1:knon),yus0(1:knon),yvs0(1:knon),yts_w(1:knon), &
                 yu_w(1:knon,:),yv_w(1:knon,:),yt_w(1:knon,:),yq_w(1:knon,:),ypplay(1:knon,:),ypaprs(1:knon,:),      &
-                ytke_w(1:knon,:),ycoefm_w(1:knon,:), ycoefh_w(1:knon,:))
+                ytke_w(1:knon,:),yeps_w(1:knon,:),ycoefm_w(1:knon,:),ycoefh_w(1:knon,:))
 
         ELSE
-
 
         CALL coef_diff_turb(dtime, nsrf, knon, ni,  &
             ypaprs, ypplay, yu_w, yv_w, yq_w, yt_w, yts_w, yqsurf_w, ycdragm_w, &
-            ycoefm_w, ycoefh_w, ytke_w,y_treedrg)
+            ycoefm_w, ycoefh_w, ytke_w,yeps_w,y_treedrg)
 !            ycoefm_w, ycoefh_w, ytke_w)
        IF (iflag_pbl>=20.AND.iflag_pbl<30) THEN
@@ -3332 +3336 @@
 !!              tke(i,k,is_ave) = tke(i,k,is_ave) + ytke(j,k)*ypct(j)
               tke_x(i,k,nsrf)    = ytke(j,k)
-              tke_x(i,k,is_ave) = tke_x(i,k,is_ave) + ytke(j,k)*ypct(j)
-
+              tke_x(i,k,is_ave)  = tke_x(i,k,is_ave) + ytke(j,k)*ypct(j)
+              eps_x(i,k,nsrf)    = yeps(j,k)
+              eps_x(i,k,is_ave)  = eps_x(i,k,is_ave) + yeps(j,k)*ypct(j)
 !>jyg
            ENDDO
@@ -3347 +3352 @@
 !!            tke(i,k,is_ave) = tke(i,k,is_ave) + tke(i,k,nsrf)*ypct(j)
             tke_x(i,k,nsrf)   = ytke_x(j,k)
-            tke_x(i,k,is_ave)   = tke_x(i,k,is_ave) + tke_x(i,k,nsrf)*ypct(j)        
+            tke_x(i,k,is_ave)   = tke_x(i,k,is_ave) + tke_x(i,k,nsrf)*ypct(j)       
+            eps_x(i,k,nsrf)   = yeps_x(j,k)
+            eps_x(i,k,is_ave)   = eps_x(i,k,is_ave) + eps_x(i,k,nsrf)*ypct(j) 
             wake_dltke(i,k,is_ave)   = wake_dltke(i,k,is_ave) + wake_dltke(i,k,nsrf)*ypct(j)