/[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 201 by guez, Mon Jun 6 17:42:15 2016 UTC revision 205 by guez, Tue Jun 21 15:16:03 2016 UTC
# Line 6  contains Line 6  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, Ma1, upwd1, dnwd1, &         fv1, precip1, VPrecip1, sig1, w01, icb1, inb1, Ma1, upwd1, dnwd1, &
9         dnwd01, qcondc1, cape1, da1, phi1, mp1)         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 94  contains Line 94  contains
94      ! total upward mass flux (adiabatic + mixed)      ! total upward mass flux (adiabatic + mixed)
95    
96      real, intent(out):: dnwd1(klon, klev) ! saturated downward mass flux (mixed)      real, intent(out):: dnwd1(klon, klev) ! saturated downward mass flux (mixed)
     real, intent(out):: dnwd01(klon, klev) ! unsaturated downward mass flux  
97    
98      real, intent(out):: qcondc1(klon, klev)      real, intent(out):: qcondc1(klon, klev)
99      ! in-cloud mixing ratio of condensed water      ! in-cloud mixing ratio of condensed water
100    
101      real, intent(out):: cape1(klon)      real, intent(out):: cape1(klon)
102      real, intent(inout):: da1(klon, klev), phi1(klon, klev, klev)      real, intent(out):: da1(:, :) ! (klon, klev)
103      real, intent(inout):: mp1(klon, klev)      real, intent(out):: phi1(:, :, :) ! (klon, klev, klev)
104    
105        real, intent(out):: mp1(:, :) ! (klon, klev) Mass flux of the
106        ! unsaturated downdraft, defined positive downward, in kg m-2
107        ! s-1. M_p in Emanuel (1991 928).
108    
109      ! Local:      ! Local:
110    
111      real da(klon, klev), phi(klon, klev, klev)      real da(klon, klev), phi(klon, klev, klev)
112      real, allocatable:: mp(:, :) ! (ncum, nl)  
113        real, allocatable:: mp(:, :) ! (ncum, nl) Mass flux of the
114        ! unsaturated downdraft, defined positive downward, in kg m-2
115        ! s-1. M_p in Emanuel (1991 928).
116    
117      integer i, k, il      integer i, k, il
118      integer icbs1(klon)      integer icbs1(klon)
119      real plcl1(klon)      real plcl1(klon)
# Line 174  contains Line 181  contains
181      real, allocatable:: b(:, :) ! (ncum, nl - 1)      real, allocatable:: b(:, :) ! (ncum, nl - 1)
182      real ft(klon, klev), fq(klon, klev)      real ft(klon, klev), fq(klon, klev)
183      real fu(klon, klev), fv(klon, klev)      real fu(klon, klev), fv(klon, klev)
184      real upwd(klon, klev), dnwd(klon, klev), dnwd0(klon, klev)      real upwd(klon, klev), dnwd(klon, klev)
185      real Ma(klon, klev), mike(klon, klev), tls(klon, klev)      real Ma(klon, klev), mike(klon, klev), tls(klon, klev)
186      real tps(klon, klev)      real tps(klon, klev)
187      real precip(klon)      real precip(klon)
# Line 188  contains Line 195  contains
195    
196      ! INITIALIZE OUTPUT ARRAYS AND PARAMETERS      ! INITIALIZE OUTPUT ARRAYS AND PARAMETERS
197    
198        da1 = 0.
199        mp1 = 0.
200        phi1 = 0.
201    
202      do k = 1, klev      do k = 1, klev
203         do i = 1, klon         do i = 1, klon
204            ft1(i, k) = 0.            ft1(i, k) = 0.
# Line 203  contains Line 214  contains
214            Ma1(i, k) = 0.            Ma1(i, k) = 0.
215            upwd1(i, k) = 0.            upwd1(i, k) = 0.
216            dnwd1(i, k) = 0.            dnwd1(i, k) = 0.
           dnwd01(i, k) = 0.  
217            qcondc1(i, k) = 0.            qcondc1(i, k) = 0.
218         end do         end do
219      end do      end do
# Line 256  contains Line 266  contains
266              lv, cpn, th, ep, clw, m, tp, mp, qp, up, vp(:ncum, 2:nl), &              lv, cpn, th, ep, clw, m, tp, mp, qp, up, vp(:ncum, 2:nl), &
267              wt(:ncum, :nl - 1), water(:ncum, :nl), evap, b, ment, qent, uent, &              wt(:ncum, :nl - 1), water(:ncum, :nl), evap, b, ment, qent, uent, &
268              vent, nent, elij, sig, tv, tvp, iflag, precip, VPrecip, ft, fq, &              vent, nent, elij, sig, tv, tvp, iflag, precip, VPrecip, ft, fq, &
269              fu, fv, upwd, dnwd, dnwd0, ma, mike, tls, tps, qcondc)              fu, fv, upwd, dnwd, ma, mike, tls, tps, qcondc)
270         CALL cv30_tracer(klon, ncum, klev, ment, sij, da, phi)         CALL cv30_tracer(klon, ncum, klev, ment, sij, da, phi)
271         CALL cv30_uncompress(idcum, iflag, precip, VPrecip, sig, w0, ft, fq, &         CALL cv30_uncompress(idcum, iflag, precip, VPrecip, sig, w0, ft, fq, &
272              fu, fv, inb, Ma, upwd, dnwd, dnwd0, qcondc, cape, da, phi, mp, &              fu, fv, inb, Ma, upwd, dnwd, qcondc, cape, da, phi, mp, iflag1, &
273              iflag1, precip1, VPrecip1, sig1, w01, ft1, fq1, fu1, fv1, inb1, &              precip1, VPrecip1, sig1, w01, ft1, fq1, fu1, fv1, inb1, Ma1, &
274              Ma1, upwd1, dnwd1, dnwd01, qcondc1, cape1, da1, phi1, mp1)              upwd1, dnwd1, qcondc1, cape1, da1, phi1, mp1)
275      ENDIF      ENDIF
276    
277    end SUBROUTINE cv_driver    end SUBROUTINE cv_driver
Legend:
Removed from v.201  
changed lines
  Added in v.205
  ViewVC Help
Powered by ViewVC 1.1.21