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

Diff of /trunk/Sources/phylmd/cv_driver.f

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

trunk/phylmd/cv_driver.f revision 99 by guez, Wed Jul 2 18:39:15 2014 UTC trunk/Sources/phylmd/cv_driver.f revision 180 by guez, Tue Mar 15 17:07:47 2016 UTC
# Line 4  module cv_driver_m Line 4  module cv_driver_m
4    
5  contains  contains
6    
7    SUBROUTINE cv_driver(len, nd, t1, q1, qs1, u1, v1, p1, ph1, iflag1, ft1, &    SUBROUTINE cv_driver(t1, q1, qs1, u1, v1, p1, ph1, iflag1, ft1, &
8         fq1, fu1, fv1, precip1, VPrecip1, cbmf1, sig1, w01, icb1, inb1, delt, &         fq1, fu1, fv1, precip1, VPrecip1, cbmf1, sig1, w01, icb1, inb1, delt, &
9         Ma1, upwd1, dnwd1, dnwd01, qcondc1, wd1, cape1, da1, phi1, mp1)         Ma1, upwd1, dnwd1, dnwd01, qcondc1, wd1, cape1, da1, phi1, mp1)
10    
# Line 20  contains Line 20  contains
20    
21      use clesphys2, only: iflag_con      use clesphys2, only: iflag_con
22      use cv3_compress_m, only: cv3_compress      use cv3_compress_m, only: cv3_compress
23        use cv3_feed_m, only: cv3_feed
24      use cv3_mixing_m, only: cv3_mixing      use cv3_mixing_m, only: cv3_mixing
25      use cv3_param_m, only: cv3_param      use cv3_param_m, only: cv3_param
26      use cv3_prelim_m, only: cv3_prelim      use cv3_prelim_m, only: cv3_prelim
# Line 27  contains Line 28  contains
28      use cv3_uncompress_m, only: cv3_uncompress      use cv3_uncompress_m, only: cv3_uncompress
29      use cv3_unsat_m, only: cv3_unsat      use cv3_unsat_m, only: cv3_unsat
30      use cv3_yield_m, only: cv3_yield      use cv3_yield_m, only: cv3_yield
31        use cv_feed_m, only: cv_feed
32      use cv_uncompress_m, only: cv_uncompress      use cv_uncompress_m, only: cv_uncompress
33      USE dimphy, ONLY: klev, klon      USE dimphy, ONLY: klev, klon
34    
35      integer, intent(in):: len ! first dimension      real, intent(in):: t1(klon, klev) ! temperature
36      integer, intent(in):: nd ! vertical dimension      real, intent(in):: q1(klon, klev) ! specific hum
37      real, intent(in):: t1(len, nd) ! temperature      real, intent(in):: qs1(klon, klev) ! sat specific hum
38      real q1(len, nd) !           Input        specific hum      real, intent(in):: u1(klon, klev) ! u-wind
39      real qs1(len, nd)      real, intent(in):: v1(klon, klev) ! v-wind
40      !      qs1           Real           Input        sat specific hum      real, intent(in):: p1(klon, klev) ! full level pressure
41      real, intent(in):: u1(len, nd)      real, intent(in):: ph1(klon, klev + 1) ! half level pressure
42      !      u1            Real           Input        u-wind      integer, intent(out):: iflag1(klon) ! flag for Emanuel conditions
43      real, intent(in):: v1(len, nd)      real, intent(out):: ft1(klon, klev) ! temp tend
44      !      v1            Real           Input        v-wind      real, intent(out):: fq1(klon, klev) ! spec hum tend
45      real p1(len, nd)      real, intent(out):: fu1(klon, klev) ! u-wind tend
46      !      p1            Real           Input        full level pressure      real, intent(out):: fv1(klon, klev) ! v-wind tend
47      real ph1(len, nd + 1)      real, intent(out):: precip1(klon) ! precipitation
48      !      ph1           Real           Input        half level pressure  
49      integer iflag1(len)      real, intent(out):: VPrecip1(klon, klev + 1)
50      !      iflag1        Integer        Output       flag for Emanuel conditions      ! vertical profile of precipitation
51      real ft1(len, nd)  
52      !      ft1           Real           Output       temp tend      real, intent(inout):: cbmf1(klon) ! cloud base mass flux
     real fq1(len, nd)  
     !      fq1           Real           Output       spec hum tend  
     real fu1(len, nd)  
     !      fu1           Real           Output       u-wind tend  
     real fv1(len, nd)  
     !      fv1           Real           Output       v-wind tend  
     real precip1(len)  
     !      precip1       Real           Output       precipitation  
     real VPrecip1(len, nd+1)  
     !      VPrecip1      Real           Output       vertical profile of precipitations  
     real cbmf1(len)  
     !      cbmf1         Real           Output       cloud base mass flux  
