/[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

trunk/Sources/phylmd/CV30_routines/cv3_compress.f revision 184 by guez, Wed Mar 16 14:50:46 2016 UTC trunk/Sources/phylmd/CV30_routines/cv30_compress.f revision 196 by guez, Mon May 23 13:50:39 2016 UTC
# Line 1  Line 1 
1  module cv3_compress_m  module cv30_compress_m
2    
3    implicit none    implicit none
4    
5  contains  contains
6    
7    SUBROUTINE cv3_compress(len, nloc, ncum, nd, iflag1, nk1, icb1, icbs1, &    SUBROUTINE cv30_compress(iflag1, nk1, icb1, icbs1, plcl1, tnk1, qnk1, &
8         plcl1, tnk1, qnk1, gznk1, pbase1, buoybase1, t1, q1, qs1, u1, v1, gz1, &         gznk1, pbase1, buoybase1, t1, q1, qs1, u1, v1, gz1, th1, h1, lv1, cpn1, &
9         th1, h1, lv1, cpn1, p1, ph1, tv1, tp1, tvp1, clw1, sig1, w01, iflag, &         p1, ph1, tv1, tp1, tvp1, clw1, sig1, w01, nk, icb, icbs, plcl, tnk, &
10         nk, icb, icbs, plcl, tnk, qnk, gznk, pbase, buoybase, t, q, qs, u, v, &         qnk, gznk, pbase, buoybase, t, q, qs, u, v, gz, th, h, lv, cpn, p, ph, &
11         gz, th, h, lv, cpn, p, ph, tv, tp, tvp, clw, sig, w0)         tv, tp, tvp, clw, sig, w0)
12    
13      use cv3_param_m      ! Compress the fields (vectorization over convective gridpoints).
14    
15        use cv30_param_m, only: nl
16        USE dimphy, ONLY: klev, klon
17    
18      ! inputs:      ! inputs:
19      integer, intent(in):: len, ncum, nd, nloc      integer, intent(in):: iflag1(:), nk1(:), icb1(:), icbs1(:) ! (klon)
20      integer iflag1(len), nk1(len), icb1(len), icbs1(len)      real, intent(in):: plcl1(klon), tnk1(klon), qnk1(klon), gznk1(klon)
21      real plcl1(len), tnk1(len), qnk1(len), gznk1(len)      real pbase1(klon), buoybase1(klon)
22      real pbase1(len), buoybase1(len)      real, intent(in):: t1(klon, klev)
23      real, intent(in):: t1(len, nd)      real, intent(in):: q1(klon, klev), qs1(klon, klev)
24      real, intent(in):: q1(len, nd), qs1(len, nd), u1(len, nd), v1(len, nd)      real, intent(in):: u1(klon, klev), v1(klon, klev)
25      real gz1(len, nd), h1(len, nd), lv1(len, nd), cpn1(len, nd)      real gz1(klon, klev), h1(klon, klev), lv1(klon, klev), cpn1(klon, klev)
26      real, intent(in):: p1(len, nd), ph1(len, nd+1), tv1(len, nd), tp1(len, nd)      real, intent(in):: p1(klon, klev), ph1(klon, klev + 1)
27      real tvp1(len, nd), clw1(len, nd)      real, intent(in):: tv1(klon, klev), tp1(klon, klev)
28      real th1(len, nd)      real tvp1(klon, klev), clw1(klon, klev)
29      real sig1(len, nd), w01(len, nd)      real th1(klon, klev)
30        real sig1(klon, klev), w01(klon, klev)
31    
32      ! outputs:      ! outputs:
33      ! en fait, on a nloc=len pour l'instant (cf cv_driver)      integer nk(:) ! (klon)
34      integer iflag(nloc), nk(nloc), icb(nloc), icbs(nloc)  
35      real plcl(nloc), tnk(nloc), qnk(nloc), gznk(nloc)      integer, intent(out):: icb(:) ! (ncum)
36      real pbase(nloc), buoybase(nloc)      ! {2 <= icb <= nl - 3}
37      real t(nloc, nd), q(nloc, nd), qs(nloc, nd), u(nloc, nd), v(nloc, nd)      ! {ph(i, icb(i) + 1) < plcl(i) <= ph(i, icb(i))}
38      real gz(nloc, nd), h(nloc, nd), lv(nloc, nd), cpn(nloc, nd)  
39      real p(nloc, nd), ph(nloc, nd+1), tv(nloc, nd), tp(nloc, nd)      integer icbs(klon)
40      real tvp(nloc, nd), clw(nloc, nd)      real, intent(out):: plcl(:) ! (ncum)
41      real th(nloc, nd)      real tnk(:), qnk(:), gznk(:) ! (klon)
42      real sig(nloc, nd), w0(nloc, nd)      real pbase(klon), buoybase(klon)
43        real t(klon, klev), q(klon, klev), qs(klon, klev)
44      ! local variables:      real u(klon, klev), v(klon, klev)
45      integer i, k, nn      real gz(klon, klev), h(klon, klev), lv(klon, klev), cpn(klon, klev)
46        real p(klon, klev)
47        real ph(:, :) ! (klon, klev + 1)
48      do  k=1, nl+1      real tv(klon, klev), tp(klon, klev)
49         nn=0      real tvp(klon, klev), clw(klon, klev)
50         do  i=1, len      real th(klon, klev)
51            if(iflag1(i).eq.0)then      real sig(klon, klev), w0(klon, klev)
52               nn=nn+1  
53               sig(nn, k)=sig1(i, k)      ! Local:
54               w0(nn, k)=w01(i, k)      integer i, k, nn, ncum
55               t(nn, k)=t1(i, k)  
56               q(nn, k)=q1(i, k)      !---------------------------------------------------------------
57               qs(nn, k)=qs1(i, k)  
58               u(nn, k)=u1(i, k)      ncum = size(icb)
59               v(nn, k)=v1(i, k)  
60               gz(nn, k)=gz1(i, k)      do k = 1, nl + 1
61               h(nn, k)=h1(i, k)         nn = 0
62               lv(nn, k)=lv1(i, k)         do i = 1, klon
63               cpn(nn, k)=cpn1(i, k)            if (iflag1(i) == 0) then
64               p(nn, k)=p1(i, k)               nn = nn + 1
65               ph(nn, k)=ph1(i, k)               sig(nn, k) = sig1(i, k)
66               tv(nn, k)=tv1(i, k)               w0(nn, k) = w01(i, k)
67               tp(nn, k)=tp1(i, k)               t(nn, k) = t1(i, k)
68               tvp(nn, k)=tvp1(i, k)               q(nn, k) = q1(i, k)
69               clw(nn, k)=clw1(i, k)               qs(nn, k) = qs1(i, k)
70               th(nn, k)=th1(i, k)               u(nn, k) = u1(i, k)
71                 v(nn, k) = v1(i, k)
72                 gz(nn, k) = gz1(i, k)
73                 h(nn, k) = h1(i, k)
74                 lv(nn, k) = lv1(i, k)
75                 cpn(nn, k) = cpn1(i, k)
76                 p(nn, k) = p1(i, k)
77                 ph(nn, k) = ph1(i, k)
78                 tv(nn, k) = tv1(i, k)
79                 tp(nn, k) = tp1(i, k)
80                 tvp(nn, k) = tvp1(i, k)
81                 clw(nn, k) = clw1(i, k)
82                 th(nn, k) = th1(i, k)
83            endif            endif
84         end do         end do
85      end do      end do
86    
87      if (nn.ne.ncum) then      nn = 0
88         print*, 'strange! nn not equal to ncum: ', nn, ncum  
89         stop      do i = 1, klon
90      endif         if (iflag1(i) == 0) then
91              nn = nn + 1
92      nn=0            pbase(nn) = pbase1(i)
93      do  i=1, len            buoybase(nn) = buoybase1(i)
94         if(iflag1(i).eq.0)then            plcl(nn) = plcl1(i)
95            nn=nn+1            tnk(nn) = tnk1(i)
96            pbase(nn)=pbase1(i)            qnk(nn) = qnk1(i)
97            buoybase(nn)=buoybase1(i)            gznk(nn) = gznk1(i)
98            plcl(nn)=plcl1(i)            nk(nn) = nk1(i)
99            tnk(nn)=tnk1(i)            icb(nn) = icb1(i)
100            qnk(nn)=qnk1(i)            icbs(nn) = icbs1(i)
           gznk(nn)=gznk1(i)  
           nk(nn)=nk1(i)  
           icb(nn)=icb1(i)  
           icbs(nn)=icbs1(i)  
           iflag(nn)=iflag1(i)  
101         endif         endif
102      end do      end do
103    
104    end SUBROUTINE cv3_compress    end SUBROUTINE cv30_compress
105    
106  end module cv3_compress_m  end module cv30_compress_m
Legend:
Removed from v.184  
changed lines
  Added in v.196
  ViewVC Help
Powered by ViewVC 1.1.21