Changeset 1816 for LMDZ5/trunk/libf/phylmd/hbtm.F

Timestamp:

Jul 25, 2013, 10:33:44 AM (12 years ago)

Author:

idelkadi

Message:

Inclusion of the convective scale velocity w* for tracers
Concerns : hbtm.F, pbl_surface_mod.F90,
physiq.F, phys_output_ctrlout_mod.F90, phys_output_write_mod.F90

File:

• LMDZ5/trunk/libf/phylmd/hbtm.F (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/hbtm.F

@@ -4 +4 @@
 
       SUBROUTINE HBTM(knon, paprs, pplay,
-     .                t2m,t10m,q2m,q10m,ustar,
+     .                t2m,t10m,q2m,q10m,ustar,wstar,
      .                flux_t,flux_q,u,v,t,q,
      .                pblh,cape,EauLiq,ctei,pblT,
@@ -54 +54 @@
       REAL q2m(klon), q10m(klon) ! q a 2 et 10m
       REAL ustar(klon)
+      REAL wstar(klon)       ! w*, convective velocity scale
       REAL paprs(klon,klev+1) ! pression a inter-couche (Pa)
       REAL pplay(klon,klev)   ! pression au milieu de couche (Pa)
@@ -158 +159 @@
       REAL fak1(klon)       ! k*ustar*pblh
       REAL fak2(klon)       ! k*wm*pblh
-      REAL fak3(klon)       ! fakn*wstr/wm
+      REAL fak3(klon)       ! fakn*wstar/wm
       REAL pblk(klon)       ! level eddy diffusivity for momentum
       REAL pr(klon)         ! Prandtl number for eddy diffusivities
@@ -164 +165 @@
       REAL zh(klon)         ! zmzp / pblh
       REAL zzh(klon)        ! (1-(zmzp/pblh))**2
-      REAL wstr(klon)       ! w*, convective velocity scale
       REAL zm(klon)         ! current level height
       REAL zp(klon)         ! current level height + one level up
@@ -625 +625 @@
           wm(i)      = ustar(i)*phiminv(i)
           fak2(i)    = wm(i)*pblh(i)*vk
-          wstr(i)    = (heatv(i)*RG*pblh(i)/zxt)**onet
-          fak3(i)    = fakn*wstr(i)/wm(i)
+          wstar(i)    = (heatv(i)*RG*pblh(i)/zxt)**onet
+          fak3(i)    = fakn*wstar(i)/wm(i)
         ENDIF
 c Computes Theta_e for thermal (all cases : to be modified)