53      real, intent(inout):: sig1(klon, klev) ! section adiabatic updraft      real, intent(inout):: sig1(klon, klev) ! section adiabatic updraft
54    
55      real, intent(inout):: w01(klon, klev)      real, intent(inout):: w01(klon, klev)
56      ! vertical velocity within adiabatic updraft      ! vertical velocity within adiabatic updraft
57    
58      integer icb1(klon)      integer, intent(out):: icb1(klon)
59      integer inb1(klon)      integer, intent(inout):: inb1(klon)
60      real, intent(in):: delt      real, intent(in):: delt ! time step
61      !      delt          Real           Input        time step      real Ma1(klon, klev)
62      real Ma1(len, nd)      ! Ma1 Real Output mass flux adiabatic updraft
63      !      Ma1           Real           Output       mass flux adiabatic updraft  
64      real, intent(out):: upwd1(len, nd) ! total upward mass flux (adiab+mixed)      real, intent(out):: upwd1(klon, klev)
65      real, intent(out):: dnwd1(len, nd) ! saturated downward mass flux (mixed)      ! total upward mass flux (adiab + mixed)
66      real, intent(out):: dnwd01(len, nd) ! unsaturated downward mass flux  
67        real, intent(out):: dnwd1(klon, klev) ! saturated downward mass flux (mixed)
68      real qcondc1(len, nd)     ! cld      real, intent(out):: dnwd01(klon, klev) ! unsaturated downward mass flux
69      !      qcondc1       Real           Output       in-cld mixing ratio of condensed water  
70      real wd1(len)            ! gust      real qcondc1(klon, klev) ! cld
71      !      wd1           Real           Output       downdraft velocity scale for sfc fluxes      ! qcondc1 Real Output in-cld mixing ratio of condensed water
72      real cape1(len)      real wd1(klon) ! gust
73      !      cape1         Real           Output       CAPE      ! wd1 Real Output downdraft velocity scale for sfc fluxes
74        real cape1(klon)
75      real, intent(inout):: da1(len, nd), phi1(len, nd, nd), mp1(len, nd)      ! cape1 Real Output CAPE
76    
77      !-------------------------------------------------------------------      real, intent(inout):: da1(klon, klev), phi1(klon, klev, klev)
78      ! --- ARGUMENTS      real, intent(inout):: mp1(klon, klev)
79      !-------------------------------------------------------------------  
80      ! --- On input:      ! ARGUMENTS
81    
82      !  t:   Array of absolute temperature (K) of dimension ND, with first      ! On input:
83      !       index corresponding to lowest model level. Note that this array  
84      !       will be altered by the subroutine if dry convective adjustment      ! t: Array of absolute temperature (K) of dimension KLEV, with first
85      !       occurs and if IPBL is not equal to 0.      ! index corresponding to lowest model level. Note that this array
86        ! will be altered by the subroutine if dry convective adjustment
87      !  q:   Array of specific humidity (gm/gm) of dimension ND, with first      ! occurs and if IPBL is not equal to 0.
88      !       index corresponding to lowest model level. Must be defined  
89      !       at same grid levels as T. Note that this array will be altered      ! q: Array of specific humidity (gm/gm) of dimension KLEV, with first
90      !       if dry convective adjustment occurs and if IPBL is not equal to 0.      ! index corresponding to lowest model level. Must be defined
91        ! at same grid levels as T. Note that this array will be altered
92      !  qs:  Array of saturation specific humidity of dimension ND, with first      ! if dry convective adjustment occurs and if IPBL is not equal to 0.
93      !       index corresponding to lowest model level. Must be defined  
94      !       at same grid levels as T. Note that this array will be altered      ! qs: Array of saturation specific humidity of dimension KLEV, with first
95      !       if dry convective adjustment occurs and if IPBL is not equal to 0.      ! index corresponding to lowest model level. Must be defined
96        ! at same grid levels as T. Note that this array will be altered
97      !  u:   Array of zonal wind velocity (m/s) of dimension ND, witth first      ! if dry convective adjustment occurs and if IPBL is not equal to 0.
98      !       index corresponding with the lowest model level. Defined at  
99      !       same levels as T. Note that this array will be altered if      ! u: Array of zonal wind velocity (m/s) of dimension KLEV, witth first
100      !       dry convective adjustment occurs and if IPBL is not equal to 0.      ! index corresponding with the lowest model level. Defined at
101        ! same levels as T. Note that this array will be altered if
102      !  v:   Same as u but for meridional velocity.      ! dry convective adjustment occurs and if IPBL is not equal to 0.
103    
104      !  p:   Array of pressure (mb) of dimension ND, with first      ! v: Same as u but for meridional velocity.
105      !       index corresponding to lowest model level. Must be defined  
106      !       at same grid levels as T.      ! p: Array of pressure (mb) of dimension KLEV, with first
107        ! index corresponding to lowest model level. Must be defined
108      !  ph:  Array of pressure (mb) of dimension ND+1, with first index      ! at same grid levels as T.
109      !       corresponding to lowest level. These pressures are defined at  
110      !       levels intermediate between those of P, T, Q and QS. The first      ! ph: Array of pressure (mb) of dimension KLEV + 1, with first index
111      !       value of PH should be greater than (i.e. at a lower level than)      ! corresponding to lowest level. These pressures are defined at
112      !       the first value of the array P.      ! levels intermediate between those of P, T, Q and QS. The first
113        ! value of PH should be greater than (i.e. at a lower level than)
114      !  nl:  The maximum number of levels to which convection can penetrate, plus 1.      ! the first value of the array P.
115      !       NL MUST be less than or equal to ND-1.  
116        ! nl: The maximum number of levels to which convection can penetrate, plus 1
117      !  delt: The model time step (sec) between calls to CONVECT      ! NL MUST be less than or equal to KLEV-1.
118    
119      !----------------------------------------------------------------------------      ! delt: The model time step (sec) between calls to CONVECT
120      ! ---   On Output:  
121        ! On Output:
122      !  iflag: An output integer whose value denotes the following:  
123      !       VALUE   INTERPRETATION      ! iflag: An output integer whose value denotes the following:
124      !       -----   --------------      ! VALUE INTERPRETATION
125      !         0     Moist convection occurs.      ! ----- --------------
126      !         1     Moist convection occurs, but a CFL condition      ! 0 Moist convection occurs.
127      !               on the subsidence warming is violated. This      ! 1 Moist convection occurs, but a CFL condition
128      !               does not cause the scheme to terminate.      ! on the subsidence warming is violated. This
129      !         2     Moist convection, but no precip because ep(inb) lt 0.0001      ! does not cause the scheme to terminate.
130      !         3     No moist convection because new cbmf is 0 and old cbmf is 0.      ! 2 Moist convection, but no precip because ep(inb) lt 0.0001
131      !         4     No moist convection; atmosphere is not      ! 3 No moist convection because new cbmf is 0 and old cbmf is 0.
132      !               unstable      ! 4 No moist convection; atmosphere is not
133      !         6     No moist convection because ihmin le minorig.      ! unstable
134      !         7     No moist convection because unreasonable      ! 6 No moist convection because ihmin le minorig.
135      !               parcel level temperature or specific humidity.      ! 7 No moist convection because unreasonable
136      !         8     No moist convection: lifted condensation      ! parcel level temperature or specific humidity.
137      !               level is above the 200 mb level.      ! 8 No moist convection: lifted condensation
138      !         9     No moist convection: cloud base is higher      ! level is above the 200 mb level.
139      !               then the level NL-1.      ! 9 No moist convection: cloud base is higher
140        ! then the level NL-1.
141      !  ft:   Array of temperature tendency (K/s) of dimension ND, defined at same  
142      !        grid levels as T, Q, QS and P.      ! ft: Array of temperature tendency (K/s) of dimension KLEV, defined at same
143        ! grid levels as T, Q, QS and P.
144      !  fq:   Array of specific humidity tendencies ((gm/gm)/s) of dimension ND,  
145      !        defined at same grid levels as T, Q, QS and P.      ! fq: Array of specific humidity tendencies ((gm/gm)/s) of dimension KLEV,
146        ! defined at same grid levels as T, Q, QS and P.
147      !  fu:   Array of forcing of zonal velocity (m/s^2) of dimension ND,  
148      !        defined at same grid levels as T.      ! fu: Array of forcing of zonal velocity (m/s^2) of dimension KLEV,
149        ! defined at same grid levels as T.
150      !  fv:   Same as FU, but for forcing of meridional velocity.  
151        ! fv: Same as FU, but for forcing of meridional velocity.
152      !  precip: Scalar convective precipitation rate (mm/day).  
153        ! precip: Scalar convective precipitation rate (mm/day).
154      !  VPrecip: Vertical profile of convective precipitation (kg/m2/s).  
155        ! VPrecip: Vertical profile of convective precipitation (kg/m2/s).
156      !  wd:   A convective downdraft velocity scale. For use in surface  
157      !        flux parameterizations. See convect.ps file for details.      ! wd: A convective downdraft velocity scale. For use in surface
158        ! flux parameterizations. See convect.ps file for details.
159      !  tprime: A convective downdraft temperature perturbation scale (K).  
160      !          For use in surface flux parameterizations. See convect.ps      ! tprime: A convective downdraft temperature perturbation scale (K).
161      !          file for details.      ! For use in surface flux parameterizations. See convect.ps
162        ! file for details.
163      !  qprime: A convective downdraft specific humidity  
164      !          perturbation scale (gm/gm).      ! qprime: A convective downdraft specific humidity
165      !          For use in surface flux parameterizations. See convect.ps      ! perturbation scale (gm/gm).
166      !          file for details.      ! For use in surface flux parameterizations. See convect.ps
167        ! file for details.
168      !  cbmf: The cloud base mass flux ((kg/m**2)/s). THIS SCALAR VALUE MUST  
169      !        BE STORED BY THE CALLING PROGRAM AND RETURNED TO CONVECT AT      ! cbmf: The cloud base mass flux ((kg/m**2)/s). THIS SCALAR VALUE MUST
170      !        ITS NEXT CALL. That is, the value of CBMF must be "remembered"      ! BE STORED BY THE CALLING PROGRAM AND RETURNED TO CONVECT AT
171      !        by the calling program between calls to CONVECT.      ! ITS NEXT CALL. That is, the value of CBMF must be "remembered"
172        ! by the calling program between calls to CONVECT.
173    
174      !  det:   Array of detrainment mass flux of dimension ND.      ! det: Array of detrainment mass flux of dimension KLEV.
175    
176      !-------------------------------------------------------------------      ! Local arrays
   
     !  Local arrays  
