Sources/phylmd/cv_driver.f

module cv_driver_m

  implicit none

contains

  SUBROUTINE cv_driver(len, nd, t1, q1, qs1, u1, v1, p1, ph1, iflag1, ft1, &
       fq1, fu1, fv1, precip1, VPrecip1, cbmf1, sig1, w01, icb1, inb1, delt, &
       Ma1, upwd1, dnwd1, dnwd01, qcondc1, wd1, cape1, da1, phi1, mp1)

    ! From LMDZ4/libf/phylmd/cv_driver.F, version 1.3, 2005/04/15 12:36:17
    ! Main driver for convection
    ! Author: S. Bony, March 2002

    ! Several modules corresponding to different physical processes

    ! Several versions of convect may be used:
    ! - iflag_con = 3: version lmd
    ! - iflag_con = 4: version 4.3b

    use clesphys2, only: iflag_con
    use cv3_compress_m, only: cv3_compress
    use cv3_mixing_m, only: cv3_mixing
    use cv3_param_m, only: cv3_param
    use cv3_prelim_m, only: cv3_prelim
    use cv3_tracer_m, only: cv3_tracer
    use cv3_uncompress_m, only: cv3_uncompress
    use cv3_unsat_m, only: cv3_unsat
    use cv3_yield_m, only: cv3_yield
    use cv_uncompress_m, only: cv_uncompress
    USE dimphy, ONLY: klev, klon

    integer, intent(in):: len ! first dimension
    integer, intent(in):: nd ! vertical dimension
    real, intent(in):: t1(len, nd) ! temperature
    real q1(len, nd) !           Input        specific hum
    real qs1(len, nd)
    !      qs1           Real           Input        sat specific hum
    real, intent(in):: u1(len, nd)
    !      u1            Real           Input        u-wind
    real, intent(in):: v1(len, nd)
    !      v1            Real           Input        v-wind
    real p1(len, nd)
    !      p1            Real           Input        full level pressure
    real ph1(len, nd + 1)
    !      ph1           Real           Input        half level pressure
    integer iflag1(len)
    !      iflag1        Integer        Output       flag for Emanuel conditions
    real ft1(len, nd)
    !      ft1           Real           Output       temp tend
    real fq1(len, nd)
    !      fq1           Real           Output       spec hum tend
    real fu1(len, nd)
    !      fu1           Real           Output       u-wind tend
    real fv1(len, nd)
    !      fv1           Real           Output       v-wind tend
    real precip1(len)
    !      precip1       Real           Output       precipitation
    real VPrecip1(len, nd+1)
    !      VPrecip1      Real           Output       vertical profile of precipitations
    real cbmf1(len)
    !      cbmf1         Real           Output       cloud base mass flux
    real, intent(inout):: sig1(klon, klev) ! section adiabatic updraft

    real, intent(inout):: w01(klon, klev) 
    ! vertical velocity within adiabatic updraft

    integer icb1(klon)
    integer inb1(klon)
    real, intent(in):: delt
    !      delt          Real           Input        time step
    real Ma1(len, nd)
    !      Ma1           Real           Output       mass flux adiabatic updraft
    real, intent(out):: upwd1(len, nd) ! total upward mass flux (adiab+mixed)
    real, intent(out):: dnwd1(len, nd) ! saturated downward mass flux (mixed)
    real, intent(out):: dnwd01(len, nd) ! unsaturated downward mass flux

    real qcondc1(len, nd)     ! cld
    !      qcondc1       Real           Output       in-cld mixing ratio of condensed water
    real wd1(len)            ! gust
    !      wd1           Real           Output       downdraft velocity scale for sfc fluxes
    real cape1(len)
    !      cape1         Real           Output       CAPE

    real, intent(inout):: da1(len, nd), phi1(len, nd, nd), mp1(len, nd)

    !-------------------------------------------------------------------
    ! --- ARGUMENTS
    !-------------------------------------------------------------------
    ! --- On input:

    !  t:   Array of absolute temperature (K) of dimension ND, with first
    !       index corresponding to lowest model level. Note that this array
    !       will be altered by the subroutine if dry convective adjustment
    !       occurs and if IPBL is not equal to 0.

    !  q:   Array of specific humidity (gm/gm) of dimension ND, with first
    !       index corresponding to lowest model level. Must be defined
    !       at same grid levels as T. Note that this array will be altered
    !       if dry convective adjustment occurs and if IPBL is not equal to 0.

    !  qs:  Array of saturation specific humidity of dimension ND, with first
    !       index corresponding to lowest model level. Must be defined
    !       at same grid levels as T. Note that this array will be altered
    !       if dry convective adjustment occurs and if IPBL is not equal to 0.

    !  u:   Array of zonal wind velocity (m/s) of dimension ND, witth first
    !       index corresponding with the lowest model level. Defined at
    !       same levels as T. Note that this array will be altered if
    !       dry convective adjustment occurs and if IPBL is not equal to 0.

    !  v:   Same as u but for meridional velocity.

    !  p:   Array of pressure (mb) of dimension ND, with first
    !       index corresponding to lowest model level. Must be defined
    !       at same grid levels as T.

    !  ph:  Array of pressure (mb) of dimension ND+1, with first index
    !       corresponding to lowest level. These pressures are defined at
    !       levels intermediate between those of P, T, Q and QS. The first
    !       value of PH should be greater than (i.e. at a lower level than)
    !       the first value of the array P.

    !  nl:  The maximum number of levels to which convection can penetrate, plus 1.
    !       NL MUST be less than or equal to ND-1.

    !  delt: The model time step (sec) between calls to CONVECT

    !----------------------------------------------------------------------------
    ! ---   On Output:

    !  iflag: An output integer whose value denotes the following:
    !       VALUE   INTERPRETATION
    !       -----   --------------
    !         0     Moist convection occurs.
    !         1     Moist convection occurs, but a CFL condition
    !               on the subsidence warming is violated. This
    !               does not cause the scheme to terminate.
    !         2     Moist convection, but no precip because ep(inb) lt 0.0001
    !         3     No moist convection because new cbmf is 0 and old cbmf is 0.
    !         4     No moist convection; atmosphere is not
    !               unstable
    !         6     No moist convection because ihmin le minorig.
    !         7     No moist convection because unreasonable
    !               parcel level temperature or specific humidity.
    !         8     No moist convection: lifted condensation
    !               level is above the 200 mb level.
    !         9     No moist convection: cloud base is higher
    !               then the level NL-1.

    !  ft:   Array of temperature tendency (K/s) of dimension ND, defined at same
    !        grid levels as T, Q, QS and P.

    !  fq:   Array of specific humidity tendencies ((gm/gm)/s) of dimension ND,
    !        defined at same grid levels as T, Q, QS and P.

    !  fu:   Array of forcing of zonal velocity (m/s^2) of dimension ND,
    !        defined at same grid levels as T.

    !  fv:   Same as FU, but for forcing of meridional velocity.

    !  precip: Scalar convective precipitation rate (mm/day).

    !  VPrecip: Vertical profile of convective precipitation (kg/m2/s).

    !  wd:   A convective downdraft velocity scale. For use in surface
    !        flux parameterizations. See convect.ps file for details.

    !  tprime: A convective downdraft temperature perturbation scale (K).
    !          For use in surface flux parameterizations. See convect.ps
    !          file for details.

    !  qprime: A convective downdraft specific humidity
    !          perturbation scale (gm/gm).
    !          For use in surface flux parameterizations. See convect.ps
    !          file for details.

    !  cbmf: The cloud base mass flux ((kg/m**2)/s). THIS SCALAR VALUE MUST
    !        BE STORED BY THE CALLING PROGRAM AND RETURNED TO CONVECT AT
    !        ITS NEXT CALL. That is, the value of CBMF must be "remembered"
    !        by the calling program between calls to CONVECT.

    !  det:   Array of detrainment mass flux of dimension ND.

    !-------------------------------------------------------------------

    !  Local arrays

    real da(len, nd), phi(len, nd, nd), mp(len, nd)

    integer i, k, il
    integer icbmax
    integer nk1(klon)
    integer icbs1(klon)

    real plcl1(klon)
    real tnk1(klon)
    real qnk1(klon)
    real gznk1(klon)
    real pbase1(klon)
    real buoybase1(klon)

    real lv1(klon, klev)
    real cpn1(klon, klev)
    real tv1(klon, klev)
    real gz1(klon, klev)
    real hm1(klon, klev)
    real h1(klon, klev)
    real tp1(klon, klev)
    real tvp1(klon, klev)
    real clw1(klon, klev)
    real th1(klon, klev)

    integer ncum

    ! (local) compressed fields:

    integer nloc
    parameter (nloc = klon) ! pour l'instant

    integer idcum(nloc)
    integer iflag(nloc), nk(nloc), icb(nloc)
    integer nent(nloc, klev)
    integer icbs(nloc)
    integer inb(nloc), inbis(nloc)

    real cbmf(nloc), plcl(nloc), tnk(nloc), qnk(nloc), gznk(nloc)
    real t(nloc, klev), q(nloc, klev), qs(nloc, klev)
    real u(nloc, klev), v(nloc, klev)
    real gz(nloc, klev), h(nloc, klev), lv(nloc, klev), cpn(nloc, klev)
    real p(nloc, klev), ph(nloc, klev+1), tv(nloc, klev), tp(nloc, klev)
    real clw(nloc, klev)
    real dph(nloc, klev)
    real pbase(nloc), buoybase(nloc), th(nloc, klev)
    real tvp(nloc, klev)
    real sig(nloc, klev), w0(nloc, klev)
    real hp(nloc, klev), ep(nloc, klev), sigp(nloc, klev)
    real frac(nloc), buoy(nloc, klev)
    real cape(nloc)
    real m(nloc, klev), ment(nloc, klev, klev), qent(nloc, klev, klev)
    real uent(nloc, klev, klev), vent(nloc, klev, klev)
    real ments(nloc, klev, klev), qents(nloc, klev, klev)
    real sij(nloc, klev, klev), elij(nloc, klev, klev)
    real qp(nloc, klev), up(nloc, klev), vp(nloc, klev)
    real wt(nloc, klev), water(nloc, klev), evap(nloc, klev)
    real b(nloc, klev), ft(nloc, klev), fq(nloc, klev)
    real fu(nloc, klev), fv(nloc, klev)
    real upwd(nloc, klev), dnwd(nloc, klev), dnwd0(nloc, klev)
    real Ma(nloc, klev), mike(nloc, klev), tls(nloc, klev)
    real tps(nloc, klev), qprime(nloc), tprime(nloc)
    real precip(nloc)
    real VPrecip(nloc, klev+1)
    real qcondc(nloc, klev)  ! cld
    real wd(nloc)           ! gust

    !-------------------------------------------------------------------
    ! --- SET CONSTANTS AND PARAMETERS
    !-------------------------------------------------------------------

    ! -- set simulation flags:
    !   (common cvflag)

    CALL cv_flag

    ! -- set thermodynamical constants:
    !     (common cvthermo)

    CALL cv_thermo

    ! -- set convect parameters

    !     includes microphysical parameters and parameters that
    !     control the rate of approach to quasi-equilibrium)
    !     (common cvparam)

    if (iflag_con.eq.3) then
       CALL cv3_param(nd, delt)
    endif

    if (iflag_con.eq.4) then
       CALL cv_param(nd)
    endif

    !---------------------------------------------------------------------
    ! --- INITIALIZE OUTPUT ARRAYS AND PARAMETERS
    !---------------------------------------------------------------------

    do k = 1, nd
       do  i = 1, len
          ft1(i, k) = 0.0
          fq1(i, k) = 0.0
          fu1(i, k) = 0.0
          fv1(i, k) = 0.0
          tvp1(i, k) = 0.0
          tp1(i, k) = 0.0
          clw1(i, k) = 0.0
          !ym
          clw(i, k) = 0.0
          gz1(i, k)  =  0.
          VPrecip1(i, k) = 0.
          Ma1(i, k) = 0.0
          upwd1(i, k) = 0.0
          dnwd1(i, k) = 0.0
          dnwd01(i, k) = 0.0
          qcondc1(i, k) = 0.0
       end do
    end do

    do  i = 1, len
       precip1(i) = 0.0
       iflag1(i) = 0
       wd1(i) = 0.0
       cape1(i) = 0.0
       VPrecip1(i, nd+1) = 0.0
    end do

    if (iflag_con.eq.3) then
       do il = 1, len
          sig1(il, nd) = sig1(il, nd) + 1.
          sig1(il, nd)  =  min(sig1(il, nd), 12.1)
       enddo
    endif

    !--------------------------------------------------------------------
    ! --- CALCULATE ARRAYS OF GEOPOTENTIAL, HEAT CAPACITY & STATIC ENERGY
    !--------------------------------------------------------------------

    if (iflag_con.eq.3) then
       CALL cv3_prelim(len, nd, nd + 1, t1, q1, p1, ph1, lv1, cpn1, tv1, gz1, &
            h1, hm1, th1)
    endif

    if (iflag_con.eq.4) then
       CALL cv_prelim(len, nd, nd + 1, t1, q1, p1, ph1 &
            , lv1, cpn1, tv1, gz1, h1, hm1)
    endif

    !--------------------------------------------------------------------
    ! --- CONVECTIVE FEED
    !--------------------------------------------------------------------

    if (iflag_con.eq.3) then
       CALL cv3_feed(len, nd, t1, q1, qs1, p1, ph1, hm1, gz1            &
            , nk1, icb1, icbmax, iflag1, tnk1, qnk1, gznk1, plcl1) ! nd->na
    endif

    if (iflag_con.eq.4) then
       CALL cv_feed(len, nd, t1, q1, qs1, p1, hm1, gz1 &
            , nk1, icb1, icbmax, iflag1, tnk1, qnk1, gznk1, plcl1)
    endif

    !--------------------------------------------------------------------
    ! --- UNDILUTE (ADIABATIC) UPDRAFT / 1st part
    ! (up through ICB for convect4, up through ICB+1 for convect3)
    !     Calculates the lifted parcel virtual temperature at nk, the
    !     actual temperature, and the adiabatic liquid water content.
    !--------------------------------------------------------------------

    if (iflag_con.eq.3) then
       CALL cv3_undilute1(len, nd, t1, q1, qs1, gz1, plcl1, p1, nk1, icb1   &
            , tp1, tvp1, clw1, icbs1) ! nd->na
    endif

    if (iflag_con.eq.4) then
       CALL cv_undilute1(len, nd, t1, q1, qs1, gz1, p1, nk1, icb1, icbmax &
            , tp1, tvp1, clw1)
    endif

    !-------------------------------------------------------------------
    ! --- TRIGGERING
    !-------------------------------------------------------------------

    if (iflag_con.eq.3) then
       CALL cv3_trigger(len, nd, icb1, plcl1, p1, th1, tv1, tvp1, pbase1, &
            buoybase1, iflag1, sig1, w01) ! nd->na
    endif

    if (iflag_con.eq.4) then
       CALL cv_trigger(len, nd, icb1, cbmf1, tv1, tvp1, iflag1)
    endif

    ! --- IF THIS POINT IS REACHED, MOIST CONVECTIVE ADJUSTMENT IS NECESSARY

    ncum = 0
    do  i = 1, len
       if(iflag1(i).eq.0)then
          ncum = ncum+1
          idcum(ncum) = i
       endif
    end do

    IF (ncum.gt.0) THEN

       !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
       ! --- COMPRESS THE FIELDS
       !        (-> vectorization over convective gridpoints)
       !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

       if (iflag_con.eq.3) then
          CALL cv3_compress(len, nloc, ncum, nd, iflag1, nk1, icb1, icbs1, &
               plcl1, tnk1, qnk1, gznk1, pbase1, buoybase1, t1, q1, qs1, u1, &
               v1, gz1, th1, h1, lv1, cpn1, p1, ph1, tv1, tp1, tvp1, clw1, &
               sig1, w01, iflag, nk, icb, icbs, plcl, tnk, qnk, gznk, pbase, &
               buoybase, t, q, qs, u, v, gz, th, h, lv, cpn, p, ph, tv, tp, &
               tvp, clw, sig, w0)
       endif

       if (iflag_con.eq.4) then
          CALL cv_compress( len, nloc, ncum, nd &
               , iflag1, nk1, icb1 &
               , cbmf1, plcl1, tnk1, qnk1, gznk1 &
               , t1, q1, qs1, u1, v1, gz1 &
               , h1, lv1, cpn1, p1, ph1, tv1, tp1, tvp1, clw1 &
               , iflag, nk, icb &
               , cbmf, plcl, tnk, qnk, gznk &
               , t, q, qs, u, v, gz, h, lv, cpn, p, ph, tv, tp, tvp, clw  &
               , dph )
       endif

       !-------------------------------------------------------------------
       ! --- UNDILUTE (ADIABATIC) UPDRAFT / second part :
       ! ---   FIND THE REST OF THE LIFTED PARCEL TEMPERATURES
       ! ---   &
       ! ---   COMPUTE THE PRECIPITATION EFFICIENCIES AND THE
       ! ---   FRACTION OF PRECIPITATION FALLING OUTSIDE OF CLOUD
       ! ---   &
       ! ---   FIND THE LEVEL OF NEUTRAL BUOYANCY
       !-------------------------------------------------------------------

       if (iflag_con.eq.3) then
          CALL cv3_undilute2(nloc, ncum, nd, icb, icbs, nk         &
               , tnk, qnk, gznk, t, q, qs, gz &
               , p, h, tv, lv, pbase, buoybase, plcl &
               , inb, tp, tvp, clw, hp, ep, sigp, buoy) !na->nd
       endif

       if (iflag_con.eq.4) then
          CALL cv_undilute2(nloc, ncum, nd, icb, nk &
               , tnk, qnk, gznk, t, q, qs, gz &
               , p, dph, h, tv, lv &
               , inb, inbis, tp, tvp, clw, hp, ep, sigp, frac)
       endif

       !-------------------------------------------------------------------
       ! --- CLOSURE
       !-------------------------------------------------------------------

       if (iflag_con.eq.3) then
          CALL cv3_closure(nloc, ncum, nd, icb, inb               &
               , pbase, p, ph, tv, buoy &
               , sig, w0, cape, m) ! na->nd
       endif

       if (iflag_con.eq.4) then
          CALL cv_closure(nloc, ncum, nd, nk, icb &
               , tv, tvp, p, ph, dph, plcl, cpn &
               , iflag, cbmf)
       endif

       !-------------------------------------------------------------------
       ! --- MIXING
       !-------------------------------------------------------------------

       if (iflag_con.eq.3) then
          CALL cv3_mixing(nloc, ncum, nd, nd, icb, nk, inb, ph, t, q, &
               qs, u, v, h, lv, qnk, hp, tv, tvp, ep, clw, m, sig, ment, &
               qent, uent, vent, nent, sij, elij, ments, qents)
       endif

       if (iflag_con.eq.4) then
          CALL cv_mixing(nloc, ncum, nd, icb, nk, inb, inbis &
               , ph, t, q, qs, u, v, h, lv, qnk &
               , hp, tv, tvp, ep, clw, cbmf &
               , m, ment, qent, uent, vent, nent, sij, elij)
       endif

       !-------------------------------------------------------------------
       ! --- UNSATURATED (PRECIPITATING) DOWNDRAFTS
       !-------------------------------------------------------------------

       if (iflag_con.eq.3) then
          CALL cv3_unsat(nloc, ncum, nd, nd, icb, inb     &
               , t, q, qs, gz, u, v, p, ph &
               , th, tv, lv, cpn, ep, sigp, clw &
               , m, ment, elij, delt, plcl &
               , mp, qp, up, vp, wt, water, evap, b)! na->nd
       endif

       if (iflag_con.eq.4) then
          CALL cv_unsat(nloc, ncum, nd, inb, t, q, qs, gz, u, v, p, ph &
               , h, lv, ep, sigp, clw, m, ment, elij &
               , iflag, mp, qp, up, vp, wt, water, evap)
       endif

       !-------------------------------------------------------------------
       ! --- YIELD
       !     (tendencies, precipitation, variables of interface with other
       !      processes, etc)
       !-------------------------------------------------------------------

       if (iflag_con.eq.3) then
          CALL cv3_yield(nloc, ncum, nd, nd             &
               , icb, inb, delt &
               , t, q, u, v, gz, p, ph, h, hp, lv, cpn, th &
               , ep, clw, m, tp, mp, qp, up, vp &
               , wt, water, evap, b &
               , ment, qent, uent, vent, nent, elij, sig &
               , tv, tvp &
               , iflag, precip, VPrecip, ft, fq, fu, fv &
               , upwd, dnwd, dnwd0, ma, mike, tls, tps, qcondc, wd)! na->nd
       endif

       if (iflag_con.eq.4) then
          CALL cv_yield(nloc, ncum, nd, nk, icb, inb, delt &
               , t, q, u, v, gz, p, ph, h, hp, lv, cpn &
               , ep, clw, frac, m, mp, qp, up, vp &
               , wt, water, evap &
               , ment, qent, uent, vent, nent, elij &
               , tv, tvp &
               , iflag, wd, qprime, tprime &
               , precip, cbmf, ft, fq, fu, fv, Ma, qcondc)
       endif

       !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
       ! --- passive tracers
       !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

       if (iflag_con.eq.3) then
          CALL cv3_tracer(nloc, len, ncum, nd, nd, &
               ment, sij, da, phi)
       endif

       !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
       ! --- UNCOMPRESS THE FIELDS
       !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
       ! set iflag1  = 42 for non convective points
       do  i = 1, len
          iflag1(i) = 42
       end do

       if (iflag_con.eq.3) then
          CALL cv3_uncompress(nloc, len, ncum, nd, idcum, iflag, precip, &
               VPrecip, sig, w0, ft, fq, fu, fv, inb, Ma, upwd, dnwd, dnwd0, &
               qcondc, wd, cape, da, phi, mp, iflag1, precip1, VPrecip1, &
               sig1, w01, ft1, fq1, fu1, fv1, inb1, Ma1, upwd1, dnwd1, &
               dnwd01, qcondc1, wd1, cape1, da1, phi1, mp1)
       endif

       if (iflag_con.eq.4) then
          CALL cv_uncompress(nloc, len, ncum, nd, idcum &
               , iflag &
               , precip, cbmf &
               , ft, fq, fu, fv &
               , Ma, qcondc             &
               , iflag1 &
               , precip1, cbmf1 &
               , ft1, fq1, fu1, fv1 &
               , Ma1, qcondc1 )
       endif
    ENDIF ! ncum>0

  end SUBROUTINE cv_driver

