My Project
 All Classes Files Functions Variables Macros
ini_bilKP_ins.h
Go to the documentation of this file.
1 c
2 c $Id: ini_bilKP_ins.h 1403 2010-07-01 09:02:53Z fairhead $
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