/[lmdze]/trunk/phylmd/cv_driver.f
ViewVC logotype

Diff of /trunk/phylmd/cv_driver.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 185 by guez, Wed Mar 16 15:04:46 2016 UTC revision 196 by guez, Mon May 23 13:50:39 2016 UTC
# Line 5  module cv_driver_m Line 5  module cv_driver_m
5  contains  contains
6    
7    SUBROUTINE cv_driver(t1, q1, qs1, u1, v1, p1, ph1, iflag1, ft1, fq1, fu1, &    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, &         fv1, precip1, VPrecip1, sig1, w01, icb1, inb1, Ma1, upwd1, dnwd1, &
9         dnwd1, dnwd01, qcondc1, wd1, cape1, da1, phi1, mp1)         dnwd01, qcondc1, cape1, da1, phi1, mp1)
10    
11      ! From LMDZ4/libf/phylmd/cv_driver.F, version 1.3, 2005/04/15 12:36:17      ! From LMDZ4/libf/phylmd/cv_driver.F, version 1.3, 2005/04/15 12:36:17
12      ! Main driver for convection      ! Main driver for convection
# Line 14  contains Line 14  contains
14    
15      ! Several modules corresponding to different physical processes      ! Several modules corresponding to different physical processes
16    
17        use comconst, only: dtphys
18        use cv30_closure_m, only: cv30_closure
19      use cv30_compress_m, only: cv30_compress      use cv30_compress_m, only: cv30_compress
20      use cv30_feed_m, only: cv30_feed      use cv30_feed_m, only: cv30_feed
21      use cv30_mixing_m, only: cv30_mixing      use cv30_mixing_m, only: cv30_mixing
22      use cv30_param_m, only: cv30_param      use cv30_param_m, only: cv30_param, nl
23      use cv30_prelim_m, only: cv30_prelim      use cv30_prelim_m, only: cv30_prelim
24      use cv30_tracer_m, only: cv30_tracer      use cv30_tracer_m, only: cv30_tracer
25        use cv30_trigger_m, only: cv30_trigger
26      use cv30_uncompress_m, only: cv30_uncompress      use cv30_uncompress_m, only: cv30_uncompress
27        use cv30_undilute1_m, only: cv30_undilute1
28      use cv30_undilute2_m, only: cv30_undilute2      use cv30_undilute2_m, only: cv30_undilute2
29      use cv30_unsat_m, only: cv30_unsat      use cv30_unsat_m, only: cv30_unsat
30      use cv30_yield_m, only: cv30_yield      use cv30_yield_m, only: cv30_yield
31        use cv_thermo_m, only: cv_thermo
32      USE dimphy, ONLY: klev, klon      USE dimphy, ONLY: klev, klon
33    
34      real, intent(in):: t1(klon, klev) ! temperature      real, intent(in):: t1(klon, klev) ! temperature (K)
35      real, intent(in):: q1(klon, klev) ! specific hum      real, intent(in):: q1(klon, klev) ! specific humidity
36      real, intent(in):: qs1(klon, klev) ! sat specific hum      real, intent(in):: qs1(klon, klev) ! saturation specific humidity
37      real, intent(in):: u1(klon, klev) ! u-wind  
38      real, intent(in):: v1(klon, klev) ! v-wind      real, intent(in):: u1(klon, klev), v1(klon, klev)
39      real, intent(in):: p1(klon, klev) ! full level pressure      ! zonal wind and meridional velocity (m/s)
40      real, intent(in):: ph1(klon, klev + 1) ! half level pressure  
41      integer, intent(out):: iflag1(klon) ! flag for Emanuel conditions      real, intent(in):: p1(klon, klev) ! full level pressure (hPa)
42      real, intent(out):: ft1(klon, klev) ! temp tend  
43      real, intent(out):: fq1(klon, klev) ! spec hum tend      real, intent(in):: ph1(klon, klev + 1)
44      real, intent(out):: fu1(klon, klev) ! u-wind tend      ! Half level pressure (hPa). These pressures are defined at levels
45      real, intent(out):: fv1(klon, klev) ! v-wind tend      ! intermediate between those of P1, T1, Q1 and QS1. The first
46      real, intent(out):: precip1(klon) ! precipitation      ! value of PH should be greater than (i.e. at a lower level than)
47        ! the first value of the array P1.
48    
49      real, intent(out):: VPrecip1(klon, klev + 1)      integer, intent(out):: iflag1(:) ! (klon)
50      ! vertical profile of precipitation      ! Flag for Emanuel conditions.
51    
52      real, intent(inout):: sig1(klon, klev) ! section adiabatic updraft      ! 0: Moist convection occurs.
53    
54      real, intent(inout):: w01(klon, klev)      ! 1: Moist convection occurs, but a CFL condition on the
55      ! vertical velocity within adiabatic updraft      ! subsidence warming is violated. This does not cause the scheme
56        ! to terminate.
57    
58      integer, intent(out):: icb1(klon)      ! 2: Moist convection, but no precipitation because ep(inb) < 1e-4
     integer, intent(inout):: inb1(klon)  
     real, intent(in):: delt ! time step  
     real Ma1(klon, klev)  
     ! Ma1 Real Output mass flux adiabatic updraft  
