Changeset 1761 for LMDZ5/trunk/libf/phylmd/physiq.F

Timestamp:

Jun 4, 2013, 12:11:38 PM (11 years ago)

Author:

Laurent Fairhead

Message:

New version of Mellor et Yamada pronostic TKE

... based on energy transfer from the mean state.

The new version is yamada_c.
It must be called after vertical diffusion rather than just before
since the source terms u_z w'u' + v_z w'theta' and g/theta w'theta'
are diagnosed from the vertical diffusion (as energy loss from the mean
state) rather than computed as K (u_z^2+v_z2) or g/\theta K theta_z
(where _z means vertical derivative).
The call to this version is controled by iflag_pbl.

iflag_pbl = 20 : with TKE computed at interfaces between layers
iflag_pbl = 25 : with TKE computed within the layer
In both cases, the dissipation of turbulence is translated into heat, and
passed to the physics as dtdiss (temperature tendency due to dissipation

of TKE).

The diffusion coefficient being computed after dissipation, it must be
kept for diffusion at the next time step, and thus be stored in the
restartphy file.
This coefficient must be computed and stored for each sub-surface.

A new way of managing subsurface variables is introduced.
For arrays of the form X(:,nbsrf) are extented to X(:,nbsrf+1), where
is_ave=nbsrf+1, is an additional sub-surface which contains the averaged values.

coef_diff_turb_mod.F90 : change of flags.
ener_conserv.F90 : energy conservation must not be applied in those

cases

indicesol.h : definition of is_ave
pbl_surface_mod.F90 : call to yamada_c and changes in the management of

coefh/coefm

physiq.F : Change in the initialisation of pmflxr/s and

modified calls to pbl_surface_mod (introduction
of dtdiss, initialisation
of pbl_tke and coefh in 1D).

phys_local_var_mod.F90 : declaration of d_t_diss
phys_output_mod.F90 : new outputs (bils_tke and bils_diss) and

coefh->coefh(:,:,is_ave)

phys_state_var_mod.F90 : modified declaration for coefh and coefm

(nbsrf -> nbsrf+1)

FH

File:

• LMDZ5/trunk/libf/phylmd/physiq.F (modified) (15 diffs)

LMDZ5/trunk/libf/phylmd/physiq.F

@@ -655 +655 @@
       real w0(klon)                                          ! Vitesse des thermiques au LCL
       real w_conv(klon)                                      ! Vitesse verticale de grande \'echelle au LCL
+      real tke0(klon,klev+1)                              ! TKE au début du pas de temps
       real therm_tke_max0(klon)                              ! TKE dans les thermiques au LCL 
       real env_tke_max0(klon)                                ! TKE dans l'environnement au LCL 
@@ -694 +695 @@
 cAA
 cAA  Pour phytrac 
-cAA
-      REAL coefh(klon,klev)     ! coef d'echange pour phytrac, valable pour 2<=k<=klev
-      REAL coefm(klon,klev)     ! coef d'echange pour U, V
       REAL u1(klon)             ! vents dans la premiere couche U
       REAL v1(klon)             ! vents dans la premiere couche V
@@ -1288 +1286 @@
           pbase=0
 cIM 180608
-c         pmflxr=0.
-c         pmflxs=0.
 
         itau_con=0
@@ -1396 +1392 @@
 
          CALL phyetat0 ("startphy.nc",clesphy0,tabcntr0)
+         IF (klon_glo==1) THEN
+         coefh=0. ; coefm=0. ; pbl_tke=0.
+         coefh(:,2,:)=1.e-2 ; coefm(:,2,:)=1.e-2 ; pbl_tke(:,2,:)=1.e-2
+         PRINT*,'FH WARNING : lignes a supprimer'
+         ENDIF
 cIM begin
           print*,'physiq: clwcon rnebcon ratqs',clwcon(1,1),rnebcon(1,1)
@@ -1770 +1771 @@
       d1a(:,:)=0.
       dam(:,:)=0.
+          pmflxr=0.
+          pmflxs=0.
 ! RomP <<<
 
@@ -1785 +1788 @@
       ENDDO
       ENDDO
+      tke0(:,:)=pbl_tke(:,:,is_ave)
       IF (nqtot.GE.3) THEN
       DO iq = 3, nqtot
@@ -2073 +2077 @@
      s     albsol1,   albsol2,   sens,    evap,  
      s     zxtsol,    zxfluxlat, zt2m,    qsat2m, 
-     s     d_t_vdf,   d_q_vdf,   d_u_vdf, d_v_vdf,
+     s     d_t_vdf,   d_q_vdf,   d_u_vdf, d_v_vdf, d_t_diss,
      s     coefh,     coefm,     slab_wfbils,                
      d     qsol,      zq2m,      s_pblh,  s_lcl,
@@ -2089 +2093 @@
 !-----------------------------------------------------------------------------------------
 ! ajout des tendances de la diffusion turbulente
-      CALL add_phys_tend(d_u_vdf,d_v_vdf,d_t_vdf,d_q_vdf,dql0,'vdf')
+      CALL add_phys_tend
+     s     (d_u_vdf,d_v_vdf,d_t_vdf+d_t_diss,d_q_vdf,dql0,'vdf')
 !-----------------------------------------------------------------------------------------
 
@@ -3175 +3180 @@
      $                          paprs,
      $                          pplay,
-     $                          coefh,
+     $                          coefh(:,:,is_ave),
      $                          pphi,
      $                          t_seri,
@@ -3696 +3701 @@
      I     paprs,    pplay,     pmfu,     pmfd,
      I     pen_u,    pde_u,     pen_d,    pde_d,
-     I     cdragh,   coefh,     fm_therm, entr_therm,
+     I     cdragh,   coefh(:,:,is_ave),     fm_therm, entr_therm,
      I     u1,       v1,        ftsol,    pctsrf,
      I     ustar,     u10m,      v10m,
@@ -3723 +3728 @@
      I                   t,pmfu, pmfd, pen_u, pde_u, pen_d, pde_d,
      I                   fm_therm,entr_therm,
-     I                   cdragh,coefh,u1,v1,ftsol,pctsrf,
+     I                   cdragh,coefh(:,:,is_ave),u1,v1,ftsol,pctsrf,
      I                   frac_impa, frac_nucl,
      I                   pphis,airephy,dtime,itap,
@@ -3758 +3763 @@
       forall (k=1: llm) exner(:, k) = (pplay(:, k)/paprs(:,1))**RKAPPA
       CALL ener_conserv(klon,klev,pdtphys,u,v,t,qx(:,:,ivap),
-     s        u_seri,v_seri,t_seri,q_seri,
+     s        u_seri,v_seri,t_seri,q_seri,pbl_tke(:,:,is_ave)-tke0(:,:),
      s        zmasse,exner,d_t_ec)
       t_seri(:,:)=t_seri(:,:)+d_t_ec(:,:)
@@ -3963 +3968 @@
 
 
+
 c=============================================================
 ! Separation entre thermiques et non thermiques dans les sorties
@@ -3991 +3997 @@
       endif
 
- 
 #include "phys_output_write.h"
 
@@ -4070 +4075 @@
          ENDIF !if callstats
      
-
       IF (lafin) THEN
          itau_phy = itau_phy + itap