Sources/phylmd/cv_driver.f

module cv_driver_m

  implicit none

contains

  SUBROUTINE cv_driver(t1, q1, qs1, u1, v1, p1, ph1, iflag1, ft1, fq1, fu1, &
       fv1, precip1, VPrecip1, 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

    use cv30_closure_m, only: cv30_closure
    use cv30_compress_m, only: cv30_compress
    use cv30_feed_m, only: cv30_feed
    use cv30_mixing_m, only: cv30_mixing
    use cv30_param_m, only: cv30_param, nl
    use cv30_prelim_m, only: cv30_prelim
    use cv30_tracer_m, only: cv30_tracer
    use cv30_uncompress_m, only: cv30_uncompress
    use cv30_undilute2_m, only: cv30_undilute2
    use cv30_unsat_m, only: cv30_unsat
    use cv30_yield_m, only: cv30_yield
    USE dimphy, ONLY: klev, klon

    real, intent(in):: t1(klon, klev)
    ! temperature (K), with first index corresponding to lowest model
    ! level

    real, intent(in):: q1(klon, klev)
    ! Specific humidity, with first index corresponding to lowest
    ! model level. Must be defined at same grid levels as T1.

    real, intent(in):: qs1(klon, klev)
    ! Saturation specific humidity, with first index corresponding to
    ! lowest model level. Must be defined at same grid levels as
    ! T1.

    real, intent(in):: u1(klon, klev), v1(klon, klev)
    ! Zonal wind and meridional velocity (m/s), witth first index
    ! corresponding with the lowest model level. Defined at same
    ! levels as T1.

    real, intent(in):: p1(klon, klev)
    ! Full level pressure (mb) of dimension KLEV, with first index
    ! corresponding to lowest model level. Must be defined at same
    ! grid levels as T1.

    real, intent(in):: ph1(klon, klev + 1) 
    ! Half level pressure (mb), with first index corresponding to
    ! lowest level. These pressures are defined at levels intermediate
    ! between those of P1, T1, Q1 and QS1. The first value of PH
    ! should be greater than (i.e. at a lower level than) the first
    ! value of the array P1.

    integer, intent(out):: iflag1(klon)
    ! Flag for Emanuel conditions.

    ! 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 precipitation because ep(inb) < 1e-4

    ! 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.

    real, intent(out):: ft1(klon, klev)
    ! Temperature tendency (K/s), defined at same grid levels as T1,
    ! Q1, QS1 and P1.

    real, intent(out):: fq1(klon, klev)
    ! Specific humidity tendencies (s-1), defined at same grid levels
    ! as T1, Q1, QS1 and P1.

    real, intent(out):: fu1(klon, klev), fv1(klon, klev)
    ! Forcing (tendency) of zonal and meridional velocity (m/s^2),
    ! defined at same grid levels as T1.

    real, intent(out):: precip1(klon) ! convective precipitation rate (mm/day)

    real, intent(out):: VPrecip1(klon, klev + 1)
    ! vertical profile of convective precipitation (kg/m2/s)

    real, intent(inout):: sig1(klon, klev) ! section adiabatic updraft

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

    integer, intent(out):: icb1(klon)
    integer, intent(inout):: inb1(klon)
    real, intent(in):: delt ! the model time step (sec) between calls

    real Ma1(klon, klev) ! Output mass flux adiabatic updraft

    real, intent(out):: upwd1(klon, klev) 
    ! total upward mass flux (adiab + mixed)

    real, intent(out):: dnwd1(klon, klev) ! saturated downward mass flux (mixed)
    real, intent(out):: dnwd01(klon, klev) ! unsaturated downward mass flux

    real qcondc1(klon, klev) ! Output in-cld mixing ratio of condensed water

    real wd1(klon) ! gust
    ! Output downdraft velocity scale for surface fluxes
    ! A convective downdraft velocity scale. For use in surface
    ! flux parameterizations. See convect.ps file for details.

    real cape1(klon) ! Output
    real, intent(inout):: da1(klon, klev), phi1(klon, klev, klev)
    real, intent(inout):: mp1(klon, klev)

    ! Local:

    real da(klon, klev), phi(klon, klev, klev), mp(klon, klev)

    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

    ! Compressed fields:

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

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

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

    ! SET CONSTANTS AND PARAMETERS

    ! 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)
    CALL cv30_param(delt)

    ! INITIALIZE OUTPUT ARRAYS AND PARAMETERS

    do k = 1, klev
       do i = 1, klon
          ft1(i, k) = 0.
          fq1(i, k) = 0.
          fu1(i, k) = 0.
          fv1(i, k) = 0.
          tvp1(i, k) = 0.
          tp1(i, k) = 0.
          clw1(i, k) = 0.
          clw(i, k) = 0.
          gz1(i, k) = 0.
          VPrecip1(i, k) = 0.
          Ma1(i, k) = 0.
          upwd1(i, k) = 0.
          dnwd1(i, k) = 0.
          dnwd01(i, k) = 0.
          qcondc1(i, k) = 0.
       end do
    end do

    do i = 1, klon
       precip1(i) = 0.
       iflag1(i) = 0
       wd1(i) = 0.
       cape1(i) = 0.
       VPrecip1(i, klev + 1) = 0.
    end do

    do il = 1, klon
       sig1(il, klev) = sig1(il, klev) + 1.
       sig1(il, klev) = min(sig1(il, klev), 12.1)
    enddo

    ! CALCULATE ARRAYS OF GEOPOTENTIAL, HEAT CAPACITY & STATIC ENERGY
    CALL cv30_prelim(klon, klev, klev + 1, t1, q1, p1, ph1, lv1, cpn1, tv1, &
         gz1, h1, hm1, th1)

    ! CONVECTIVE FEED
    CALL cv30_feed(klon, klev, t1, q1, qs1, p1, ph1, gz1, nk1, icb1, &
         icbmax, iflag1, tnk1, qnk1, gznk1, plcl1) ! klev->na

    ! 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.
    CALL cv30_undilute1(klon, klev, t1, q1, qs1, gz1, plcl1, p1, nk1, icb1, &
         tp1, tvp1, clw1, icbs1) ! klev->na

    ! TRIGGERING
    CALL cv30_trigger(klon, klev, icb1, plcl1, p1, th1, tv1, tvp1, pbase1, &
         buoybase1, iflag1, sig1, w01) ! klev->na

    ! Moist convective adjustment is necessary

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

    IF (ncum > 0) THEN
       allocate(b(ncum, nl))

       ! COMPRESS THE FIELDS
       ! (-> vectorization over convective gridpoints)
       CALL cv30_compress(klon, klon, ncum, klev, 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)

       CALL cv30_undilute2(ncum, icb, icbs, nk, tnk, qnk, gznk, t, qs, gz, p, &
            h, tv, lv, pbase, buoybase, plcl, inb(:ncum), tp, tvp, clw, hp, &
            ep, sigp, buoy)

       ! CLOSURE
       CALL cv30_closure(klon, ncum, klev, icb, inb, pbase, p, ph, tv, &
            buoy, sig, w0, cape, m) ! na->klev

       ! MIXING
       CALL cv30_mixing(klon, ncum, klev, klev, icb, nk, inb, t, q, qs, u, &
            v, h, lv, hp, ep, clw, m, sig, ment, qent, uent, vent, nent, &
            sij, elij, ments, qents)

       ! Unsaturated (precipitating) downdrafts
       CALL cv30_unsat(icb(:ncum), inb(:ncum), 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)

       ! Yield (tendencies, precipitation, variables of interface with
       ! other processes, etc)
       CALL cv30_yield(icb(:ncum), inb(:ncum), delt, t, q, u, v, gz, p, ph, &
            h, hp, lv, cpn, th, ep, clw, m, tp, mp, qp, up, vp, wt, &
            water(:ncum, :nl), evap(:ncum, :nl), 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)

       ! passive tracers
       CALL cv30_tracer(klon, ncum, klev, ment, sij, da, phi)

       ! UNCOMPRESS THE FIELDS

       ! set iflag1 = 42 for non convective points
       iflag1 = 42

       CALL cv30_uncompress(idcum(:ncum), 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

  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	183	SUBROUTINE cv_driver(t1, q1, qs1, u1, v1, p1, ph1, iflag1, ft1, fq1, fu1, &
8			fv1, precip1, VPrecip1, sig1, w01, icb1, inb1, delt, Ma1, upwd1, &
9			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	guez	187	use cv30_closure_m, only: cv30_closure
18	guez	185	use cv30_compress_m, only: cv30_compress
19			use cv30_feed_m, only: cv30_feed
20			use cv30_mixing_m, only: cv30_mixing
21	guez	188	use cv30_param_m, only: cv30_param, nl
22	guez	185	use cv30_prelim_m, only: cv30_prelim
23			use cv30_tracer_m, only: cv30_tracer
24			use cv30_uncompress_m, only: cv30_uncompress
25			use cv30_undilute2_m, only: cv30_undilute2
26			use cv30_unsat_m, only: cv30_unsat
27			use cv30_yield_m, only: cv30_yield
28	guez	62	USE dimphy, ONLY: klev, klon
29
30	guez	187	real, intent(in):: t1(klon, klev)
31			! temperature (K), with first index corresponding to lowest model
32			! level
33	guez	103
34	guez	187	real, intent(in):: q1(klon, klev)
35			! Specific humidity, with first index corresponding to lowest
36			! model level. Must be defined at same grid levels as T1.
37	guez	103
38	guez	187	real, intent(in):: qs1(klon, klev)
39			! Saturation specific humidity, with first index corresponding to
40			! lowest model level. Must be defined at same grid levels as
41			! T1.
42	guez	72
43	guez	187	real, intent(in):: u1(klon, klev), v1(klon, klev)
44			! Zonal wind and meridional velocity (m/s), witth first index
45			! corresponding with the lowest model level. Defined at same
46			! levels as T1.
47	guez	72
48	guez	187	real, intent(in):: p1(klon, klev)
49			! Full level pressure (mb) of dimension KLEV, with first index
50			! corresponding to lowest model level. Must be defined at same
51			! grid levels as T1.
52	guez	180
53	guez	187	real, intent(in):: ph1(klon, klev + 1)
54			! Half level pressure (mb), with first index corresponding to
55			! lowest level. These pressures are defined at levels intermediate
56			! between those of P1, T1, Q1 and QS1. The first value of PH
57			! should be greater than (i.e. at a lower level than) the first
58			! value of the array P1.
59	guez	180
60	guez	187	integer, intent(out):: iflag1(klon)
61			! Flag for Emanuel conditions.
62	guez	3
63	guez	187	! 0: Moist convection occurs.
64	guez	3
65	guez	187	! 1: Moist convection occurs, but a CFL condition on the
66			! subsidence warming is violated. This does not cause the scheme
67			! to terminate.
68	guez	3
69	guez	187	! 2: Moist convection, but no precipitation because ep(inb) < 1e-4
70	guez	103
71	guez	187	! 3: No moist convection because new cbmf is 0 and old cbmf is 0.
72	guez	62
73	guez	187	! 4: No moist convection; atmosphere is not unstable
74	guez	62
75	guez	187	! 6: No moist convection because ihmin le minorig.
76	guez	62
77	guez	187	! 7: No moist convection because unreasonable parcel level
78			! temperature or specific humidity.
79	guez	62
80	guez	187	! 8: No moist convection: lifted condensation level is above the
81			! 200 mb level.
82	guez	62
83	guez	187	! 9: No moist convection: cloud base is higher then the level NL-1.
84	guez	62
85	guez	187	real, intent(out):: ft1(klon, klev)
86			! Temperature tendency (K/s), defined at same grid levels as T1,
87			! Q1, QS1 and P1.
88	guez	62
89	guez	187	real, intent(out):: fq1(klon, klev)
90			! Specific humidity tendencies (s-1), defined at same grid levels
91			! as T1, Q1, QS1 and P1.
92	guez	62
93	guez	187	real, intent(out):: fu1(klon, klev), fv1(klon, klev)
94			! Forcing (tendency) of zonal and meridional velocity (m/s^2),
95			! defined at same grid levels as T1.
96	guez	62
97	guez	187	real, intent(out):: precip1(klon) ! convective precipitation rate (mm/day)
98	guez	62
99	guez	187	real, intent(out):: VPrecip1(klon, klev + 1)
100			! vertical profile of convective precipitation (kg/m2/s)
101	guez	62
102	guez	187	real, intent(inout):: sig1(klon, klev) ! section adiabatic updraft
103	guez	62
104	guez	187	real, intent(inout):: w01(klon, klev)
105			! vertical velocity within adiabatic updraft
106	guez	62
107	guez	187	integer, intent(out):: icb1(klon)
108			integer, intent(inout):: inb1(klon)
109			real, intent(in):: delt ! the model time step (sec) between calls
110	guez	62
111	guez	187	real Ma1(klon, klev) ! Output mass flux adiabatic updraft
112	guez	62
113	guez	187	real, intent(out):: upwd1(klon, klev)
114			! total upward mass flux (adiab + mixed)
115	guez	62
116	guez	187	real, intent(out):: dnwd1(klon, klev) ! saturated downward mass flux (mixed)
117			real, intent(out):: dnwd01(klon, klev) ! unsaturated downward mass flux
118	guez	62
119	guez	187	real qcondc1(klon, klev) ! Output in-cld mixing ratio of condensed water
120	guez	62
121	guez	187	real wd1(klon) ! gust
122			! Output downdraft velocity scale for surface fluxes
123			! A convective downdraft velocity scale. For use in surface
124	guez	103	! flux parameterizations. See convect.ps file for details.
125	guez	62
126	guez	187	real cape1(klon) ! Output
127			real, intent(inout):: da1(klon, klev), phi1(klon, klev, klev)
128			real, intent(inout):: mp1(klon, klev)
129	guez	62
130	guez	187	! Local:
131	guez	62
132	guez	103	real da(klon, klev), phi(klon, klev, klev), mp(klon, klev)
133	guez	3
134	guez	97	integer i, k, il
135	guez	52	integer icbmax
136			integer nk1(klon)
137			integer icbs1(klon)
138	guez	3
139	guez	52	real plcl1(klon)
140			real tnk1(klon)
141			real qnk1(klon)
142			real gznk1(klon)
143			real pbase1(klon)
144			real buoybase1(klon)
145	guez	3
146	guez	52	real lv1(klon, klev)
147			real cpn1(klon, klev)
148			real tv1(klon, klev)
149			real gz1(klon, klev)
150			real hm1(klon, klev)
151			real h1(klon, klev)
152			real tp1(klon, klev)
153			real tvp1(klon, klev)
154			real clw1(klon, klev)
155			real th1(klon, klev)
156	guez	62
157	guez	52	integer ncum
158	guez	62
159	guez	187	! Compressed fields:
160	guez	62
161	guez	103	integer idcum(klon)
162			integer iflag(klon), nk(klon), icb(klon)
163			integer nent(klon, klev)
164			integer icbs(klon)
165	guez	183	integer inb(klon)
166	guez	3
167	guez	182	real plcl(klon), tnk(klon), qnk(klon), gznk(klon)
168	guez	103	real t(klon, klev), q(klon, klev), qs(klon, klev)
169			real u(klon, klev), v(klon, klev)
170			real gz(klon, klev), h(klon, klev), lv(klon, klev), cpn(klon, klev)
171	guez	180	real p(klon, klev), ph(klon, klev + 1), tv(klon, klev), tp(klon, klev)
172	guez	103	real clw(klon, klev)
173			real pbase(klon), buoybase(klon), th(klon, klev)
174			real tvp(klon, klev)
175			real sig(klon, klev), w0(klon, klev)
176			real hp(klon, klev), ep(klon, klev), sigp(klon, klev)
177	guez	183	real buoy(klon, klev)
178	guez	103	real cape(klon)
179			real m(klon, klev), ment(klon, klev, klev), qent(klon, klev, klev)
180			real uent(klon, klev, klev), vent(klon, klev, klev)
181			real ments(klon, klev, klev), qents(klon, klev, klev)
182			real sij(klon, klev, klev), elij(klon, klev, klev)
183			real qp(klon, klev), up(klon, klev), vp(klon, klev)
184			real wt(klon, klev), water(klon, klev), evap(klon, klev)
185	guez	188	real, allocatable:: b(:, :) ! (ncum, nl)
186			real ft(klon, klev), fq(klon, klev)
187	guez	103	real fu(klon, klev), fv(klon, klev)
188			real upwd(klon, klev), dnwd(klon, klev), dnwd0(klon, klev)
189			real Ma(klon, klev), mike(klon, klev), tls(klon, klev)
190	guez	183	real tps(klon, klev)
191	guez	103	real precip(klon)
192	guez	180	real VPrecip(klon, klev + 1)
193	guez	103	real qcondc(klon, klev) ! cld
194			real wd(klon) ! gust
195	guez	3
196	guez	52	!-------------------------------------------------------------------
197	guez	3
198	guez	180	! SET CONSTANTS AND PARAMETERS
199
200			! set thermodynamical constants:
201	guez	103	! (common cvthermo)
202	guez	69	CALL cv_thermo
203	guez	52
204	guez	180	! set convect parameters
205	guez	103	! includes microphysical parameters and parameters that
206			! control the rate of approach to quasi-equilibrium)
207			! (common cvparam)
208	guez	186	CALL cv30_param(delt)
209	guez	52
210	guez	180	! INITIALIZE OUTPUT ARRAYS AND PARAMETERS
211	guez	3
212	guez	103	do k = 1, klev
213			do i = 1, klon
214	guez	187	ft1(i, k) = 0.
215			fq1(i, k) = 0.
216			fu1(i, k) = 0.
217			fv1(i, k) = 0.
218			tvp1(i, k) = 0.
219			tp1(i, k) = 0.
220			clw1(i, k) = 0.
221			clw(i, k) = 0.
222	guez	103	gz1(i, k) = 0.
223	guez	52	VPrecip1(i, k) = 0.
224	guez	187	Ma1(i, k) = 0.
225			upwd1(i, k) = 0.
226			dnwd1(i, k) = 0.
227			dnwd01(i, k) = 0.
228			qcondc1(i, k) = 0.
229	guez	52	end do
230			end do
231	guez	3
232	guez	103	do i = 1, klon
233	guez	187	precip1(i) = 0.
234	guez	91	iflag1(i) = 0
235	guez	187	wd1(i) = 0.
236			cape1(i) = 0.
237			VPrecip1(i, klev + 1) = 0.
238	guez	52	end do
239	guez	3
240	guez	181	do il = 1, klon
241			sig1(il, klev) = sig1(il, klev) + 1.
242			sig1(il, klev) = min(sig1(il, klev), 12.1)
243			enddo
244	guez	3
245	guez	180	! CALCULATE ARRAYS OF GEOPOTENTIAL, HEAT CAPACITY & STATIC ENERGY
246	guez	185	CALL cv30_prelim(klon, klev, klev + 1, t1, q1, p1, ph1, lv1, cpn1, tv1, &
247	guez	181	gz1, h1, hm1, th1)
248	guez	3
249	guez	180	! CONVECTIVE FEED
250	guez	185	CALL cv30_feed(klon, klev, t1, q1, qs1, p1, ph1, gz1, nk1, icb1, &
251	guez	181	icbmax, iflag1, tnk1, qnk1, gznk1, plcl1) ! klev->na
252	guez	3
253	guez	180	! UNDILUTE (ADIABATIC) UPDRAFT / 1st part
254			! (up through ICB for convect4, up through ICB + 1 for convect3)
255	guez	103	! Calculates the lifted parcel virtual temperature at nk, the
256			! actual temperature, and the adiabatic liquid water content.
257	guez	185	CALL cv30_undilute1(klon, klev, t1, q1, qs1, gz1, plcl1, p1, nk1, icb1, &
258	guez	181	tp1, tvp1, clw1, icbs1) ! klev->na
259	guez	3
260	guez	180	! TRIGGERING
261	guez	185	CALL cv30_trigger(klon, klev, icb1, plcl1, p1, th1, tv1, tvp1, pbase1, &
262	guez	181	buoybase1, iflag1, sig1, w01) ! klev->na
263	guez	3
264	guez	180	! Moist convective adjustment is necessary
265	guez	3
266	guez	91	ncum = 0
267	guez	103	do i = 1, klon
268	guez	180	if (iflag1(i) == 0) then
269			ncum = ncum + 1
270	guez	91	idcum(ncum) = i
271	guez	52	endif
272			end do
273	guez	3
274	guez	103	IF (ncum > 0) THEN
275	guez	188	allocate(b(ncum, nl))
276
277	guez	180	! COMPRESS THE FIELDS
278	guez	103	! (-> vectorization over convective gridpoints)
279	guez	185	CALL cv30_compress(klon, klon, ncum, klev, iflag1, nk1, icb1, icbs1, &
280	guez	181	plcl1, tnk1, qnk1, gznk1, pbase1, buoybase1, t1, q1, qs1, u1, &
281			v1, gz1, th1, h1, lv1, cpn1, p1, ph1, tv1, tp1, tvp1, clw1, &
282			sig1, w01, iflag, nk, icb, icbs, plcl, tnk, qnk, gznk, pbase, &
283			buoybase, t, q, qs, u, v, gz, th, h, lv, cpn, p, ph, tv, tp, &
284			tvp, clw, sig, w0)
285	guez	3
286	guez	187	CALL cv30_undilute2(ncum, icb, icbs, nk, tnk, qnk, gznk, t, qs, gz, p, &
287			h, tv, lv, pbase, buoybase, plcl, inb(:ncum), tp, tvp, clw, hp, &
288			ep, sigp, buoy)
289	guez	3
290	guez	180	! CLOSURE
291	guez	185	CALL cv30_closure(klon, ncum, klev, icb, inb, pbase, p, ph, tv, &
292	guez	181	buoy, sig, w0, cape, m) ! na->klev
293	guez	3
294	guez	180	! MIXING
295	guez	185	CALL cv30_mixing(klon, ncum, klev, klev, icb, nk, inb, t, q, qs, u, &
296	guez	181	v, h, lv, hp, ep, clw, m, sig, ment, qent, uent, vent, nent, &
297			sij, elij, ments, qents)
298	guez	3
299	guez	186	! Unsaturated (precipitating) downdrafts
300	guez	188	CALL cv30_unsat(icb(:ncum), inb(:ncum), t, q, qs, gz, u, v, p, ph, th, &
301			tv, lv, cpn, ep, sigp, clw, m, ment, elij, delt, plcl, mp, qp, &
302			up, vp, wt, water, evap, b)
303	guez	3
304	guez	186	! Yield (tendencies, precipitation, variables of interface with
305			! other processes, etc)
306	guez	188	CALL cv30_yield(icb(:ncum), inb(:ncum), delt, t, q, u, v, gz, p, ph, &
307			h, hp, lv, cpn, th, ep, clw, m, tp, mp, qp, up, vp, wt, &
308			water(:ncum, :nl), evap(:ncum, :nl), b, ment, qent, uent, vent, &
309			nent, elij, sig, tv, tvp, iflag, precip, VPrecip, ft, fq, fu, fv, &
310			upwd, dnwd, dnwd0, ma, mike, tls, tps, qcondc, wd)
311	guez	3
312	guez	180	! passive tracers
313	guez	185	CALL cv30_tracer(klon, ncum, klev, ment, sij, da, phi)
314	guez	3
315	guez	180	! UNCOMPRESS THE FIELDS
316	guez	103
317			! set iflag1 = 42 for non convective points
318	guez	186	iflag1 = 42
319	guez	62
320	guez	185	CALL cv30_uncompress(idcum(:ncum), iflag, precip, VPrecip, sig, w0, &
321	guez	181	ft, fq, fu, fv, inb, Ma, upwd, dnwd, dnwd0, qcondc, wd, cape, &
322			da, phi, mp, iflag1, precip1, VPrecip1, sig1, w01, ft1, fq1, &
323			fu1, fv1, inb1, Ma1, upwd1, dnwd1, dnwd01, qcondc1, wd1, &
324			cape1, da1, phi1, mp1)
325	guez	183	ENDIF
326	guez	3
327	guez	52	end SUBROUTINE cv_driver
328	guez	3
329	guez	52	end module cv_driver_m