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