177    
178      real da(len, nd), phi(len, nd, nd), mp(len, nd)      real da(klon, klev), phi(klon, klev, klev), mp(klon, klev)
179    
180      integer i, k, il      integer i, k, il
181      integer icbmax      integer icbmax
# Line 215  contains Line 204  contains
204    
205      ! (local) compressed fields:      ! (local) compressed fields:
206    
207      integer nloc      integer idcum(klon)
208      parameter (nloc = klon) ! pour l'instant      integer iflag(klon), nk(klon), icb(klon)
209        integer nent(klon, klev)
210        integer icbs(klon)
211        integer inb(klon), inbis(klon)
212    
213        real cbmf(klon), plcl(klon), tnk(klon), qnk(klon), gznk(klon)
214        real t(klon, klev), q(klon, klev), qs(klon, klev)
215        real u(klon, klev), v(klon, klev)
216        real gz(klon, klev), h(klon, klev), lv(klon, klev), cpn(klon, klev)
217        real p(klon, klev), ph(klon, klev + 1), tv(klon, klev), tp(klon, klev)
218        real clw(klon, klev)
219        real dph(klon, klev)
220        real pbase(klon), buoybase(klon), th(klon, klev)
221        real tvp(klon, klev)
222        real sig(klon, klev), w0(klon, klev)
223        real hp(klon, klev), ep(klon, klev), sigp(klon, klev)
224        real frac(klon), buoy(klon, klev)
225        real cape(klon)
226        real m(klon, klev), ment(klon, klev, klev), qent(klon, klev, klev)
227        real uent(klon, klev, klev), vent(klon, klev, klev)
228        real ments(klon, klev, klev), qents(klon, klev, klev)
229        real sij(klon, klev, klev), elij(klon, klev, klev)
230        real qp(klon, klev), up(klon, klev), vp(klon, klev)
231        real wt(klon, klev), water(klon, klev), evap(klon, klev)
232        real b(klon, klev), ft(klon, klev), fq(klon, klev)
233        real fu(klon, klev), fv(klon, klev)
234        real upwd(klon, klev), dnwd(klon, klev), dnwd0(klon, klev)
235        real Ma(klon, klev), mike(klon, klev), tls(klon, klev)
236        real tps(klon, klev), qprime(klon), tprime(klon)
237        real precip(klon)
238        real VPrecip(klon, klev + 1)
239        real qcondc(klon, klev) ! cld
240        real wd(klon) ! gust
241    
     integer idcum(nloc)  
     integer iflag(nloc), nk(nloc), icb(nloc)  
     integer nent(nloc, klev)  
     integer icbs(nloc)  
     integer inb(nloc), inbis(nloc)  
   
     real cbmf(nloc), plcl(nloc), tnk(nloc), qnk(nloc), gznk(nloc)  
     real t(nloc, klev), q(nloc, klev), qs(nloc, klev)  
     real u(nloc, klev), v(nloc, klev)  
     real gz(nloc, klev), h(nloc, klev), lv(nloc, klev), cpn(nloc, klev)  
     real p(nloc, klev), ph(nloc, klev+1), tv(nloc, klev), tp(nloc, klev)  
     real clw(nloc, klev)  
     real dph(nloc, klev)  
     real pbase(nloc), buoybase(nloc), th(nloc, klev)  
     real tvp(nloc, klev)  
     real sig(nloc, klev), w0(nloc, klev)  
     real hp(nloc, klev), ep(nloc, klev), sigp(nloc, klev)  
     real frac(nloc), buoy(nloc, klev)  
     real cape(nloc)  
     real m(nloc, klev), ment(nloc, klev, klev), qent(nloc, klev, klev)  
     real uent(nloc, klev, klev), vent(nloc, klev, klev)  
     real ments(nloc, klev, klev), qents(nloc, klev, klev)  
     real sij(nloc, klev, klev), elij(nloc, klev, klev)  
     real qp(nloc, klev), up(nloc, klev), vp(nloc, klev)  
     real wt(nloc, klev), water(nloc, klev), evap(nloc, klev)  
     real b(nloc, klev), ft(nloc, klev), fq(nloc, klev)  
     real fu(nloc, klev), fv(nloc, klev)  
     real upwd(nloc, klev), dnwd(nloc, klev), dnwd0(nloc, klev)  
     real Ma(nloc, klev), mike(nloc, klev), tls(nloc, klev)  
     real tps(nloc, klev), qprime(nloc), tprime(nloc)  
     real precip(nloc)  
     real VPrecip(nloc, klev+1)  
     real qcondc(nloc, klev)  ! cld  
     real wd(nloc)           ! gust  
   
     !-------------------------------------------------------------------  
     ! --- SET CONSTANTS AND PARAMETERS  
