Changeset 2824


Ignore:
Timestamp:
Mar 15, 2017, 12:15:55 PM (7 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
Files:
12 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/DefLists/field_def_lmdz.xml

    r2823 r2824  
    208208        <field id="plfc"    long_name="Level of Free Convection"    unit="hPa" />
    209209        <field id="wbeff"    long_name="Conv. updraft velocity at LFC (inf100)"    unit="m/s" />
     210        <field id="convoccur"    long_name="Convective occurence"    unit="-" />
    210211        <field id="prw"    long_name="Precipitable water"    unit="kg/m2" />
    211212        <field id="prlw"   long_name="Precipitable liquid water"  unit="kg/m2" />
  • LMDZ5/trunk/DefLists/file_def_histLES_lmdz.xml

    r2823 r2824  
    195195                <field field_ref="plfc" level="10" />
    196196                <field field_ref="wbeff" level="10" />
     197                <field field_ref="convoccur" level="10" />
    197198                <field field_ref="prw" level="10" />
    198199                <field field_ref="prlw" level="10" />
  • LMDZ5/trunk/DefLists/file_def_histday_lmdz.xml

    r2823 r2824  
    195195                <field field_ref="plfc" level="10" />
    196196                <field field_ref="wbeff" level="10" />
     197                <field field_ref="convoccur" level="10" />
    197198                <field field_ref="prw" level="1" />
    198199                <field field_ref="prlw" level="1" />
  • LMDZ5/trunk/DefLists/file_def_histhf_lmdz.xml

    r2823 r2824  
    195195                <field field_ref="plfc" level="10" />
    196196                <field field_ref="wbeff" level="10" />
     197                <field field_ref="convoccur" level="10" />
    197198                <field field_ref="prw" level="10" />
    198199                <field field_ref="prlw" level="10" />
  • LMDZ5/trunk/DefLists/file_def_histins_lmdz.xml

    r2823 r2824  
    195195                <field field_ref="plfc" level="10" />
    196196                <field field_ref="wbeff" level="10" />
     197                <field field_ref="convoccur" level="10" />
    197198                <field field_ref="prw" level="10" />
    198199                <field field_ref="prlw" level="10" />
  • LMDZ5/trunk/DefLists/file_def_histmth_lmdz.xml

    r2823 r2824  
    215215                <field field_ref="plfc" level="1" />
    216216                <field field_ref="wbeff" level="1" />
     217                <field field_ref="convoccur" level="10" />
    217218                <field field_ref="prw" level="1" />
    218219                <field field_ref="prlw" level="1" />
  • LMDZ5/trunk/DefLists/file_def_histstn_lmdz.xml

    r2823 r2824  
    195195                <field field_ref="plfc" level="10" />
    196196                <field field_ref="wbeff" level="10" />
     197                <field field_ref="convoccur" level="10" />
    197198                <field field_ref="prw" level="10" />
    198199                <field field_ref="prlw" level="10" />
  • LMDZ5/trunk/libf/phylmd/concvl.F90

    r2481 r2824  
    55                  d_t, d_q, d_u, d_v, d_tra, &
    66                  rain, snow, kbas, ktop, sigd, &
    7                   cbmf, plcl, plfc, wbeff, upwd, dnwd, dnwdbis, &
     7                  cbmf, plcl, plfc, wbeff, convoccur, &
     8                  upwd, dnwd, dnwdbis, &
    89                  Ma, mip, Vprecip, &
    910                  cape, cin, tvp, Tconv, iflag, &
     
    157158  REAL qs(klon, klev), qs_wake(klon, klev)
    158159  REAL cbmf(klon), plcl(klon), plfc(klon), wbeff(klon)
     160  REAL convoccur(klon)
    159161!LF          SAVE cbmf
    160162!IM/JYG      REAL, SAVE, ALLOCATABLE :: cbmf(:)
     
    228230    ALLOCATE (t1(klon,klev))
    229231    ALLOCATE (q1(klon,klev))
     232!
     233    convoccur(:) = 0.
     234!
    230235    itap = 0
    231236    igout = klon/2 + 1/klon
     
    450455  END DO
    451456
     457  IF (iflag_con==3) THEN
     458    DO i = 1,klon
     459      IF (wbeff(i) > 100. .OR. wbeff(i) == 0 .OR. iflag(i) > 3) THEN
     460        wbeff(i) = 0.
     461        convoccur(i) = 0. 
     462      ELSE
     463        convoccur(i) = 1.
     464      ENDIF
     465    ENDDO
     466  ENDIF
     467
    452468  IF (iflag_con==30) THEN
    453469    DO itra = 1, ntra
  • LMDZ5/trunk/libf/phylmd/phys_local_var_mod.F90

    r2823 r2824  
    348348!>jyg+nrlmd
    349349  !
    350       REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: wbeff, zmax_th, zq2m, zt2m
    351 !$OMP THREADPRIVATE(wbeff, zmax_th, zq2m, zt2m)
     350      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: wbeff, convoccur, zmax_th, zq2m, zt2m
     351!$OMP THREADPRIVATE(wbeff, convoccur, zmax_th, zq2m, zt2m)
    352352      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zt2m_min_mon, zt2m_max_mon
    353353!$OMP THREADPRIVATE(zt2m_min_mon, zt2m_max_mon)
     
    673673      ALLOCATE(kh(klon), kh_x(klon), kh_w(klon))
    674674!
    675       ALLOCATE(wbeff(klon), zmax_th(klon))
     675      ALLOCATE(wbeff(klon), convoccur(klon), zmax_th(klon))
    676676      ALLOCATE(zq2m(klon), zt2m(klon), weak_inversion(klon))
    677677      ALLOCATE(zt2m_min_mon(klon), zt2m_max_mon(klon))
     
    930930      DEALLOCATE(kh, kh_x, kh_w)
    931931!
    932       DEALLOCATE(wbeff, zmax_th)
     932      DEALLOCATE(wbeff, convoccur, zmax_th)
    933933      DEALLOCATE(zq2m, zt2m, weak_inversion)
    934934      DEALLOCATE(zt2m_min_mon, zt2m_max_mon)
  • LMDZ5/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

    r2823 r2824  
    669669  TYPE(ctrl_out), SAVE :: o_wbeff = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
    670670    'wbeff', 'Conv. updraft velocity at LFC (<100)', 'm/s', (/ ('', i=1, 10) /))
     671  TYPE(ctrl_out), SAVE :: o_convoccur = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     672    'convoccur', 'Convective occurence', '', (/ ('', i=1, 10) /))
    671673  TYPE(ctrl_out), SAVE :: o_prw = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11, 11/), &
    672674    'prw', 'Precipitable water', 'kg/m2', (/ ('', i=1, 10) /))
  • LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

    r2823 r2824  
    6868         o_ue, o_ve, o_uq, o_vq, o_cape, o_pbase, &
    6969         o_ptop, o_fbase, o_plcl, o_plfc, &
    70          o_wbeff, o_cape_max, o_upwd, o_ep,o_epmax_diag, o_Ma, &
     70         o_wbeff, o_convoccur, o_cape_max, o_upwd, o_ep,o_epmax_diag, o_Ma, &
    7171         o_dnwd, o_dnwd0, o_ftime_con, o_mc, &
    7272         o_prw, o_prlw, o_prsw, o_s_pblh, o_s_pblt, o_s_lcl, &
     
    236236         cldh, cldt, JrNt, cldljn, cldmjn, cldhjn, &
    237237         cldtjn, cldq, flwp, fiwp, ue, ve, uq, vq, &
    238          plcl, plfc, wbeff, upwd, dnwd, dnwd0, prw, prlw, prsw, &
     238         plcl, plfc, wbeff, convoccur, upwd, dnwd, dnwd0, prw, prlw, prsw, &
    239239         s_pblh, s_pblt, s_lcl, s_therm, uwriteSTD, &
    240240         vwriteSTD, wwriteSTD, phiwriteSTD, qwriteSTD, &
     
    887887             CALL histwrite_phy(o_plfc, plfc)
    888888             CALL histwrite_phy(o_wbeff, wbeff)
     889             CALL histwrite_phy(o_convoccur, convoccur)
    889890          ENDIF
    890891
  • LMDZ5/trunk/libf/phylmd/physiq_mod.F90

    r2820 r2824  
    155155!!!       d_s_the, d_dens_the, &            ! due to thermals
    156156       !                                 
    157        wbeff, zmax_th, &
     157       wbeff, convoccur, zmax_th, &
    158158       sens, flwp, fiwp,  &
    159159       ale_bl_stat,alp_bl_conv,alp_bl_det,  &
     
    24842484               d_t_con,d_q_con,d_u_con,d_v_con,d_tr, &
    24852485               rain_con, snow_con, ibas_con, itop_con, sigd, &
    2486                ema_cbmf,plcl,plfc,wbeff,upwd,dnwd,dnwd0, &
     2486               ema_cbmf,plcl,plfc,wbeff,convoccur,upwd,dnwd,dnwd0, &
    24872487               Ma,mip,Vprecip,cape,cin,tvp,Tconv,iflagctrl, &
    24882488               pbase,bbase,dtvpdt1,dtvpdq1,dplcldt,dplcldr,qcondc,wd, &
Note: See TracChangeset for help on using the changeset viewer.