/[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 187 by guez, Mon Mar 21 18:01:02 2016 UTC revision 192 by guez, Thu May 12 13:00:07 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, delt, 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 18  contains Line 18  contains
18      use cv30_compress_m, only: cv30_compress      use cv30_compress_m, only: cv30_compress
19      use cv30_feed_m, only: cv30_feed      use cv30_feed_m, only: cv30_feed
20      use cv30_mixing_m, only: cv30_mixing      use cv30_mixing_m, only: cv30_mixing
21      use cv30_param_m, only: cv30_param      use cv30_param_m, only: cv30_param, nl
22      use cv30_prelim_m, only: cv30_prelim      use cv30_prelim_m, only: cv30_prelim
23      use cv30_tracer_m, only: cv30_tracer      use cv30_tracer_m, only: cv30_tracer
24        use cv30_trigger_m, only: cv30_trigger
25      use cv30_uncompress_m, only: cv30_uncompress      use cv30_uncompress_m, only: cv30_uncompress
26      use cv30_undilute2_m, only: cv30_undilute2      use cv30_undilute2_m, only: cv30_undilute2
27      use cv30_unsat_m, only: cv30_unsat      use cv30_unsat_m, only: cv30_unsat
28      use cv30_yield_m, only: cv30_yield      use cv30_yield_m, only: cv30_yield
29        use cv_thermo_m, only: cv_thermo
30      USE dimphy, ONLY: klev, klon      USE dimphy, ONLY: klev, klon
31    
32      real, intent(in):: t1(klon, klev)      real, intent(in):: t1(klon, klev) ! temperature (K)
33      ! temperature (K), with first index corresponding to lowest model      real, intent(in):: q1(klon, klev) ! specific humidity
34      ! 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.  
35    
36      real, intent(in):: u1(klon, klev), v1(klon, klev)      real, intent(in):: u1(klon, klev), v1(klon, klev)
37      ! Zonal wind and meridional velocity (m/s), witth first index      ! zonal wind and meridional velocity (m/s)
38      ! corresponding with the lowest model level. Defined at same  
39      ! 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.  
40    
41      real, intent(in):: ph1(klon, klev + 1)      real, intent(in):: ph1(klon, klev + 1)
42      ! Half level pressure (mb), with first index corresponding to      ! Half level pressure (hPa). These pressures are defined at levels
43      ! lowest level. These pressures are defined at levels intermediate      ! intermediate between those of P1, T1, Q1 and QS1. The first
44      ! 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)
45      ! 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.  
46    
47      integer, intent(out):: iflag1(klon)      integer, intent(out):: iflag1(klon)
48      ! Flag for Emanuel conditions.      ! Flag for Emanuel conditions.
# Line 82  contains Line 69  contains
69    
70      ! 9: No moist convection: cloud base is higher then the level NL-1.      ! 9: No moist convection: cloud base is higher then the level NL-1.
71    
72      real, intent(out):: ft1(klon, klev)      real, intent(out):: ft1(klon, klev) ! temperature tendency (K/s)
73      ! Temperature tendency (K/s), defined at same grid levels as T1,      real, intent(out):: fq1(klon, klev) ! specific humidity tendency (s-1)
     ! Q1, QS1 and P1.  
   
     real, intent(out):: fq1(klon, klev)  
     ! Specific humidity tendencies (s-1), defined at same grid levels  
     ! as T1, Q1, QS1 and P1.  
74    
75      real, intent(out):: fu1(klon, klev), fv1(klon, klev)      real, intent(out):: fu1(klon, klev), fv1(klon, klev)
76      ! 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.  
77    
78      real, intent(out):: precip1(klon) ! convective precipitation rate (mm/day)      real, intent(out):: precip1(klon) ! convective precipitation rate (mm/day)
79    
80      real, intent(out):: VPrecip1(klon, klev + 1)      real, intent(out):: VPrecip1(klon, klev + 1)
81      ! vertical profile of convective precipitation (kg/m2/s)      ! vertical profile of convective precipitation (kg/m2/s)
82    
83      real, intent(inout):: sig1(klon, klev) ! section adiabatic updraft      real, intent(inout):: sig1(klon, klev) ! section of adiabatic updraft
84    
85      real, intent(inout):: w01(klon, klev)      real, intent(inout):: w01(klon, klev)
86      ! vertical velocity within adiabatic updraft      ! vertical velocity within adiabatic updraft
# Line 108  contains Line 89  contains
89      integer, intent(inout):: inb1(klon)      integer, intent(inout):: inb1(klon)
90      real, intent(in):: delt ! the model time step (sec) between calls      real, intent(in):: delt ! the model time step (sec) between calls
91    
92      real Ma1(klon, klev) ! Output mass flux adiabatic updraft      real, intent(out):: Ma1(klon, klev) ! mass flux of 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
     real wd1(klon) ! gust  
     ! Output downdraft velocity scale for surface fluxes  
     ! A convective downdraft velocity scale. For use in surface  
     ! flux parameterizations. See convect.ps file for details.  