242      !-------------------------------------------------------------------      !-------------------------------------------------------------------
243    
244      ! -- set simulation flags:      ! SET CONSTANTS AND PARAMETERS
245      !   (common cvflag)  
246        ! set simulation flags:
247        ! (common cvflag)
248    
249      CALL cv_flag      CALL cv_flag
250    
251      ! -- set thermodynamical constants:      ! set thermodynamical constants:
252      !     (common cvthermo)      ! (common cvthermo)
253    
254      CALL cv_thermo      CALL cv_thermo
255    
256      ! -- set convect parameters      ! set convect parameters
   
     !     includes microphysical parameters and parameters that  
     !     control the rate of approach to quasi-equilibrium)  
     !     (common cvparam)  
257    
258      if (iflag_con.eq.3) then      ! includes microphysical parameters and parameters that
259         CALL cv3_param(nd, delt)      ! control the rate of approach to quasi-equilibrium)
260      endif      ! (common cvparam)
261    
262      if (iflag_con.eq.4) then      if (iflag_con == 3) CALL cv3_param(klev, delt)
        CALL cv_param(nd)  
     endif  
263    
264      !---------------------------------------------------------------------      ! INITIALIZE OUTPUT ARRAYS AND PARAMETERS
     ! --- INITIALIZE OUTPUT ARRAYS AND PARAMETERS  
     !---------------------------------------------------------------------  
