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

LMDZ6/trunk/libf/phylmd/lmdz_call_atke.F90

@@ -10 +10 @@
 subroutine call_atke(dtime,ngrid,nlay,cdrag_uv,cdrag_t,u_surf,v_surf,temp_surf, &
                         wind_u,wind_v,temp,qvap,play,pinterf, &
-                        tke,Km_out,Kh_out)
+                        tke,eps,Km_out,Kh_out)
 
 
@@ -43 +43 @@
 
 
-REAL, DIMENSION(ngrid,nlay+1), INTENT(INOUT)  :: tke  ! turbulent kinetic energy at interface between layers
+REAL, DIMENSION(ngrid,nlay+1), INTENT(INOUT)  :: tke      ! turbulent kinetic energy at interface between layers
 
+REAL, DIMENSION(ngrid,nlay+1), INTENT(OUT)    :: eps      ! output: tke dissipation rate at interface between layers
 REAL, DIMENSION(ngrid,nlay), INTENT(OUT)      :: Km_out   ! output: Exchange coefficient for momentum at interface between layers
 REAL, DIMENSION(ngrid,nlay), INTENT(OUT)      :: Kh_out   ! output: Exchange coefficient for heat flux at interface between layers
@@ -56 +57 @@
 call atke_compute_km_kh(ngrid,nlay,dtime,&
                         wind_u,wind_v,temp,qvap,play,pinterf,cdrag_uv,&
-                        tke,Km_out,Kh_out)
+                        tke,eps,Km_out,Kh_out)
 
 
@@ -75 +76 @@
    call atke_compute_km_kh(ngrid,nlay,dtime,&
                         wind_u_predict,wind_v_predict,temp,qvap,play,pinterf,cdrag_uv, &
-                        tke,Km_out,Kh_out)
+                        tke,eps,Km_out,Kh_out)
 
 end if