/[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

revision 188 by guez, Tue Mar 22 16:31:39 2016 UTC revision 195 by guez, Wed May 18 17:56:44 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      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
# Line 21  contains Line 22  contains
22      use cv30_param_m, only: cv30_param, nl      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)      real, intent(in):: t1(klon, klev) ! temperature (K)
35      ! temperature (K), with first index corresponding to lowest model      real, intent(in):: q1(klon, klev) ! specific humidity
36      ! level      real, intent(in):: qs1(klon, klev) ! saturation specific humidity
   
     real, intent(in):: q1(klon, klev)  
     ! Specific humidity, with first index corresponding to lowest  
     ! model level. Must be defined at same grid levels as T1.  
   
     real, intent(in):: qs1(klon, klev)  
     ! Saturation specific humidity, with first index corresponding to  
     ! lowest model level. Must be defined at same grid levels as  
     ! T1.  
37    
38      real, intent(in):: u1(klon, klev), v1(klon, klev)      real, intent(in):: u1(klon, klev), v1(klon, klev)
39      ! Zonal wind and meridional velocity (m/s), witth first index      ! zonal wind and meridional velocity (m/s)
40      ! corresponding with the lowest model level. Defined at same  
41      ! levels as T1.      real, intent(in):: p1(klon, klev) ! full level pressure (hPa)
   
     real, intent(in):: p1(klon, klev)  
     ! Full level pressure (mb) of dimension KLEV, with first index  
     ! corresponding to lowest model level. Must be defined at same  
     ! grid levels as T1.  
42    
43      real, intent(in):: ph1(klon, klev + 1)      real, intent(in):: ph1(klon, klev + 1)
44      ! Half level pressure (mb), with first index corresponding to      ! Half level pressure (hPa). These pressures are defined at levels
45      ! lowest level. These pressures are defined at levels intermediate      ! intermediate between those of P1, T1, Q1 and QS1. The first
46      ! between those of P1, T1, Q1 and QS1. The first value of PH      ! value of PH should be greater than (i.e. at a lower level than)
47      ! should be greater than (i.e. at a lower level than) the first      ! the first value of the array P1.
     ! value of the array P1.  
48    
49      integer, intent(out):: iflag1(klon)      integer, intent(out):: iflag1(klon)
50      ! Flag for Emanuel conditions.      ! Flag for Emanuel conditions.
# Line 70  contains Line 59  contains
59    
60      ! 3: No moist convection because new cbmf is 0 and old cbmf is 0.      ! 3: No moist convection because new cbmf is 0 and old cbmf is 0.
61    
62      ! 4: No moist convection; atmosphere is not unstable      ! 4: No moist convection; atmosphere is not unstable.
63    
64      ! 6: No moist convection because ihmin le minorig.      ! 6: No moist convection because ihmin <= minorig.
65    
66      ! 7: No moist convection because unreasonable parcel level      ! 7: No moist convection because unreasonable parcel level
67      ! temperature or specific humidity.      ! temperature or specific humidity.
68    
69      ! 8: No moist convection: lifted condensation level is above the      ! 8: No moist convection: lifted condensation level is above the
70      ! 200 mb level.      ! 200 mbar level.
   
     ! 9: No moist convection: cloud base is higher then the level NL-1.  
71    
72      real, intent(out):: ft1(klon, klev)      ! 9: No moist convection: cloud base is higher than the level NL-1.
     ! Temperature tendency (K/s), defined at same grid levels as T1,  
     ! Q1, QS1 and P1.  
73    
74      real, intent(out):: fq1(klon, klev)      real, intent(out):: ft1(klon, klev) ! temperature tendency (K/s)
75      ! Specific humidity tendencies (s-1), defined at same grid levels      real, intent(out):: fq1(klon, klev) ! specific humidity tendency (s-1)
     ! as T1, Q1, QS1 and P1.  
76    
77      real, intent(out):: fu1(klon, klev), fv1(klon, klev)      real, intent(out):: fu1(klon, klev), fv1(klon, klev)
78      ! Forcing (tendency) of zonal and meridional velocity (m/s^2),      ! forcing (tendency) of zonal and meridional velocity (m/s^2)
     ! defined at same grid levels as T1.  