265    
266      do k = 1, nd      do k = 1, klev
267         do  i = 1, len         do i = 1, klon
268            ft1(i, k) = 0.0            ft1(i, k) = 0.0
269            fq1(i, k) = 0.0            fq1(i, k) = 0.0
270            fu1(i, k) = 0.0            fu1(i, k) = 0.0
# Line 296  contains Line 274  contains
274            clw1(i, k) = 0.0            clw1(i, k) = 0.0
275            !ym            !ym
276            clw(i, k) = 0.0            clw(i, k) = 0.0
277            gz1(i, k)  =  0.            gz1(i, k) = 0.
278            VPrecip1(i, k) = 0.            VPrecip1(i, k) = 0.
279            Ma1(i, k) = 0.0            Ma1(i, k) = 0.0
280            upwd1(i, k) = 0.0            upwd1(i, k) = 0.0
# Line 306  contains Line 284  contains
284         end do         end do
285      end do      end do
286    
287      do  i = 1, len      do i = 1, klon
288         precip1(i) = 0.0         precip1(i) = 0.0
289         iflag1(i) = 0         iflag1(i) = 0
290         wd1(i) = 0.0         wd1(i) = 0.0
291         cape1(i) = 0.0         cape1(i) = 0.0
292         VPrecip1(i, nd+1) = 0.0         VPrecip1(i, klev + 1) = 0.0
293      end do      end do
294    
295      if (iflag_con.eq.3) then      if (iflag_con == 3) then
296         do il = 1, len         do il = 1, klon
297            sig1(il, nd) = sig1(il, nd) + 1.            sig1(il, klev) = sig1(il, klev) + 1.
298            sig1(il, nd)  =  min(sig1(il, nd), 12.1)            sig1(il, klev) = min(sig1(il, klev), 12.1)
299         enddo         enddo
300      endif      endif
301    
302      !--------------------------------------------------------------------      ! CALCULATE ARRAYS OF GEOPOTENTIAL, HEAT CAPACITY & STATIC ENERGY
     ! --- CALCULATE ARRAYS OF GEOPOTENTIAL, HEAT CAPACITY & STATIC ENERGY  
     !--------------------------------------------------------------------  
   
     if (iflag_con.eq.3) then  
        CALL cv3_prelim(len, nd, nd + 1, t1, q1, p1, ph1, lv1, cpn1, tv1, gz1, &  
             h1, hm1, th1)  
     endif  
   
     if (iflag_con.eq.4) then  
        CALL cv_prelim(len, nd, nd + 1, t1, q1, p1, ph1 &  
             , lv1, cpn1, tv1, gz1, h1, hm1)  
     endif  
   
     !--------------------------------------------------------------------  
     ! --- CONVECTIVE FEED  
     !--------------------------------------------------------------------  
   
     if (iflag_con.eq.3) then  
        CALL cv3_feed(len, nd, t1, q1, qs1, p1, ph1, hm1, gz1            &  
             , nk1, icb1, icbmax, iflag1, tnk1, qnk1, gznk1, plcl1) ! nd->na  
     endif  
   
     if (iflag_con.eq.4) then  
        CALL cv_feed(len, nd, t1, q1, qs1, p1, hm1, gz1 &  
             , nk1, icb1, icbmax, iflag1, tnk1, qnk1, gznk1, plcl1)  
     endif  
   
     !--------------------------------------------------------------------  
     ! --- 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.  
     !--------------------------------------------------------------------  
   
     if (iflag_con.eq.3) then  
        CALL cv3_undilute1(len, nd, t1, q1, qs1, gz1, plcl1, p1, nk1, icb1   &  
             , tp1, tvp1, clw1, icbs1) ! nd->na  
     endif  
   
     if (iflag_con.eq.4) then  
        CALL cv_undilute1(len, nd, t1, q1, qs1, gz1, p1, nk1, icb1, icbmax &  
             , tp1, tvp1, clw1)  
     endif  
   
     !-------------------------------------------------------------------  
     ! --- TRIGGERING  
     !-------------------------------------------------------------------  
   
     if (iflag_con.eq.3) then  
        CALL cv3_trigger(len, nd, icb1, plcl1, p1, th1, tv1, tvp1, pbase1, &  
             buoybase1, iflag1, sig1, w01) ! nd->na  
     endif  
