/[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 189 by guez, Tue Mar 29 15:20:23 2016 UTC revision 198 by guez, Tue May 31 16:17:35 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, dnwd1, &         fv1, precip1, VPrecip1, sig1, w01, icb1, inb1, Ma1, upwd1, dnwd1, &
9         dnwd01, qcondc1, 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
# 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 23  contains Line 24  contains
24      use cv30_tracer_m, only: cv30_tracer      use cv30_tracer_m, only: cv30_tracer
25      use cv30_trigger_m, only: cv30_trigger      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 (K)      real, intent(in):: t1(klon, klev) ! temperature (K)
# Line 43  contains Line 46  contains
46      ! value of PH should be greater than (i.e. at a lower level than)      ! value of PH should be greater than (i.e. at a lower level than)
47      ! the first value of the array P1.      ! the first 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.
51    
52      ! 0: Moist convection occurs.      ! 0: Moist convection occurs.
# Line 56  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.
71    
72      ! 9: No moist convection: cloud base is higher then the level NL-1.      ! 9: No moist convection: cloud base is higher than the level NL-1.
73    
74      real, intent(out):: ft1(klon, klev) ! temperature tendency (K/s)      real, intent(out):: ft1(klon, klev) ! temperature tendency (K/s)
75      real, intent(out):: fq1(klon, klev) ! specific humidity tendency (s-1)      real, intent(out):: fq1(klon, klev) ! specific humidity tendency (s-1)
# Line 86  contains Line 89  contains
89    
90      integer, intent(out):: icb1(klon)      integer, intent(out):: icb1(klon)
91      integer, intent(inout):: inb1(klon)      integer, intent(inout):: inb1(klon)
     real, intent(in):: delt ! the model time step (sec) between calls  
   
92      real, intent(out):: Ma1(klon, klev) ! mass flux of 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)
# Line 107  contains Line 108  contains
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  
     integer nk1(klon)  
111      integer icbs1(klon)      integer icbs1(klon)
112      real plcl1(klon)      real plcl1(klon)
113      real tnk1(klon)      real tnk1(klon)
# Line 129  contains Line 128  contains
128      integer ncum      integer ncum
129    
130      ! Compressed fields:      ! Compressed fields:
131      integer idcum(klon)      integer, allocatable:: idcum(:), iflag(:) ! (ncum)
132      integer iflag(klon), nk(klon), icb(klon)      integer, allocatable:: icb(:) ! (ncum)
133      integer nent(klon, klev)      integer nent(klon, klev)
134      integer icbs(klon)      integer icbs(klon)
135      integer inb(klon)      integer inb(klon)
136      real plcl(klon), tnk(klon), qnk(klon), gznk(klon)      real, allocatable:: plcl(:) ! (ncum)
137        real tnk(klon), qnk(klon), gznk(klon)
138      real t(klon, klev), q(klon, klev), qs(klon, klev)      real t(klon, klev), q(klon, klev), qs(klon, klev)
139      real u(klon, klev), v(klon, klev)      real u(klon, klev), v(klon, klev)
140      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)
141      real p(klon, klev), ph(klon, klev + 1), tv(klon, klev), tp(klon, klev)      real p(klon, klev) ! pressure at full level, in hPa
142        real ph(klon, klev + 1), tv(klon, klev), tp(klon, klev)
143      real clw(klon, klev)      real clw(klon, klev)
144      real pbase(klon), buoybase(klon), th(klon, klev)      real pbase(klon), buoybase(klon), th(klon, klev)
145      real tvp(klon, klev)      real tvp(klon, klev)
146      real sig(klon, klev), w0(klon, klev)      real sig(klon, klev), w0(klon, klev)
147      real hp(klon, klev), ep(klon, klev), sigp(klon, klev)      real hp(klon, klev), ep(klon, klev)
148      real buoy(klon, klev)      real buoy(klon, klev)
149      real cape(klon)      real cape(klon)
150      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 151  contains Line 152  contains
152      real ments(klon, klev, klev), qents(klon, klev, klev)      real ments(klon, klev, klev), qents(klon, klev, klev)
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)
156        real, allocatable:: evap(:, :) ! (ncum, nl)
157      real, allocatable:: b(:, :) ! (ncum, nl - 1)      real, allocatable:: b(:, :) ! (ncum, nl - 1)
158      real ft(klon, klev), fq(klon, klev)      real ft(klon, klev), fq(klon, klev)
159      real fu(klon, klev), fv(klon, klev)      real fu(klon, klev), fv(klon, klev)
# Line 165  contains Line 167  contains
167      !-------------------------------------------------------------------      !-------------------------------------------------------------------
168    
169      ! SET CONSTANTS AND PARAMETERS      ! SET CONSTANTS AND PARAMETERS
   
     ! set thermodynamical constants:  
     ! (common cvthermo)  
170      CALL cv_thermo      CALL cv_thermo
171        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(delt)  
172    
173      ! INITIALIZE OUTPUT ARRAYS AND PARAMETERS      ! INITIALIZE OUTPUT ARRAYS AND PARAMETERS
174    
# Line 198  contains Line 192  contains
192         end do         end do
193      end do      end do
194    
195      do i = 1, klon      precip1 = 0.
196         precip1(i) = 0.      cape1 = 0.
197         iflag1(i) = 0      VPrecip1(:, klev + 1) = 0.
        cape1(i) = 0.  
        VPrecip1(i, klev + 1) = 0.  
     end do  
198    
199      do il = 1, klon      do il = 1, klon
200         sig1(il, klev) = sig1(il, klev) + 1.         sig1(il, klev) = sig1(il, klev) + 1.
201         sig1(il, klev) = min(sig1(il, klev), 12.1)         sig1(il, klev) = min(sig1(il, klev), 12.1)
202      enddo      enddo
203    
204      ! CALCULATE ARRAYS OF GEOPOTENTIAL, HEAT CAPACITY & STATIC ENERGY      CALL cv30_prelim(t1, q1, p1, ph1, lv1, cpn1, tv1, gz1, h1, hm1, th1)
205      CALL cv30_prelim(klon, klev, klev + 1, t1, q1, p1, ph1, lv1, cpn1, tv1, &      CALL cv30_feed(t1, q1, qs1, p1, ph1, gz1, icb1, iflag1, tnk1, qnk1, &
206           gz1, h1, hm1, th1)           gznk1, plcl1)
207        CALL cv30_undilute1(t1, q1, qs1, gz1, plcl1, p1, icb1, tp1, tvp1, clw1, &
208      ! CONVECTIVE FEED           icbs1)
209      CALL cv30_feed(klon, klev, t1, q1, qs1, p1, ph1, gz1, nk1, icb1, &      CALL cv30_trigger(icb1, plcl1, p1, th1, tv1, tvp1, pbase1, buoybase1, &
210           icbmax, iflag1, tnk1, qnk1, gznk1, plcl1) ! klev->na           iflag1, sig1, w01)
   
     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  
211    
212      IF (ncum > 0) THEN      ncum = count(iflag1 == 0)
        allocate(b(ncum, nl - 1))  
        CALL cv30_compress(ncum, iflag1, nk1, icb1, icbs1, plcl1, tnk1, qnk1, &  
             gznk1, pbase1, buoybase1, t1, q1, qs1, u1, v1, gz1, th1, h1, lv1, &  
             cpn1, p1, ph1, tv1, tp1, tvp1, clw1, sig1, w01, iflag, nk, icb, &  
             icbs, plcl, tnk, qnk, gznk, pbase, buoybase, t, q, qs, u, v, gz, &  
             th, h, lv, cpn, p, ph, tv, tp, tvp, clw, sig, w0)  
        CALL cv30_undilute2(ncum, icb, icbs, nk, tnk, qnk, gznk, t, qs, gz, p, &  
             h, tv, lv, pbase, buoybase, plcl, inb(:ncum), tp, tvp, clw, hp, &  
             ep, sigp, buoy)  
   
        ! CLOSURE  
        CALL cv30_closure(klon, ncum, klev, icb, inb, pbase, p, ph, tv, &  
             buoy, sig, w0, cape, m) ! na->klev  
   
        ! MIXING  
        CALL cv30_mixing(klon, ncum, klev, klev, icb, nk, inb, t, q, qs, u, &  
             v, h, lv, hp, ep, clw, m, sig, ment, qent, uent, vent, nent, &  
             sij, elij, ments, qents)  
   
        ! Unsaturated (precipitating) downdrafts  
        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(:ncum, :nl), up(:ncum, :nl), vp(:ncum, :nl), 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)  