59    
60      real, intent(out):: upwd1(klon, klev)      ! 3: No moist convection because new cbmf is 0 and old cbmf is 0.
     ! total upward mass flux (adiab + mixed)  
61    
62      real, intent(out):: dnwd1(klon, klev) ! saturated downward mass flux (mixed)      ! 4: No moist convection; atmosphere is not unstable.
     real, intent(out):: dnwd01(klon, klev) ! unsaturated downward mass flux  
63    
64      real qcondc1(klon, klev) ! cld      ! 6: No moist convection because ihmin <= minorig.
     ! qcondc1 Real Output in-cld mixing ratio of condensed water  
     real wd1(klon) ! gust  
     ! wd1 Real Output downdraft velocity scale for sfc fluxes  
     real cape1(klon)  
     ! cape1 Real Output CAPE  
65    
66      real, intent(inout):: da1(klon, klev), phi1(klon, klev, klev)      ! 7: No moist convection because unreasonable parcel level
67      real, intent(inout):: mp1(klon, klev)      ! temperature or specific humidity.
68    
69      ! ARGUMENTS      ! 8: No moist convection: lifted condensation level is above the
70        ! 200 mbar level.
71    
72      ! On input:      ! 9: No moist convection: cloud base is higher than the level NL-1.
73    
74      ! t: Array of absolute temperature (K) of dimension KLEV, with first      real, intent(out):: ft1(klon, klev) ! temperature tendency (K/s)
75      ! index corresponding to lowest model level. Note that this array      real, intent(out):: fq1(klon, klev) ! specific humidity tendency (s-1)
     ! will be altered by the subroutine if dry convective adjustment  
     ! occurs and if IPBL is not equal to 0.  
   
     ! q: Array of specific humidity (gm/gm) of dimension KLEV, with first  
     ! index corresponding to lowest model level. Must be defined  
     ! at same grid levels as T. Note that this array will be altered  
     ! if dry convective adjustment occurs and if IPBL is not equal to 0.  
   
     ! qs: Array of saturation specific humidity of dimension KLEV, with first  
     ! index corresponding to lowest model level. Must be defined  
     ! at same grid levels as T. Note that this array will be altered  
     ! if dry convective adjustment occurs and if IPBL is not equal to 0.  
   
     ! u: Array of zonal wind velocity (m/s) of dimension KLEV, witth first  
     ! index corresponding with the lowest model level. Defined at  
     ! same levels as T. Note that this array will be altered if  
     ! dry convective adjustment occurs and if IPBL is not equal to 0.  
   
     ! v: Same as u but for meridional velocity.  
   
     ! p: Array of pressure (mb) of dimension KLEV, with first  
     ! index corresponding to lowest model level. Must be defined  
     ! at same grid levels as T.  
   
     ! ph: Array of pressure (mb) of dimension KLEV + 1, with first index  
     ! corresponding to lowest level. These pressures are defined at  
     ! levels intermediate between those of P, T, Q and QS. The first  
     ! value of PH should be greater than (i.e. at a lower level than)  
     ! the first value of the array P.  
