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
    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 b(klon, klev), 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
       ! 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(ncum, 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(:ncum, :))

       ! Yield (tendencies, precipitation, variables of interface with
       ! other processes, etc)
       CALL cv30_yield(klon, ncum, klev, klev, 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->klev

       ! 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			use cv30_param_m, only: cv30_param
22			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			real b(klon, klev), ft(klon, klev), fq(klon, klev)
186			real fu(klon, klev), fv(klon, klev)
187			real upwd(klon, klev), dnwd(klon, klev), dnwd0(klon, klev)
188			real Ma(klon, klev), mike(klon, klev), tls(klon, klev)
189	guez	183	real tps(klon, klev)
190	guez	103	real precip(klon)
191	guez	180	real VPrecip(klon, klev + 1)
192	guez	103	real qcondc(klon, klev) ! cld
193			real wd(klon) ! gust
194	guez	3
195	guez	52	!-------------------------------------------------------------------
196	guez	3
197	guez	180	! SET CONSTANTS AND PARAMETERS
198
199			! set thermodynamical constants:
200	guez	103	! (common cvthermo)
201	guez	69	CALL cv_thermo
202	guez	52
203	guez	180	! set convect parameters
204	guez	103	! includes microphysical parameters and parameters that
205			! control the rate of approach to quasi-equilibrium)
206			! (common cvparam)
207	guez	186	CALL cv30_param(delt)
208	guez	52
209	guez	180	! INITIALIZE OUTPUT ARRAYS AND PARAMETERS
210	guez	3
211	guez	103	do k = 1, klev
212			do i = 1, klon
213	guez	187	ft1(i, k) = 0.
214			fq1(i, k) = 0.
215			fu1(i, k) = 0.
216			fv1(i, k) = 0.
217			tvp1(i, k) = 0.
218			tp1(i, k) = 0.
219			clw1(i, k) = 0.
220			clw(i, k) = 0.
221	guez	103	gz1(i, k) = 0.
222	guez	52	VPrecip1(i, k) = 0.
223	guez	187	Ma1(i, k) = 0.
224			upwd1(i, k) = 0.
225			dnwd1(i, k) = 0.
226			dnwd01(i, k) = 0.
227			qcondc1(i, k) = 0.
228	guez	52	end do
229			end do
230	guez	3
231	guez	103	do i = 1, klon
232	guez	187	precip1(i) = 0.
233	guez	91	iflag1(i) = 0
234	guez	187	wd1(i) = 0.
235			cape1(i) = 0.
236			VPrecip1(i, klev + 1) = 0.
237	guez	52	end do
238	guez	3
239	guez	181	do il = 1, klon
240			sig1(il, klev) = sig1(il, klev) + 1.
241			sig1(il, klev) = min(sig1(il, klev), 12.1)
242			enddo
243	guez	3
244	guez	180	! CALCULATE ARRAYS OF GEOPOTENTIAL, HEAT CAPACITY & STATIC ENERGY
245	guez	185	CALL cv30_prelim(klon, klev, klev + 1, t1, q1, p1, ph1, lv1, cpn1, tv1, &
246	guez	181	gz1, h1, hm1, th1)
247	guez	3
248	guez	180	! CONVECTIVE FEED
249	guez	185	CALL cv30_feed(klon, klev, t1, q1, qs1, p1, ph1, gz1, nk1, icb1, &
250	guez	181	icbmax, iflag1, tnk1, qnk1, gznk1, plcl1) ! klev->na
251	guez	3
252	guez	180	! UNDILUTE (ADIABATIC) UPDRAFT / 1st part
253			! (up through ICB for convect4, up through ICB + 1 for convect3)
254	guez	103	! Calculates the lifted parcel virtual temperature at nk, the
255			! actual temperature, and the adiabatic liquid water content.
256	guez	185	CALL cv30_undilute1(klon, klev, t1, q1, qs1, gz1, plcl1, p1, nk1, icb1, &
257	guez	181	tp1, tvp1, clw1, icbs1) ! klev->na
258	guez	3
259	guez	180	! TRIGGERING
260	guez	185	CALL cv30_trigger(klon, klev, icb1, plcl1, p1, th1, tv1, tvp1, pbase1, &
261	guez	181	buoybase1, iflag1, sig1, w01) ! klev->na
262	guez	3
263	guez	180	! Moist convective adjustment is necessary
264	guez	3
265	guez	91	ncum = 0
266	guez	103	do i = 1, klon
267	guez	180	if (iflag1(i) == 0) then
268			ncum = ncum + 1
269	guez	91	idcum(ncum) = i
270	guez	52	endif
271			end do
272	guez	3
273	guez	103	IF (ncum > 0) THEN
274	guez	180	! COMPRESS THE FIELDS
275	guez	103	! (-> vectorization over convective gridpoints)
276	guez	185	CALL cv30_compress(klon, klon, ncum, klev, iflag1, nk1, icb1, icbs1, &
277	guez	181	plcl1, tnk1, qnk1, gznk1, pbase1, buoybase1, t1, q1, qs1, u1, &
278			v1, gz1, th1, h1, lv1, cpn1, p1, ph1, tv1, tp1, tvp1, clw1, &
279			sig1, w01, iflag, nk, icb, icbs, plcl, tnk, qnk, gznk, pbase, &
280			buoybase, t, q, qs, u, v, gz, th, h, lv, cpn, p, ph, tv, tp, &
281			tvp, clw, sig, w0)
282	guez	3
283	guez	187	CALL cv30_undilute2(ncum, icb, icbs, nk, tnk, qnk, gznk, t, qs, gz, p, &
284			h, tv, lv, pbase, buoybase, plcl, inb(:ncum), tp, tvp, clw, hp, &
285			ep, sigp, buoy)
286	guez	3
287	guez	180	! CLOSURE
288	guez	185	CALL cv30_closure(klon, ncum, klev, icb, inb, pbase, p, ph, tv, &
289	guez	181	buoy, sig, w0, cape, m) ! na->klev
290	guez	3
291	guez	180	! MIXING
292	guez	185	CALL cv30_mixing(klon, ncum, klev, klev, icb, nk, inb, t, q, qs, u, &
293	guez	181	v, h, lv, hp, ep, clw, m, sig, ment, qent, uent, vent, nent, &
294			sij, elij, ments, qents)
295	guez	3
296	guez	186	! Unsaturated (precipitating) downdrafts
297	guez	187	CALL cv30_unsat(ncum, icb(:ncum), inb(:ncum), t, q, qs, gz, u, v, p, &
298			ph, th, tv, lv, cpn, ep, sigp, clw, m, ment, elij, delt, plcl, &
299			mp, qp, up, vp, wt, water, evap, b(:ncum, :))
300	guez	3
301	guez	186	! Yield (tendencies, precipitation, variables of interface with
302			! other processes, etc)
303	guez	185	CALL cv30_yield(klon, ncum, klev, klev, icb, inb, delt, t, q, u, v, &
304	guez	181	gz, p, ph, h, hp, lv, cpn, th, ep, clw, m, tp, mp, qp, up, vp, &
305			wt, water, evap, b, ment, qent, uent, vent, nent, elij, sig, &
306			tv, tvp, iflag, precip, VPrecip, ft, fq, fu, fv, upwd, dnwd, &
307			dnwd0, ma, mike, tls, tps, qcondc, wd)! na->klev
308	guez	3
309	guez	180	! passive tracers
310	guez	185	CALL cv30_tracer(klon, ncum, klev, ment, sij, da, phi)
311	guez	3
312	guez	180	! UNCOMPRESS THE FIELDS
313	guez	103
314			! set iflag1 = 42 for non convective points
315	guez	186	iflag1 = 42
316	guez	62
317	guez	185	CALL cv30_uncompress(idcum(:ncum), iflag, precip, VPrecip, sig, w0, &
318	guez	181	ft, fq, fu, fv, inb, Ma, upwd, dnwd, dnwd0, qcondc, wd, cape, &
319			da, phi, mp, iflag1, precip1, VPrecip1, sig1, w01, ft1, fq1, &
320			fu1, fv1, inb1, Ma1, upwd1, dnwd1, dnwd01, qcondc1, wd1, &
321			cape1, da1, phi1, mp1)
322	guez	183	ENDIF
323	guez	3
324	guez	52	end SUBROUTINE cv_driver
325	guez	3
326	guez	52	end module cv_driver_m