79    
80      real, intent(out):: precip1(klon) ! convective precipitation rate (mm/day)      real, intent(out):: precip1(klon) ! convective precipitation rate (mm/day)
81    
82      real, intent(out):: VPrecip1(klon, klev + 1)      real, intent(out):: VPrecip1(klon, klev + 1)
83      ! vertical profile of convective precipitation (kg/m2/s)      ! vertical profile of convective precipitation (kg/m2/s)
84    
85      real, intent(inout):: sig1(klon, klev) ! section adiabatic updraft      real, intent(inout):: sig1(klon, klev) ! section of adiabatic updraft
86    
87      real, intent(inout):: w01(klon, klev)      real, intent(inout):: w01(klon, klev)
88      ! vertical velocity within adiabatic updraft      ! vertical velocity within adiabatic updraft
89    
90      integer, intent(out):: icb1(klon)      integer, intent(out):: icb1(klon)
91      integer, intent(inout):: inb1(klon)      integer, intent(inout):: inb1(klon)
92      real, intent(in):: delt ! the model time step (sec) between calls      real, intent(out):: Ma1(klon, klev) ! mass flux of adiabatic updraft
   
     real Ma1(klon, klev) ! Output mass flux adiabatic updraft  
93    
94      real, intent(out):: upwd1(klon, klev)      real, intent(out):: upwd1(klon, klev)
95      ! total upward mass flux (adiab + mixed)      ! total upward mass flux (adiabatic + mixed)
96    
97      real, intent(out):: dnwd1(klon, klev) ! saturated downward mass flux (mixed)      real, intent(out):: dnwd1(klon, klev) ! saturated downward mass flux (mixed)
98      real, intent(out):: dnwd01(klon, klev) ! unsaturated downward mass flux      real, intent(out):: dnwd01(klon, klev) ! unsaturated downward mass flux
99    
100      real qcondc1(klon, klev) ! Output in-cld mixing ratio of condensed water      real, intent(out):: qcondc1(klon, klev)
101        ! in-cloud mixing ratio of condensed water
102    
103      real wd1(klon) ! gust      real, intent(out):: cape1(klon)
     ! Output downdraft velocity scale for surface fluxes  
     ! A convective downdraft velocity scale. For use in surface  
     ! flux parameterizations. See convect.ps file for details.  
   
     real cape1(klon) ! Output  
104      real, intent(inout):: da1(klon, klev), phi1(klon, klev, klev)      real, intent(inout):: da1(klon, klev), phi1(klon, klev, klev)
105      real, intent(inout):: mp1(klon, klev)      real, intent(inout):: mp1(klon, klev)
106    
107      ! Local:      ! 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
111      integer icbmax      integer icbmax
112      integer nk1(klon)      integer nk1(klon)
113      integer icbs1(klon)      integer icbs1(klon)
   
114      real plcl1(klon)      real plcl1(klon)
115      real tnk1(klon)      real tnk1(klon)
116      real qnk1(klon)      real qnk1(klon)
117      real gznk1(klon)      real gznk1(klon)
118      real pbase1(klon)      real pbase1(klon)
119      real buoybase1(klon)      real buoybase1(klon)
   
120      real lv1(klon, klev)      real lv1(klon, klev)
121      real cpn1(klon, klev)      real cpn1(klon, klev)
122      real tv1(klon, klev)      real tv1(klon, klev)
# Line 153  contains Line 127  contains
127      real tvp1(klon, klev)      real tvp1(klon, klev)
128      real clw1(klon, klev)      real clw1(klon, klev)
129      real th1(klon, klev)      real th1(klon, klev)
   
130      integer ncum      integer ncum
131    
132      ! Compressed fields:      ! Compressed fields:
   
133      integer idcum(klon)      integer idcum(klon)
134      integer iflag(klon), nk(klon), icb(klon)      integer iflag(klon), nk(klon)
135        integer, allocatable:: icb(:) ! (ncum)
136      integer nent(klon, klev)      integer nent(klon, klev)
137      integer icbs(klon)      integer icbs(klon)
138      integer inb(klon)      integer inb(klon)
   
139      real plcl(klon), tnk(klon), qnk(klon), gznk(klon)      real plcl(klon), 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 181  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, allocatable:: b(:, :) ! (ncum, nl)      real, allocatable:: evap(:, :) ! (ncum, nl)
159        real, allocatable:: b(:, :) ! (ncum, nl - 1)
160      real ft(klon, klev), fq(klon, klev)      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)
# Line 191  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
172    
173      ! set thermodynamical constants:      ! set thermodynamical constants:
     ! (common cvthermo)  