76    
77      ! nl: The maximum number of levels to which convection can penetrate, plus 1      real, intent(out):: fu1(klon, klev), fv1(klon, klev)
78      ! NL MUST be less than or equal to KLEV-1.      ! forcing (tendency) of zonal and meridional velocity (m/s^2)
79    
80      ! delt: The model time step (sec) between calls to CONVECT      real, intent(out):: precip1(klon) ! convective precipitation rate (mm/day)
81    
82      ! On Output:      real, intent(out):: VPrecip1(klon, klev + 1)
83        ! vertical profile of convective precipitation (kg/m2/s)
     ! iflag: An output integer whose value denotes the following:  
     ! VALUE INTERPRETATION  
     ! ----- --------------  
     ! 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 precip because ep(inb) lt 0.0001  
     ! 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.  
   
     ! ft: Array of temperature tendency (K/s) of dimension KLEV, defined at same  
     ! grid levels as T, Q, QS and P.  
   
     ! fq: Array of specific humidity tendencies ((gm/gm)/s) of dimension KLEV,  
     ! defined at same grid levels as T, Q, QS and P.  
   
     ! fu: Array of forcing of zonal velocity (m/s^2) of dimension KLEV,  
     ! defined at same grid levels as T.  
   
     ! fv: Same as FU, but for forcing of meridional velocity.  
84    
85      ! precip: Scalar convective precipitation rate (mm/day).      real, intent(inout):: sig1(klon, klev) ! section of adiabatic updraft
86    
87      ! VPrecip: Vertical profile of convective precipitation (kg/m2/s).      real, intent(inout):: w01(klon, klev)
88        ! vertical velocity within adiabatic updraft
89    
90      ! wd: A convective downdraft velocity scale. For use in surface      integer, intent(out):: icb1(klon)
91      ! flux parameterizations. See convect.ps file for details.      integer, intent(inout):: inb1(klon)
92        real, intent(out):: Ma1(klon, klev) ! mass flux of adiabatic updraft
93    
94      ! tprime: A convective downdraft temperature perturbation scale (K).      real, intent(out):: upwd1(klon, klev)
95      ! For use in surface flux parameterizations. See convect.ps      ! total upward mass flux (adiabatic + mixed)
     ! file for details.  
96    
97      ! qprime: A convective downdraft specific humidity      real, intent(out):: dnwd1(klon, klev) ! saturated downward mass flux (mixed)
98      ! perturbation scale (gm/gm).      real, intent(out):: dnwd01(klon, klev) ! unsaturated downward mass flux
     ! For use in surface flux parameterizations. See convect.ps  
     ! file for details.  
99    
100      ! cbmf: The cloud base mass flux ((kg/m**2)/s). THIS SCALAR VALUE MUST      real, intent(out):: qcondc1(klon, klev)
101      ! BE STORED BY THE CALLING PROGRAM AND RETURNED TO CONVECT AT      ! in-cloud mixing ratio of condensed water
     ! ITS NEXT CALL. That is, the value of CBMF must be "remembered"  
     ! by the calling program between calls to CONVECT.  
102    
103      ! det: Array of detrainment mass flux of dimension KLEV.      real, intent(out):: cape1(klon)
104        real, intent(inout):: da1(klon, klev), phi1(klon, klev, klev)
105        real, intent(inout):: mp1(klon, klev)
106    
107      ! Local arrays      ! Local:
108    
109      real da(klon, klev), phi(klon, klev, klev), mp(klon, klev)      real da(klon, klev), phi(klon, klev, klev), mp(klon, klev)
   
110      integer i, k, il      integer i, k, il
     integer icbmax  
111      integer nk1(klon)      integer nk1(klon)
112      integer icbs1(klon)      integer icbs1(klon)
   
113      real plcl1(klon)      real plcl1(klon)
114      real tnk1(klon)      real tnk1(klon)
115      real qnk1(klon)      real qnk1(klon)
116      real gznk1(klon)      real gznk1(klon)
117      real pbase1(klon)      real pbase1(klon)
118      real buoybase1(klon)      real buoybase1(klon)
   
119      real lv1(klon, klev)      real lv1(klon, klev)
120      real cpn1(klon, klev)      real cpn1(klon, klev)
121      real tv1(klon, klev)      real tv1(klon, klev)
# Line 192  contains Line 126  contains
126      real tvp1(klon, klev)      real tvp1(klon, klev)
127      real clw1(klon, klev)      real clw1(klon, klev)
128      real th1(klon, klev)      real th1(klon, klev)
   
