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, 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_trigger_m, only: cv30_trigger
    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 cv_thermo_m, only: cv_thermo
    USE dimphy, ONLY: klev, klon

    real, intent(in):: t1(klon, klev) ! temperature (K)
    real, intent(in):: q1(klon, klev) ! specific humidity
    real, intent(in):: qs1(klon, klev) ! saturation specific humidity

    real, intent(in):: u1(klon, klev), v1(klon, klev)
    ! zonal wind and meridional velocity (m/s)

    real, intent(in):: p1(klon, klev) ! full level pressure (hPa)

    real, intent(in):: ph1(klon, klev + 1) 
    ! Half level pressure (hPa). 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)
    real, intent(out):: fq1(klon, klev) ! specific humidity tendency (s-1)

    real, intent(out):: fu1(klon, klev), fv1(klon, klev)
    ! forcing (tendency) of zonal and meridional velocity (m/s^2)

    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 of 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, intent(out):: Ma1(klon, klev) ! mass flux of adiabatic updraft

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

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

    real, intent(out):: qcondc1(klon, klev)
    ! in-cloud mixing ratio of condensed water

    real, intent(out):: cape1(klon)
    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 - 1)
    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

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

    ! SET CONSTANTS AND PARAMETERS

    ! set thermodynamical constants:
    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
       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

    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 - 1))
       CALL cv30_compress(ncum, 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(:ncum, :), sigp(:ncum, :), clw(:ncum, :), &
            m(:ncum, :), ment(:ncum, :, :), elij(:ncum, :, :), delt, plcl, &
            mp, qp(:ncum, :nl), up(:ncum, :nl), vp(:ncum, :nl), 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)

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

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