source: LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_bilKP_ave.h @ 5220

Last change on this file since 5220 was 740, checked in by lmdzadmin, 18 years ago

Correction bogues: les ecrit_ sont des REALs lus dans conf_phys.F90 en
nombre de jours sauf pour ecrit_ins et ecrit_tra en secondes!
Les ecrit_ sont initialises dans conf_phys.F90 et peuvent etre modifies dans
physiq.def.
IM
  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 8.6 KB
Line 
1c
2c $Header$
3c
4      IF (ok_journe) THEN
5c
6         zsto = dtime
7         zout = ecrit_day
8         typeval=tave
9c
10         idayref = day_ref
11         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
12c
13         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
14         DO i = 1, iim
15            zx_lon(i,1) = rlon(i+1)
16            zx_lon(i,jjmp1) = rlon(i+1)
17         ENDDO
18         DO ll=1,klev
19            znivsig(ll)=float(ll)
20         ENDDO
21         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
22         write(*,*)'zx_lon = ',zx_lon(:,1)
23         write(*,*)'zx_lat = ',zx_lat(1,:)
24         CALL histbeg("histbilKP_ave", iim,zx_lon(:,1), jjmp1,
25     .                zx_lat(1,:),
26     .                1,iim,1,jjmp1, itau_phy, zjulian, dtime,
27     .                nhori, nid_bilKPave)
28         write(*,*)'Journee ', itau_phy, zjulian
29         CALL histvert(nid_bilKPave, "presnivs",
30     .                "Vertical levels","mb",
31     .                 klev, presnivs/100., nvert)
32c
33c
34c Champs 3D:
35c
36         CALL histdef(nid_bilKPave,"ue",
37     .   "Zonal energy transport","-",
38     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32, 
39     .                typeval, zsto,zout)
40c
41         CALL histdef(nid_bilKPave,"ve",
42     .   "Merid energy transport","-",
43     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32, 
44     .                typeval, zsto,zout)
45c
46         CALL histdef(nid_bilKPave,"uq",
47     .   "Zonal humidity transport","-",
48     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32, 
49     .                typeval, zsto,zout)
50c
51         CALL histdef(nid_bilKPave,"vq",
52     .   "Merid humidity transport","-",
53     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32, 
54     .                typeval, zsto,zout)
55c
56c Champs 3D:
57c
58         CALL histdef(nid_bilKPave,"temp",
59     .   "Air temperature","K",
60     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
61     .                typeval, zsto,zout)
62c
63         CALL histdef(nid_bilKPave,"ovap",
64     .   "Specific humidity","Kg/Kg",
65     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
66     .                typeval, zsto,zout)
67c
68         CALL histdef(nid_bilKPave,"geop",
69     .   "Geopotential height","m",
70     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
71     .                typeval, zsto,zout)
72c
73         CALL histdef(nid_bilKPave,"vitu",
74     .   "Zonal wind","m/s",
75     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
76     .                typeval, zsto,zout)
77c
78         CALL histdef(nid_bilKPave,"vitv",
79     .   "Meridional wind","m/s",
80     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
81     .                typeval, zsto,zout)
82c
83         CALL histdef(nid_bilKPave, "vitw", 
84     .   "Vertical wind", "m/s",
85     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
86     .                typeval, zsto,zout)
87c
88         CALL histdef(nid_bilKPave, "pres", 
89     .   "Inter-Layer Air pressure",
90     .                "Pa",
91     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
92     .                typeval, zsto,zout)
93c
94         CALL histdef(nid_bilKPave, "play", 
95     .   "Mean-Layer Air pressure",
96     .                "Pa",
97     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
98     .                typeval, zsto,zout)
99c
100         CALL histdef(nid_bilKPave, "oliq", 
101     .   "Liquid water content", 
102     .                "kg/kg",
103     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
104     .                typeval, zsto,zout)
105c
106         CALL histdef(nid_bilKPave, "dtdyn", 
107     .   "Dynamics dT", "K/s",
108     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
109     .                typeval, zsto,zout)
110c
111         CALL histdef(nid_bilKPave, "dqdyn", 
112     .   "Dynamics dQ", "Kg/Kg/s",
113     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
114     .                typeval, zsto,zout)
115c
116         CALL histdef(nid_bilKPave, "dtcon", 
117     .   "Convection dT", "K/s",
118     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
119     .                typeval, zsto,zout)
120c
121         CALL histdef(nid_bilKPave, "ducon", 
122     .   "Convection du", "m/s2",
123     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
124     .                typeval, zsto,zout)
125c
126         CALL histdef(nid_bilKPave, "dvcon", 
127     .   "Convection dv", "m/s2",
128     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
129     .                typeval, zsto,zout)
130c
131         CALL histdef(nid_bilKPave,"dqcon",
132     .   "Convection dQ","Kg/Kg/s",
133     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
134     .                typeval, zsto,zout)
135c
136         CALL histdef(nid_bilKPave, "dtlsc", 
137     .   "Condensation dT", "K/s",
138     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
139     .                typeval, zsto,zout)
140c
141         CALL histdef(nid_bilKPave,"dqlsc",
142     .   "Condensation dQ","Kg/Kg/s",
143     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
144     .                typeval, zsto,zout)
145c
146         CALL histdef(nid_bilKPave,"dtvdf",
147     .   "Boundary-layer dT","K/s",
148     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
149     .                typeval, zsto,zout)
150c
151         CALL histdef(nid_bilKPave, "dqvdf", 
152     .   "Boundary-layer dQ", 
153     .               "Kg/Kg/s",
154     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
155     .                typeval, zsto,zout)
156c
157         CALL histdef(nid_bilKPave,"dtajs",
158     .   "Ajustement sec dT","K/s",
159     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
160     .                typeval, zsto,zout)
161c
162         CALL histdef(nid_bilKPave, "dqajs",
163     .   "Ajustement sec dQ", 
164     .               "Kg/Kg/s",
165     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
166     .                typeval, zsto,zout)
167c
168         CALL histdef(nid_bilKPave, "dteva", 
169     .   "Reevaporation dT", "K/s",
170     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
171     .                typeval, zsto,zout)
172c
173         CALL histdef(nid_bilKPave,"dqeva",
174     .   "Reevaporation dQ",
175     .                "Kg/Kg/s",
176     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
177     .                typeval, zsto,zout)
178
179c
180         CALL histdef(nid_bilKPave, "dtswr", 
181     .   "SW radiation dT", "K/s",
182     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
183     .                typeval, zsto,zout)
184c
185         CALL histdef(nid_bilKPave, "dtsw0", 
186     .   "SW radiation dT", "K/s",
187     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
188     .                typeval, zsto,zout)
189c
190         CALL histdef(nid_bilKPave, "dtlwr", 
191     .   "LW radiation dT", "K/s",
192     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
193     .                typeval, zsto,zout)
194c
195         CALL histdef(nid_bilKPave, "dtlw0", 
196     .   "LW radiation dT", "K/s",
197     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
198     .                typeval, zsto,zout)
199c
200         CALL histdef(nid_bilKPave,"duvdf",
201     .   "Boundary-layer dU","m/s2",
202     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
203     .                typeval, zsto,zout)
204c
205         CALL histdef(nid_bilKPave,"dvvdf",
206     .   "Boundary-layer dV","m/s2",
207     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
208     .                typeval, zsto,zout)
209c
210         IF (ok_orodr) THEN
211         IF (ok_orolf) THEN
212         CALL histdef(nid_bilKPave, "duoli",
213     .   "Orography dU", "m/s2",
214     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
215     .                typeval, zsto,zout)
216c
217         CALL histdef(nid_bilKPave, "dvoli", 
218     .   "Orography dV", "m/s2",
219     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
220     .                typeval, zsto,zout)
221c
222         ENDIF
223         ENDIF
224C
225         CALL histdef(nid_bilKPave, "duphy",
226     .   "Physiq dU","-",
227     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
228     .                typeval, zsto,zout)
229C
230         CALL histdef(nid_bilKPave, "dvphy",
231     .   "Physiq dV","-",
232     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
233     .                typeval, zsto,zout)
234C
235         CALL histdef(nid_bilKPave, "dtphy",
236     .   "Physiq dT","-",
237     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
238     .                typeval, zsto,zout)
239C
240         CALL histdef(nid_bilKPave, "dqphy",
241     .   "Physiq dQ","-",
242     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
243     .                typeval, zsto,zout)
244C
245         CALL histdef(nid_bilKPave, "dqlphy",
246     .   "Physiq dQl","-",
247     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
248     .                typeval, zsto,zout)
249C
250C
251         CALL histend(nid_bilKPave)
252c
253         ndex2d = 0
254         ndex3d = 0
255c
256      ENDIF ! fin de test sur ok_journe
Note: See TracBrowser for help on using the repository browser.