303    
304      if (iflag_con.eq.4) then      if (iflag_con == 3) then
305         CALL cv_trigger(len, nd, icb1, cbmf1, tv1, tvp1, iflag1)         CALL cv3_prelim(klon, klev, klev + 1, t1, q1, p1, ph1, lv1, cpn1, tv1, &
306      endif              gz1, h1, hm1, th1)
307        else
308           ! iflag_con == 4
309           CALL cv_prelim(klon, klev, klev + 1, t1, q1, p1, ph1, lv1, cpn1, tv1, &
310                gz1, h1, hm1)
311        endif
312    
313        ! CONVECTIVE FEED
314    
315        if (iflag_con == 3) then
316           CALL cv3_feed(klon, klev, t1, q1, qs1, p1, ph1, gz1, nk1, icb1, &
317                icbmax, iflag1, tnk1, qnk1, gznk1, plcl1) ! klev->na
318        else
319           ! iflag_con == 4
320           CALL cv_feed(klon, klev, t1, q1, qs1, p1, hm1, gz1, nk1, icb1, icbmax, &
321                iflag1, tnk1, qnk1, gznk1, plcl1)
322        endif
323    
324        ! UNDILUTE (ADIABATIC) UPDRAFT / 1st part
325        ! (up through ICB for convect4, up through ICB + 1 for convect3)
326        ! Calculates the lifted parcel virtual temperature at nk, the
327        ! actual temperature, and the adiabatic liquid water content.
328    
329        if (iflag_con == 3) then
330           CALL cv3_undilute1(klon, klev, t1, q1, qs1, gz1, plcl1, p1, nk1, icb1, &
331                tp1, tvp1, clw1, icbs1) ! klev->na
332        else
333           ! iflag_con == 4
334           CALL cv_undilute1(klon, klev, t1, q1, qs1, gz1, p1, nk1, icb1, icbmax, &
335                tp1, tvp1, clw1)
336        endif
337    
338        ! TRIGGERING
339    
340        if (iflag_con == 3) then
341           CALL cv3_trigger(klon, klev, icb1, plcl1, p1, th1, tv1, tvp1, pbase1, &
342                buoybase1, iflag1, sig1, w01) ! klev->na
343        else
344           ! iflag_con == 4
345           CALL cv_trigger(klon, klev, icb1, cbmf1, tv1, tvp1, iflag1)
346        end if
347    
348      ! --- IF THIS POINT IS REACHED, MOIST CONVECTIVE ADJUSTMENT IS NECESSARY      ! Moist convective adjustment is necessary
349    
350      ncum = 0      ncum = 0
351      do  i = 1, len      do i = 1, klon
352         if(iflag1(i).eq.0)then         if (iflag1(i) == 0) then
353            ncum = ncum+1            ncum = ncum + 1
354            idcum(ncum) = i            idcum(ncum) = i
355         endif         endif
356      end do      end do
357    
358      IF (ncum.gt.0) THEN      IF (ncum > 0) THEN
359           ! COMPRESS THE FIELDS
360           ! (-> vectorization over convective gridpoints)
361    
362         !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^         if (iflag_con == 3) then
363         ! --- COMPRESS THE FIELDS            CALL cv3_compress(klon, klon, ncum, klev, iflag1, nk1, icb1, icbs1, &
        !        (-> vectorization over convective gridpoints)  
        !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^  
   
        if (iflag_con.eq.3) then  
           CALL cv3_compress(len, nloc, ncum, nd, iflag1, nk1, icb1, icbs1, &  
364                 plcl1, tnk1, qnk1, gznk1, pbase1, buoybase1, t1, q1, qs1, u1, &                 plcl1, tnk1, qnk1, gznk1, pbase1, buoybase1, t1, q1, qs1, u1, &
365                 v1, gz1, th1, h1, lv1, cpn1, p1, ph1, tv1, tp1, tvp1, clw1, &                 v1, gz1, th1, h1, lv1, cpn1, p1, ph1, tv1, tp1, tvp1, clw1, &
366                 sig1, w01, iflag, nk, icb, icbs, plcl, tnk, qnk, gznk, pbase, &                 sig1, w01, iflag, nk, icb, icbs, plcl, tnk, qnk, gznk, pbase, &
367                 buoybase, t, q, qs, u, v, gz, th, h, lv, cpn, p, ph, tv, tp, &                 buoybase, t, q, qs, u, v, gz, th, h, lv, cpn, p, ph, tv, tp, &
368                 tvp, clw, sig, w0)                 tvp, clw, sig, w0)
369           else
370              ! iflag_con == 4
371              CALL cv_compress(klon, klon, ncum, klev, iflag1, nk1, icb1, cbmf1, &
372                   plcl1, tnk1, qnk1, gznk1, t1, q1, qs1, u1, v1, gz1, h1, lv1, &
373                   cpn1, p1, ph1, tv1, tp1, tvp1, clw1, iflag, nk, icb, cbmf, &
374                   plcl, tnk, qnk, gznk, t, q, qs, u, v, gz, h, lv, cpn, p, ph, &
375                   tv, tp, tvp, clw, dph)
376           endif
377    
378           ! UNDILUTE (ADIABATIC) UPDRAFT / second part :
379           ! FIND THE REST OF THE LIFTED PARCEL TEMPERATURES
380           ! &
381           ! COMPUTE THE PRECIPITATION EFFICIENCIES AND THE
382           ! FRACTION OF PRECIPITATION FALLING OUTSIDE OF CLOUD
383           ! &
384           ! FIND THE LEVEL OF NEUTRAL BUOYANCY
385    
386           if (iflag_con == 3) then
387              CALL cv3_undilute2(klon, ncum, klev, icb, icbs, nk, tnk, qnk, gznk, &
388                   t, qs, gz, p, h, tv, lv, pbase, buoybase, plcl, inb, tp, &
389                   tvp, clw, hp, ep, sigp, buoy) !na->klev
390           else
391              ! iflag_con == 4
392              CALL cv_undilute2(klon, ncum, klev, icb, nk, tnk, qnk, gznk, t, &
393                   qs, gz, p, dph, h, tv, lv, inb, inbis, tp, tvp, clw, hp, ep, &
394                   sigp, frac)
395           endif
396    
397           ! CLOSURE
398    
399           if (iflag_con == 3) then
400              CALL cv3_closure(klon, ncum, klev, icb, inb, pbase, p, ph, tv, &
401                   buoy, sig, w0, cape, m) ! na->klev
402           else
403              ! iflag_con == 4
404              CALL cv_closure(klon, ncum, klev, nk, icb, tv, tvp, p, ph, dph, &
405                   plcl, cpn, iflag, cbmf)
406           endif
407    
408           ! MIXING
409    
410           if (iflag_con == 3) then
411              CALL cv3_mixing(klon, ncum, klev, klev, icb, nk, inb, t, q, qs, u, &
412                   v, h, lv, hp, ep, clw, m, sig, ment, qent, uent, vent, nent, &
413                   sij, elij, ments, qents)
414           else
415              ! iflag_con == 4
416              CALL cv_mixing(klon, ncum, klev, icb, nk, inb, inbis, ph, t, q, qs, &
417                   u, v, h, lv, qnk, hp, tv, tvp, ep, clw, cbmf, m, ment, qent, &
418                   uent, vent, nent, sij, elij)
419           endif
420    
421           ! UNSATURATED (PRECIPITATING) DOWNDRAFTS
422    
423           if (iflag_con == 3) then
424              CALL cv3_unsat(klon, ncum, klev, klev, icb, inb, t, q, qs, gz, u, &
425                   v, p, ph, th, tv, lv, cpn, ep, sigp, clw, m, ment, elij, delt, &
426                   plcl, mp, qp, up, vp, wt, water, evap, b)! na->klev
427           else
428              ! iflag_con == 4
429              CALL cv_unsat(klon, ncum, klev, inb, t, q, qs, gz, u, v, p, ph, h, &
430                   lv, ep, sigp, clw, m, ment, elij, iflag, mp, qp, up, vp, wt, &
431                   water, evap)
432           endif
433    
434           ! YIELD
435           ! (tendencies, precipitation, variables of interface with other
436           ! processes, etc)
437    
438           if (iflag_con == 3) then
439              CALL cv3_yield(klon, ncum, klev, klev, icb, inb, delt, t, q, u, v, &
440                   gz, p, ph, h, hp, lv, cpn, th, ep, clw, m, tp, mp, qp, up, vp, &
441                   wt, water, evap, b, ment, qent, uent, vent, nent, elij, sig, &
442                   tv, tvp, iflag, precip, VPrecip, ft, fq, fu, fv, upwd, dnwd, &
443                   dnwd0, ma, mike, tls, tps, qcondc, wd)! na->klev
444           else
445              ! iflag_con == 4
446              CALL cv_yield(klon, ncum, klev, nk, icb, inb, delt, t, q, u, v, gz, &
447                   p, ph, h, hp, lv, cpn, ep, clw, frac, m, mp, qp, up, vp, wt, &
448                   water, evap, ment, qent, uent, vent, nent, elij, tv, tvp, &
449                   iflag, wd, qprime, tprime, precip, cbmf, ft, fq, fu, fv, Ma, &
450                   qcondc)
451         endif         endif
452    
453         if (iflag_con.eq.4) then         ! passive tracers
           CALL cv_compress( len, nloc, ncum, nd &  
                , iflag1, nk1, icb1 &  
                , cbmf1, plcl1, tnk1, qnk1, gznk1 &  
                , t1, q1, qs1, u1, v1, gz1 &  
                , h1, lv1, cpn1, p1, ph1, tv1, tp1, tvp1, clw1 &  
                , iflag, nk, icb &  
                , cbmf, plcl, tnk, qnk, gznk &  
                , t, q, qs, u, v, gz, h, lv, cpn, p, ph, tv, tp, tvp, clw  &  
                , dph )  
        endif  
   
        !-------------------------------------------------------------------  
        ! --- UNDILUTE (ADIABATIC) UPDRAFT / second part :  
        ! ---   FIND THE REST OF THE LIFTED PARCEL TEMPERATURES  
        ! ---   &  
        ! ---   COMPUTE THE PRECIPITATION EFFICIENCIES AND THE  
        ! ---   FRACTION OF PRECIPITATION FALLING OUTSIDE OF CLOUD  
        ! ---   &  
        ! ---   FIND THE LEVEL OF NEUTRAL BUOYANCY  
        !-------------------------------------------------------------------  
   
        if (iflag_con.eq.3) then  
           CALL cv3_undilute2(nloc, ncum, nd, icb, icbs, nk         &  
                , tnk, qnk, gznk, t, q, qs, gz &  
                , p, h, tv, lv, pbase, buoybase, plcl &  
                , inb, tp, tvp, clw, hp, ep, sigp, buoy) !na->nd  
        endif  
   
        if (iflag_con.eq.4) then  
           CALL cv_undilute2(nloc, ncum, nd, icb, nk &  
                , tnk, qnk, gznk, t, q, qs, gz &  
                , p, dph, h, tv, lv &  
                , inb, inbis, tp, tvp, clw, hp, ep, sigp, frac)  
        endif  
   
        !-------------------------------------------------------------------  
        ! --- CLOSURE  
        !-------------------------------------------------------------------  
   
        if (iflag_con.eq.3) then  
           CALL cv3_closure(nloc, ncum, nd, icb, inb               &  
                , pbase, p, ph, tv, buoy &  
                , sig, w0, cape, m) ! na->nd  
        endif  
   
        if (iflag_con.eq.4) then  
           CALL cv_closure(nloc, ncum, nd, nk, icb &  
                , tv, tvp, p, ph, dph, plcl, cpn &  
                , iflag, cbmf)  
        endif  
   
        !-------------------------------------------------------------------  
        ! --- MIXING  
        !-------------------------------------------------------------------  
   
        if (iflag_con.eq.3) then  
           CALL cv3_mixing(nloc, ncum, nd, nd, icb, nk, inb, ph, t, q, &  
                qs, u, v, h, lv, qnk, hp, tv, tvp, ep, clw, m, sig, ment, &  
                qent, uent, vent, nent, sij, elij, ments, qents)  
        endif  
   
        if (iflag_con.eq.4) then  
           CALL cv_mixing(nloc, ncum, nd, icb, nk, inb, inbis &  
                , ph, t, q, qs, u, v, h, lv, qnk &  
                , hp, tv, tvp, ep, clw, cbmf &  
                , m, ment, qent, uent, vent, nent, sij, elij)  
        endif  
   
        !-------------------------------------------------------------------  
        ! --- UNSATURATED (PRECIPITATING) DOWNDRAFTS  
        !-------------------------------------------------------------------  
   
        if (iflag_con.eq.3) then  
           CALL cv3_unsat(nloc, ncum, nd, nd, 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->nd  
        endif  
   
        if (iflag_con.eq.4) then  
           CALL cv_unsat(nloc, ncum, nd, inb, t, q, qs, gz, u, v, p, ph &  
                , h, lv, ep, sigp, clw, m, ment, elij &  
                , iflag, mp, qp, up, vp, wt, water, evap)  
        endif  
454    
455         !-------------------------------------------------------------------         if (iflag_con == 3) CALL cv3_tracer(klon, ncum, klev, ment, sij, da, phi)
        ! --- YIELD  
        !     (tendencies, precipitation, variables of interface with other  
        !      processes, etc)  
        !-------------------------------------------------------------------  
   
        if (iflag_con.eq.3) then  
           CALL cv3_yield(nloc, ncum, nd, nd             &  
                , 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->nd  
        endif  
456    
457         if (iflag_con.eq.4) then         ! UNCOMPRESS THE FIELDS
           CALL cv_yield(nloc, ncum, nd, nk, icb, inb, delt &  
                , t, q, u, v, gz, p, ph, h, hp, lv, cpn &  
                , ep, clw, frac, m, mp, qp, up, vp &  
                , wt, water, evap &  
                , ment, qent, uent, vent, nent, elij &  
                , tv, tvp &  
                , iflag, wd, qprime, tprime &  
                , precip, cbmf, ft, fq, fu, fv, Ma, qcondc)  
        endif  
458    
459         !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^         ! set iflag1 = 42 for non convective points
460         ! --- passive tracers         do i = 1, klon
        !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^  
   
        if (iflag_con.eq.3) then  
           CALL cv3_tracer(nloc, len, ncum, nd, nd, &  
                ment, sij, da, phi)  
        endif  
   
        !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^  
        ! --- UNCOMPRESS THE FIELDS  
        !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^  
        ! set iflag1  = 42 for non convective points  
        do  i = 1, len  
461            iflag1(i) = 42            iflag1(i) = 42
462         end do         end do
463    
464         if (iflag_con.eq.3) then         if (iflag_con == 3) then
465            CALL cv3_uncompress(nloc, len, ncum, nd, idcum, iflag, precip, &            CALL cv3_uncompress(idcum(:ncum), iflag, precip, VPrecip, sig, w0, &
466                 VPrecip, sig, w0, ft, fq, fu, fv, inb, Ma, upwd, dnwd, dnwd0, &                 ft, fq, fu, fv, inb, Ma, upwd, dnwd, dnwd0, qcondc, wd, cape, &
467                 qcondc, wd, cape, da, phi, mp, iflag1, precip1, VPrecip1, &                 da, phi, mp, iflag1, precip1, VPrecip1, sig1, w01, ft1, fq1, &
468                 sig1, w01, ft1, fq1, fu1, fv1, inb1, Ma1, upwd1, dnwd1, &                 fu1, fv1, inb1, Ma1, upwd1, dnwd1, dnwd01, qcondc1, wd1, &
469                 dnwd01, qcondc1, wd1, cape1, da1, phi1, mp1)                 cape1, da1, phi1, mp1)
470         endif         else
471              ! iflag_con == 4
472         if (iflag_con.eq.4) then            CALL cv_uncompress(idcum(:ncum), iflag, precip, cbmf, ft, fq, fu, &
473            CALL cv_uncompress(nloc, len, ncum, nd, idcum &                 fv, Ma, qcondc, iflag1, precip1, cbmf1, ft1, fq1, fu1, fv1, &
474                 , iflag &                 Ma1, qcondc1)
                , precip, cbmf &  
                , ft, fq, fu, fv &  
                , Ma, qcondc             &  
                , iflag1 &  
                , precip1, cbmf1 &  
                , ft1, fq1, fu1, fv1 &  
                , Ma1, qcondc1 )  
475         endif         endif
476      ENDIF ! ncum>0      ENDIF ! ncum>0
477    
Legend:
Removed from v.99  
changed lines
  Added in v.180
  ViewVC Help
Powered by ViewVC 1.1.21