Changeset 2824 for LMDZ5/trunk/libf/phylmd

Timestamp:

Mar 15, 2017, 12:15:55 PM (8 years ago)

Author:

Laurent Fairhead

Message:

Introduction of convoccur variable and modification of the wbeff variable
in the convection so that wbeff/convoccur is the mean value of wbeff when
convection occurs
Jyg

Location:

LMDZ5/trunk/libf/phylmd

Files:

• concvl.F90 (modified) (4 diffs)
• phys_local_var_mod.F90 (modified) (3 diffs)
• phys_output_ctrlout_mod.F90 (modified) (1 diff)
• phys_output_write_mod.F90 (modified) (3 diffs)
• physiq_mod.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/concvl.F90

@@ -5 +5 @@
                   d_t, d_q, d_u, d_v, d_tra, &
                   rain, snow, kbas, ktop, sigd, &
-                  cbmf, plcl, plfc, wbeff, upwd, dnwd, dnwdbis, &
+                  cbmf, plcl, plfc, wbeff, convoccur, &
+                  upwd, dnwd, dnwdbis, &
                   Ma, mip, Vprecip, &
                   cape, cin, tvp, Tconv, iflag, &
@@ -157 +158 @@
   REAL qs(klon, klev), qs_wake(klon, klev)
   REAL cbmf(klon), plcl(klon), plfc(klon), wbeff(klon)
+  REAL convoccur(klon)
 !LF          SAVE cbmf
 !IM/JYG      REAL, SAVE, ALLOCATABLE :: cbmf(:)
@@ -228 +230 @@
     ALLOCATE (t1(klon,klev))
     ALLOCATE (q1(klon,klev))
+!
+    convoccur(:) = 0.
+!
     itap = 0
     igout = klon/2 + 1/klon
@@ -450 +455 @@
   END DO
 
+  IF (iflag_con==3) THEN
+    DO i = 1,klon
+      IF (wbeff(i) > 100. .OR. wbeff(i) == 0 .OR. iflag(i) > 3) THEN
+        wbeff(i) = 0.
+        convoccur(i) = 0.  
+      ELSE
+        convoccur(i) = 1.
+      ENDIF
+    ENDDO
+  ENDIF
+
   IF (iflag_con==30) THEN
     DO itra = 1, ntra

LMDZ5/trunk/libf/phylmd/phys_local_var_mod.F90

@@ -348 +348 @@
 !>jyg+nrlmd
   !
-      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: wbeff, zmax_th, zq2m, zt2m
-!$OMP THREADPRIVATE(wbeff, zmax_th, zq2m, zt2m)
+      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: wbeff, convoccur, zmax_th, zq2m, zt2m
+!$OMP THREADPRIVATE(wbeff, convoccur, zmax_th, zq2m, zt2m)
       REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zt2m_min_mon, zt2m_max_mon
 !$OMP THREADPRIVATE(zt2m_min_mon, zt2m_max_mon)
@@ -673 +673 @@
       ALLOCATE(kh(klon), kh_x(klon), kh_w(klon))
 !
-      ALLOCATE(wbeff(klon), zmax_th(klon))
+      ALLOCATE(wbeff(klon), convoccur(klon), zmax_th(klon))
       ALLOCATE(zq2m(klon), zt2m(klon), weak_inversion(klon))
       ALLOCATE(zt2m_min_mon(klon), zt2m_max_mon(klon))
@@ -930 +930 @@
       DEALLOCATE(kh, kh_x, kh_w)
 !
-      DEALLOCATE(wbeff, zmax_th)
+      DEALLOCATE(wbeff, convoccur, zmax_th)
       DEALLOCATE(zq2m, zt2m, weak_inversion)
       DEALLOCATE(zt2m_min_mon, zt2m_max_mon)

LMDZ5/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -669 +669 @@
   TYPE(ctrl_out), SAVE :: o_wbeff = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'wbeff', 'Conv. updraft velocity at LFC (<100)', 'm/s', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_convoccur = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'convoccur', 'Convective occurence', '', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_prw = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'prw', 'Precipitable water', 'kg/m2', (/ ('', i=1, 10) /))

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -68 +68 @@
          o_ue, o_ve, o_uq, o_vq, o_cape, o_pbase, &
          o_ptop, o_fbase, o_plcl, o_plfc, &
-         o_wbeff, o_cape_max, o_upwd, o_ep,o_epmax_diag, o_Ma, &
+         o_wbeff, o_convoccur, o_cape_max, o_upwd, o_ep,o_epmax_diag, o_Ma, &
          o_dnwd, o_dnwd0, o_ftime_con, o_mc, &
          o_prw, o_prlw, o_prsw, o_s_pblh, o_s_pblt, o_s_lcl, &
@@ -236 +236 @@
          cldh, cldt, JrNt, cldljn, cldmjn, cldhjn, &
          cldtjn, cldq, flwp, fiwp, ue, ve, uq, vq, &
-         plcl, plfc, wbeff, upwd, dnwd, dnwd0, prw, prlw, prsw, &
+         plcl, plfc, wbeff, convoccur, upwd, dnwd, dnwd0, prw, prlw, prsw, &
          s_pblh, s_pblt, s_lcl, s_therm, uwriteSTD, &
          vwriteSTD, wwriteSTD, phiwriteSTD, qwriteSTD, &
@@ -887 +887 @@
              CALL histwrite_phy(o_plfc, plfc)
              CALL histwrite_phy(o_wbeff, wbeff)
+             CALL histwrite_phy(o_convoccur, convoccur)
           ENDIF
 

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -155 +155 @@
 !!!       d_s_the, d_dens_the, &            ! due to thermals
        !                                  
-       wbeff, zmax_th, &
+       wbeff, convoccur, zmax_th, &
        sens, flwp, fiwp,  &
        ale_bl_stat,alp_bl_conv,alp_bl_det,  &
@@ -2484 +2484 @@
                d_t_con,d_q_con,d_u_con,d_v_con,d_tr, &
                rain_con, snow_con, ibas_con, itop_con, sigd, &
-               ema_cbmf,plcl,plfc,wbeff,upwd,dnwd,dnwd0, &
+               ema_cbmf,plcl,plfc,wbeff,convoccur,upwd,dnwd,dnwd0, &
                Ma,mip,Vprecip,cape,cin,tvp,Tconv,iflagctrl, &
                pbase,bbase,dtvpdt1,dtvpdq1,dplcldt,dplcldr,qcondc,wd, &