129      integer ncum      integer ncum
130    
131      ! (local) compressed fields:      ! Compressed fields:
132        integer, allocatable:: idcum(:), iflag(:) ! (ncum)
133      integer idcum(klon)      integer nk(klon)
134      integer iflag(klon), nk(klon), icb(klon)      integer, allocatable:: icb(:) ! (ncum)
135      integer nent(klon, klev)      integer nent(klon, klev)
136      integer icbs(klon)      integer icbs(klon)
137      integer inb(klon)      integer inb(klon)
138        real, allocatable:: plcl(:) ! (ncum)
139      real plcl(klon), tnk(klon), qnk(klon), gznk(klon)      real tnk(klon), qnk(klon), gznk(klon)
140      real t(klon, klev), q(klon, klev), qs(klon, klev)      real t(klon, klev), q(klon, klev), qs(klon, klev)
141      real u(klon, klev), v(klon, klev)      real u(klon, klev), v(klon, klev)
142      real gz(klon, klev), h(klon, klev), lv(klon, klev), cpn(klon, klev)      real gz(klon, klev), h(klon, klev), lv(klon, klev), cpn(klon, klev)
143      real p(klon, klev), ph(klon, klev + 1), tv(klon, klev), tp(klon, klev)      real p(klon, klev) ! pressure at full level, in hPa
144        real ph(klon, klev + 1), tv(klon, klev), tp(klon, klev)
145      real clw(klon, klev)      real clw(klon, klev)
146      real pbase(klon), buoybase(klon), th(klon, klev)      real pbase(klon), buoybase(klon), th(klon, klev)
147      real tvp(klon, klev)      real tvp(klon, klev)
148      real sig(klon, klev), w0(klon, klev)      real sig(klon, klev), w0(klon, klev)
149      real hp(klon, klev), ep(klon, klev), sigp(klon, klev)      real hp(klon, klev), ep(klon, klev)
150      real buoy(klon, klev)      real buoy(klon, klev)
151      real cape(klon)      real cape(klon)
152      real m(klon, klev), ment(klon, klev, klev), qent(klon, klev, klev)      real m(klon, klev), ment(klon, klev, klev), qent(klon, klev, klev)
# Line 220  contains Line 154  contains
154      real ments(klon, klev, klev), qents(klon, klev, klev)      real ments(klon, klev, klev), qents(klon, klev, klev)
155      real sij(klon, klev, klev), elij(klon, klev, klev)      real sij(klon, klev, klev), elij(klon, klev, klev)
156      real qp(klon, klev), up(klon, klev), vp(klon, klev)      real qp(klon, klev), up(klon, klev), vp(klon, klev)
157      real wt(klon, klev), water(klon, klev), evap(klon, klev)      real wt(klon, klev), water(klon, klev)
158      real b(klon, klev), ft(klon, klev), fq(klon, klev)      real, allocatable:: evap(:, :) ! (ncum, nl)
159        real, allocatable:: b(:, :) ! (ncum, nl - 1)
160        real ft(klon, klev), fq(klon, klev)
161      real fu(klon, klev), fv(klon, klev)      real fu(klon, klev), fv(klon, klev)
162      real upwd(klon, klev), dnwd(klon, klev), dnwd0(klon, klev)      real upwd(klon, klev), dnwd(klon, klev), dnwd0(klon, klev)
163      real Ma(klon, klev), mike(klon, klev), tls(klon, klev)      real Ma(klon, klev), mike(klon, klev), tls(klon, klev)
# Line 229  contains Line 165  contains
165      real precip(klon)      real precip(klon)
166      real VPrecip(klon, klev + 1)      real VPrecip(klon, klev + 1)
167      real qcondc(klon, klev) ! cld      real qcondc(klon, klev) ! cld
     real wd(klon) ! gust  
168    
169      !-------------------------------------------------------------------      !-------------------------------------------------------------------
170    
171      ! SET CONSTANTS AND PARAMETERS      ! SET CONSTANTS AND PARAMETERS
   
     ! set simulation flags:  
     ! (common cvflag)  
     CALL cv_flag  
   
     ! set thermodynamical constants:  
     ! (common cvthermo)  
172      CALL cv_thermo      CALL cv_thermo
173        CALL cv30_param
     ! set convect parameters  
     ! includes microphysical parameters and parameters that  
     ! control the rate of approach to quasi-equilibrium)  
     ! (common cvparam)  
   
     CALL cv30_param(klev, delt)  
