1 | c |
---|
2 | c $Id: ini_bilKP_ins.h 1403 2010-07-01 09:02:53Z acaubel $ |
---|
3 | c |
---|
4 | IF (ok_journe) THEN |
---|
5 | c |
---|
6 | zsto = dtime |
---|
7 | zout = dtime |
---|
8 | typeval=tinst |
---|
9 | c |
---|
10 | idayref = day_ref |
---|
11 | CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) |
---|
12 | c |
---|
13 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) |
---|
14 | cym DO i = 1, iim |
---|
15 | cym zx_lon(i,1) = rlon(i+1) |
---|
16 | cym zx_lon(i,jjmp1) = rlon(i+1) |
---|
17 | cym ENDDO |
---|
18 | DO ll=1,klev |
---|
19 | znivsig(ll)=REAL(ll) |
---|
20 | ENDDO |
---|
21 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) |
---|
22 | cym write(*,*)'zx_lon = ',zx_lon(:,1) |
---|
23 | cym write(*,*)'zx_lat = ',zx_lat(1,:) |
---|
24 | c |
---|
25 | cIM 280405 BEG |
---|
26 | c |
---|
27 | cIM cf. AM 081204 BEG region |
---|
28 | imin_ins=1 |
---|
29 | imax_ins=iim |
---|
30 | jmin_ins=1 |
---|
31 | jmax_ins=jjmp1 |
---|
32 | cym do i=1,iim-1 |
---|
33 | cym if(zx_lon(i,1).lt.lonmin_ins) imin_ins=i |
---|
34 | cym if(zx_lon(i,1).le.lonmax_ins) imax_ins=i+1 |
---|
35 | cym enddo |
---|
36 | cym do j=1,jjmp1 |
---|
37 | cym if(zx_lat(1,j).ge.latmin_ins) jmax_ins=j |
---|
38 | cym if(zx_lat(1,j).gt.latmax_ins) jmin_ins=j |
---|
39 | cym enddo |
---|
40 | c |
---|
41 | print*,'On stoke le fichier bilKP instantanne sur ', |
---|
42 | s imin_ins,imax_ins,jmin_ins,jmax_ins |
---|
43 | print*,'On stoke le fichier bilKP instantanne sur ', |
---|
44 | s zx_lon(imin_ins,1),zx_lon(imax_ins,1), |
---|
45 | s zx_lat(1,jmin_ins),zx_lat(1,jmax_ins) |
---|
46 | cIM cf. AM 081204 END region |
---|
47 | c |
---|
48 | cIM 280405 END |
---|
49 | c |
---|
50 | cym IF(1.EQ.0) THEN |
---|
51 | cym CALL histbeg("histbilKP_ins", iim,zx_lon(:,1), jjmp1, |
---|
52 | cym . zx_lat(1,:), |
---|
53 | cym . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, |
---|
54 | cym . nhori, nid_bilKPins) |
---|
55 | ENDIF |
---|
56 | c |
---|
57 | cIM 280405 BEG |
---|
58 | c |
---|
59 | cIM cf. AM 081204 BEG region |
---|
60 | cym CALL histbeg("histbilKP_ins", iim,zx_lon(:,1), |
---|
61 | cym . jjmp1,zx_lat(1,:), |
---|
62 | cym . imin_ins,imax_ins-imin_ins+1, |
---|
63 | cym . jmin_ins,jmax_ins-jmin_ins+1, |
---|
64 | cym . itau_phy, zjulian, dtime, |
---|
65 | cym . nhori, nid_bilKPins) |
---|
66 | CALL histbeg_phy("histbilKP_ins", itau_phy, zjulian, dtime, |
---|
67 | . nhori, nid_bilKPins) |
---|
68 | cIM 081204 END |
---|
69 | c |
---|
70 | cIM 280405 END |
---|
71 | c |
---|
72 | write(*,*)'Journee ', itau_phy, zjulian |
---|
73 | CALL histvert(nid_bilKPins, "presnivs", |
---|
74 | . "Vertical levels","mb", |
---|
75 | . klev, presnivs/100., nvert) |
---|
76 | c |
---|
77 | c Champs 3D: |
---|
78 | c |
---|
79 | CALL histdef(nid_bilKPins,"ue", |
---|
80 | . "Zonal energy transport","-", |
---|
81 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
82 | . typeval, zsto,zout) |
---|
83 | c |
---|
84 | CALL histdef(nid_bilKPins,"ve", |
---|
85 | . "Merid energy transport","-", |
---|
86 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
87 | . typeval, zsto,zout) |
---|
88 | c |
---|
89 | CALL histdef(nid_bilKPins,"uq", |
---|
90 | . "Zonal humidity transport","-", |
---|
91 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
92 | . typeval, zsto,zout) |
---|
93 | c |
---|
94 | CALL histdef(nid_bilKPins,"vq", |
---|
95 | . "Merid humidity transport","-", |
---|
96 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
97 | . typeval, zsto,zout) |
---|
98 | c |
---|
99 | c Champs 3D: |
---|
100 | c |
---|
101 | CALL histdef(nid_bilKPins, "temp", |
---|
102 | . "Air temperature", "K", |
---|
103 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
104 | . typeval, zsto,zout) |
---|
105 | c |
---|
106 | CALL histdef(nid_bilKPins,"ovap", |
---|
107 | . "Specific humidity","Kg/Kg", |
---|
108 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
109 | . typeval, zsto,zout) |
---|
110 | c |
---|
111 | CALL histdef(nid_bilKPins,"geop", |
---|
112 | . "Geopotential height", "m", |
---|
113 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
114 | . typeval, zsto,zout) |
---|
115 | c |
---|
116 | CALL histdef(nid_bilKPins,"vitu", |
---|
117 | . "Zonal wind", "m/s", |
---|
118 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
119 | . typeval, zsto,zout) |
---|
120 | c |
---|
121 | CALL histdef(nid_bilKPins,"vitv", |
---|
122 | . "Meridional wind", "m/s", |
---|
123 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
124 | . typeval, zsto,zout) |
---|
125 | c |
---|
126 | CALL histdef(nid_bilKPins, "vitw", |
---|
127 | . "Vertical wind", "m/s", |
---|
128 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
129 | . typeval, zsto,zout) |
---|
130 | c |
---|
131 | CALL histdef(nid_bilKPins, "pres", |
---|
132 | . "Inter-Layer Air pressure", |
---|
133 | . "Pa", |
---|
134 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
135 | . typeval, zsto,zout) |
---|
136 | c |
---|
137 | CALL histdef(nid_bilKPins, "play", |
---|
138 | . "Mean-Layer Air pressure", |
---|
139 | . "Pa", |
---|
140 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
141 | . typeval, zsto,zout) |
---|
142 | c |
---|
143 | CALL histdef(nid_bilKPins, "oliq", |
---|
144 | . "Liquid water content", |
---|
145 | . "kg/kg", |
---|
146 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
147 | . typeval, zsto,zout) |
---|
148 | c |
---|
149 | CALL histdef(nid_bilKPins, "dtdyn", |
---|
150 | . "Dynamics dT", "K/s", |
---|
151 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
152 | . typeval, zsto,zout) |
---|
153 | c |
---|
154 | CALL histdef(nid_bilKPins, "dqdyn", |
---|
155 | . "Dynamics dQ", "Kg/Kg/s", |
---|
156 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
157 | . typeval, zsto,zout) |
---|
158 | c |
---|
159 | CALL histdef(nid_bilKPins, "dtcon", |
---|
160 | . "Convection dT", "K/s", |
---|
161 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
162 | . typeval, zsto,zout) |
---|
163 | c |
---|
164 | CALL histdef(nid_bilKPins, "ducon", |
---|
165 | . "Convection du", "m/s2", |
---|
166 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
167 | . typeval, zsto,zout) |
---|
168 | c |
---|
169 | CALL histdef(nid_bilKPins, "dvcon", |
---|
170 | . "Convection dv", "m/s2", |
---|
171 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
172 | . typeval, zsto,zout) |
---|
173 | c |
---|
174 | CALL histdef(nid_bilKPins,"dqcon", |
---|
175 | . "Convection dQ","Kg/Kg/s", |
---|
176 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
177 | . typeval, zsto,zout) |
---|
178 | c |
---|
179 | CALL histdef(nid_bilKPins, "dtlsc", |
---|
180 | . "Condensation dT", "K/s", |
---|
181 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
182 | . typeval, zsto,zout) |
---|
183 | c |
---|
184 | CALL histdef(nid_bilKPins,"dqlsc", |
---|
185 | . "Condensation dQ","Kg/Kg/s", |
---|
186 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
187 | . typeval, zsto,zout) |
---|
188 | c |
---|
189 | CALL histdef(nid_bilKPins,"dtvdf", |
---|
190 | . "Boundary-layer dT","K/s", |
---|
191 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
192 | . typeval, zsto,zout) |
---|
193 | c |
---|
194 | CALL histdef(nid_bilKPins, "dqvdf", |
---|
195 | . "Boundary-layer dQ", |
---|
196 | . "Kg/Kg/s", |
---|
197 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
198 | . typeval, zsto,zout) |
---|
199 | c |
---|
200 | CALL histdef(nid_bilKPins,"dtajs", |
---|
201 | . "Ajustement sec dT","K/s", |
---|
202 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
203 | . typeval, zsto,zout) |
---|
204 | c |
---|
205 | CALL histdef(nid_bilKPins,"dqajs", |
---|
206 | . "Ajustement sec dQ", |
---|
207 | . "Kg/Kg/s", |
---|
208 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
209 | . typeval, zsto,zout) |
---|
210 | c |
---|
211 | CALL histdef(nid_bilKPins,"dteva", |
---|
212 | . "Reevaporation dT","K/s", |
---|
213 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
214 | . typeval, zsto,zout) |
---|
215 | c |
---|
216 | CALL histdef(nid_bilKPins,"dqeva", |
---|
217 | . "Reevaporation dQ", |
---|
218 | . "Kg/Kg/s", |
---|
219 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
220 | . typeval, zsto,zout) |
---|
221 | |
---|
222 | c |
---|
223 | CALL histdef(nid_bilKPins, "dtswr", |
---|
224 | . "SW radiation dT", "K/s", |
---|
225 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
226 | . typeval, zsto,zout) |
---|
227 | c |
---|
228 | CALL histdef(nid_bilKPins, "dtsw0", |
---|
229 | . "SW radiation dT", "K/s", |
---|
230 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
231 | . typeval, zsto,zout) |
---|
232 | c |
---|
233 | CALL histdef(nid_bilKPins, "dtlwr", |
---|
234 | . "LW radiation dT", "K/s", |
---|
235 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
236 | . typeval, zsto,zout) |
---|
237 | c |
---|
238 | CALL histdef(nid_bilKPins, "dtlw0", |
---|
239 | . "LW radiation dT", "K/s", |
---|
240 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
241 | . typeval, zsto,zout) |
---|
242 | c |
---|
243 | CALL histdef(nid_bilKPins,"duvdf", |
---|
244 | . "Boundary-layer dU","m/s2", |
---|
245 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
246 | . typeval, zsto,zout) |
---|
247 | c |
---|
248 | CALL histdef(nid_bilKPins,"dvvdf", |
---|
249 | . "Boundary-layer dV","m/s2", |
---|
250 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
251 | . typeval, zsto,zout) |
---|
252 | c |
---|
253 | IF (ok_orodr) THEN |
---|
254 | IF (ok_orolf) THEN |
---|
255 | CALL histdef(nid_bilKPins, "duoli", |
---|
256 | . "Orography dU", "m/s2", |
---|
257 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
258 | . typeval, zsto,zout) |
---|
259 | c |
---|
260 | CALL histdef(nid_bilKPins, "dvoli", |
---|
261 | . "Orography dV", "m/s2", |
---|
262 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
263 | . typeval, zsto,zout) |
---|
264 | c |
---|
265 | ENDIF |
---|
266 | ENDIF |
---|
267 | C |
---|
268 | CALL histdef(nid_bilKPins, "duphy", |
---|
269 | . "Physiq dU","-", |
---|
270 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
271 | . typeval, zsto,zout) |
---|
272 | C |
---|
273 | CALL histdef(nid_bilKPins, "dvphy", |
---|
274 | . "Physiq dV","-", |
---|
275 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
276 | . typeval, zsto,zout) |
---|
277 | C |
---|
278 | CALL histdef(nid_bilKPins, "dtphy", |
---|
279 | . "Physiq dT","-", |
---|
280 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
281 | . typeval, zsto,zout) |
---|
282 | C |
---|
283 | CALL histdef(nid_bilKPins, "dqphy", |
---|
284 | . "Physiq dQ","-", |
---|
285 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
286 | . typeval, zsto,zout) |
---|
287 | C |
---|
288 | CALL histdef(nid_bilKPins, "dqlphy", |
---|
289 | . "Physiq dQl","-", |
---|
290 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
291 | . typeval, zsto,zout) |
---|
292 | cIM 280405 BEG |
---|
293 | c |
---|
294 | c Champs 2D: |
---|
295 | c |
---|
296 | c u850, v850 |
---|
297 | c DO k=1, nlevSTD |
---|
298 | DO k=1, 12 |
---|
299 | c |
---|
300 | IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k) |
---|
301 | IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k) |
---|
302 | c |
---|
303 | IF(bb2.EQ."850") THEN |
---|
304 | c |
---|
305 | CALL histdef(nid_bilKPins, "u"//bb2, |
---|
306 | . "Zonal wind "//bb2//"mb","m/s", |
---|
307 | . iim,jjphy_nb,nhori, 1,1,1, -99, 32, |
---|
308 | . typeval, zsto,zout) |
---|
309 | c |
---|
310 | CALL histdef(nid_bilKPins, "v"//bb2, |
---|
311 | . "Meridional wind "//bb2//"mb","m/s", |
---|
312 | . iim,jjphy_nb,nhori, 1,1,1, -99, 32, |
---|
313 | . typeval, zsto,zout) |
---|
314 | c |
---|
315 | ENDIF !(bb2.EQ."850") |
---|
316 | c |
---|
317 | ENDDO !k=1, 12 |
---|
318 | c |
---|
319 | cIM 280405 END |
---|
320 | c |
---|
321 | CALL histend(nid_bilKPins) |
---|
322 | c |
---|
323 | ndex2d = 0 |
---|
324 | ndex3d = 0 |
---|
325 | c |
---|
326 | ENDIF ! fin de test sur ok_journe |
---|