174      CALL cv_thermo      CALL cv_thermo
175    
176      ! set convect parameters      CALL cv30_param
     ! includes microphysical parameters and parameters that  
     ! control the rate of approach to quasi-equilibrium)  
     ! (common cvparam)  
     CALL cv30_param(delt)  
177    
178      ! INITIALIZE OUTPUT ARRAYS AND PARAMETERS      ! INITIALIZE OUTPUT ARRAYS AND PARAMETERS
179    
# Line 229  contains Line 197  contains
197         end do         end do
198      end do      end do
199    
200      do i = 1, klon      precip1 = 0.
201         precip1(i) = 0.      iflag1 = 0
202         iflag1(i) = 0      cape1 = 0.
203         wd1(i) = 0.      VPrecip1(:, klev + 1) = 0.
        cape1(i) = 0.  
        VPrecip1(i, klev + 1) = 0.  
     end do  
204    
205      do il = 1, klon      do il = 1, klon
206         sig1(il, klev) = sig1(il, klev) + 1.         sig1(il, klev) = sig1(il, klev) + 1.
207         sig1(il, klev) = min(sig1(il, klev), 12.1)         sig1(il, klev) = min(sig1(il, klev), 12.1)
208      enddo      enddo
209    
     ! CALCULATE ARRAYS OF GEOPOTENTIAL, HEAT CAPACITY & STATIC ENERGY  