102    
103      real cape1(klon) ! Output      real, intent(out):: cape1(klon)
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), icb(klon)
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 plcl(klon), tnk(klon), qnk(klon), gznk(klon)      real plcl(klon), tnk(klon), qnk(klon), gznk(klon)
139      real t(klon, klev), q(klon, klev), qs(klon, klev)      real t(klon, klev), q(klon, klev), qs(klon, klev)
140      real u(klon, klev), v(klon, klev)      real u(klon, klev), v(klon, klev)
# Line 182  contains Line 153  contains
153      real sij(klon, klev, klev), elij(klon, klev, klev)      real sij(klon, klev, klev), elij(klon, klev, klev)
154      real qp(klon, klev), up(klon, klev), vp(klon, klev)      real qp(klon, klev), up(klon, klev), vp(klon, klev)
155      real wt(klon, klev), water(klon, klev), evap(klon, klev)      real wt(klon, klev), water(klon, klev), evap(klon, klev)
156      real b(klon, klev), ft(klon, klev), fq(klon, klev)      real, allocatable:: b(:, :) ! (ncum, nl - 1)
157        real ft(klon, klev), fq(klon, klev)
158      real fu(klon, klev), fv(klon, klev)      real fu(klon, klev), fv(klon, klev)
159      real upwd(klon, klev), dnwd(klon, klev), dnwd0(klon, klev)      real upwd(klon, klev), dnwd(klon, klev), dnwd0(klon, klev)
160      real Ma(klon, klev), mike(klon, klev), tls(klon, klev)      real Ma(klon, klev), mike(klon, klev), tls(klon, klev)
# Line 190  contains Line 162  contains
162      real precip(klon)      real precip(klon)
163      real VPrecip(klon, klev + 1)      real VPrecip(klon, klev + 1)
164      real qcondc(klon, klev) ! cld      real qcondc(klon, klev) ! cld
     real wd(klon) ! gust  
165    
166      !-------------------------------------------------------------------      !-------------------------------------------------------------------
167    
168      ! SET CONSTANTS AND PARAMETERS      ! SET CONSTANTS AND PARAMETERS
169    
170      ! set thermodynamical constants:      ! set thermodynamical constants:
     ! (common cvthermo)  
171      CALL cv_thermo      CALL cv_thermo
172    
173      ! set convect parameters      ! set convect parameters
# Line 231  contains Line 201  contains
201      do i = 1, klon      do i = 1, klon
202         precip1(i) = 0.         precip1(i) = 0.
203         iflag1(i) = 0         iflag1(i) = 0
        wd1(i) = 0.  
204         cape1(i) = 0.         cape1(i) = 0.
205         VPrecip1(i, klev + 1) = 0.         VPrecip1(i, klev + 1) = 0.
206      end do      end do
# Line 249  contains Line 218  contains
218      CALL cv30_feed(klon, klev, t1, q1, qs1, p1, ph1, gz1, nk1, icb1, &      CALL cv30_feed(klon, klev, t1, q1, qs1, p1, ph1, gz1, nk1, icb1, &
219           icbmax, iflag1, tnk1, qnk1, gznk1, plcl1) ! klev->na           icbmax, iflag1, tnk1, qnk1, gznk1, plcl1) ! klev->na
220    
     ! 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.  
221      CALL cv30_undilute1(klon, klev, t1, q1, qs1, gz1, plcl1, p1, nk1, icb1, &      CALL cv30_undilute1(klon, klev, t1, q1, qs1, gz1, plcl1, p1, nk1, icb1, &
222           tp1, tvp1, clw1, icbs1) ! klev->na           tp1, tvp1, clw1, icbs1) ! klev->na
223    
# Line 271  contains Line 236  contains
236      end do      end do
237    
238      IF (ncum > 0) THEN      IF (ncum > 0) THEN
239         ! COMPRESS THE FIELDS         allocate(b(ncum, nl - 1))
240         ! (-> vectorization over convective gridpoints)         CALL cv30_compress(ncum, iflag1, nk1, icb1, icbs1, plcl1, tnk1, qnk1, &
241         CALL cv30_compress(klon, klon, ncum, klev, iflag1, nk1, icb1, icbs1, &              gznk1, pbase1, buoybase1, t1, q1, qs1, u1, v1, gz1, th1, h1, lv1, &
242              plcl1, tnk1, qnk1, gznk1, pbase1, buoybase1, t1, q1, qs1, u1, &              cpn1, p1, ph1, tv1, tp1, tvp1, clw1, sig1, w01, iflag, nk, icb, &
243              v1, gz1, th1, h1, lv1, cpn1, p1, ph1, tv1, tp1, tvp1, clw1, &              icbs, plcl, tnk, qnk, gznk, pbase, buoybase, t, q, qs, u, v, gz, &
244              sig1, w01, iflag, nk, icb, icbs, plcl, tnk, qnk, gznk, pbase, &              th, h, lv, cpn, p, ph, tv, tp, tvp, clw, sig, w0)
             buoybase, t, q, qs, u, v, gz, th, h, lv, cpn, p, ph, tv, tp, &  
             tvp, clw, sig, w0)  
   
