Context Navigation

source: LMDZ4/trunk/libf/phylmd/ini_histhf.h @ 648

Last change on this file since 648 was 644, checked in by Laurent Fairhead, 20 years ago
Synchronisation avec tous les diagnostiques de Ionela IM Inclusion du slab ocean IM LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 12.6 KB

Line
1	!
2	! $Header$
3	!
4
5	IF (ok_hf) THEN
6	c
7	zstohf = dtime * FLOAT(ecrit_hf)
8	zstophy = dtime
9	zstoday = dtime * FLOAT(ecrit_day)
10	zout = dtime * FLOAT(ecrit_hf)
11	c
12	idayref = day_ref
13	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
14
15	c
16	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
17	DO i = 1, iim
18	zx_lon(i,1) = rlon(i+1)
19	zx_lon(i,jjmp1) = rlon(i+1)
20	ENDDO
21
22	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
23	CALL histbeg("histhf", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
24	. 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
25	. nhori, nid_hf)
26
27	CALL histvert(nid_hf, "presnivs", "Vertical levels", "mb",
28	. klev, presnivs/100., nvert)
29	c
30
31	IF(type_run.EQ."CLIM".OR.type_run.EQ."ENSP") THEN
32	c
33	IF(lev_histhf.GE.1) THEN
34	c
35	CALL histdef(nid_hf, "aireTER","Grid area CONT","-",
36	. iim,jjmp1,nhori, 1,1,1, -99, 32,
37	. "once", zstohf,zout)
38	c
39	CALL histdef(nid_hf, "contfracATM","% sfce ter+lic ","-",
40	. iim,jjmp1,nhori, 1,1,1, -99, 32,
41	. "once", zstohf,zout)
42	c
43	CALL histdef(nid_hf, "contfracOR","% sfce terre OR", "-",
44	. iim,jjmp1,nhori, 1,1,1, -99, 32,
45	. "once", zstohf,zout)
46	c
47	c Champs 2D:
48	c
49	CALL histdef(nid_hf, "t2m", "Temperature 2m", "K",
50	. iim,jjmp1,nhori, 1,1,1, -99, 32,
51	. "inst(X)", zstohf,zout)
52	c
53	CALL histdef(nid_hf, "q2m", "Specific humidity", "kg/kg",
54	. iim,jjmp1,nhori, 1,1,1, -99, 32,
55	. "inst(X)", zstohf,zout)
56	c
57	CALL histdef(nid_hf, "psol", "Surface Pressure", "Pa",
58	. iim,jjmp1,nhori, 1,1,1, -99, 32,
59	. "inst(X)", zstohf,zout)
60	c
61	CALL histdef(nid_hf, "precip", "Precipitation", "kg/(s*m2)",
62	. iim,jjmp1,nhori, 1,1,1, -99, 32,
63	. "ave(X)", zstophy,zout)
64	c
65	cIM ENSEMBLES BEG
66	CALL histdef(nid_hf, "tsol", "Surface Temperature", "K",
67	. iim,jjmp1,nhori, 1,1,1, -99, 32,
68	. "inst(X)", zstohf,zout)
69	c
70	CALL histdef(nid_hf, "slp", "Sea Level Pressure", "Pa",
71	. iim,jjmp1,nhori, 1,1,1, -99, 32,
72	. "inst(X)", zstohf,zout)
73	c
74	CALL histdef(nid_hf, "u10m", "Vent zonal 10m", "m/s",
75	. iim,jjmp1,nhori, 1,1,1, -99, 32,
76	. "inst(X)", zstohf,zout)
77	c
78	CALL histdef(nid_hf, "v10m", "Vent meridien 10m", "m/s",
79	. iim,jjmp1,nhori, 1,1,1, -99, 32,
80	. "inst(X)", zstohf,zout)
81	c
82	CALL histdef(nid_hf, "wind10m","10-m wind speed","m/s",
83	. iim,jjmp1,nhori, 1,1,1, -99, 32,
84	. "inst(X)", zstohf,zout)
85	c
86	DO k=1, nlevSTD
87	c
88	c bb=clevSTD(k)
89	c
90	IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
91	IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
92	c
93	c IF(k.GE.2) THEN
94	c aa=clevSTD(k)
95	c bb=aa(1:lnblnk1(aa))
96	c ENDIF
97	c
98	IF(bb2.EQ."500") THEN
99	CALL histdef(nid_hf, "phi"//bb2,
100	. "Geopotential "//bb2//"mb", "m",
101	. iim,jjmp1,nhori, 1,1,1, -99, 32,
102	. "inst(X)", zstohf,zout)
103	ENDIF !bb2.EQ."500"
104	ENDDO !k=1, nlevSTD
105	c
106	ENDIF !lev_histhf.GE.1
107	c
108	IF(lev_histhf.GE.2) THEN
109	c
110	CALL histdef(nid_hf, "cldt", "Total cloudiness", "%",
111	. iim,jjmp1,nhori, 1,1,1, -99, 32,
112	. "ave(X)", zstophy,zout)
113	c
114	CALL histdef(nid_hf, "SWnetOR",
115	. "Sfce net SW radiation OR", "W/m2",
116	. iim,jjmp1,nhori, 1,1,1, -99, 32,
117	. "ave(X)", zstophy,zout)
118	c
119	CALL histdef(nid_hf, "SWdownOR",
120	. "Sfce incident SW radiation OR", "W/m2",
121	. iim,jjmp1,nhori, 1,1,1, -99, 32,
122	. "ave(X)", zstophy,zout)
123	c
124	CALL histdef(nid_hf, "LWdownOR",
125	. "Sfce incident LW radiation OR", "W/m2",
126	. iim,jjmp1,nhori, 1,1,1, -99, 32,
127	. "ave(X)", zstophy,zout)
128	ENDIF !lev_histhf.GE.2
129	c
130	cIM ENSEMBLES END
131	c
132	IF(lev_histhf.GE.3) THEN
133	c
134	DO k=1, nlevSTD
135	c
136	c bb=clevSTD(k)
137	c
138	c IF(k.GE.2) THEN
139	IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
140	IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
141	c aa=clevSTD(k)
142	c bb=aa(1:lnblnk1(aa))
143	c ENDIF
144	c
145	IF(bb2.EQ."850".or.bb2.EQ."700".OR.
146	$ bb2.EQ."500".or.bb2.EQ."200") THEN
147	c
148	CALL histdef(nid_hf, "t"//bb2,
149	. "Temperature "//bb2//"mb","K",
150	. iim,jjmp1,nhori, 1,1,1, -99, 32,
151	. "inst(X)", zstohf,zout)
152	c
153	IF(bb2.NE."500") THEN
154	CALL histdef(nid_hf, "phi"//bb2,
155	. "Geopotential "//bb2//"mb", "m",
156	. iim,jjmp1,nhori, 1,1,1, -99, 32,
157	. "inst(X)", zstohf,zout)
158	ENDIF
159	c
160	CALL histdef(nid_hf, "q"//bb2,
161	. "Specific humidity "//bb2//"mb","kg/kg",
162	. iim,jjmp1,nhori, 1,1,1, -99, 32,
163	. "inst(X)", zstohf,zout)
164	c
165	CALL histdef(nid_hf, "u"//bb2,
166	. "Zonal wind "//bb2//"mb","m/s",
167	. iim,jjmp1,nhori, 1,1,1, -99, 32,
168	. "inst(X)", zstohf,zout)
169	c
170	CALL histdef(nid_hf, "v"//bb2,
171	. "Meridional wind "//bb2//"mb","m/s",
172	. iim,jjmp1,nhori, 1,1,1, -99, 32,
173	. "inst(X)", zstohf,zout)
174	c
175	ENDIF !IF(bb2.EQ."850".or.bb2.EQ."700".OR.
176	c
177	ENDDO !nlevSTD
178	c
179	ENDIF !lev_histhf.GE.3
180	c
181	IF(lev_histhf.GE.4) THEN
182	c
183	cIM sorties hf 3D
184	c
185	#undef histhf3d
186	#define histhf3d
187	#ifdef histhf3d
188	#include "ini_histhf3d.h"
189	#endif
190	ENDIF !lev_histhf.GE.4
191	c
192	ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
193	c
194	IF(lev_histhf.GE.1) THEN
195	c
196	DO k=1, nlevSTD
197	c
198	c
199	c bb=clevSTD(k)
200	c
201	c IF(k.GE.2) THEN
202	IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
203	IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
204	c aa=clevSTD(k)
205	c bb=aa(1:lnblnk1(aa))
206	c ENDIF
207	c
208	IF(bb2.EQ."850".or.bb2.EQ."200") THEN
209	c
210	CALL histdef(nid_hf, "u"//bb2,
211	. "Zonal wind "//bb2//"mb","m/s",
212	. iim,jjmp1,nhori, 1,1,1, -99, 32,
213	. "inst(X)", zstohf,zstohf)
214	c
215	CALL histdef(nid_hf, "v"//bb2,
216	. "Meridional wind "//bb2//"mb","m/s",
217	. iim,jjmp1,nhori, 1,1,1, -99, 32,
218	. "inst(X)", zstohf,zstohf)
219	c
220	ENDIF !bb2.EQ."850".or.bb2.EQ."200"
221	ENDDO !k=1, nlevSTD
222	c
223	CALL histdef(nid_hf, "topl",
224	. "Outgoing longwave radiation", "W/m2",
225	. iim,jjmp1,nhori, 1,1,1, -99, 32,
226	. "ave(X)", zstophy,zout)
227	c
228	CALL histdef(nid_hf, "precip","Total precipitation rate",
229	. "kg/(s*m2)",iim,jjmp1,nhori, 1,1,1, -99, 32,
230	. "ave(X)", zstophy,zout)
231	c
232	CALL histdef(nid_hf, "slp",
233	. "Mean sea-level pressure", "Pa",
234	. iim,jjmp1,nhori, 1,1,1, -99, 32,
235	. "inst(X)", zstohf,zout)
236	c
237	cIM 280504 BEG
238	ENDIF !(lev_histhf.GE.1) THEN
239	c
240	IF(lev_histhf.GE.2) THEN
241	c
242	DO k=1, nlevSTD
243	c
244	c bb=clevSTD(k)
245	c
246	c IF(k.GE.2) THEN
247	IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
248	IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
249	c aa=clevSTD(k)
250	c bb=aa(1:lnblnk1(aa))
251	c ENDIF
252	c
253	IF(bb2.EQ."850".OR.bb2.EQ."500") THEN
254	c
255	CALL histdef(nid_hf, "t"//bb2,
256	. "Temperature "//bb2//"mb","K",
257	. iim,jjmp1,nhori, 1,1,1, -99, 32,
258	. "inst(X)", zstohf,zstohf)
259	c
260	CALL histdef(nid_hf, "q"//bb2,
261	. "Specific humidity "//bb2//"mb","kg/kg",
262	. iim,jjmp1,nhori, 1,1,1, -99, 32,
263	. "inst(X)", zstohf,zstohf)
264	c
265	ENDIF !bb2
266	c
267	IF(bb2.EQ."500") THEN
268	c
269	CALL histdef(nid_hf, "phi"//bb2,
270	. "Geopotential "//bb2//"mb", "m",
271	. iim,jjmp1,nhori, 1,1,1, -99, 32,
272	. "inst(X)", zstohf,zstohf)
273	c
274	CALL histdef(nid_hf, "w"//bb2,
275	. "Vertical motion "//bb2//"mb","Pa/s",
276	. iim,jjmp1,nhori, 1,1,1, -99, 32,
277	. "inst(X)", zstohf,zstohf)
278	ENDIF !bb2.EQ."500"
279	c
280	IF(bb3.EQ."50") THEN
281	c
282	CALL histdef(nid_hf, "t"//bb3,
283	. "Temperature "//bb3//"mb","K",
284	. iim,jjmp1,nhori, 1,1,1, -99, 32,
285	. "inst(X)", zstohf,zstohf)
286	c
287	CALL histdef(nid_hf, "u"//bb3,
288	. "Zonal wind "//bb3//"mb","m/s",
289	. iim,jjmp1,nhori, 1,1,1, -99, 32,
290	. "inst(X)", zstohf,zstohf)
291	c
292	CALL histdef(nid_hf, "v"//bb3,
293	. "Meridional wind "//bb3//"mb","m/s",
294	. iim,jjmp1,nhori, 1,1,1, -99, 32,
295	. "inst(X)", zstohf,zstohf)
296	c
297	ENDIF !bb3
298	c
299	ENDDO !k=1, nlevSTD
300	c
301	CALL histdef(nid_hf, "u10m", "Vent zonal 10m", "m/s",
302	. iim,jjmp1,nhori, 1,1,1, -99, 32,
303	. "inst(X)", zstohf,zout)
304	c
305	CALL histdef(nid_hf, "v10m", "Vent meridien 10m", "m/s",
306	. iim,jjmp1,nhori, 1,1,1, -99, 32,
307	. "inst(X)", zstohf,zout)
308	c
309	CALL histdef(nid_hf, "t2m", "Temperature 2m", "K",
310	. iim,jjmp1,nhori, 1,1,1, -99, 32,
311	. "inst(X)", zstohf,zout)
312	c
313	CALL histdef(nid_hf, "q2m", "Specific humidity", "kg/kg",
314	. iim,jjmp1,nhori, 1,1,1, -99, 32,
315	. "inst(X)", zstohf,zout)
316	c
317	CALL histdef(nid_hf, "psol", "Surface Pressure", "Pa",
318	. iim,jjmp1,nhori, 1,1,1, -99, 32,
319	. "inst(X)", zstohf,zout)
320	c
321	CALL histdef(nid_hf, "cldt", "Total cloudiness", "%",
322	. iim,jjmp1,nhori, 1,1,1, -99, 32,
323	. "ave(X)", zstohf,zout)
324	c
325	CALL histdef(nid_hf, "taux",
326	$ "Zonal wind stress","Pa",
327	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
328	$ "ave(X)", zstohf,zout)
329	c
330	CALL histdef(nid_hf, "tauy",
331	$ "Meridional wind stress ", "Pa",
332	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
333	$ "ave(X)", zstohf,zout)
334	c
335	CALL histdef(nid_hf, "prw", "Precipitable water", "kg/m2",
336	. iim,jjmp1,nhori, 1,1,1, -99, 32,
337	cIM 090904 . "ave(X)", zstohf,zout)
338	. "ave(X)", zstophy,zout)
339	c
340	CALL histdef(nid_hf, "sens", "Sensible heat flux", "W/m2",
341	. iim,jjmp1,nhori, 1,1,1, -99, 32,
342	cIM 090904 . "ave(X)", zstohf,zout)
343	. "ave(X)", zstophy,zout)
344	c
345	CALL histdef(nid_hf, "flat", "Latent heat flux", "W/m2",
346	. iim,jjmp1,nhori, 1,1,1, -99, 32,
347	. "ave(X)", zstohf,zout)
348	c
349	c CALL histdef(nid_hf, "snow_mass", "Snow Mass", "kg/m2",
350	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
351	c . "ave(X)", zstohf,zout)
352	c
353	CALL histdef(nid_hf, "SWdnSFC", "SWdn at surface","W/m2",
354	. iim,jjmp1,nhori, 1,1,1,-99,
355	. 32, "ave(X)", zstohf,zout)
356	c
357	CALL histdef(nid_hf, "SWupSFC", "SWup at surface","W/m2",
358	. iim,jjmp1,nhori, 1,1,1,-99,
359	. 32, "ave(X)", zstohf,zout)
360	c
361	CALL histdef(nid_hf,"LWdnSFC","Down. IR rad. at surface",
362	. "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
363	. "ave(X)", zstohf,zout)
364	c
365	CALL histdef(nid_hf,"LWupSFC","Upwd. IR rad. at surface",
366	. "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
367	. "ave(X)", zstohf,zout)
368	c
369	CALL histdef(nid_hf, "SWdnTOA", "SWdn at TOA","W/m2",
370	. iim,jjmp1,nhori, 1,1,1,-99,
371	. 32, "ave(X)", zstohf,zout)
372	c
373	CALL histdef(nid_hf, "SWupTOA", "SWup at TOA","W/m2",
374	. iim,jjmp1,nhori, 1,1,1,-99,
375	. 32, "ave(X)", zstohf,zout)
376	c
377	DO k=1, nbteta
378	CALL histdef(nid_hf, "PV"//ctetaSTD(k),
379	$ "Vorticite potentielle "//ctetaSTD(k)//"K",
380	$ "1/(Pa*s)",iim,jjmp1,nhori, 1,1,1, -99,
381	$ 32, "inst(X)", zstohf,zout)
382	ENDDO !k=1, nbteta
383	c
384	c planetary boundary layer height
385	c
386	c ENDIF !1.EQ.0
387	c
388	ENDIF !(lev_histhf.GE.2) THEN
389	cIM 280504 END
390	c
391	ENDIF !type_run
392	c
393	CALL histend(nid_hf)
394	c
395	ENDIF ! ok_hf

Note: See TracBrowser for help on using the repository browser.

Download in other formats: