/[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 192 by guez, Thu May 12 13:00:07 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, 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
# Line 57  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 87  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 131  contains Line 131  contains
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)
# Line 139  contains Line 140  contains
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 152  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:: evap(:, :) ! (ncum, nl)
159      real, allocatable:: b(:, :) ! (ncum, nl - 1)      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)
# Line 170  contains Line 173  contains
173      ! set thermodynamical constants:      ! set thermodynamical constants:
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 198  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         cape1(i) = 0.      VPrecip1(:, klev + 1) = 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      CALL cv30_undilute1(klon, klev, t1, q1, qs1, gz1, plcl1, p1, nk1, icb1, &           iflag1, sig1, w01)
          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 236  contains Line 227  contains
227      end do      end do
228    
229      IF (ncum > 0) THEN      IF (ncum > 0) THEN
230         allocate(b(ncum, nl - 1))         allocate(b(ncum, nl - 1), evap(ncum, nl), icb(ncum))
231         CALL cv30_compress(ncum, iflag1, nk1, icb1, icbs1, plcl1, tnk1, qnk1, &         CALL cv30_compress(ncum, iflag1, nk1, icb1, icbs1, plcl1, tnk1, qnk1, &
232              gznk1, pbase1, buoybase1, t1, q1, qs1, u1, v1, gz1, th1, h1, lv1, &              gznk1, pbase1, buoybase1, t1, q1, qs1, u1, v1, gz1, th1, h1, lv1, &
233              cpn1, p1, ph1, tv1, tp1, tvp1, clw1, sig1, w01, iflag, nk, icb, &              cpn1, p1, ph1, tv1, tp1, tvp1, clw1, sig1, w01, iflag, nk, icb, &
234              icbs, plcl, tnk, qnk, gznk, pbase, buoybase, t, q, qs, u, v, gz, &              icbs, plcl, tnk, qnk, gznk, pbase, buoybase, t, q, qs, u, v, gz, &
235              th, h, lv, cpn, p, ph, tv, tp, tvp, clw, sig, w0)              th, h, lv, cpn, p, ph, tv, tp, tvp, clw, sig, w0)
236         CALL cv30_undilute2(ncum, icb, icbs, nk, tnk, qnk, gznk, t, qs, gz, p, &         CALL cv30_undilute2(icb, icbs(:ncum), nk, tnk, qnk, gznk, t, qs, gz, &
237              h, tv, lv, pbase, buoybase, plcl, inb(:ncum), tp, tvp, clw, hp, &              p, h, tv, lv, pbase, buoybase, plcl, inb(:ncum), tp, tvp, clw, &
238              ep, sigp, buoy)              hp, ep, buoy)
239           CALL cv30_closure(icb, inb(:ncum), pbase, p, ph, tv, buoy, sig, w0, &
240         ! CLOSURE              cape, m)
241         CALL cv30_closure(klon, ncum, klev, icb, inb, pbase, p, ph, tv, &         CALL cv30_mixing(icb, nk(:ncum), inb(:ncum), t, q, qs, u, v, h, lv, &
242              buoy, sig, w0, cape, m) ! na->klev              hp, ep, clw, m, sig, ment, qent, uent, vent, nent, sij, elij, &
243                ments, qents)
244         ! MIXING         CALL cv30_unsat(icb, inb(:ncum), t(:ncum, :nl), q(:ncum, :nl), &
245         CALL cv30_mixing(klon, ncum, klev, klev, icb, nk, inb, t, q, qs, u, &              qs(:ncum, :nl), gz, u, v, p, ph, th(:ncum, :nl - 1), tv, lv, cpn, &
246              v, h, lv, hp, ep, clw, m, sig, ment, qent, uent, vent, nent, &              ep(:ncum, :), clw(:ncum, :), m(:ncum, :), ment(:ncum, :, :), &
247              sij, elij, ments, qents)              elij(:ncum, :, :), dtphys, plcl, mp, qp(:ncum, :nl), &
248                up(:ncum, :nl), vp(:ncum, :nl), wt(:ncum, :nl), &
249         ! Unsaturated (precipitating) downdrafts              water(:ncum, :nl), evap, b)
250         CALL cv30_unsat(icb(:ncum), inb(:ncum), t, q, qs, gz, u, v, p, ph, th, &         CALL cv30_yield(icb, inb(:ncum), dtphys, t, q, u, v, gz, p, ph, h, hp, &
251              tv, lv, cpn, ep(:ncum, :), sigp(:ncum, :), clw(:ncum, :), &              lv, cpn, th, ep, clw, m, tp, mp, qp, up, vp(:ncum, 2:nl), &
252              m(:ncum, :), ment(:ncum, :, :), elij(:ncum, :, :), delt, plcl, &              wt(:ncum, :nl - 1), water(:ncum, :nl), evap, b, ment, qent, uent, &
253              mp, qp(:ncum, :nl), up(:ncum, :nl), vp(:ncum, :nl), wt, water, &              vent, nent, elij, sig, tv, tvp, iflag, precip, VPrecip, ft, fq, &
254              evap, b)              fu, fv, upwd, dnwd, dnwd0, ma, mike, tls, tps, qcondc)
   
        ! 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)  
   
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
Legend:
Removed from v.192  
changed lines
  Added in v.195
  ViewVC Help
Powered by ViewVC 1.1.21