source: LMDZ4/trunk/libf/phylmd/cv3a_uncompress.F @ 5446

Last change on this file since 5446 was 1398, checked in by musat, 15 years ago

Last corrections for CMIP5:

  • Add O3 at standard level files histmthNMC.nc
  • Add positive attribute "down" for vertical axes for all output files
  • Replace "inst" by "ave" for hist*NMC.nc files to have time_counter and bounds for time axis (Marie-Alice's hint)
  • Correct units for vertical axes : mb instead of hPa
  • Add mass flux at the bottom of clouds
  • Comment non initialized variables (s_capCL, s_oliqCL, s_cteiCL, s_trmb1, s_trmb2, s_trmb3) for the output files
  • Geopotential field phy850, phi700, phi500, etc are modified to "geopotential height and are called z850, z700, z500, etc
  • Meaning of specific humidity outputs - ovapinit and ovap - were interchanged
  • Fields albs, albslw become alb1, alb2 in output files
  • Correct title for rugs_* fields
  • Correct units for pbase and ptop are Pa (not mb)
  • Correct ndayrain field

FH/JLD/JYG/MAF/IM

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 4.2 KB
RevLine 
[879]1      SUBROUTINE cv3a_uncompress(nloc,len,ncum,nd,ntra,idcum
2     :         ,iflag,kbas,ktop
[1398]3     :         ,precip,cbmf,sig,w0,ptop2
[879]4     :         ,ft,fq,fu,fv,ftra
5     :         ,Ma,mip,Vprecip,upwd,dnwd,dnwd0
6     :         ,qcondc,wd,cape,cin
7     :         ,tvp
8     :         ,ftd,fqd
9     :         ,Plim1,Plim2,asupmax,supmax0
10     :         ,asupmaxmin
11     o         ,iflag1,kbas1,ktop1
[1398]12     :         ,precip1,cbmf1,sig1,w01,ptop21
[879]13     :         ,ft1,fq1,fu1,fv1,ftra1
14     :         ,Ma1,mip1,Vprecip1,upwd1,dnwd1,dnwd01
15     :         ,qcondc1,wd1,cape1,cin1
16     :         ,tvp1
17     :         ,ftd1,fqd1
18     :         ,Plim11,Plim21,asupmax1,supmax01
19     :         ,asupmaxmin1     )
20***************************************************************
21*                                                             *
22* CV3A_UNCOMPRESS                                             *
23*                                                             *
24*                                                             *
25* written by   : Sandrine Bony-Lena , 17/05/2003, 11.22.15    *
26* modified by  : Jean-Yves Grandpeix, 23/06/2003, 10.36.17    *
27***************************************************************
28*
29      implicit none
30
31#include "cv3param.h"
32
33c inputs:
34      integer nloc, len, ncum, nd, ntra
35      integer idcum(nloc)
36      integer iflag(nloc),kbas(nloc),ktop(nloc)
[1398]37      real precip(nloc),cbmf(nloc)
[879]38      real sig(nloc,nd), w0(nloc,nd),ptop2(nloc)
39      real ft(nloc,nd), fq(nloc,nd), fu(nloc,nd), fv(nloc,nd)
40      real ftra(nloc,nd,ntra)
[1334]41      real Ma(nloc,nd),mip(nloc,nd),Vprecip(nloc,nd+1)
[879]42      real upwd(nloc,nd),dnwd(nloc,nd),dnwd0(nloc,nd)
43      real qcondc(nloc,nd)
44      real wd(nloc),cape(nloc),cin(nloc)
45      real tvp(nloc,nd)
46      real ftd(nloc,nd), fqd(nloc,nd)
47      real Plim1(nloc),Plim2(nloc)
48      real asupmax(nloc,nd),supmax0(nloc)
49      real asupmaxmin(nloc)
50
51c outputs:
52      integer iflag1(len),kbas1(len),ktop1(len)
[1398]53      real precip1(len),cbmf1(len)
[879]54      real sig1(len,nd), w01(len,nd),ptop21(len)
55      real ft1(len,nd), fq1(len,nd), fu1(len,nd), fv1(len,nd)
56      real ftra1(len,nd,ntra)
[1334]57      real Ma1(len,nd),mip1(len,nd),Vprecip1(len,nd+1)
[879]58      real upwd1(len,nd),dnwd1(len,nd),dnwd01(len,nd)
59      real qcondc1(len,nd)
60      real wd1(len),cape1(len),cin1(len)
61      real tvp1(len,nd)
62      real ftd1(len,nd), fqd1(len,nd)
63      real Plim11(len),Plim21(len)
64      real asupmax1(len,nd),supmax01(len)
65      real asupmaxmin1(len)
66c
67c local variables:
68      integer i,k,j,k1,k2
69
70        do 2000 i=1,ncum
71         ptop21(idcum(i))=ptop2(i)
72         precip1(idcum(i))=precip(i)
[1398]73         cbmf1(idcum(i))=cbmf(i)
[879]74         iflag1(idcum(i))=iflag(i)
75         kbas1(idcum(i))=kbas(i)
76         ktop1(idcum(i))=ktop(i)
77         wd1(idcum(i))=wd(i)
78         cape1(idcum(i))=cape(i)
79         cin1(idcum(i))=cin(i)
80         Plim11(idcum(i))=Plim1(i)
81         Plim21(idcum(i))=Plim2(i)
82         supmax01(idcum(i))=supmax0(i)
83         asupmaxmin1(idcum(i))=asupmaxmin(i)
84 2000   continue
85
86        do 2020 k=1,nd
87          do 2010 i=1,ncum
88            sig1(idcum(i),k)=sig(i,k)
89            w01(idcum(i),k)=w0(i,k)
90            ft1(idcum(i),k)=ft(i,k)
91            fq1(idcum(i),k)=fq(i,k)
92            fu1(idcum(i),k)=fu(i,k)
93            fv1(idcum(i),k)=fv(i,k)
94            Ma1(idcum(i),k)=Ma(i,k)
95            mip1(idcum(i),k)=mip(i,k)
96            Vprecip1(idcum(i),k)=Vprecip(i,k)
97            upwd1(idcum(i),k)=upwd(i,k)
98            dnwd1(idcum(i),k)=dnwd(i,k)
99            dnwd01(idcum(i),k)=dnwd0(i,k)
100            qcondc1(idcum(i),k)=qcondc(i,k)
101            tvp1(idcum(i),k)=tvp(i,k)
102            ftd1(idcum(i),k)=ftd(i,k)
103            fqd1(idcum(i),k)=fqd(i,k)
104            asupmax1(idcum(i),k)=asupmax(i,k)
105 2010     continue
106 2020   continue
107
108        do 2040 i=1,ncum
109          sig1(idcum(i),nd)=sig(i,nd)
1102040    continue
111
112
113        do 2100 j=1,ntra
114c oct3         do 2110 k=1,nl
115         do 2110 k=1,nd ! oct3
116          do 2120 i=1,ncum
117            ftra1(idcum(i),k,j)=ftra(i,k,j)
118 2120     continue
119 2110    continue
120 2100   continue
121c
122c        do 2220 k2=1,nd
123c         do 2210 k1=1,nd
124c          do 2200 i=1,ncum
125c            ment1(idcum(i),k1,k2) = ment(i,k1,k2)
126c            sij1(idcum(i),k1,k2) = sij(i,k1,k2)
127c2200      enddo
128c2210     enddo
129c2220    enddo
130
131       RETURN
132      END
133
Note: See TracBrowser for help on using the repository browser.