245         CALL cv30_undilute2(ncum, icb, icbs, nk, tnk, qnk, gznk, t, qs, gz, p, &         CALL cv30_undilute2(ncum, icb, icbs, nk, tnk, qnk, gznk, t, qs, gz, p, &
246              h, tv, lv, pbase, buoybase, plcl, inb(:ncum), tp, tvp, clw, hp, &              h, tv, lv, pbase, buoybase, plcl, inb(:ncum), tp, tvp, clw, hp, &
247              ep, sigp, buoy)              ep, sigp, buoy)
# Line 294  contains Line 256  contains
256              sij, elij, ments, qents)              sij, elij, ments, qents)
257    
258         ! Unsaturated (precipitating) downdrafts         ! Unsaturated (precipitating) downdrafts
259         CALL cv30_unsat(ncum, icb(:ncum), inb(:ncum), t, q, qs, gz, u, v, p, &         CALL cv30_unsat(icb(:ncum), inb(:ncum), t, q, qs, gz, u, v, p, ph, th, &
260              ph, th, tv, lv, cpn, ep, sigp, clw, m, ment, elij, delt, plcl, &              tv, lv, cpn, ep(:ncum, :), sigp(:ncum, :), clw(:ncum, :), &
261              mp, qp, up, vp, wt, water, evap, b(:ncum, :))              m(:ncum, :), ment(:ncum, :, :), elij(:ncum, :, :), delt, plcl, &
262                mp, qp(:ncum, :nl), up(:ncum, :nl), vp(:ncum, :nl), wt, water, &
263                evap, b)
264    
265         ! Yield (tendencies, precipitation, variables of interface with         ! Yield (tendencies, precipitation, variables of interface with
266         ! other processes, etc)         ! other processes, etc)
267         CALL cv30_yield(klon, ncum, klev, klev, icb, inb, delt, t, q, u, v, &         CALL cv30_yield(icb(:ncum), inb(:ncum), delt, t, q, u, v, gz, p, ph, &
268              gz, p, ph, h, hp, lv, cpn, th, ep, clw, m, tp, mp, qp, up, vp, &              h, hp, lv, cpn, th, ep, clw, m, tp, mp, qp, up, vp, wt, &
269              wt, water, evap, b, ment, qent, uent, vent, nent, elij, sig, &              water(:ncum, :nl), evap(:ncum, :nl), b, ment, qent, uent, vent, &
270              tv, tvp, iflag, precip, VPrecip, ft, fq, fu, fv, upwd, dnwd, &              nent, elij, sig, tv, tvp, iflag, precip, VPrecip, ft, fq, fu, fv, &
271              dnwd0, ma, mike, tls, tps, qcondc, wd)! na->klev              upwd, dnwd, dnwd0, ma, mike, tls, tps, qcondc)
272    
        ! passive tracers  
273         CALL cv30_tracer(klon, ncum, klev, ment, sij, da, phi)         CALL cv30_tracer(klon, ncum, klev, ment, sij, da, phi)
274    
275         ! UNCOMPRESS THE FIELDS         ! UNCOMPRESS THE FIELDS
276           iflag1 = 42 ! for non convective points
        ! set iflag1 = 42 for non convective points  
        iflag1 = 42  
   
277         CALL cv30_uncompress(idcum(:ncum), iflag, precip, VPrecip, sig, w0, &         CALL cv30_uncompress(idcum(:ncum), iflag, precip, VPrecip, sig, w0, &
278              ft, fq, fu, fv, inb, Ma, upwd, dnwd, dnwd0, qcondc, wd, cape, &              ft, fq, fu, fv, inb, Ma, upwd, dnwd, dnwd0, qcondc, cape, &
279              da, phi, mp, iflag1, precip1, VPrecip1, sig1, w01, ft1, fq1, &              da, phi, mp, iflag1, precip1, VPrecip1, sig1, w01, ft1, fq1, &
280              fu1, fv1, inb1, Ma1, upwd1, dnwd1, dnwd01, qcondc1, wd1, &              fu1, fv1, inb1, Ma1, upwd1, dnwd1, dnwd01, qcondc1, cape1, da1, &
281              cape1, da1, phi1, mp1)              phi1, mp1)
282      ENDIF      ENDIF
283    
284    end SUBROUTINE cv_driver    end SUBROUTINE cv_driver
Legend:
Removed from v.187  
changed lines
  Added in v.192
  ViewVC Help
Powered by ViewVC 1.1.21