210      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, &
211           gz1, h1, hm1, th1)           gz1, h1, hm1, th1)
212        CALL cv30_feed(t1, q1, qs1, p1, ph1, gz1, nk1, icb1, icbmax, iflag1, &
213      ! CONVECTIVE FEED           tnk1, qnk1, gznk1, plcl1)
214      CALL cv30_feed(klon, klev, t1, q1, qs1, p1, ph1, gz1, nk1, icb1, &      CALL cv30_undilute1(t1, q1, qs1, gz1, plcl1, p1, nk1, icb1, tp1, tvp1, &
215           icbmax, iflag1, tnk1, qnk1, gznk1, plcl1) ! klev->na           clw1, icbs1)
216        CALL cv30_trigger(icb1, plcl1, p1, th1, tv1, tvp1, pbase1, buoybase1, &
217      ! UNDILUTE (ADIABATIC) UPDRAFT / 1st part           iflag1, sig1, w01)
     ! (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  
218    
219      ! Moist convective adjustment is necessary      ! Moist convective adjustment is necessary
220    
# Line 272  contains Line 227  contains
227      end do      end do
228    
229      IF (ncum > 0) THEN      IF (ncum > 0) THEN
230         allocate(b(ncum, nl))         allocate(b(ncum, nl - 1), evap(ncum, nl), icb(ncum))
231           CALL cv30_compress(ncum, iflag1, nk1, icb1, icbs1, plcl1, tnk1, qnk1, &
232         ! COMPRESS THE FIELDS              gznk1, pbase1, buoybase1, t1, q1, qs1, u1, v1, gz1, th1, h1, lv1, &
233         ! (-> vectorization over convective gridpoints)              cpn1, p1, ph1, tv1, tp1, tvp1, clw1, sig1, w01, iflag, nk, icb, &
234         CALL cv30_compress(klon, klon, ncum, klev, iflag1, nk1, icb1, icbs1, &              icbs, plcl, tnk, qnk, gznk, pbase, buoybase, t, q, qs, u, v, gz, &
235              plcl1, tnk1, qnk1, gznk1, pbase1, buoybase1, t1, q1, qs1, u1, &              th, h, lv, cpn, p, ph, tv, tp, tvp, clw, sig, w0)
236              v1, gz1, th1, h1, lv1, cpn1, p1, ph1, tv1, tp1, tvp1, clw1, &         CALL cv30_undilute2(icb, icbs(:ncum), nk, tnk, qnk, gznk, t, qs, gz, &
237              sig1, w01, iflag, nk, icb, icbs, plcl, tnk, qnk, gznk, pbase, &              p, h, tv, lv, pbase, buoybase, plcl, inb(:ncum), tp, tvp, clw, &
238              buoybase, t, q, qs, u, v, gz, th, h, lv, cpn, p, ph, tv, tp, &              hp, ep, buoy)
239              tvp, clw, sig, w0)         CALL cv30_closure(icb, inb(:ncum), pbase, p, ph, tv, buoy, sig, w0, &
240                cape, m)
241         CALL cv30_undilute2(ncum, icb, icbs, nk, tnk, qnk, gznk, t, qs, gz, p, &         CALL cv30_mixing(icb, nk(:ncum), inb(:ncum), t, q, qs, u, v, h, lv, &
242              h, tv, lv, pbase, buoybase, plcl, inb(:ncum), tp, tvp, clw, hp, &              hp, ep, clw, m, sig, ment, qent, uent, vent, nent, sij, elij, &
243              ep, sigp, buoy)              ments, qents)
244           CALL cv30_unsat(icb, inb(:ncum), t(:ncum, :nl), q(:ncum, :nl), &
245         ! CLOSURE              qs(:ncum, :nl), gz, u, v, p, ph, th(:ncum, :nl - 1), tv, lv, cpn, &
246         CALL cv30_closure(klon, ncum, klev, icb, inb, pbase, p, ph, tv, &              ep(:ncum, :), clw(:ncum, :), m(:ncum, :), ment(:ncum, :, :), &
247              buoy, sig, w0, cape, m) ! na->klev              elij(:ncum, :, :), dtphys, plcl, mp, qp(:ncum, :nl), &
248                up(:ncum, :nl), vp(:ncum, :nl), wt(:ncum, :nl), &
249         ! MIXING              water(:ncum, :nl), evap, b)
250         CALL cv30_mixing(klon, ncum, klev, klev, icb, nk, inb, t, q, qs, u, &         CALL cv30_yield(icb, inb(:ncum), dtphys, t, q, u, v, gz, p, ph, h, hp, &
251              v, h, lv, hp, ep, clw, m, sig, ment, qent, uent, vent, nent, &              lv, cpn, th, ep, clw, m, tp, mp, qp, up, vp(:ncum, 2:nl), &
252              sij, elij, ments, qents)              wt(:ncum, :nl - 1), water(:ncum, :nl), evap, b, ment, qent, uent, &
253                vent, nent, elij, sig, tv, tvp, iflag, precip, VPrecip, ft, fq, &
254         ! Unsaturated (precipitating) downdrafts              fu, fv, upwd, dnwd, dnwd0, ma, mike, tls, tps, qcondc)
        CALL cv30_unsat(icb(:ncum), inb(:ncum), 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)  
   
        ! Yield (tendencies, precipitation, variables of interface with  
        ! other processes, etc)  
        CALL cv30_yield(icb(:ncum), inb(:ncum), delt, t, q, u, v, gz, p, ph, &  
             h, hp, lv, cpn, th, ep, clw, m, tp, mp, qp, up, vp, wt, &  
             water(:ncum, :nl), evap(:ncum, :nl), b, ment, qent, uent, vent, &  
             nent, elij, sig, tv, tvp, iflag, precip, VPrecip, ft, fq, fu, fv, &  
             upwd, dnwd, dnwd0, ma, mike, tls, tps, qcondc, wd)  
   
        ! passive tracers  
255         CALL cv30_tracer(klon, ncum, klev, ment, sij, da, phi)         CALL cv30_tracer(klon, ncum, klev, ment, sij, da, phi)
256    
257         ! UNCOMPRESS THE FIELDS         ! UNCOMPRESS THE FIELDS
258           iflag1 = 42 ! for non convective points
        ! set iflag1 = 42 for non convective points  
        iflag1 = 42  
   
259         CALL cv30_uncompress(idcum(:ncum), iflag, precip, VPrecip, sig, w0, &         CALL cv30_uncompress(idcum(:ncum), iflag, precip, VPrecip, sig, w0, &
260              ft, fq, fu, fv, inb, Ma, upwd, dnwd, dnwd0, qcondc, wd, cape, &              ft, fq, fu, fv, inb, Ma, upwd, dnwd, dnwd0, qcondc, cape, &
261              da, phi, mp, iflag1, precip1, VPrecip1, sig1, w01, ft1, fq1, &              da, phi, mp, iflag1, precip1, VPrecip1, sig1, w01, ft1, fq1, &
262              fu1, fv1, inb1, Ma1, upwd1, dnwd1, dnwd01, qcondc1, wd1, &              fu1, fv1, inb1, Ma1, upwd1, dnwd1, dnwd01, qcondc1, cape1, da1, &
263              cape1, da1, phi1, mp1)              phi1, mp1)
264      ENDIF      ENDIF
265    
266    end SUBROUTINE cv_driver    end SUBROUTINE cv_driver
Legend:
Removed from v.188  
changed lines
  Added in v.195
  ViewVC Help
Powered by ViewVC 1.1.21