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

Diff of /trunk/Sources/phylmd/CV30_routines/cv30_compress.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 195 by guez, Wed May 18 17:56:44 2016 UTC
# Line 14  contains Line 14  contains
14    
15      use cv30_param_m, only: nl      use cv30_param_m, only: nl
16      USE dimphy, ONLY: klev, klon      USE dimphy, ONLY: klev, klon
17        use nr_util, only: assert
18    
19      ! inputs:      ! inputs:
20      integer, intent(in):: ncum      integer, intent(in):: ncum
21      integer iflag1(klon), nk1(klon), icb1(klon), icbs1(klon)      integer, intent(in):: iflag1(klon), nk1(klon), icb1(klon), icbs1(klon)
22      real plcl1(klon), tnk1(klon), qnk1(klon), gznk1(klon)      real, intent(in):: plcl1(klon), tnk1(klon), qnk1(klon), gznk1(klon)
23      real pbase1(klon), buoybase1(klon)      real pbase1(klon), buoybase1(klon)
24      real, intent(in):: t1(klon, klev)      real, intent(in):: t1(klon, klev)
25      real, intent(in):: q1(klon, klev), qs1(klon, klev)      real, intent(in):: q1(klon, klev), qs1(klon, klev)
# Line 31  contains Line 32  contains
32      real sig1(klon, klev), w01(klon, klev)      real sig1(klon, klev), w01(klon, klev)
33    
34      ! outputs:      ! outputs:
35      integer iflag(klon), nk(klon), icb(klon), icbs(klon)      integer iflag(klon), nk(klon)
36        integer, intent(out):: icb(:) ! (ncum)
37        integer icbs(klon)
38      real plcl(klon), tnk(klon), qnk(klon), gznk(klon)      real plcl(klon), tnk(klon), qnk(klon), gznk(klon)
39      real pbase(klon), buoybase(klon)      real pbase(klon), buoybase(klon)
40      real t(klon, klev), q(klon, klev), qs(klon, klev)      real t(klon, klev), q(klon, klev), qs(klon, klev)
# Line 74  contains Line 77  contains
77         end do         end do
78      end do      end do
79    
80      if (nn /= ncum) then      call assert(nn == ncum, "cv30_compress")
        print*, 'strange! nn not equal to ncum: ', nn, ncum  
        stop 1  
     endif  
   
81      nn=0      nn=0
82    
83      do i=1, klon      do i=1, klon
84         if (iflag1(i) == 0) then         if (iflag1(i) == 0) then
85            nn=nn+1            nn=nn+1
Legend:
Removed from v.189  
changed lines
  Added in v.195
  ViewVC Help
Powered by ViewVC 1.1.21