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, 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 comconst, only: dtphys
    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_undilute1_m, only: cv30_undilute1
    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, in 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, in hPa

    real, intent(in):: ph1(klon, klev + 1) 
    ! Half level pressure, in 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 <= minorig.

    ! 7: No moist convection because unreasonable parcel level
    ! temperature or specific humidity.

    ! 8: No moist convection: lifted condensation level is above the
    ! 200 mbar level.

    ! 9: No moist convection: cloud base is higher than 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(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)
    real, allocatable:: mp(:, :) ! (ncum, nl)
    integer i, k, il
    integer icbs1(klon)
    real plcl1(klon)
    real tnk1(klon)
    real qnk1(klon)
    real gznk1(klon)
    real pbase1(klon)
    real buoybase1(klon)

    real lv1(klon, nl)
    ! specific latent heat of vaporization of water, in J kg-1

    real cpn1(klon, nl)
    ! specific heat capacity at constant pressure of humid air, in J K-1 kg-1

    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, nl) ! potential temperature, in K
    integer ncum

    ! Compressed fields:
    integer, allocatable:: idcum(:), iflag(:) ! (ncum)
    integer, allocatable:: icb(:) ! (ncum)
    integer nent(klon, klev)
    integer icbs(klon)

    integer, allocatable:: inb(:) ! (ncum)
    ! first model level above the level of neutral buoyancy of the
    ! parcel (1 <= inb <= nl - 1)

    real, allocatable:: plcl(:) ! (ncum)
    real 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)

    real, allocatable:: lv(:, :) ! (ncum, nl)
    ! specific latent heat of vaporization of water, in J kg-1

    real, allocatable:: cpn(:, :) ! (ncum, nl)
    ! specific heat capacity at constant pressure of humid air, in J K-1 kg-1

    real p(klon, klev) ! pressure at full level, in hPa
    real ph(klon, klev + 1), tv(klon, klev), tp(klon, klev)
    real clw(klon, klev)
    real pbase(klon), buoybase(klon)
    real, allocatable:: th(:, :) ! (ncum, nl)
    real tvp(klon, klev)
    real sig(klon, klev), w0(klon, klev)
    real hp(klon, klev), ep(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)
    real, allocatable:: evap(:, :) ! (ncum, nl)
    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
    CALL cv30_param

    ! 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

    precip1 = 0.
    cape1 = 0.
    VPrecip1(:, klev + 1) = 0.

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

    CALL cv30_prelim(t1, q1, p1, ph1, lv1, cpn1, tv1, gz1, h1, hm1, th1)
    CALL cv30_feed(t1, q1, qs1, p1, ph1, gz1, icb1, iflag1, tnk1, qnk1, &
         gznk1, plcl1)
    CALL cv30_undilute1(t1, q1, qs1, gz1, plcl1, p1, icb1, tp1, tvp1, clw1, &
         icbs1)
    CALL cv30_trigger(icb1, plcl1, p1, th1, tv1, tvp1, pbase1, buoybase1, &
         iflag1, sig1, w01)

    ncum = count(iflag1 == 0)

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