end module cv_driver_m
1	guez	52	module cv_driver_m
2	guez	3
3	guez	52	implicit none
4	guez	3
5	guez	52	contains
6	guez	3
7	guez	97	SUBROUTINE cv_driver(len, nd, t1, q1, qs1, u1, v1, p1, ph1, iflag1, ft1, &
8			fq1, fu1, fv1, precip1, VPrecip1, cbmf1, sig1, w01, icb1, inb1, delt, &
9			Ma1, upwd1, dnwd1, dnwd01, qcondc1, wd1, cape1, da1, phi1, mp1)
10	guez	3
11	guez	91	! From LMDZ4/libf/phylmd/cv_driver.F, version 1.3, 2005/04/15 12:36:17
12	guez	69	! Main driver for convection
13	guez	97	! Author: S. Bony, March 2002
14	guez	69
15	guez	72	! Several modules corresponding to different physical processes
16
17			! Several versions of convect may be used:
18	guez	97	! - iflag_con = 3: version lmd
19			! - iflag_con = 4: version 4.3b
20	guez	72
21	guez	69	use clesphys2, only: iflag_con
22	guez	91	use cv3_compress_m, only: cv3_compress
23	guez	97	use cv3_mixing_m, only: cv3_mixing
24	guez	69	use cv3_param_m, only: cv3_param
25	guez	97	use cv3_prelim_m, only: cv3_prelim
26			use cv3_tracer_m, only: cv3_tracer
27			use cv3_uncompress_m, only: cv3_uncompress
28			use cv3_unsat_m, only: cv3_unsat
29			use cv3_yield_m, only: cv3_yield
30			use cv_uncompress_m, only: cv_uncompress
31	guez	62	USE dimphy, ONLY: klev, klon
32
33	guez	97	integer, intent(in):: len ! first dimension
34			integer, intent(in):: nd ! vertical dimension
35			real, intent(in):: t1(len, nd) ! temperature
36			real q1(len, nd) ! Input specific hum
37			real qs1(len, nd)
38	guez	52	! qs1 Real Input sat specific hum
39	guez	97	real, intent(in):: u1(len, nd)
40	guez	52	! u1 Real Input u-wind
41	guez	97	real, intent(in):: v1(len, nd)
42	guez	52	! v1 Real Input v-wind
43	guez	97	real p1(len, nd)
44	guez	52	! p1 Real Input full level pressure
45	guez	97	real ph1(len, nd + 1)
46	guez	52	! ph1 Real Input half level pressure
47	guez	97	integer iflag1(len)
48	guez	52	! iflag1 Integer Output flag for Emanuel conditions
49	guez	97	real ft1(len, nd)
50	guez	52	! ft1 Real Output temp tend
51	guez	97	real fq1(len, nd)
52	guez	52	! fq1 Real Output spec hum tend
53	guez	97	real fu1(len, nd)
54	guez	52	! fu1 Real Output u-wind tend
55	guez	97	real fv1(len, nd)
56	guez	52	! fv1 Real Output v-wind tend
57	guez	97	real precip1(len)
58	guez	52	! precip1 Real Output precipitation
59	guez	97	real VPrecip1(len, nd+1)
60	guez	52	! VPrecip1 Real Output vertical profile of precipitations
61	guez	97	real cbmf1(len)
62	guez	52	! cbmf1 Real Output cloud base mass flux
63	guez	72	real, intent(inout):: sig1(klon, klev) ! section adiabatic updraft
64
65			real, intent(inout):: w01(klon, klev)
66			! vertical velocity within adiabatic updraft
67
68			integer icb1(klon)
69			integer inb1(klon)
70			real, intent(in):: delt
71	guez	97	! delt Real Input time step
72	guez	52	real Ma1(len, nd)
73	guez	97	! Ma1 Real Output mass flux adiabatic updraft
74	guez	62	real, intent(out):: upwd1(len, nd) ! total upward mass flux (adiab+mixed)
75			real, intent(out):: dnwd1(len, nd) ! saturated downward mass flux (mixed)
76			real, intent(out):: dnwd01(len, nd) ! unsaturated downward mass flux
77	guez	3
78	guez	52	real qcondc1(len, nd) ! cld
79	guez	97	! qcondc1 Real Output in-cld mixing ratio of condensed water
80	guez	52	real wd1(len) ! gust
81	guez	97	! wd1 Real Output downdraft velocity scale for sfc fluxes
82	guez	52	real cape1(len)
83	guez	97	! cape1 Real Output CAPE
84	guez	3
85	guez	99	real, intent(inout):: da1(len, nd), phi1(len, nd, nd), mp1(len, nd)
86	guez	3
87	guez	52	!-------------------------------------------------------------------
88			! --- ARGUMENTS
89			!-------------------------------------------------------------------
90			! --- On input:
91	guez	62
92	guez	52	! t: Array of absolute temperature (K) of dimension ND, with first
93			! index corresponding to lowest model level. Note that this array
94			! will be altered by the subroutine if dry convective adjustment
95			! occurs and if IPBL is not equal to 0.
96	guez	62
97	guez	52	! q: Array of specific humidity (gm/gm) of dimension ND, with first
98			! index corresponding to lowest model level. Must be defined
99			! at same grid levels as T. Note that this array will be altered
100			! if dry convective adjustment occurs and if IPBL is not equal to 0.
101	guez	62
102	guez	52	! qs: Array of saturation specific humidity of dimension ND, with first
103			! index corresponding to lowest model level. Must be defined
104			! at same grid levels as T. Note that this array will be altered
105			! if dry convective adjustment occurs and if IPBL is not equal to 0.
106	guez	62
107	guez	52	! u: Array of zonal wind velocity (m/s) of dimension ND, witth first
108			! index corresponding with the lowest model level. Defined at
109			! same levels as T. Note that this array will be altered if
110			! dry convective adjustment occurs and if IPBL is not equal to 0.
111	guez	62
112	guez	52	! v: Same as u but for meridional velocity.
113	guez	62
114	guez	52	! p: Array of pressure (mb) of dimension ND, with first
115			! index corresponding to lowest model level. Must be defined
116			! at same grid levels as T.
117	guez	62
118	guez	52	! ph: Array of pressure (mb) of dimension ND+1, with first index
119			! corresponding to lowest level. These pressures are defined at
120			! levels intermediate between those of P, T, Q and QS. The first
121			! value of PH should be greater than (i.e. at a lower level than)
122			! the first value of the array P.
123	guez	62
124	guez	52	! nl: The maximum number of levels to which convection can penetrate, plus 1.
125			! NL MUST be less than or equal to ND-1.
126	guez	62
127	guez	52	! delt: The model time step (sec) between calls to CONVECT
128	guez	62
129	guez	52	!----------------------------------------------------------------------------
130			! --- On Output:
131	guez	62
132	guez	52	! iflag: An output integer whose value denotes the following:
133			! VALUE INTERPRETATION
134			! ----- --------------
135			! 0 Moist convection occurs.
136			! 1 Moist convection occurs, but a CFL condition
137			! on the subsidence warming is violated. This
138			! does not cause the scheme to terminate.
139			! 2 Moist convection, but no precip because ep(inb) lt 0.0001
140			! 3 No moist convection because new cbmf is 0 and old cbmf is 0.
141			! 4 No moist convection; atmosphere is not
142			! unstable
143			! 6 No moist convection because ihmin le minorig.
144			! 7 No moist convection because unreasonable
145			! parcel level temperature or specific humidity.
146			! 8 No moist convection: lifted condensation
147			! level is above the 200 mb level.
148			! 9 No moist convection: cloud base is higher
149			! then the level NL-1.
150	guez	62
151	guez	52	! ft: Array of temperature tendency (K/s) of dimension ND, defined at same
152			! grid levels as T, Q, QS and P.
153	guez	62
154	guez	52	! fq: Array of specific humidity tendencies ((gm/gm)/s) of dimension ND,
155			! defined at same grid levels as T, Q, QS and P.
156	guez	62
157	guez	52	! fu: Array of forcing of zonal velocity (m/s^2) of dimension ND,
158			! defined at same grid levels as T.
159	guez	62
160	guez	52	! fv: Same as FU, but for forcing of meridional velocity.
161	guez	62
162	guez	52	! precip: Scalar convective precipitation rate (mm/day).
163	guez	62
164	guez	52	! VPrecip: Vertical profile of convective precipitation (kg/m2/s).
165	guez	62
166	guez	52	! wd: A convective downdraft velocity scale. For use in surface
167			! flux parameterizations. See convect.ps file for details.
168	guez	62
169	guez	52	! tprime: A convective downdraft temperature perturbation scale (K).
170			! For use in surface flux parameterizations. See convect.ps
171			! file for details.
172	guez	62
173	guez	52	! qprime: A convective downdraft specific humidity
174			! perturbation scale (gm/gm).
175			! For use in surface flux parameterizations. See convect.ps
176			! file for details.
177	guez	62
178	guez	52	! cbmf: The cloud base mass flux ((kg/m**2)/s). THIS SCALAR VALUE MUST
179			! BE STORED BY THE CALLING PROGRAM AND RETURNED TO CONVECT AT
180			! ITS NEXT CALL. That is, the value of CBMF must be "remembered"
181			! by the calling program between calls to CONVECT.
182	guez	62
183	guez	52	! det: Array of detrainment mass flux of dimension ND.
184	guez	62
185	guez	52	!-------------------------------------------------------------------
186	guez	62
187	guez	52	! Local arrays
188	guez	3
189	guez	72	real da(len, nd), phi(len, nd, nd), mp(len, nd)
190
191	guez	97	integer i, k, il
192	guez	52	integer icbmax
193			integer nk1(klon)
194			integer icbs1(klon)
195	guez	3
196	guez	52	real plcl1(klon)
197			real tnk1(klon)
198			real qnk1(klon)
199			real gznk1(klon)
200			real pbase1(klon)
201			real buoybase1(klon)
202	guez	3
203	guez	52	real lv1(klon, klev)
204			real cpn1(klon, klev)
205			real tv1(klon, klev)
206			real gz1(klon, klev)
207			real hm1(klon, klev)
208			real h1(klon, klev)
209			real tp1(klon, klev)
210			real tvp1(klon, klev)
211			real clw1(klon, klev)
212			real th1(klon, klev)
213	guez	62
214	guez	52	integer ncum
215	guez	62
216	guez	52	! (local) compressed fields:
217	guez	62
218	guez	52	integer nloc
219	guez	91	parameter (nloc = klon) ! pour l'instant
220	guez	3
221	guez	52	integer idcum(nloc)
222			integer iflag(nloc), nk(nloc), icb(nloc)
223			integer nent(nloc, klev)
224			integer icbs(nloc)
225			integer inb(nloc), inbis(nloc)
226	guez	3
227	guez	52	real cbmf(nloc), plcl(nloc), tnk(nloc), qnk(nloc), gznk(nloc)
228			real t(nloc, klev), q(nloc, klev), qs(nloc, klev)
229			real u(nloc, klev), v(nloc, klev)
230			real gz(nloc, klev), h(nloc, klev), lv(nloc, klev), cpn(nloc, klev)
231			real p(nloc, klev), ph(nloc, klev+1), tv(nloc, klev), tp(nloc, klev)
232			real clw(nloc, klev)
233			real dph(nloc, klev)
234			real pbase(nloc), buoybase(nloc), th(nloc, klev)
235			real tvp(nloc, klev)
236			real sig(nloc, klev), w0(nloc, klev)
237			real hp(nloc, klev), ep(nloc, klev), sigp(nloc, klev)
238			real frac(nloc), buoy(nloc, klev)
239			real cape(nloc)
240			real m(nloc, klev), ment(nloc, klev, klev), qent(nloc, klev, klev)
241			real uent(nloc, klev, klev), vent(nloc, klev, klev)
242			real ments(nloc, klev, klev), qents(nloc, klev, klev)
243			real sij(nloc, klev, klev), elij(nloc, klev, klev)
244			real qp(nloc, klev), up(nloc, klev), vp(nloc, klev)
245			real wt(nloc, klev), water(nloc, klev), evap(nloc, klev)
246			real b(nloc, klev), ft(nloc, klev), fq(nloc, klev)
247			real fu(nloc, klev), fv(nloc, klev)
248			real upwd(nloc, klev), dnwd(nloc, klev), dnwd0(nloc, klev)
249			real Ma(nloc, klev), mike(nloc, klev), tls(nloc, klev)
250			real tps(nloc, klev), qprime(nloc), tprime(nloc)
251			real precip(nloc)
252			real VPrecip(nloc, klev+1)
253			real qcondc(nloc, klev) ! cld
254			real wd(nloc) ! gust
255	guez	3
256	guez	52	!-------------------------------------------------------------------
257			! --- SET CONSTANTS AND PARAMETERS
258			!-------------------------------------------------------------------
259	guez	3
260	guez	52	! -- set simulation flags:
261			! (common cvflag)
262	guez	3
263	guez	52	CALL cv_flag
264
265			! -- set thermodynamical constants:
266			! (common cvthermo)
267
268	guez	69	CALL cv_thermo
269	guez	52
270			! -- set convect parameters
271	guez	62
272	guez	52	! includes microphysical parameters and parameters that
273			! control the rate of approach to quasi-equilibrium)
274			! (common cvparam)
275
276			if (iflag_con.eq.3) then
277			CALL cv3_param(nd, delt)
278			endif
279
280			if (iflag_con.eq.4) then
281	guez	3	CALL cv_param(nd)
282	guez	52	endif
283	guez	3
284	guez	52	!---------------------------------------------------------------------
285			! --- INITIALIZE OUTPUT ARRAYS AND PARAMETERS
286			!---------------------------------------------------------------------
287	guez	3
288	guez	91	do k = 1, nd
289			do i = 1, len
290			ft1(i, k) = 0.0
291			fq1(i, k) = 0.0
292			fu1(i, k) = 0.0
293			fv1(i, k) = 0.0
294			tvp1(i, k) = 0.0
295			tp1(i, k) = 0.0
296			clw1(i, k) = 0.0
297	guez	52	!ym
298	guez	91	clw(i, k) = 0.0
299			gz1(i, k) = 0.
300	guez	52	VPrecip1(i, k) = 0.
301	guez	91	Ma1(i, k) = 0.0
302			upwd1(i, k) = 0.0
303			dnwd1(i, k) = 0.0
304			dnwd01(i, k) = 0.0
305			qcondc1(i, k) = 0.0
306	guez	52	end do
307			end do
308	guez	3
309	guez	91	do i = 1, len
310			precip1(i) = 0.0
311			iflag1(i) = 0
312			wd1(i) = 0.0
313			cape1(i) = 0.0
314			VPrecip1(i, nd+1) = 0.0
315	guez	52	end do
316	guez	3
317	guez	52	if (iflag_con.eq.3) then
318	guez	91	do il = 1, len
319			sig1(il, nd) = sig1(il, nd) + 1.
320			sig1(il, nd) = min(sig1(il, nd), 12.1)
321	guez	52	enddo
322			endif
323	guez	3
324	guez	52	!--------------------------------------------------------------------
325			! --- CALCULATE ARRAYS OF GEOPOTENTIAL, HEAT CAPACITY & STATIC ENERGY
326			!--------------------------------------------------------------------
327	guez	3
328	guez	52	if (iflag_con.eq.3) then
329	guez	97	CALL cv3_prelim(len, nd, nd + 1, t1, q1, p1, ph1, lv1, cpn1, tv1, gz1, &
330			h1, hm1, th1)
331	guez	52	endif
332	guez	3
333	guez	52	if (iflag_con.eq.4) then
334	guez	97	CALL cv_prelim(len, nd, nd + 1, t1, q1, p1, ph1 &
335	guez	52	, lv1, cpn1, tv1, gz1, h1, hm1)
336			endif
337	guez	3
338	guez	52	!--------------------------------------------------------------------
339			! --- CONVECTIVE FEED
340			!--------------------------------------------------------------------
341	guez	3
342	guez	52	if (iflag_con.eq.3) then
343			CALL cv3_feed(len, nd, t1, q1, qs1, p1, ph1, hm1, gz1 &
344			, nk1, icb1, icbmax, iflag1, tnk1, qnk1, gznk1, plcl1) ! nd->na
345			endif
346	guez	3
347	guez	52	if (iflag_con.eq.4) then
348			CALL cv_feed(len, nd, t1, q1, qs1, p1, hm1, gz1 &
349			, nk1, icb1, icbmax, iflag1, tnk1, qnk1, gznk1, plcl1)
350			endif
351	guez	3
352	guez	52	!--------------------------------------------------------------------
353			! --- UNDILUTE (ADIABATIC) UPDRAFT / 1st part
354			! (up through ICB for convect4, up through ICB+1 for convect3)
355			! Calculates the lifted parcel virtual temperature at nk, the
356			! actual temperature, and the adiabatic liquid water content.
357			!--------------------------------------------------------------------
358	guez	3
359	guez	52	if (iflag_con.eq.3) then
360			CALL cv3_undilute1(len, nd, t1, q1, qs1, gz1, plcl1, p1, nk1, icb1 &
361			, tp1, tvp1, clw1, icbs1) ! nd->na
362			endif
363	guez	3
364	guez	52	if (iflag_con.eq.4) then
365			CALL cv_undilute1(len, nd, t1, q1, qs1, gz1, p1, nk1, icb1, icbmax &
366			, tp1, tvp1, clw1)
367			endif
368	guez	3
369	guez	52	!-------------------------------------------------------------------
370			! --- TRIGGERING
371			!-------------------------------------------------------------------
372	guez	3
373	guez	52	if (iflag_con.eq.3) then
374	guez	72	CALL cv3_trigger(len, nd, icb1, plcl1, p1, th1, tv1, tvp1, pbase1, &
375			buoybase1, iflag1, sig1, w01) ! nd->na
376	guez	52	endif
377	guez	3
378	guez	52	if (iflag_con.eq.4) then
379			CALL cv_trigger(len, nd, icb1, cbmf1, tv1, tvp1, iflag1)
380			endif
381	guez	3
382	guez	52	! --- IF THIS POINT IS REACHED, MOIST CONVECTIVE ADJUSTMENT IS NECESSARY
383	guez	3
384	guez	91	ncum = 0
385			do i = 1, len
386	guez	52	if(iflag1(i).eq.0)then
387	guez	91	ncum = ncum+1
388			idcum(ncum) = i
389	guez	52	endif
390			end do
391	guez	3
392	guez	52	IF (ncum.gt.0) THEN
393	guez	3
394	guez	52	!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
395			! --- COMPRESS THE FIELDS
396			! (-> vectorization over convective gridpoints)
397			!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
398	guez	3
399	guez	52	if (iflag_con.eq.3) then
400	guez	97	CALL cv3_compress(len, nloc, ncum, nd, iflag1, nk1, icb1, icbs1, &
401			plcl1, tnk1, qnk1, gznk1, pbase1, buoybase1, t1, q1, qs1, u1, &
402			v1, gz1, th1, h1, lv1, cpn1, p1, ph1, tv1, tp1, tvp1, clw1, &
403			sig1, w01, iflag, nk, icb, icbs, plcl, tnk, qnk, gznk, pbase, &
404			buoybase, t, q, qs, u, v, gz, th, h, lv, cpn, p, ph, tv, tp, &
405			tvp, clw, sig, w0)
406	guez	52	endif
407	guez	3
408	guez	52	if (iflag_con.eq.4) then
409			CALL cv_compress( len, nloc, ncum, nd &
410			, iflag1, nk1, icb1 &
411			, cbmf1, plcl1, tnk1, qnk1, gznk1 &
412			, t1, q1, qs1, u1, v1, gz1 &
413			, h1, lv1, cpn1, p1, ph1, tv1, tp1, tvp1, clw1 &
414			, iflag, nk, icb &
415			, cbmf, plcl, tnk, qnk, gznk &
416			, t, q, qs, u, v, gz, h, lv, cpn, p, ph, tv, tp, tvp, clw &
417			, dph )
418			endif
419	guez	3
420	guez	52	!-------------------------------------------------------------------
421			! --- UNDILUTE (ADIABATIC) UPDRAFT / second part :
422			! --- FIND THE REST OF THE LIFTED PARCEL TEMPERATURES
423			! --- &
424			! --- COMPUTE THE PRECIPITATION EFFICIENCIES AND THE
425			! --- FRACTION OF PRECIPITATION FALLING OUTSIDE OF CLOUD
426			! --- &
427			! --- FIND THE LEVEL OF NEUTRAL BUOYANCY
428			!-------------------------------------------------------------------
429	guez	3
430	guez	52	if (iflag_con.eq.3) then
431			CALL cv3_undilute2(nloc, ncum, nd, icb, icbs, nk &
432			, tnk, qnk, gznk, t, q, qs, gz &
433			, p, h, tv, lv, pbase, buoybase, plcl &
434			, inb, tp, tvp, clw, hp, ep, sigp, buoy) !na->nd
435			endif
436	guez	3
437	guez	52	if (iflag_con.eq.4) then
438			CALL cv_undilute2(nloc, ncum, nd, icb, nk &
439			, tnk, qnk, gznk, t, q, qs, gz &
440			, p, dph, h, tv, lv &
441			, inb, inbis, tp, tvp, clw, hp, ep, sigp, frac)
442			endif
443	guez	3
444	guez	52	!-------------------------------------------------------------------
445			! --- CLOSURE
446			!-------------------------------------------------------------------
447	guez	3
448	guez	52	if (iflag_con.eq.3) then
449			CALL cv3_closure(nloc, ncum, nd, icb, inb &
450			, pbase, p, ph, tv, buoy &
451			, sig, w0, cape, m) ! na->nd
452			endif
453	guez	3
454	guez	52	if (iflag_con.eq.4) then
455			CALL cv_closure(nloc, ncum, nd, nk, icb &
456			, tv, tvp, p, ph, dph, plcl, cpn &
457			, iflag, cbmf)
458			endif
459	guez	3
460	guez	52	!-------------------------------------------------------------------
461			! --- MIXING
462			!-------------------------------------------------------------------
463	guez	3
464	guez	52	if (iflag_con.eq.3) then
465	guez	97	CALL cv3_mixing(nloc, ncum, nd, nd, icb, nk, inb, ph, t, q, &
466			qs, u, v, h, lv, qnk, hp, tv, tvp, ep, clw, m, sig, ment, &
467			qent, uent, vent, nent, sij, elij, ments, qents)
468	guez	52	endif
469	guez	3
470	guez	52	if (iflag_con.eq.4) then
471			CALL cv_mixing(nloc, ncum, nd, icb, nk, inb, inbis &
472			, ph, t, q, qs, u, v, h, lv, qnk &
473			, hp, tv, tvp, ep, clw, cbmf &
474			, m, ment, qent, uent, vent, nent, sij, elij)
475			endif
476	guez	3
477	guez	52	!-------------------------------------------------------------------
478			! --- UNSATURATED (PRECIPITATING) DOWNDRAFTS
479			!-------------------------------------------------------------------
480	guez	3
481	guez	52	if (iflag_con.eq.3) then
482	guez	97	CALL cv3_unsat(nloc, ncum, nd, nd, icb, inb &
483			, t, q, qs, gz, u, v, p, ph &
484	guez	52	, th, tv, lv, cpn, ep, sigp, clw &
485			, m, ment, elij, delt, plcl &
486	guez	97	, mp, qp, up, vp, wt, water, evap, b)! na->nd
487	guez	52	endif
488	guez	3
489	guez	52	if (iflag_con.eq.4) then
490			CALL cv_unsat(nloc, ncum, nd, inb, t, q, qs, gz, u, v, p, ph &
491			, h, lv, ep, sigp, clw, m, ment, elij &
492			, iflag, mp, qp, up, vp, wt, water, evap)
493			endif
494	guez	3
495	guez	52	!-------------------------------------------------------------------
496			! --- YIELD
497			! (tendencies, precipitation, variables of interface with other
498			! processes, etc)
499			!-------------------------------------------------------------------
500	guez	3
501	guez	52	if (iflag_con.eq.3) then
502	guez	97	CALL cv3_yield(nloc, ncum, nd, nd &
503	guez	52	, icb, inb, delt &
504	guez	97	, t, q, u, v, gz, p, ph, h, hp, lv, cpn, th &
505			, ep, clw, m, tp, mp, qp, up, vp &
506	guez	52	, wt, water, evap, b &
507	guez	97	, ment, qent, uent, vent, nent, elij, sig &
508	guez	52	, tv, tvp &
509	guez	97	, iflag, precip, VPrecip, ft, fq, fu, fv &
510	guez	52	, upwd, dnwd, dnwd0, ma, mike, tls, tps, qcondc, wd)! na->nd
511			endif
512	guez	3
513	guez	52	if (iflag_con.eq.4) then
514			CALL cv_yield(nloc, ncum, nd, nk, icb, inb, delt &
515			, t, q, u, v, gz, p, ph, h, hp, lv, cpn &
516			, ep, clw, frac, m, mp, qp, up, vp &
517			, wt, water, evap &
518			, ment, qent, uent, vent, nent, elij &
519			, tv, tvp &
520			, iflag, wd, qprime, tprime &
521			, precip, cbmf, ft, fq, fu, fv, Ma, qcondc)
522			endif
523	guez	3
524	guez	52	!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
525			! --- passive tracers
526			!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
527	guez	3
528	guez	52	if (iflag_con.eq.3) then
529			CALL cv3_tracer(nloc, len, ncum, nd, nd, &
530			ment, sij, da, phi)
531			endif
532	guez	3
533	guez	52	!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
534			! --- UNCOMPRESS THE FIELDS
535			!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
536	guez	91	! set iflag1 = 42 for non convective points
537			do i = 1, len
538			iflag1(i) = 42
539	guez	52	end do
540	guez	62
541	guez	52	if (iflag_con.eq.3) then
542	guez	99	CALL cv3_uncompress(nloc, len, ncum, nd, idcum, iflag, precip, &
543			VPrecip, sig, w0, ft, fq, fu, fv, inb, Ma, upwd, dnwd, dnwd0, &
544			qcondc, wd, cape, da, phi, mp, iflag1, precip1, VPrecip1, &
545			sig1, w01, ft1, fq1, fu1, fv1, inb1, Ma1, upwd1, dnwd1, &
546			dnwd01, qcondc1, wd1, cape1, da1, phi1, mp1)
547	guez	52	endif
548	guez	3
549	guez	52	if (iflag_con.eq.4) then
550			CALL cv_uncompress(nloc, len, ncum, nd, idcum &
551			, iflag &
552			, precip, cbmf &
553			, ft, fq, fu, fv &
554			, Ma, qcondc &
555			, iflag1 &
556			, precip1, cbmf1 &
557			, ft1, fq1, fu1, fv1 &
558			, Ma1, qcondc1 )
559			endif
560			ENDIF ! ncum>0
561	guez	3
562	guez	52	end SUBROUTINE cv_driver
563	guez	3
564	guez	52	end module cv_driver_m