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

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

Parent Directory Parent Directory | Revision Log Revision Log

Revision 189 - (hide annotations)
Tue Mar 29 15:20:23 2016 UTC (8 years, 1 month ago) by guez
File size: 2678 byte(s) 
There was a function gr_phy_write_3d in dyn3d and a function
gr_phy_write_2d in module grid_change. Moved them into a new module
gr_phy_write_m under a generic interface gr_phy_write. Replaced calls
to gr_fi_ecrit by calls to gr_phy_write.

Removed arguments len, nloc and nd of cv30_compress.

Removed arguments wd and wd1 of cv30_uncompress, wd of cv30_yield, wd
of concvl, wd1 of cv_driver. Was just filled with 0. Removed option
ok_gust in physiq, never used.

In cv30_unsat, cv30_yield and cv_driver, we only need to define b to
level nl - 1.
1 guez 185 module cv30_uncompress_m
2 guez 47
3     implicit none
4    
5 guez 97 contains
6 guez 47
7 guez 185 SUBROUTINE cv30_uncompress(idcum, iflag, precip, VPrecip, sig, w0, ft, fq, &
8 guez 189 fu, fv, inb, Ma, upwd, dnwd, dnwd0, qcondc, cape, da, phi, mp, iflag1, &
9     precip1, VPrecip1, sig1, w01, ft1, fq1, fu1, fv1, inb1, Ma1, upwd1, &
10     dnwd1, dnwd01, qcondc1, cape1, da1, phi1, mp1)
11 guez 47
12 guez 185 USE cv30_param_m, ONLY: nl
13 guez 103 use dimphy, only: klon, klev
14 guez 47
15 guez 103 integer, intent(in):: idcum(:) ! (ncum)
16     integer, intent(in):: iflag(klon)
17     real, intent(in):: precip(klon)
18     real, intent(in):: VPrecip(klon, klev+1)
19     real, intent(in):: sig(klon, klev), w0(klon, klev)
20     real, intent(in), dimension(klon, klev):: ft, fq, fu, fv
21     integer, intent(in):: inb(klon)
22     real, intent(in):: Ma(klon, klev)
23     real, intent(in):: upwd(klon, klev), dnwd(klon, klev), dnwd0(klon, klev)
24     real, intent(in):: qcondc(klon, klev)
25 guez 189 real, intent(in):: cape(klon)
26 guez 103 real, intent(in):: da(klon, klev), phi(klon, klev, klev), mp(klon, klev)
27 guez 47
28 guez 97 ! outputs:
29 guez 103 integer iflag1(klon)
30     real precip1(klon)
31     real VPrecip1(klon, klev+1)
32     real sig1(klon, klev), w01(klon, klev)
33     real ft1(klon, klev), fq1(klon, klev), fu1(klon, klev), fv1(klon, klev)
34     integer, intent(inout):: inb1(klon)
35     real Ma1(klon, klev)
36     real upwd1(klon, klev), dnwd1(klon, klev), dnwd01(klon, klev)
37     real qcondc1(klon, klev)
38 guez 189 real cape1(klon)
39 guez 103 real, intent(inout):: da1(klon, klev), phi1(klon, klev, klev)
40     real, intent(inout):: mp1(klon, klev)
41 guez 47
42 guez 180 ! Local:
43 guez 103 integer ncum, i, k, j
44 guez 47
45 guez 97 !-------------------------------------------------------------------
46 guez 47
47 guez 103 ncum = size(idcum)
48    
49 guez 97 do i=1, ncum
50     precip1(idcum(i))=precip(i)
51     iflag1(idcum(i))=iflag(i)
52     inb1(idcum(i))=inb(i)
53     cape1(idcum(i))=cape(i)
54     end do
55 guez 47
56 guez 97 do k=1, nl
57     do i=1, ncum
58     VPrecip1(idcum(i), k)=VPrecip(i, k)
59     sig1(idcum(i), k)=sig(i, k)
60     w01(idcum(i), k)=w0(i, k)
61     ft1(idcum(i), k)=ft(i, k)
62     fq1(idcum(i), k)=fq(i, k)
63     fu1(idcum(i), k)=fu(i, k)
64     fv1(idcum(i), k)=fv(i, k)
65     Ma1(idcum(i), k)=Ma(i, k)
66     upwd1(idcum(i), k)=upwd(i, k)
67     dnwd1(idcum(i), k)=dnwd(i, k)
68     dnwd01(idcum(i), k)=dnwd0(i, k)
69     qcondc1(idcum(i), k)=qcondc(i, k)
70     da1(idcum(i), k)=da(i, k)
71     mp1(idcum(i), k)=mp(i, k)
72     end do
73     end do
74    
75     do i=1, ncum
76 guez 103 sig1(idcum(i), klev)=sig(i, klev)
77 guez 97 end do
78    
79 guez 103 do j=1, klev
80     do k=1, klev
81 guez 97 do i=1, ncum
82     phi1(idcum(i), k, j)=phi(i, k, j)
83     end do
84     end do
85     end do
86    
87 guez 185 end SUBROUTINE cv30_uncompress
88 guez 97
89 guez 185 end module cv30_uncompress_m
  ViewVC Help
Powered by ViewVC 1.1.21