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

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -282 +282 @@
        zxfluxlat_x, zxfluxlat_w, &
        !
-       pbl_tke_input, tke_dissip, l_mix, wprime,&
+       pbl_tke_input, pbl_eps, l_mix, wprime,&
        t_therm, q_therm, u_therm, v_therm, &
        cdragh_x, cdragh_w, &
@@ -1257 +1257 @@
     integer iostat
 
-    REAL, dimension(klon,klev+1) :: tke_dissip_ave, l_mix_ave, wprime_ave
+    REAL, dimension(klon,klev+1) :: l_mix_ave, wprime_ave
     REAL zzz
     !albedo SB >>>
@@ -2905 +2905 @@
             fluxt,   fluxu,  fluxv, &
             dsens,     devap,     zxsnow, &
-            zxfluxt,   zxfluxq,  zxfluxqbs,  q2m, fluxq, fluxqbs, pbl_tke, &
+            zxfluxt,   zxfluxq,  zxfluxqbs,  q2m, fluxq, fluxqbs, pbl_tke, pbl_eps,  &
                                 !nrlmd+jyg<
             wake_delta_pbl_TKE, &
@@ -3849 +3849 @@
     ! water distribution
 
-    tke_dissip_ave(:,:)=0.
     l_mix_ave(:,:)=0.
     wprime_ave(:,:)=0.
@@ -3855 +3854 @@
     DO nsrf = 1, nbsrf
        DO i = 1, klon
-          tke_dissip_ave(i,:) = tke_dissip_ave(i,:) + tke_dissip(i,:,nsrf)*pctsrf(i,nsrf)
           l_mix_ave(i,:) = l_mix_ave(i,:) + l_mix(i,:,nsrf)*pctsrf(i,nsrf)
           wprime_ave(i,:) = wprime_ave(i,:) + wprime(i,:,nsrf)*pctsrf(i,nsrf)
@@ -3869 +3867 @@
          paprs,pplay,t_seri,q_seri, &
          qtc_cv, sigt_cv,detrain_cv,fm_cv,fqd,fqcomp,sigd,zqsat, &
-         omega,pbl_tke(:,:,is_ave),tke_dissip_ave,l_mix_ave,wprime_ave, &
+         omega,pbl_tke(:,:,is_ave),pbl_eps(:,:,is_ave),l_mix_ave,wprime_ave, &
          t2m,q2m,fm_therm,entr_therm,detr_therm,cell_area, &
          ratqs,ratqsc,ratqs_inter_)