174    
175      ! INITIALIZE OUTPUT ARRAYS AND PARAMETERS      ! INITIALIZE OUTPUT ARRAYS AND PARAMETERS
176    
177      do k = 1, klev      do k = 1, klev
178         do i = 1, klon         do i = 1, klon
179            ft1(i, k) = 0.0            ft1(i, k) = 0.
180            fq1(i, k) = 0.0            fq1(i, k) = 0.
181            fu1(i, k) = 0.0            fu1(i, k) = 0.
182            fv1(i, k) = 0.0            fv1(i, k) = 0.
183            tvp1(i, k) = 0.0            tvp1(i, k) = 0.
184            tp1(i, k) = 0.0            tp1(i, k) = 0.
185            clw1(i, k) = 0.0            clw1(i, k) = 0.
186            clw(i, k) = 0.0            clw(i, k) = 0.
187            gz1(i, k) = 0.            gz1(i, k) = 0.
188            VPrecip1(i, k) = 0.            VPrecip1(i, k) = 0.
189            Ma1(i, k) = 0.0            Ma1(i, k) = 0.
190            upwd1(i, k) = 0.0            upwd1(i, k) = 0.
191            dnwd1(i, k) = 0.0            dnwd1(i, k) = 0.
192            dnwd01(i, k) = 0.0            dnwd01(i, k) = 0.
193            qcondc1(i, k) = 0.0            qcondc1(i, k) = 0.
194         end do         end do
195      end do      end do
196    
197      do i = 1, klon      precip1 = 0.
198         precip1(i) = 0.0      cape1 = 0.
199         iflag1(i) = 0      VPrecip1(:, klev + 1) = 0.
        wd1(i) = 0.0  
        cape1(i) = 0.0  
        VPrecip1(i, klev + 1) = 0.0  
     end do  
200    
201      do il = 1, klon      do il = 1, klon
202         sig1(il, klev) = sig1(il, klev) + 1.         sig1(il, klev) = sig1(il, klev) + 1.
203         sig1(il, klev) = min(sig1(il, klev), 12.1)         sig1(il, klev) = min(sig1(il, klev), 12.1)
204      enddo      enddo
205    
     ! CALCULATE ARRAYS OF GEOPOTENTIAL, HEAT CAPACITY & STATIC ENERGY  
206      CALL cv30_prelim(klon, klev, klev + 1, t1, q1, p1, ph1, lv1, cpn1, tv1, &      CALL cv30_prelim(klon, klev, klev + 1, t1, q1, p1, ph1, lv1, cpn1, tv1, &
207           gz1, h1, hm1, th1)           gz1, h1, hm1, th1)
208        CALL cv30_feed(t1, q1, qs1, p1, ph1, gz1, nk1, icb1, iflag1, tnk1, qnk1, &
209             gznk1, plcl1)
210        CALL cv30_undilute1(t1, q1, qs1, gz1, plcl1, p1, nk1, icb1, tp1, tvp1, &
211             clw1, icbs1)
212        CALL cv30_trigger(icb1, plcl1, p1, th1, tv1, tvp1, pbase1, buoybase1, &
213             iflag1, sig1, w01)
214    
215      ! CONVECTIVE FEED      ncum = count(iflag1 == 0)
     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  
