/[lmdze]/trunk/Sources/phylmd/CV30_routines/cv30_undilute2.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/CV30_routines/cv30_undilute2.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 4  module cv30_undilute2_m Line 4  module cv30_undilute2_m
4    
5  contains  contains
6    
7    SUBROUTINE cv30_undilute2(ncum, icb, icbs, nk, tnk, qnk, gznk, t, qs, gz, &    SUBROUTINE cv30_undilute2(icb, icbs, nk, tnk, qnk, gznk, t, qs, gz, p, h, &
8         p, h, tv, lv, pbase, buoybase, plcl, inb, tp, tvp, clw, hp, ep, sigp, &         tv, lv, pbase, buoybase, plcl, inb, tp, tvp, clw, hp, ep, buoy)
        buoy)  
9    
10      ! Undilute (adiabatic) updraft, second part. Purpose: find the      ! Undilute (adiabatic) updraft, second part. Purpose: find the
11      ! rest of the lifted parcel temperatures; compute the      ! rest of the lifted parcel temperatures; compute the
# Line 16  contains Line 15  contains
15      ! Vertical profile of buoyancy computed here (use of buoybase).      ! Vertical profile of buoyancy computed here (use of buoybase).
16    
17      use conema3_m, only: epmax      use conema3_m, only: epmax
18      use cv30_param_m, only: dtovsh, minorig, nl, pbcrit, ptcrit, spfac      use cv30_param_m, only: minorig, nl
19      use cv_thermo_m, only: cl, clmcpv, cpd, cpv, eps, lv0, rrv      use cv_thermo_m, only: cl, clmcpv, cpd, cpv, eps, lv0, rrv
20      USE dimphy, ONLY: klon, klev      USE dimphy, ONLY: klon, klev
21    
22      integer, intent(in):: ncum      integer, intent(in):: icb(:), icbs(:) ! (ncum)
   
     integer, intent(in):: icb(klon), icbs(klon)  
23      ! icbs is the first level above LCL (may differ from icb)      ! icbs is the first level above LCL (may differ from icb)
24    
25      integer, intent(in):: nk(klon)      integer, intent(in):: nk(klon)
# Line 39  contains Line 36  contains
36    
37      real tp(klon, klev), tvp(klon, klev), clw(klon, klev)      real tp(klon, klev), tvp(klon, klev), clw(klon, klev)
38      ! condensed water not removed from tvp      ! condensed water not removed from tvp
39      real hp(klon, klev), ep(klon, klev), sigp(klon, klev)      real hp(klon, klev), ep(klon, klev)
40      real buoy(klon, klev)      real buoy(klon, klev)
41    
42      ! Local:      ! Local:
43    
44        integer ncum
45    
46        real, parameter:: pbcrit = 150.
47        ! critical cloud depth (mbar) beneath which the precipitation
48        ! efficiency is assumed to be zero
49    
50        real, parameter:: ptcrit = 500.
51        ! cloud depth (mbar) above which the precipitation efficiency is
52        ! assumed to be unity
53    
54        real, parameter:: dtovsh = - 0.2 ! dT for overshoot
55    
56      integer i, k      integer i, k
57      real tg, qg, ahg, alv, s, tc, es, denom      real tg, qg, ahg, alv, s, tc, es, denom
58      real pden      real pden
# Line 50  contains Line 60  contains
60    
61      !---------------------------------------------------------------------      !---------------------------------------------------------------------
62    
63        ncum = size(icb)
64    
65      ! SOME INITIALIZATIONS      ! SOME INITIALIZATIONS
66    
67      do k = 1, nl      do k = 1, nl
68         do i = 1, ncum         do i = 1, ncum
69            ep(i, k) = 0.0            ep(i, k) = 0.
           sigp(i, k) = spfac  
70         end do         end do
71      end do      end do
72    
# Line 91  contains Line 102  contains
102    
103               tc = tg - 273.15               tc = tg - 273.15
104               denom = 243.5 + tc               denom = 243.5 + tc
105               denom = MAX(denom, 1.0)               denom = MAX(denom, 1.)
106    
107               es = 6.112 * exp(17.67 * tc / denom)               es = 6.112 * exp(17.67 * tc / denom)
108    
# Line 104  contains Line 115  contains
115    
116               tc = tg - 273.15               tc = tg - 273.15
117               denom = 243.5 + tc               denom = 243.5 + tc
118               denom = MAX(denom, 1.0)               denom = MAX(denom, 1.)
119    
120               es = 6.112 * exp(17.67 * tc / denom)               es = 6.112 * exp(17.67 * tc / denom)
121    
# Line 117  contains Line 128  contains
128                    / (cpd + (cl - cpd) * qnk(i))                    / (cpd + (cl - cpd) * qnk(i))
129    
130               clw(i, k) = qnk(i) - qg               clw(i, k) = qnk(i) - qg
131               clw(i, k) = max(0.0, clw(i, k))               clw(i, k) = max(0., clw(i, k))
132               ! qg utilise au lieu du vrai mixing ratio rg:               ! qg utilise au lieu du vrai mixing ratio rg:
133               tvp(i, k) = tp(i, k) * (1. + qg / eps - qnk(i)) ! whole thing               tvp(i, k) = tp(i, k) * (1. + qg / eps - qnk(i)) ! whole thing
134            endif            endif
135         end do         end do
136      end do      end do
137    
138      ! SET THE PRECIPITATION EFFICIENCIES AND THE FRACTION OF      ! SET THE PRECIPITATION EFFICIENCIES
139      ! PRECIPITATION FALLING OUTSIDE OF CLOUD      ! It MAY BE a FUNCTION OF TP(I), P(I) AND CLW(I)
     ! THESE MAY BE FUNCTIONS OF TP(I), P(I) AND CLW(I)  
140      do k = 1, nl      do k = 1, nl
141         do i = 1, ncum         do i = 1, ncum
142            pden = ptcrit - pbcrit            pden = ptcrit - pbcrit
143            ep(i, k) = (plcl(i) - p(i, k) - pbcrit) / pden * epmax            ep(i, k) = (plcl(i) - p(i, k) - pbcrit) / pden * epmax
144            ep(i, k) = max(ep(i, k), 0.0)            ep(i, k) = max(ep(i, k), 0.)
145            ep(i, k) = min(ep(i, k), epmax)            ep(i, k) = min(ep(i, k), epmax)
           sigp(i, k) = spfac  
146         end do         end do
147      end do      end do
148    
Legend:
Removed from v.192  
changed lines
  Added in v.195
  ViewVC Help
Powered by ViewVC 1.1.21