213    
214        IF (ncum > 0) THEN
215           ! Moist convective adjustment is necessary
216           allocate(idcum(ncum), plcl(ncum))
217           allocate(b(ncum, nl - 1), evap(ncum, nl), icb(ncum), iflag(ncum))
218           idcum = pack((/(i, i = 1, klon)/), iflag1 == 0)
219           CALL cv30_compress(iflag1, icb1, icbs1, plcl1, tnk1, qnk1, gznk1, &
220                pbase1, buoybase1, t1, q1, qs1, u1, v1, gz1, th1, h1, lv1, cpn1, &
221                p1, ph1, tv1, tp1, tvp1, clw1, sig1, w01, icb, icbs, plcl, tnk, &
222                qnk, gznk, pbase, buoybase, t, q, qs, u, v, gz, th, h, lv, cpn, &
223                p, ph, tv, tp, tvp, clw, sig, w0)
224           CALL cv30_undilute2(icb, icbs(:ncum), tnk, qnk, gznk, t, qs, gz, p, h, &
225                tv, lv, pbase(:ncum), buoybase(:ncum), plcl, inb(:ncum), tp, tvp, &
226                clw, hp, ep, buoy)
227           CALL cv30_closure(icb, inb(:ncum), pbase, p, ph(:ncum, :), tv, buoy, &
228                sig, w0, cape, m)
229           CALL cv30_mixing(icb, inb(:ncum), t, q, qs, u, v, h, lv, hp, ep, clw, &
230                m, sig, ment, qent, uent, vent, nent, sij, elij, ments, qents)
231           CALL cv30_unsat(icb, inb(:ncum), t(:ncum, :nl), q(:ncum, :nl), &
232                qs(:ncum, :nl), gz, u(:ncum, :nl), v(:ncum, :nl), p, &
233                ph(:ncum, :), th(:ncum, :nl - 1), tv, lv(:ncum, :), &
234                cpn(:ncum, :nl), ep(:ncum, :), clw(:ncum, :), m(:ncum, :), &
235                ment(:ncum, :, :), elij(:ncum, :, :), dtphys, plcl, mp, &
236                qp(:ncum, :nl), up(:ncum, :nl), vp(:ncum, :nl), wt(:ncum, :nl), &
237                water(:ncum, :nl), evap, b)
238           CALL cv30_yield(icb, inb(:ncum), dtphys, t, q, u, v, gz, p, ph, h, hp, &
239                lv, cpn, th, ep, clw, m, tp, mp, qp, up, vp(:ncum, 2:nl), &
240                wt(:ncum, :nl - 1), water(:ncum, :nl), evap, b, ment, qent, uent, &
241                vent, nent, elij, sig, tv, tvp, iflag, precip, VPrecip, ft, fq, &
242                fu, fv, upwd, dnwd, dnwd0, ma, mike, tls, tps, qcondc)
243         CALL cv30_tracer(klon, ncum, klev, ment, sij, da, phi)         CALL cv30_tracer(klon, ncum, klev, ment, sij, da, phi)
244           CALL cv30_uncompress(idcum, iflag, precip, VPrecip, sig, w0, ft, fq, &
245         ! UNCOMPRESS THE FIELDS              fu, fv, inb, Ma, upwd, dnwd, dnwd0, qcondc, cape, da, phi, mp, &
246         iflag1 = 42 ! for non convective points              iflag1, precip1, VPrecip1, sig1, w01, ft1, fq1, fu1, fv1, inb1, &
247         CALL cv30_uncompress(idcum(:ncum), iflag, precip, VPrecip, sig, w0, &              Ma1, upwd1, dnwd1, dnwd01, qcondc1, cape1, da1, phi1, mp1)
             ft, fq, fu, fv, inb, Ma, upwd, dnwd, dnwd0, qcondc, cape, &  
             da, phi, mp, iflag1, precip1, VPrecip1, sig1, w01, ft1, fq1, &  
             fu1, fv1, inb1, Ma1, upwd1, dnwd1, dnwd01, qcondc1, cape1, da1, &  
             phi1, mp1)  
248      ENDIF      ENDIF
249    
250    end SUBROUTINE cv_driver    end SUBROUTINE cv_driver
Legend:
Removed from v.189  
changed lines
  Added in v.198
  ViewVC Help
Powered by ViewVC 1.1.21