/[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 199 by guez, Tue May 31 16:17:35 2016 UTC revision 200 by guez, Thu Jun 2 15:40:30 2016 UTC
# Line 26  contains Line 26  contains
26    
27      integer, intent(out):: iflag1(:) ! (klon)      integer, intent(out):: iflag1(:) ! (klon)
28      real tnk1(:), qnk1(:), 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 qsnk(klon), rh(klon), chi(klon)      real rh(klon)
34      real, parameter:: A = 1669., B = 122.      real, parameter:: A = 1669., B = 122.
35    
36      !--------------------------------------------------------------------      !--------------------------------------------------------------------
37    
     iflag1 = 0  
     plcl1 = 0.  
   
     ! Check whether parcel level temperature and specific humidity  
     ! are reasonable  
     do i = 1, klon  
        if (t1(i, minorig) < 250. .or. q1(i, minorig) <= 0.) iflag1(i) = 7  
     end do  
   
38      ! Calculate lifted condensation level of air at parcel origin level      ! Calculate lifted condensation level of air at parcel origin level
39      ! (within 0.2 % of formula of Bolton, Mon. Wea. Rev., 1980)      ! (within 0.2 % of formula of Bolton, Mon. Wea. Rev., 1980)
40      do i = 1, klon      where (t1(:, minorig) >= 250. .and. q1(:, minorig) > 0.)
41         if (iflag1(i) == 0) then         ! Parcel level temperature and specific humidity are reasonable.
42            tnk1(i) = t1(i, minorig)         tnk1 = t1(:, minorig)
43            qnk1(i) = q1(i, minorig)         qnk1 = q1(:, minorig)
44            gznk1(i) = gz1(i, minorig)         gznk1 = gz1(:, minorig)
45            qsnk(i) = qs1(i, minorig)  
46           rh = qnk1 / qs1(:, minorig)
47            rh(i) = qnk1(i) / qsnk(i)         plcl1 = p1(:, minorig) * rh**(tnk1 / (A - B * rh - tnk1))
48            chi(i) = tnk1(i) / (A - B * rh(i) - tnk1(i))         iflag1 = 0
49            plcl1(i) = p1(i, minorig) * (rh(i)**chi(i))      elsewhere
50            if (plcl1(i) < 200. .or. plcl1(i) >= 2000.) iflag1(i) = 8         plcl1 = 0.
51         endif         iflag1 = 7
52      end do      end where
53    
54        where (iflag1 == 0 .and. (plcl1 < 200. .or. plcl1 >= 2000.)) iflag1 = 8
55    
56      ! Compute icb1:      ! Compute icb1:
57      do i = 1, klon      do i = 1, klon
Legend:
Removed from v.199  
changed lines
  Added in v.200
  ViewVC Help
Powered by ViewVC 1.1.21