216    
217      IF (ncum > 0) THEN      IF (ncum > 0) THEN
218         ! COMPRESS THE FIELDS         ! Moist convective adjustment is necessary
219         ! (-> vectorization over convective gridpoints)         allocate(idcum(ncum), plcl(ncum))
220         CALL cv30_compress(klon, klon, ncum, klev, iflag1, nk1, icb1, icbs1, &         allocate(b(ncum, nl - 1), evap(ncum, nl), icb(ncum), iflag(ncum))
221              plcl1, tnk1, qnk1, gznk1, pbase1, buoybase1, t1, q1, qs1, u1, &         idcum = pack((/(i, i = 1, klon)/), iflag1 == 0)
222              v1, gz1, th1, h1, lv1, cpn1, p1, ph1, tv1, tp1, tvp1, clw1, &         CALL cv30_compress(iflag1, nk1, icb1, icbs1, plcl1, tnk1, qnk1, gznk1, &
223              sig1, w01, iflag, nk, icb, icbs, plcl, tnk, qnk, gznk, pbase, &              pbase1, buoybase1, t1, q1, qs1, u1, v1, gz1, th1, h1, lv1, cpn1, &
224              buoybase, t, q, qs, u, v, gz, th, h, lv, cpn, p, ph, tv, tp, &              p1, ph1, tv1, tp1, tvp1, clw1, sig1, w01, nk, icb, icbs, plcl, &
225              tvp, clw, sig, w0)              tnk, qnk, gznk, pbase, buoybase, t, q, qs, u, v, gz, th, h, lv, &
226                cpn, p, ph, tv, tp, tvp, clw, sig, w0)
227         ! UNDILUTE (ADIABATIC) UPDRAFT / second part :         CALL cv30_undilute2(icb, icbs(:ncum), nk, tnk, qnk, gznk, t, qs, gz, &
228         ! FIND THE REST OF THE LIFTED PARCEL TEMPERATURES              p, h, tv, lv, pbase(:ncum), buoybase(:ncum), plcl, inb(:ncum), &
229         ! &              tp, tvp, clw, hp, ep, buoy)
230         ! COMPUTE THE PRECIPITATION EFFICIENCIES AND THE         CALL cv30_closure(icb, inb(:ncum), pbase, p, ph(:ncum, :), tv, buoy, &
231         ! FRACTION OF PRECIPITATION FALLING OUTSIDE OF CLOUD              sig, w0, cape, m)
232         ! &         CALL cv30_mixing(icb, nk(:ncum), inb(:ncum), t, q, qs, u, v, h, lv, &
233         ! FIND THE LEVEL OF NEUTRAL BUOYANCY              hp, ep, clw, m, sig, ment, qent, uent, vent, nent, sij, elij, &
234         CALL cv30_undilute2(klon, ncum, klev, icb, icbs, nk, tnk, qnk, gznk, &              ments, qents)
235              t, qs, gz, p, h, tv, lv, pbase, buoybase, plcl, inb, tp, &         CALL cv30_unsat(icb, inb(:ncum), t(:ncum, :nl), q(:ncum, :nl), &
236              tvp, clw, hp, ep, sigp, buoy) !na->klev              qs(:ncum, :nl), gz, u, v, p, ph(:ncum, :), th(:ncum, :nl - 1), &
237                tv, lv, cpn, ep(:ncum, :), clw(:ncum, :), m(:ncum, :), &
238         ! CLOSURE              ment(:ncum, :, :), elij(:ncum, :, :), dtphys, plcl, mp, &
239         CALL cv30_closure(klon, ncum, klev, icb, inb, pbase, p, ph, tv, &              qp(:ncum, :nl), up(:ncum, :nl), vp(:ncum, :nl), wt(:ncum, :nl), &
240              buoy, sig, w0, cape, m) ! na->klev              water(:ncum, :nl), evap, b)
241           CALL cv30_yield(icb, inb(:ncum), dtphys, t, q, u, v, gz, p, ph, h, hp, &
242         ! MIXING              lv, cpn, th, ep, clw, m, tp, mp, qp, up, vp(:ncum, 2:nl), &
243         CALL cv30_mixing(klon, ncum, klev, klev, icb, nk, inb, t, q, qs, u, &              wt(:ncum, :nl - 1), water(:ncum, :nl), evap, b, ment, qent, uent, &
244              v, h, lv, hp, ep, clw, m, sig, ment, qent, uent, vent, nent, &              vent, nent, elij, sig, tv, tvp, iflag, precip, VPrecip, ft, fq, &
245              sij, elij, ments, qents)              fu, fv, upwd, dnwd, dnwd0, ma, mike, tls, tps, qcondc)
   
        ! UNSATURATED (PRECIPITATING) DOWNDRAFTS  
        CALL cv30_unsat(klon, ncum, klev, klev, icb, inb, 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)! na->klev  
   
        ! 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  
246         CALL cv30_tracer(klon, ncum, klev, ment, sij, da, phi)         CALL cv30_tracer(klon, ncum, klev, ment, sij, da, phi)
247           CALL cv30_uncompress(idcum, iflag, precip, VPrecip, sig, w0, ft, fq, &
248         ! UNCOMPRESS THE FIELDS              fu, fv, inb, Ma, upwd, dnwd, dnwd0, qcondc, cape, da, phi, mp, &
249                iflag1, precip1, VPrecip1, sig1, w01, ft1, fq1, fu1, fv1, inb1, &
250         ! set iflag1 = 42 for non convective points              Ma1, upwd1, dnwd1, dnwd01, qcondc1, cape1, da1, phi1, mp1)
        do i = 1, klon  
           iflag1(i) = 42  
        end do  
   
        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)  
251      ENDIF      ENDIF
252    
253    end SUBROUTINE cv_driver    end SUBROUTINE cv_driver
Legend:
Removed from v.185  
changed lines
  Added in v.196
  ViewVC Help
Powered by ViewVC 1.1.21