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

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

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

trunk/phylmd/CV3_routines/cv3_uncompress.f revision 97 by guez, Fri Apr 25 14:58:31 2014 UTC trunk/Sources/phylmd/CV30_routines/cv30_uncompress.f revision 205 by guez, Tue Jun 21 15:16:03 2016 UTC
# Line 1  Line 1 
1  module cv3_uncompress_m  module cv30_uncompress_m
2    
3    implicit none    implicit none
4    
5  contains  contains
6    
7    SUBROUTINE cv3_uncompress(nloc, len, ncum, nd, idcum, iflag, precip, &    SUBROUTINE cv30_uncompress(idcum, iflag, precip, VPrecip, sig, w0, ft, fq, &
8         VPrecip, sig, w0, ft, fq, fu, fv, inb, Ma, upwd, dnwd, dnwd0, qcondc, &         fu, fv, inb, Ma, upwd, dnwd, qcondc, cape, da, phi, mp, iflag1, &
9         wd, cape, da, phi, mp, iflag1, precip1, VPrecip1, sig1, w01, ft1, fq1, &         precip1, VPrecip1, sig1, w01, ft1, fq1, fu1, fv1, inb1, Ma1, upwd1, &
10         fu1, fv1, inb1, Ma1, upwd1, dnwd1, dnwd01, qcondc1, wd1, cape1, da1, &         dnwd1, qcondc1, cape1, da1, phi1, mp1)
11         phi1, mp1)  
12        ! UNCOMPRESS THE FIELDS
13      USE cv3_param_m, ONLY: nl  
14        USE cv30_param_m, ONLY: nl
15      ! inputs:      use dimphy, only: klon, klev
16      integer, intent(in):: nloc, len, ncum, nd  
17      integer, intent(in):: idcum(nloc)      integer, intent(in):: idcum(:) ! (ncum)
18      integer, intent(in):: iflag(nloc)      integer, intent(in):: iflag(:) ! (ncum)
19      real, intent(in):: precip(nloc)      real, intent(in):: precip(klon)
20      real, intent(in):: VPrecip(nloc, nd+1)      real, intent(in):: VPrecip(klon, klev+1)
21      real, intent(in):: sig(nloc, nd), w0(nloc, nd)      real, intent(in):: sig(klon, klev), w0(klon, klev)
22      real, intent(in):: ft(nloc, nd), fq(nloc, nd), fu(nloc, nd), fv(nloc, nd)      real, intent(in), dimension(klon, klev):: ft, fq, fu, fv
23      integer, intent(in):: inb(nloc)      integer, intent(in):: inb(:) ! (ncum)
24      real, intent(in):: Ma(nloc, nd)      real, intent(in):: Ma(klon, klev)
25      real, intent(in):: upwd(nloc, nd), dnwd(nloc, nd), dnwd0(nloc, nd)      real, intent(in):: upwd(klon, klev), dnwd(klon, klev)
26      real, intent(in):: qcondc(nloc, nd)      real, intent(in):: qcondc(klon, klev)
27      real, intent(in):: wd(nloc), cape(nloc)      real, intent(in):: cape(klon)
28      real, intent(in):: da(nloc, nd), phi(nloc, nd, nd), mp(nloc, nd)      real, intent(in):: da(:, :) ! (klon, klev)
29        real, intent(in):: phi(:, :, :) ! (klon, klev, klev)
30        real, intent(in):: mp(:, :) ! (ncum, nl)
31    
32      ! outputs:      ! outputs:
33      integer iflag1(len)      integer, intent(out):: iflag1(:) ! (klon)
34      real precip1(len)      real precip1(klon)
35      real VPrecip1(len, nd+1)      real VPrecip1(klon, klev+1)
36      real sig1(len, nd), w01(len, nd)      real sig1(klon, klev), w01(klon, klev)
37      real ft1(len, nd), fq1(len, nd), fu1(len, nd), fv1(len, nd)      real ft1(klon, klev), fq1(klon, klev), fu1(klon, klev), fv1(klon, klev)
38      integer inb1(len)      integer, intent(inout):: inb1(klon)
39      real Ma1(len, nd)      real Ma1(klon, klev)
40      real upwd1(len, nd), dnwd1(len, nd), dnwd01(len, nd)      real upwd1(klon, klev), dnwd1(klon, klev)
41      real qcondc1(nloc, nd)      real qcondc1(klon, klev)
42      real wd1(nloc), cape1(nloc)      real cape1(klon)
43      real da1(nloc, nd), phi1(nloc, nd, nd), mp1(nloc, nd)      real, intent(inout):: da1(klon, klev), phi1(klon, klev, klev)
44        real, intent(inout):: mp1(klon, klev)
45    
46      ! local variables:      ! Local:
47      integer i, k, j      integer ncum, i, k, j
48    
49      !-------------------------------------------------------------------      !-------------------------------------------------------------------
50    
51      do  i=1, ncum      ncum = size(idcum)
52         precip1(idcum(i))=precip(i)      iflag1 = 42 ! for non convective points
53         iflag1(idcum(i))=iflag(i)  
54         wd1(idcum(i))=wd(i)      do  i = 1, ncum
55         inb1(idcum(i))=inb(i)         precip1(idcum(i)) = precip(i)
56         cape1(idcum(i))=cape(i)         iflag1(idcum(i)) = iflag(i)
57           inb1(idcum(i)) = inb(i)
58           cape1(idcum(i)) = cape(i)
59      end do      end do
60    
61      do  k=1, nl      do  k = 1, nl
62         do  i=1, ncum         do  i = 1, ncum
63            VPrecip1(idcum(i), k)=VPrecip(i, k)            VPrecip1(idcum(i), k) = VPrecip(i, k)
64            sig1(idcum(i), k)=sig(i, k)            sig1(idcum(i), k) = sig(i, k)
65            w01(idcum(i), k)=w0(i, k)            w01(idcum(i), k) = w0(i, k)
66            ft1(idcum(i), k)=ft(i, k)            ft1(idcum(i), k) = ft(i, k)
67            fq1(idcum(i), k)=fq(i, k)            fq1(idcum(i), k) = fq(i, k)
68            fu1(idcum(i), k)=fu(i, k)            fu1(idcum(i), k) = fu(i, k)
69            fv1(idcum(i), k)=fv(i, k)            fv1(idcum(i), k) = fv(i, k)
70            Ma1(idcum(i), k)=Ma(i, k)            Ma1(idcum(i), k) = Ma(i, k)
71            upwd1(idcum(i), k)=upwd(i, k)            upwd1(idcum(i), k) = upwd(i, k)
72            dnwd1(idcum(i), k)=dnwd(i, k)            dnwd1(idcum(i), k) = dnwd(i, k)
73            dnwd01(idcum(i), k)=dnwd0(i, k)            qcondc1(idcum(i), k) = qcondc(i, k)
74            qcondc1(idcum(i), k)=qcondc(i, k)            da1(idcum(i), k) = da(i, k)
75            da1(idcum(i), k)=da(i, k)            mp1(idcum(i), k) = mp(i, k)
           mp1(idcum(i), k)=mp(i, k)  
76         end do         end do
77      end do      end do
78    
79      do  i=1, ncum      do  i = 1, ncum
80         sig1(idcum(i), nd)=sig(i, nd)         sig1(idcum(i), klev) = sig(i, klev)
81      end do      end do
82    
83      do j=1, nd      do j = 1, klev
84         do k=1, nd         do k = 1, klev
85            do i=1, ncum            do i = 1, ncum
86               phi1(idcum(i), k, j)=phi(i, k, j)               phi1(idcum(i), k, j) = phi(i, k, j)
87            end do            end do
88         end do         end do
89      end do      end do
90    
91    end SUBROUTINE cv3_uncompress    end SUBROUTINE cv30_uncompress
92    
93  end module cv3_uncompress_m  end module cv30_uncompress_m
Legend:
Removed from v.97  
changed lines
  Added in v.205
  ViewVC Help
Powered by ViewVC 1.1.21