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

Diff of /trunk/Sources/phylmd/CV30_routines/cv30_feed.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 197 by guez, Tue May 24 12:25:29 2016 UTC revision 198 by guez, Tue May 31 16:17:35 2016 UTC
# Line 4  module cv30_feed_m Line 4  module cv30_feed_m
4    
5  contains  contains
6    
7    SUBROUTINE cv30_feed(t1, q1, qs1, p1, ph1, gz1, nk1, icb1, iflag1, tnk1, &    SUBROUTINE cv30_feed(t1, q1, qs1, p1, ph1, gz1, icb1, iflag1, tnk1, qnk1, &
8         qnk1, gznk1, plcl1)         gznk1, plcl1)
9    
10      ! Purpose: convective feed      ! Purpose: convective feed
11        ! Assuming origin level of ascending parcels is minorig.
12    
13      use cv30_param_m, only: minorig, nl      use cv30_param_m, only: minorig, nl
14      USE dimphy, ONLY: klev, klon      USE dimphy, ONLY: klev, klon
# Line 20  contains Line 21  contains
21    
22      ! outputs:      ! outputs:
23    
     integer, intent(out):: nk1(:) ! (klon)  
   
24      integer, intent(out):: icb1(:) ! (klon)      integer, intent(out):: icb1(:) ! (klon)
25      ! first level above lcl, 2 <= icb1 <= nl - 2      ! first level above lcl, 2 <= icb1 <= nl - 2
26    
27      integer, intent(out):: iflag1(:) ! (klon)      integer, intent(out):: iflag1(:) ! (klon)
28      real tnk1(klon), qnk1(klon), gznk1(klon)      real tnk1(:), qnk1(:), gznk1(:) ! (klon)
29      real, intent(out):: plcl1(klon)      real, intent(out):: plcl1(klon)
30    
31      ! Local:      ! Local:
32      integer i      integer i
33      real pnk(klon), qsnk(klon), rh(klon), chi(klon)      real qsnk(klon), rh(klon), chi(klon)
34      real, parameter:: A = 1669., B = 122.      real, parameter:: A = 1669., B = 122.
35    
36      !--------------------------------------------------------------------      !--------------------------------------------------------------------
# Line 39  contains Line 38  contains
38      iflag1 = 0      iflag1 = 0
39      plcl1 = 0.      plcl1 = 0.
40    
     ! Origin level of ascending parcels  
     nk1 = minorig  
   
41      ! Check whether parcel level temperature and specific humidity      ! Check whether parcel level temperature and specific humidity
42      ! are reasonable      ! are reasonable
   
43      do i = 1, klon      do i = 1, klon
44         if (t1(i, nk1(i)) < 250. .or. q1(i, nk1(i)) <= 0.) iflag1(i) = 7         if (t1(i, minorig) < 250. .or. q1(i, minorig) <= 0.) iflag1(i) = 7
45      end do      end do
46    
47      ! Calculate lifted condensation level of air at parcel origin level      ! Calculate lifted condensation level of air at parcel origin level
48      ! (within 0.2 % of formula of Bolton, Mon. Wea. Rev., 1980)      ! (within 0.2 % of formula of Bolton, Mon. Wea. Rev., 1980)
   
49      do i = 1, klon      do i = 1, klon
50         if (iflag1(i) == 0) then         if (iflag1(i) == 0) then
51            tnk1(i) = t1(i, nk1(i))            tnk1(i) = t1(i, minorig)
52            qnk1(i) = q1(i, nk1(i))            qnk1(i) = q1(i, minorig)
53            gznk1(i) = gz1(i, nk1(i))            gznk1(i) = gz1(i, minorig)
54            pnk(i) = p1(i, nk1(i))            qsnk(i) = qs1(i, minorig)
           qsnk(i) = qs1(i, nk1(i))  
55    
56            rh(i) = qnk1(i) / qsnk(i)            rh(i) = qnk1(i) / qsnk(i)
57            chi(i) = tnk1(i) / (A - B*rh(i) - tnk1(i))            chi(i) = tnk1(i) / (A - B * rh(i) - tnk1(i))
58            plcl1(i) = pnk(i)*(rh(i)**chi(i))            plcl1(i) = p1(i, minorig) * (rh(i)**chi(i))
59            if (plcl1(i) < 200. .or. plcl1(i) >= 2000.) iflag1(i) = 8            if (plcl1(i) < 200. .or. plcl1(i) >= 2000.) iflag1(i) = 8
60         endif         endif
61      end do      end do
Legend:
Removed from v.197  
changed lines
  Added in v.198
  ViewVC Help
Powered by ViewVC 1.1.21