Context Navigation

source: LMDZ4/trunk/libf/phylmd/ini_histmth.h @ 900

Last change on this file since 900 was 895, checked in by Laurent Fairhead, 17 years ago
Un peu de menage sur les fichiers de sortie apres les tests aquaplanete FH LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 74.8 KB

Line
1	c
2	c $Header$
3	c
4	IF (ok_mensuel) THEN
5	c$OMP MASTER
6	c
7	zstophy = dtime
8	zstomth = ecrit_mth !valeurs mensuelles "instantannees"
9	zstorad = dtime * FLOAT(radpas) !variables stockees toutes les 2h
10	zstohf = ecrit_hf !variables stockees toutes les 6h
11	zout = ecrit_mth !frequence des sorties : 1 val. par mois
12	c
13	idayref = day_ref
14	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
15	c
16	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
17	cym DO i = 1, iim
18	cym zx_lon(i,1) = rlon(i+1)
19	cym zx_lon(i,jjmp1) = rlon(i+1)
20	cym ENDDO
21	DO ll=1,klev
22	znivsig(ll)=float(ll)
23	ENDDO
24	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
25	cym CALL histbeg("histmth.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
26	cym . 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
27	cym . nhori, nid_mth)
28
29	CALL histbeg_phy("histmth",itau_phy, zjulian, dtime,
30	. nhori, nid_mth)
31
32	write(,)'Mensuel ', itau_phy, zjulian
33	CALL histvert(nid_mth, "presnivs", "Vertical levels", "mb",
34	. klev, presnivs/100., nvert)
35	c
36	! CALL histdef(nid_mth, "dtthe", "Dry adjust. dT", "K/s",
37	! . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
38	! . "ave(X)", zstophy,zout)
39	!
40	! CALL histdef(nid_mth, "dtcon", "Convection dT", "K/s",
41	! . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
42	! . "ave(X)", zstophy,zout)
43	c
44	IF(type_run.EQ."CLIM".OR.type_run.EQ."ENSP") THEN
45	c
46	IF(lev_histmth.GE.1) THEN
47	c
48	CALL histdef(nid_mth, "phis",
49	. "Surface geop. height", "m2/s2",
50	. iim,jj_nb,nhori, 1,1,1, -99, 32,
51	. "once", zstophy,zout)
52	c
53	CALL histdef(nid_mth, "aire", "Grid area", "-",
54	. iim,jj_nb,nhori, 1,1,1, -99, 32,
55	. "once", zstophy,zout)
56	c
57	CALL histdef(nid_mth, "pourc_"//clnsurf(is_ter),
58	. "% "//clnsurf(is_ter), "%",
59	. iim,jj_nb,nhori, 1,1,1, -99, 32,
60	. "once", zstophy,zout)
61	c
62	CALL histdef(nid_mth, "fract_"//clnsurf(is_ter),
63	. "Fraction "//clnsurf(is_ter), "1",
64	. iim,jj_nb,nhori, 1,1,1, -99, 32,
65	. "once", zstophy,zout)
66	c
67	c Champs 2D:
68	c
69	CALL histdef(nid_mth, "slp", "Sea Level Pressure", "Pa",
70	. iim,jj_nb,nhori, 1,1,1, -99, 32,
71	. "ave(X)", zstophy,zout)
72	c
73	CALL histdef(nid_mth, "tsol", "Surface Temperature", "K",
74	. iim,jj_nb,nhori, 1,1,1, -99, 32,
75	. "ave(X)", zstophy,zout)
76	c
77	CALL histdef(nid_mth, "t2m", "Temperature 2m", "K",
78	. iim,jj_nb,nhori, 1,1,1, -99, 32,
79	. "ave(X)", zstophy,zout)
80	c
81	c ENSEMBLES BEG
82	c CALL histdef(nid_mth, "t2m_min", "Temp. 2m min.",
83	c . "K",
84	c . iim,jj_nb,nhori, 1,1,1, -99, 32,
85	c . t2mincels, zstophy,zout)
86	c
87	c CALL histdef(nid_mth, "t2m_max", "Temp. 2m max.",
88	c . "K",
89	c . iim,jj_nb,nhori, 1,1,1, -99, 32,
90	c . t2maxcels, zstophy,zout)
91	c
92	c CALL histdef(nid_mth, "tsoil", "Sfce soil Temperature",
93	c . "K",
94	c . iim,jj_nb,nhori, 1,1,1, -99, 32,
95	c . "ave(X)", zstophy,zout)
96	c
97	CALL histdef(nid_mth, "wind10m","10-m wind speed","m/s",
98	. iim,jj_nb,nhori, 1,1,1, -99, 32,
99	. "ave(X)", zstophy,zout)
100	c
101	CALL histdef(nid_mth, "sicf", "Sea-ice fraction", "-",
102	. iim,jj_nb,nhori, 1,1,1, -99, 32,
103	. "ave(X)", zstophy,zout)
104	c
105	c ENSEMBLES END
106	CALL histdef(nid_mth, "q2m", "Specific humidity 2m", "kg/kg",
107	. iim,jj_nb,nhori, 1,1,1, -99, 32,
108	. "ave(X)", zstophy,zout)
109	c
110	CALL histdef(nid_mth, "u10m", "Vent zonal 10m", "m/s",
111	. iim,jj_nb,nhori, 1,1,1, -99, 32,
112	. "ave(X)", zstophy,zout)
113	c
114	CALL histdef(nid_mth, "v10m", "Vent meridien 10m", "m/s",
115	. iim,jj_nb,nhori, 1,1,1, -99, 32,
116	. "ave(X)", zstophy,zout)
117	c
118	CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa",
119	. iim,jj_nb,nhori, 1,1,1, -99, 32,
120	. "ave(X)", zstophy,zout)
121	c
122	CALL histdef(nid_mth, "qsurf", "Surface Air humidity", "kg/kg",
123	. iim,jj_nb,nhori, 1,1,1, -99, 32,
124	. "ave(X)", zstophy,zout)
125	c
126	if (.not. ok_veget) then
127	CALL histdef(nid_mth, "qsol", "Soil watter content", "mm",
128	. iim,jj_nb,nhori, 1,1,1, -99, 32,
129	. "ave(X)", zstophy,zout)
130	endif
131	c
132	c CALL histdef(nid_mth, "ndayrain",
133	c . "Number of day with rain (liq+sol)", "-",
134	c . iim,jj_nb,nhori, 1,1,1, -99, 32,
135	c . "inst(X)", zstomth,zout)
136	c
137	CALL histdef(nid_mth, "precip", "Precipitation Totale liq+sol",
138	. "kg/(s*m2)",
139	. iim,jj_nb,nhori, 1,1,1, -99, 32,
140	. "ave(X)", zstophy,zout)
141	c
142	CALL histdef(nid_mth, "plul", "Large-scale Precip.",
143	. "kg/(s*m2)",
144	. iim,jj_nb,nhori, 1,1,1, -99, 32,
145	. "ave(X)", zstophy,zout)
146	c
147	CALL histdef(nid_mth, "pluc", "Convective Precip.",
148	. "kg/(s*m2)",
149	. iim,jj_nb,nhori, 1,1,1, -99, 32,
150	. "ave(X)", zstophy,zout)
151	c
152	CALL histdef(nid_mth, "snow", "Snow fall", "kg/(s*m2)",
153	. iim,jj_nb,nhori, 1,1,1, -99, 32,
154	. "ave(X)", zstophy,zout)
155	c
156	c CALL histdef(nid_mth, "snow_mass", "Snow Mass", "kg/m2",
157	c . iim,jj_nb,nhori, 1,1,1, -99, 32,
158	c . "ave(X)", zstophy,zout)
159	c
160	CALL histdef(nid_mth, "evap", "Evaporation", "kg/(s*m2)",
161	. iim,jj_nb,nhori, 1,1,1, -99, 32,
162	. "ave(X)", zstophy,zout)
163	c
164	CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2",
165	. iim,jj_nb,nhori, 1,1,1, -99, 32,
166	. "ave(X)", zstophy,zout)
167	c
168	CALL histdef(nid_mth, "tops0", "CS Solar rad. at TOA", "W/m2",
169	. iim,jj_nb,nhori, 1,1,1, -99, 32,
170	. "ave(X)", zstophy,zout)
171	c
172	CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2",
173	. iim,jj_nb,nhori, 1,1,1, -99, 32,
174	. "ave(X)", zstophy,zout)
175	c
176	CALL histdef(nid_mth, "topl0", "IR rad. at TOA", "W/m2",
177	. iim,jj_nb,nhori, 1,1,1, -99, 32,
178	. "ave(X)", zstophy,zout)
179	c
180	CALL histdef(nid_mth, "SWupTOA", "SWup at TOA","W/m2",
181	. iim,jj_nb,nhori, 1,1,1,-99,
182	. 32, "ave(X)", zstophy,zout)
183	c
184	CALL histdef(nid_mth, "SWupTOAclr",
185	. "SWup clear sky at TOA","W/m2",
186	. iim,jj_nb,nhori, 1,1,1,-99,
187	. 32, "ave(X)", zstophy,zout)
188	c
189	CALL histdef(nid_mth, "SWdnTOA", "SWdn at TOA","W/m2",
190	. iim,jj_nb,nhori, 1,1,1,-99,
191	. 32, "ave(X)", zstophy,zout)
192	c
193	CALL histdef(nid_mth, "SWdnTOAclr",
194	. "SWdn clear sky at TOA","W/m2",
195	. iim,jj_nb,nhori, 1,1,1,-99,
196	. 32, "ave(X)", zstophy,zout)
197	c
198	CALL histdef(nid_mth, "SWup200", "SWup at 200mb","W/m2",
199	. iim,jj_nb,nhori, 1,1,1,-99,
200	. 32, "ave(X)", zstophy,zout)
201	c
202	CALL histdef(nid_mth, "SWup200clr",
203	. "SWup clear sky at 200mb","W/m2",
204	. iim,jj_nb,nhori, 1,1,1,-99,
205	. 32, "ave(X)", zstophy,zout)
206	c
207	CALL histdef(nid_mth, "SWdn200", "SWdn at 200mb","W/m2",
208	. iim,jj_nb,nhori, 1,1,1,-99,
209	. 32, "ave(X)", zstophy,zout)
210	c
211	CALL histdef(nid_mth, "SWdn200clr",
212	. "SWdn clear sky at 200mb","W/m2",
213	. iim,jj_nb,nhori, 1,1,1,-99,
214	. 32, "ave(X)", zstophy,zout)
215	c
216	CALL histdef(nid_mth, "LWup200", "LWup at 200mb","W/m2",
217	. iim,jj_nb,nhori, 1,1,1,-99,
218	. 32, "ave(X)", zstophy,zout)
219	c
220	CALL histdef(nid_mth, "LWup200clr",
221	. "LWup clear sky at 200mb","W/m2",
222	. iim,jj_nb,nhori, 1,1,1,-99,
223	. 32, "ave(X)", zstophy,zout)
224	c
225	CALL histdef(nid_mth, "LWdn200", "LWdn at 200mb","W/m2",
226	. iim,jj_nb,nhori, 1,1,1,-99,
227	. 32, "ave(X)", zstophy,zout)
228	c
229	CALL histdef(nid_mth, "LWdn200clr",
230	. "LWdn clear sky at 200mb","W/m2",
231	. iim,jj_nb,nhori, 1,1,1,-99,
232	. 32, "ave(X)", zstophy,zout)
233	c
234	CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2",
235	. iim,jj_nb,nhori, 1,1,1, -99, 32,
236	. "ave(X)", zstophy,zout)
237	c
238	CALL histdef(nid_mth, "sols0", "Solar rad. at surf.", "W/m2",
239	. iim,jj_nb,nhori, 1,1,1, -99, 32,
240	. "ave(X)", zstophy,zout)
241	c
242	CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2",
243	. iim,jj_nb,nhori, 1,1,1, -99, 32,
244	. "ave(X)", zstophy,zout)
245	c
246	CALL histdef(nid_mth, "radsol", "Rayonnement au sol", "W/m2",
247	. iim,jj_nb,nhori, 1,1,1, -99, 32,
248	. "ave(X)", zstophy,zout)
249	c
250	CALL histdef(nid_mth, "soll0", "IR rad. at surface", "W/m2",
251	. iim,jj_nb,nhori, 1,1,1, -99, 32,
252	. "ave(X)", zstophy,zout)
253	c
254	CALL histdef(nid_mth, "SWupSFC", "SWup at surface","W/m2",
255	. iim,jj_nb,nhori, 1,1,1,-99,
256	. 32, "ave(X)", zstophy,zout)
257	c
258	CALL histdef(nid_mth, "SWupSFCclr",
259	. "SWup clear sky at surface","W/m2",
260	. iim,jj_nb,nhori, 1,1,1,-99,
261	. 32, "ave(X)", zstophy,zout)
262	c
263	CALL histdef(nid_mth, "SWdnSFC", "SWdn at surface","W/m2",
264	. iim,jj_nb,nhori, 1,1,1,-99,
265	. 32, "ave(X)", zstophy,zout)
266	c
267	CALL histdef(nid_mth, "SWdnSFCclr",
268	. "SWdn clear sky at surface","W/m2",
269	. iim,jj_nb,nhori, 1,1,1,-99,
270	. 32, "ave(X)", zstophy,zout)
271	c
272	CALL histdef(nid_mth,"LWupSFC","Upwd. IR rad. at surface",
273	. "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32,
274	. "ave(X)", zstophy,zout)
275	c
276	CALL histdef(nid_mth,"LWdnSFC","Down. IR rad. at surface",
277	. "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32,
278	. "ave(X)", zstophy,zout)
279	c
280	CALL histdef(nid_mth,"LWupSFCclr",
281	. "CS Upwd. IR rad. at surface",
282	. "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32,
283	. "ave(X)", zstophy,zout)
284	c
285	CALL histdef(nid_mth,"LWdnSFCclr",
286	. "Down. CS IR rad. at surface",
287	. "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32,
288	. "ave(X)", zstophy,zout)
289	c
290	CALL histdef(nid_mth, "bils", "Surf. total heat flux", "W/m2",
291	. iim,jj_nb,nhori, 1,1,1, -99, 32,
292	. "ave(X)", zstophy,zout)
293	c
294	CALL histdef(nid_mth, "sens", "Sensible heat flux", "W/m2",
295	. iim,jj_nb,nhori, 1,1,1, -99, 32,
296	. "ave(X)", zstophy,zout)
297	c
298	CALL histdef(nid_mth, "fder", "Heat flux derivation", "W/m2",
299	. iim,jj_nb,nhori, 1,1,1, -99, 32,
300	. "ave(X)", zstophy,zout)
301	c
302	c CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa",
303	c . iim,jj_nb,nhori, 1,1,1, -99, 32,
304	c . "ave(X)", zstophy,zout)
305	c
306	c CALL histdef(nid_mth, "frtv", "Meridional wind stress", "Pa",
307	c . iim,jj_nb,nhori, 1,1,1, -99, 32,
308	c . "ave(X)", zstophy,zout)
309	c
310	CALL histdef(nid_mth, "ffonte","Thermal flux for snow melting",
311	. "W/m2",iim,jj_nb,nhori, 1,1,1, -99, 32,
312	. "ave(X)", zstophy,zout)
313	c
314	CALL histdef(nid_mth, "fqcalving","Ice Calving",
315	. "kg/m2/s",iim,jj_nb,nhori, 1,1,1, -99, 32,
316	. "ave(X)", zstophy,zout)
317	c
318	CALL histdef(nid_mth, "fqfonte","Land ice melt",
319	. "kg/m2/s",iim,jj_nb,nhori, 1,1,1, -99, 32,
320	. "ave(X)", zstophy,zout)
321
322	DO nsrf = 1, nbsrf
323	call histdef(nid_mth, "taux_"//clnsurf(nsrf),
324	$ "Zonal wind stress"//clnsurf(nsrf), "Pa",
325	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
326	$ "ave(X)", zstophy,zout)
327
328	call histdef(nid_mth, "tauy_"//clnsurf(nsrf),
329	$ "Meridional wind stress "//clnsurf(nsrf), "Pa",
330	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
331	$ "ave(X)", zstophy,zout)
332	ENDDO
333	c
334	DO nsrf = 1, nbsrf
335	C
336	IF(nsrf.GE.2) THEN
337	call histdef(nid_mth, "pourc_"//clnsurf(nsrf),
338	$ "% "//clnsurf(nsrf), "%",
339	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
340	$ "ave(X)", zstophy,zout)
341	c
342	call histdef(nid_mth, "fract_"//clnsurf(nsrf),
343	$ "Fraction "//clnsurf(nsrf), "1",
344	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
345	$ "ave(X)", zstophy,zout)
346	ENDIF !nsrf.GT.2
347	C
348	call histdef(nid_mth, "tsol_"//clnsurf(nsrf),
349	$ "Temperature "//clnsurf(nsrf), "K",
350	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
351	$ "ave(X)", zstophy,zout)
352	C
353	call histdef(nid_mth, "sens_"//clnsurf(nsrf),
354	$ "Sensible heat flux "//clnsurf(nsrf), "W/m2",
355	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
356	$ "ave(X)", zstophy,zout)
357	c
358	call histdef(nid_mth, "lat_"//clnsurf(nsrf),
359	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
360	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
361	$ "ave(X)", zstophy,zout)
362	C
363	call histdef(nid_mth, "flw_"//clnsurf(nsrf),
364	$ "LW "//clnsurf(nsrf), "W/m2",
365	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
366	$ "ave(X)", zstophy,zout)
367	c
368	call histdef(nid_mth, "fsw_"//clnsurf(nsrf),
369	$ "SW "//clnsurf(nsrf), "W/m2",
370	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
371	$ "ave(X)", zstophy,zout)
372	C
373	call histdef(nid_mth, "wbils_"//clnsurf(nsrf),
374	$ "Bilan sol "//clnsurf(nsrf), "W/m2",
375	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
376	$ "ave(X)", zstophy,zout)
377	C
378	call histdef(nid_mth, "wbilo_"//clnsurf(nsrf),
379	$ "Bilan eau "//clnsurf(nsrf), "kg/(m2*s)",
380	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
381	$ "ave(X)", zstophy,zout)
382	C
383	cIM 250105 BEG
384	IF(1.EQ.0) THEN
385	C
386	call histdef(nid_mth, "evapour_"//clnsurf(nsrf),
387	$ "Evappourc "//clnsurf(nsrf), "kg/(sm2)",
388	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
389	$ "ave(X)", zstophy,zout)
390	C
391	call histdef(nid_mth, "prepour_"//clnsurf(nsrf),
392	$ "Precippourc "//clnsurf(nsrf), "kg/(sm2)",
393	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
394	$ "ave(X)", zstophy,zout)
395	C
396	ENDIF !(1.EQ.0) THEN
397	cIM 250105 END
398	C
399	END DO
400	c
401	CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
402	. iim,jj_nb,nhori, 1,1,1, -99, 32,
403	. "ave(X)", zstophy,zout)
404	c
405	CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
406	. iim,jj_nb,nhori, 1,1,1, -99, 32,
407	. "ave(X)", zstophy,zout)
408	c
409	CALL histdef(nid_mth, "cldl", "Low-level cloudiness", "-",
410	. iim,jj_nb,nhori, 1,1,1, -99, 32,
411	. "ave(X)", zstophy,zout)
412	c
413	CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", "-",
414	. iim,jj_nb,nhori, 1,1,1, -99, 32,
415	. "ave(X)", zstophy,zout)
416	c
417	CALL histdef(nid_mth, "cldh", "High-level cloudiness", "-",
418	. iim,jj_nb,nhori, 1,1,1, -99, 32,
419	. "ave(X)", zstophy,zout)
420	c
421	CALL histdef(nid_mth, "cldt", "Total cloudiness", "-",
422	. iim,jj_nb,nhori, 1,1,1, -99, 32,
423	. "ave(X)", zstophy,zout)
424	c
425	CALL histdef(nid_mth,"cldq","Cloud liquid water path","kg/m2",
426	. iim,jj_nb,nhori, 1,1,1, -99, 32,
427	. "ave(X)", zstophy,zout)
428	c
429	CALL histdef(nid_mth,"lwp","Cloud water path","kg/m2",
430	. iim,jj_nb,nhori, 1,1,1, -99, 32,
431	. "ave(X)", zstophy,zout)
432	c
433	CALL histdef(nid_mth,"iwp","Cloud ice water path","kg/m2",
434	. iim,jj_nb,nhori, 1,1,1, -99, 32,
435	. "ave(X)", zstophy,zout)
436	c
437	CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
438	. iim,jj_nb,nhori, 1,1,1, -99, 32,
439	. "ave(X)", zstophy,zout)
440	c
441	CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
442	. iim,jj_nb,nhori, 1,1,1, -99, 32,
443	. "ave(X)", zstophy,zout)
444	c
445	CALL histdef(nid_mth, "uq", "Zonal humidity transport", "-",
446	. iim,jj_nb,nhori, 1,1,1, -99, 32,
447	. "ave(X)", zstophy,zout)
448	c
449	CALL histdef(nid_mth, "vq", "Merid humidity transport", "-",
450	. iim,jj_nb,nhori, 1,1,1, -99, 32,
451	. "ave(X)", zstophy,zout)
452	cKE43
453	IF(iflag_con.GE.3) THEN ! sb
454	c
455	CALL histdef(nid_mth, "cape", "Conv avlbl pot ener", "J/kg",
456	. iim,jj_nb,nhori, 1,1,1, -99, 32,
457	. "ave(X)", zstophy,zout)
458	c
459	CALL histdef(nid_mth, "pbase", "Cld base pressure", "mb",
460	. iim,jj_nb,nhori, 1,1,1, -99, 32,
461	. "ave(X)", zstophy,zout)
462	c
463	CALL histdef(nid_mth, "ptop", "Cld top pressure", "mb",
464	. iim,jj_nb,nhori, 1,1,1, -99, 32,
465	. "ave(X)", zstophy,zout)
466	c
467	CALL histdef(nid_mth, "fbase", "Cld base mass flux", "kg/m2/s",
468	. iim,jj_nb,nhori, 1,1,1, -99, 32,
469	. "ave(X)", zstophy,zout)
470	c
471	CALL histdef(nid_mth, "prw", "Precipitable water", "kg/m2",
472	. iim,jj_nb,nhori, 1,1,1, -99, 32,
473	. "ave(X)", zstophy,zout)
474	c
475	ENDIF !iflag_con .GE. 3
476	c34EK
477	cIM cf. AM 081204 BEG
478	c34EK
479	c HBTM2
480	CALL histdef(nid_mth, "s_pblh", "Boundary Layer Height", "m",
481	. iim,jj_nb,nhori, 1,1,1, -99, 32,
482	. "ave(X)", zstophy,zout)
483	CALL histdef(nid_mth, "s_pblt", "t at Boundary Layer Height",
484	. "K",
485	. iim,jj_nb,nhori, 1,1,1, -99, 32,
486	. "ave(X)", zstophy,zout)
487	CALL histdef(nid_mth, "s_lcl", "Condensation level", "m",
488	. iim,jj_nb,nhori, 1,1,1, -99, 32,
489	. "ave(X)", zstophy,zout)
490	CALL histdef(nid_mth, "s_capCL", "Conv avlbl pot ener
491	. for ABL", "J/m2",
492	. iim,jj_nb,nhori, 1,1,1, -99, 32,
493	. "ave(X)", zstophy,zout)
494	CALL histdef(nid_mth, "s_oliqCL", "Liq Water in BL", "kg/m2",
495	. iim,jj_nb,nhori, 1,1,1, -99, 32,
496	. "ave(X)", zstophy,zout)
497	CALL histdef(nid_mth, "s_cteiCL", "Instability criteria
498	. (ABL)", "K",
499	. iim,jj_nb,nhori, 1,1,1, -99, 32,
500	. "ave(X)", zstophy,zout)
501	CALL histdef(nid_mth, "s_therm", "Exces du thermique", "K",
502	. iim,jj_nb,nhori, 1,1,1, -99, 32,
503	. "ave(X)", zstophy,zout)
504	CALL histdef(nid_mth, "s_trmb1", "deep_cape(HBTM2)", "J/m2",
505	. iim,jj_nb,nhori, 1,1,1, -99, 32,
506	. "ave(X)", zstophy,zout)
507	CALL histdef(nid_mth, "s_trmb2", "inhibition (HBTM2)", "J/m2",
508	. iim,jj_nb,nhori, 1,1,1, -99, 32,
509	. "ave(X)", zstophy,zout)
510	CALL histdef(nid_mth, "s_trmb3", "Point Omega (HBTM2)", "m",
511	. iim,jj_nb,nhori, 1,1,1, -99, 32,
512	. "ave(X)", zstophy,zout)
513	c
514	cIM cf. AM 081204 END
515	c
516	c Champs interpolles sur des niveaux de pression
517	DO k=1, nlevSTD
518	c
519	c bb=clevSTD(k)
520	c
521	c IF(k.GE.2) THEN
522	IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
523	IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
524	c aa=clevSTD(k)
525	c bb=aa(1:lnblnk1(aa))
526	c ENDIF
527	c
528	IF(bb2.EQ."850".OR.bb2.EQ."700".OR.
529	$ bb2.EQ."500".OR.bb2.EQ."200") THEN
530	c
531	CALL histdef(nid_mth, "u"//bb2,
532	. "Zonal wind "//bb2//"mb","m/s",
533	. iim,jj_nb,nhori, 1,1,1, -99, 32,
534	. "inst(X)", zout,zout)
535	c
536	CALL histdef(nid_mth, "v"//bb2,
537	. "Meridional wind "//bb2//"mb","m/s",
538	. iim,jj_nb,nhori, 1,1,1, -99, 32,
539	. "inst(X)", zout,zout)
540	c
541	CALL histdef(nid_mth, "w"//bb2,
542	. "Vertical wind "//bb2//"mb","m/s",
543	. iim,jj_nb,nhori, 1,1,1, -99, 32,
544	. "inst(X)", zout,zout)
545	c
546	CALL histdef(nid_mth, "phi"//bb2,
547	. "Geopotential "//bb2//"mb","m",
548	. iim,jj_nb,nhori, 1,1,1, -99, 32,
549	. "inst(X)", zout,zout)
550	c
551	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.
552	ENDDO
553	c
554	cIM diagnostiques flux ocean-atm ou ocean-glace de mer
555	cIM pour utilisation dans un modele de "slab" ocean
556	c
557	CALL histdef(nid_mth, "fluxo",
558	. "Flux turbulents ocean-atmosphere", "W/m2",
559	. iim,jj_nb,nhori, 1,1,1, -99, 32,
560	. "ave(X)", zstophy,zout)
561	c
562	CALL histdef(nid_mth, "fluxg",
563	. "Flux turbulents ocean-glace de mer","W/m2",
564	. iim,jj_nb,nhori, 1,1,1, -99, 32,
565	. "ave(X)", zstophy,zout)
566	c
567	CALL histdef(nid_mth, "t_oce_sic",
568	. "Temp. mixte oce-sic","K",
569	. iim,jj_nb,nhori, 1,1,1, -99, 32,
570	. "ave(X)", zstophy,zout)
571	c
572	IF (OCEAN.EQ.'force ') THEN
573	c
574	CALL histdef(nid_mth, "lmt_bils",
575	. "Bilan au sol atmosphere forcee", "W/m2",
576	. iim,jj_nb,nhori, 1,1,1, -99, 32,
577	. "ave(X)", zstophy,zout)
578	c
579	ELSE IF (OCEAN.EQ.'slab ') THEN
580	c
581	CALL histdef(nid_mth, "slab_bils",
582	. "Bilan au sol Slab", "W/m2",
583	. iim,jj_nb,nhori, 1,1,1, -99, 32,
584	. "ave(X)", zstophy,zout)
585	c
586	CALL histdef(nid_mth, "tslab", "Slab SST ", "K",
587	. iim,jj_nb,nhori, 1,1,1, -99, 32,
588	. "ave(X)", zstophy,zout)
589	c
590	CALL histdef(nid_mth, "seaice", "Slab seaice", "kg/m2",
591	. iim,jj_nb,nhori, 1,1,1, -99, 32,
592	. "ave(X)", zstophy,zout)
593	c
594	CALL histdef(nid_mth, "siceh", "Slab seaice height", "m",
595	. iim,jj_nb,nhori, 1,1,1, -99, 32,
596	. "ave(X)", zstophy,zout)
597	c
598	ENDIF
599	c
600	ENDIF !lev_histmth.GE.1
601	c
602	IF(lev_histmth.GE.2) THEN
603	c
604	c Champs 3D:
605	c
606	c
607	CALL histdef(nid_mth,"lwcon","Cloud liquid water content",
608	. "kg/kg",
609	. iim,jj_nb,nhori, klev,1,klev, nvert, 32,
610	. "ave(X)", zstophy,zout)
611	c
612	CALL histdef(nid_mth,"iwcon","Cloud ice water content","kg/kg",
613	. iim,jj_nb,nhori, klev,1,klev, nvert, 32,
614	. "ave(X)", zstophy,zout)
615	c
616	CALL histdef(nid_mth, "temp", "Air temperature", "K",
617	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
618	. "ave(X)", zstophy,zout)
619	c
620	CALL histdef(nid_mth, "ovap", "Specific humidity", "kg/kg",
621	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
622	. "ave(X)", zstophy,zout)
623	c
624	c CALL histdef(nid_mth,"wvap","Water vapor mixing ratio","kg/kg",
625	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
626	c . "ave(X)", zstophy,zout)
627	c
628	CALL histdef(nid_mth, "geop", "Geopotential height", "m2/s2",
629	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
630	. "ave(X)", zstophy,zout)
631	c
632	CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s",
633	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
634	. "ave(X)", zstophy,zout)
635	c
636	CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s",
637	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
638	. "ave(X)", zstophy,zout)
639	c
640	CALL histdef(nid_mth, "vitw", "Vertical wind", "Pa/s",
641	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
642	. "ave(X)", zstophy,zout)
643	c
644	CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
645	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
646	. "ave(X)", zstophy,zout)
647	c
648	CALL histdef(nid_mth, "rneb", "Cloud fraction", "-",
649	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
650	. "ave(X)", zstophy,zout)
651	c
652	CALL histdef(nid_mth, "rnebcon", "Convective Cloud Fraction"
653	. , "-",
654	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
655	. "ave(X)", zstophy,zout)
656	c
657	CALL histdef(nid_mth, "rhum", "Relative humidity", "-",
658	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
659	. "ave(X)", zstophy,zout)
660	c
661	CALL histdef(nid_mth, "ozone", "Ozone concentration", "ppmv",
662	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
663	. "ave(X)", zstophy,zout)
664	c
665	CALL histdef(nid_mth, "upwd", "saturated updraft", "kg/m2/s",
666	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
667	. "ave(X)", zstophy,zout)
668	c
669	CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
670	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
671	. "ave(X)", zstophy,zout)
672	c
673	CALL histdef(nid_mth, "dqphy", "Physics dQ", "(kg/kg)/s",
674	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
675	. "ave(X)", zstophy,zout)
676	c
677	ENDIF !lev_histmth.GE.2
678	c
679	IF(lev_histmth.GE.3) THEN
680	c
681	DO nsrf=1, nbsrf
682	c
683	call histdef(nid_mth, "albe_"//clnsurf(nsrf),
684	$ "Albedo surf. "//clnsurf(nsrf), "-",
685	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
686	$ "ave(X)", zstophy,zout)
687	c
688	call histdef(nid_mth, "rugs_"//clnsurf(nsrf),
689	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
690	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
691	$ "ave(X)", zstophy,zout)
692	c
693	CALL histdef(nid_mth, "ages_"//clnsurf(nsrf), "Snow age","day",
694	. iim,jj_nb,nhori, 1,1,1, -99, 32,
695	. "ave(X)", zstophy,zout)
696	c
697	ENDDO !nsrf=1, nbsrf
698	c
699	CALL histdef(nid_mth, "albs", "Surface albedo", "-",
700	. iim,jj_nb,nhori, 1,1,1, -99, 32,
701	. "ave(X)", zstophy,zout)
702	CALL histdef(nid_mth, "albslw", "Surface albedo LW", "-",
703	. iim,jj_nb,nhori, 1,1,1, -99, 32,
704	. "ave(X)", zstophy,zout)
705	c
706	ENDIF !lev_histmth.GE.3
707	c
708	IF(lev_histmth.GE.4) THEN
709	c
710	!FH Sorties pour la couche limite
711	if (iflag_pbl>1) then
712	CALL histdef(nid_mth, "tke","TKE","m2/s2",
713	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
714	. "ave(X)", zstophy,zout)
715	c
716	CALL histdef(nid_mth, "tke_max","TKE max","m2/s2",
717	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
718	. "t_max(X)", zstophy,zout)
719	endif
720	c
721	CALL histdef(nid_mth, "kz","Kz melange","m2/s",
722	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
723	. "ave(X)", zstophy,zout)
724	c
725	CALL histdef(nid_mth, "kz_max","Kz melange max","m2/s",
726	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
727	. "t_max(X)", zstophy,zout)
728
729
730	CALL histdef(nid_mth, "clwcon",
731	. "Convective Cloud Liquid water content"
732	. , "kg/kg",
733	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
734	. "ave(X)", zstophy,zout)
735	c
736	CALL histdef(nid_mth,"Ma","undilute adiab updraft","kg/m2/s",
737	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
738	. "ave(X)", zstophy,zout)
739	c
740	CALL histdef(nid_mth, "dnwd", "saturated downdraft","kg/m2/s",
741	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
742	. "ave(X)", zstophy,zout)
743	c
744	CALL histdef(nid_mth, "dnwd0", "unsat. downdraft", "kg/m2/s",
745	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
746	. "ave(X)", zstophy,zout)
747	c
748	CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
749	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
750	. "ave(X)", zstophy,zout)
751	c
752	CALL histdef(nid_mth, "dqdyn", "Dynamics dQ", "(kg/kg)/s",
753	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
754	. "ave(X)", zstophy,zout)
755	c
756	CALL histdef(nid_mth, "dtcon", "Convection dT", "K/s",
757	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
758	. "ave(X)", zstophy,zout)
759	c
760
761
762	c
763	c CALL histdef(nid_mth, "ducon", "Convection du", "m/s2",
764	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
765	c . "ave(X)", zstophy,zout)
766	c
767	CALL histdef(nid_mth, "dqcon", "Convection dQ", "(kg/kg)/s",
768	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
769	. "ave(X)", zstophy,zout)
770	c
771	CALL histdef(nid_mth, "dtlsc", "Condensation dT", "K/s",
772	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
773	. "ave(X)", zstophy,zout)
774	c
775	CALL histdef(nid_mth, "dtlschr",
776	$ "Large-scale condensational heating rate", "K/s",iim,jj_nb
777	$ ,nhori, klev,1,klev,nvert, 32,"ave(X)", zstophy,zout)
778	c
779	CALL histdef(nid_mth, "dqlsc", "Condensation dQ", "(kg/kg)/s",
780	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
781	. "ave(X)", zstophy,zout)
782	c
783	CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s",
784	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
785	. "ave(X)", zstophy,zout)
786	c
787	CALL histdef(nid_mth, "dqvdf","Boundary-layer dQ","(kg/kg)/s",
788	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
789	. "ave(X)", zstophy,zout)
790	c
791	CALL histdef(nid_mth, "dteva", "Reevaporation dT", "K/s",
792	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
793	. "ave(X)", zstophy,zout)
794	c
795	CALL histdef(nid_mth, "dqeva","Reevaporation dQ","(kg/kg)/s",
796	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
797	. "ave(X)", zstophy,zout)
798
799	CALL histdef(nid_mth, "ptconv", "POINTS CONVECTIFS"," ",
800	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
801	. "ave(X)", zstophy,zout)
802
803	CALL histdef(nid_mth, "ratqs", "RATQS"," ",
804	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
805	. "ave(X)", zstophy,zout)
806
807	c
808	CALL histdef(nid_mth, "dtthe", "Dry adjust. dT", "K/s",
809	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
810	. "ave(X)", zstophy,zout)
811
812	CALL histdef(nid_mth,"dqthe","Dry adjust. dQ","(kg/kg)/s",
813	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
814	. "ave(X)", zstophy,zout)
815	c
816	CALL histdef(nid_mth, "dtajs", "Dry adjust. dT", "K/s",
817	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
818	. "ave(X)", zstophy,zout)
819
820	CALL histdef(nid_mth,"dqajs","Dry adjust. dQ","(kg/kg)/s",
821	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
822	. "ave(X)", zstophy,zout)
823	c
824	CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s",
825	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
826	. "ave(X)", zstophy,zout)
827	c
828	CALL histdef(nid_mth, "dtsw0", "CS SW radiation dT", "K/s",
829	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
830	. "ave(X)", zstophy,zout)
831	c
832	CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
833	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
834	. "ave(X)", zstophy,zout)
835	c
836	CALL histdef(nid_mth,"dtlw0","CS LW radiation dT","K/s",
837	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
838	. "ave(X)", zstophy,zout)
839	c
840	CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
841	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
842	. "ave(X)", zstophy,zout)
843	c
844	CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
845	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
846	. "ave(X)", zstophy,zout)
847	c
848	CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
849	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
850	. "ave(X)", zstophy,zout)
851	c
852	IF (ok_orodr) THEN
853	CALL histdef(nid_mth, "duoro", "Orography dU", "m/s2",
854	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
855	. "ave(X)", zstophy,zout)
856	c
857	CALL histdef(nid_mth, "dvoro", "Orography dV", "m/s2",
858	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
859	. "ave(X)", zstophy,zout)
860	c
861	ENDIF
862	C
863	IF (ok_orolf) THEN
864	CALL histdef(nid_mth, "dulif", "Orography dU", "m/s2",
865	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
866	. "ave(X)", zstophy,zout)
867	c
868	CALL histdef(nid_mth, "dvlif", "Orography dV", "m/s2",
869	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
870	. "ave(X)", zstophy,zout)
871	ENDIF
872	C
873	#if defined(INCA_AER) && defined(CPP_COUPLE)
874	CALL histdef(nid_mth, "tau", "cloud opt thickness", "-",
875	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
876	. "ave(X)", zstophy,zout)
877	CALL histdef(nid_mth, "CCM1", "CCM1", "",
878	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
879	. "ave(X)", zstophy,zout)
880	CALL histdef(nid_mth, "CCM2", "CCM2", "",
881	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
882	. "ave(X)", zstophy,zout)
883	CALL histdef(nid_mth, "taupi", "cloud opt thickn. (pi)", "-",
884	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
885	. "ave(X)", zstophy,zout)
886	CALL histdef(nid_mth, "cldemi", "emissivite IR", "-",
887	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
888	. "ave(X)", zstophy,zout)
889	CALL histdef(nid_mth, "SWTOAAS_ZERO", "?", "?",
890	. iim,jj_nb,nhori, 1,1,1, -99, 32,
891	. "ave(X)", zstophy,zout)
892	CALL histdef(nid_mth, "SWTOAAS_AER", "?", "?",
893	. iim,jj_nb,nhori, 1,1,1, -99, 32,
894	. "ave(X)", zstophy,zout)
895	CALL histdef(nid_mth, "SWTOAAS_NAT", "?", "?",
896	. iim,jj_nb,nhori, 1,1,1, -99, 32,
897	. "ave(X)", zstophy,zout)
898	CALL histdef(nid_mth, "SWTOAAS_AD", "?", "?",
899	. iim,jj_nb,nhori, 1,1,1, -99, 32,
900	. "ave(X)", zstophy,zout)
901	CALL histdef(nid_mth, "SWTOAAS_AI", "?", "?",
902	. iim,jj_nb,nhori, 1,1,1, -99, 32,
903	. "ave(X)", zstophy,zout)
904	CALL histdef(nid_mth, "SWTOACS_ZERO", "?", "?",
905	. iim,jj_nb,nhori, 1,1,1, -99, 32,
906	. "ave(X)", zstophy,zout)
907	CALL histdef(nid_mth, "SWTOACS_AER", "?", "?",
908	. iim,jj_nb,nhori, 1,1,1, -99, 32,
909	. "ave(X)", zstophy,zout)
910	CALL histdef(nid_mth, "SWTOACS_AD", "?", "?",
911	. iim,jj_nb,nhori, 1,1,1, -99, 32,
912	. "ave(X)", zstophy,zout)
913	CALL histdef(nid_mth, "SWTOACS_NAT", "?", "?",
914	. iim,jj_nb,nhori, 1,1,1, -99, 32,
915	. "ave(X)", zstophy,zout)
916	#endif
917	#ifndef INCA
918	if (nqmax.GE.3) THEN
919	DO iq=3,nqmax
920	iiq=niadv(iq)
921	CALL histdef(nid_mth, tnom(iq), ttext(iiq), "-",
922	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
923	. "ave(X)", zstophy,zout)
924	ENDDO
925	ENDIF
926	#endif
927	c
928	ENDIF !lev_histmth.GE.4
929	c
930	ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
931	c
932	IF(lev_histmth.GE.1) THEN
933	c
934	CALL histdef(nid_mth, "phis", "Surface geop. height", "m",
935	. iim,jj_nb,nhori, 1,1,1, -99, 32,
936	. "once", zstophy,zout)
937	c
938	CALL histdef(nid_mth, "aire", "Grid area", "-",
939	. iim,jj_nb,nhori, 1,1,1, -99, 32,
940	. "once", zstophy,zout)
941	c
942	CALL histdef(nid_mth, "pourc_"//clnsurf(is_ter),
943	. "% "//clnsurf(is_ter), "%",
944	. iim,jj_nb,nhori, 1,1,1, -99, 32,
945	. "once", zstophy,zout)
946	c
947	CALL histdef(nid_mth, "fract_"//clnsurf(is_ter),
948	. "Fraction "//clnsurf(is_ter), "1",
949	. iim,jj_nb,nhori, 1,1,1, -99, 32,
950	. "once", zstophy,zout)
951	c
952	CALL histdef(nid_mth, "pourc_"//clnsurf(is_lic),
953	. "% "//clnsurf(is_lic), "%",
954	. iim,jj_nb,nhori, 1,1,1, -99, 32,
955	. "once", zstophy,zout)
956	c
957	CALL histdef(nid_mth, "fract_"//clnsurf(is_lic),
958	. "Fraction "//clnsurf(is_lic), "1",
959	. iim,jj_nb,nhori, 1,1,1, -99, 32,
960	. "once", zstophy,zout)
961	c
962	CALL histdef(nid_mth, "pourc_"//clnsurf(is_sic),
963	. "% "//clnsurf(is_sic), "%",
964	. iim,jj_nb,nhori, 1,1,1, -99, 32,
965	. "ave(X)", zstophy,zout)
966	c
967	CALL histdef(nid_mth, "fract_"//clnsurf(is_sic),
968	. "Fraction "//clnsurf(is_sic), "1",
969	. iim,jj_nb,nhori, 1,1,1, -99, 32,
970	. "ave(X)", zstophy,zout)
971	c
972	CALL histdef(nid_mth, "flat", "Latent heat flux", "W/m2",
973	. iim,jj_nb,nhori, 1,1,1, -99, 32,
974	. "ave(X)", zstophy,zout)
975	c
976	CALL histdef(nid_mth, "bilTOA",
977	. "Net radiation at TOA", "W/m2",
978	. iim,jj_nb,nhori, 1,1,1, -99, 32,
979	. "ave(X)", zstophy,zout)
980	c
981	c
982	c CALL histdef(nid_mth, "tops01",
983	c . "TOA reflected clear-sky shortwave
984	c . radiation (method 1)", "W/m2",
985	c . iim,jj_nb,nhori, 1,1,1, -99, 32,
986	c . "ave(X)", zstophy,zout)
987	c
988	CALL histdef(nid_mth, "dtdia","dtdia","K/s",
989	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
990	. "ave(X)", zstophy,zout)
991	c
992	CALL histdef(nid_mth, "dtsw","dtsw","K/s",
993	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
994	. "ave(X)", zstophy,zout)
995	c
996	CALL histdef(nid_mth, "dtlw","dtlw","K/s",
997	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
998	. "ave(X)", zstophy,zout)
999	c
1000	CALL histdef(nid_mth,"dqthe","Dry adjust. dQ","(kg/kg)/s",
1001	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1002	. "ave(X)", zstophy,zout)
1003	c
1004	CALL histdef(nid_mth, "dtajs","dtajs","K/s",
1005	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1006	. "ave(X)", zstophy,zout)
1007	c
1008	CALL histdef(nid_mth, "dtlsc","dtlsc","K/s",
1009	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1010	. "ave(X)", zstophy,zout)
1011	c
1012	CALL histdef(nid_mth, "dqdia","dqdia","(kg/kg)/s",
1013	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1014	. "ave(X)", zstophy,zout)
1015	c
1016	CALL histdef(nid_mth, "duoli",
1017	. "duoli","m/s2",
1018	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1019	. "ave(X)", zstophy,zout)
1020	c
1021	CALL histdef(nid_mth, "dvoli",
1022	. "dvoli","m/s2",
1023	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1024	. "ave(X)", zstophy,zout)
1025	c
1026	CALL histdef(nid_mth, "ducon",
1027	. "ducon","m/s2",
1028	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1029	. "ave(X)", zstophy,zout)
1030	c
1031	CALL histdef(nid_mth, "dvcon",
1032	. "dvcon","m/s2",
1033	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1034	. "ave(X)", zstophy,zout)
1035	c
1036	CALL histdef(nid_mth, "dutot",
1037	. "dutot","m/s2",
1038	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1039	. "ave(X)", zstophy,zout)
1040	c
1041	CALL histdef(nid_mth, "dvtot",
1042	. "dvtot","m/s2",
1043	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1044	. "ave(X)", zstophy,zout)
1045	c
1046	c CALL histdef(nid_mth, "cldtrue",
1047	c . "TRUE Cloud fraction ","%",
1048	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1049	c . "inst(X)", zout,zout)
1050	c
1051	c CALL histdef(nid_mth, "cldTOA",
1052	c . "Cloud amount (satellite view)","%",
1053	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1054	c . "ave(X)", zstophy,zout)
1055	c
1056	c CALL histdef(nid_mth, "cldSFC",
1057	c . "Cloud amount (surface view)","%",
1058	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1059	c . "ave(X)", zstophy,zout)
1060	c
1061	CALL histdef(nid_mth, "GWDtaux",
1062	. "GWD induced Ewd wind stress","Pa",
1063	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1064	. "ave(X)", zstophy,zout)
1065	c
1066	CALL histdef(nid_mth, "GWDtauy",
1067	. "GWD induced Nwd wind stress","Pa",
1068	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1069	. "ave(X)", zstophy,zout)
1070	c
1071	c Champs 2D:
1072	c
1073	CALL histdef(nid_mth, "slp", "Sea Level Pressure", "Pa",
1074	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1075	. "ave(X)", zstorad,zout)
1076	c
1077	CALL histdef(nid_mth, "tsol", "Surface Temperature", "K",
1078	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1079	. "ave(X)", zstophy,zout)
1080	c
1081	CALL histdef(nid_mth, "t2m", "Temperature 2m", "K",
1082	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1083	. "ave(X)", zstorad,zout)
1084	c
1085	c CALL histdef(nid_mth, "t2m_min", "Temp. 2m min.",
1086	c . "K",
1087	c . iim,jj_nb,nhori, 1,1,1, -99, 32,
1088	c . t2mincels, zstophy,zout)
1089	c
1090	c CALL histdef(nid_mth, "t2m_max", "Temp. 2m max.",
1091	c . "K",
1092	c . iim,jj_nb,nhori, 1,1,1, -99, 32,
1093	c . t2maxcels, zstophy,zout)
1094	c
1095	c CALL histdef(nid_mth, "tsoil", "Sfce soil Temperature",
1096	c . "K",
1097	c . iim,jj_nb,nhori, 1,1,1, -99, 32,
1098	c . "ave(X)", zstophy,zout)
1099	c
1100	CALL histdef(nid_mth, "wind10m","10-m wind speed","m/s",
1101	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1102	. "ave(X)", zstorad,zout)
1103	c
1104	CALL histdef(nid_mth, "sicf", "Sea-ice fraction", "-",
1105	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1106	. "ave(X)", zstophy,zout)
1107
1108	c
1109	CALL histdef(nid_mth, "q2m","Specific humidity 2m","kg/kg",
1110	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1111	. "ave(X)", zstorad,zout)
1112	c
1113	CALL histdef(nid_mth, "u10m", "Vent zonal 10m", "m/s",
1114	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1115	. "ave(X)", zstorad,zout)
1116	c
1117	CALL histdef(nid_mth, "v10m", "Vent meridien 10m", "m/s",
1118	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1119	. "ave(X)", zstorad,zout)
1120	c
1121	CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa",
1122	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1123	. "ave(X)", zstorad,zout)
1124	c
1125	CALL histdef(nid_mth, "qsurf","Surface Air humidity","kg/kg",
1126	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1127	. "ave(X)", zstophy,zout)
1128	c
1129	if (.not. ok_veget) then
1130	CALL histdef(nid_mth, "qsol", "Soil watter content", "mm",
1131	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1132	. "ave(X)", zstophy,zout)
1133	endif
1134	c
1135	CALL histdef(nid_mth, "ndayrain",
1136	. "Number of day with rain (liq+sol)", "-",
1137	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1138	. "inst(X)", zstomth,zout)
1139	c
1140	CALL histdef(nid_mth, "precip", "Precipitation Totale liq+sol",
1141	. "kg/(s*m2)",
1142	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1143	. "ave(X)", zstophy,zout)
1144	c
1145	CALL histdef(nid_mth, "plul", "Large-scale Precip.",
1146	. "kg/(s*m2)",
1147	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1148	. "ave(X)", zstophy,zout)
1149	c
1150	CALL histdef(nid_mth, "pluc", "Convective Precip.",
1151	. "kg/(s*m2)",
1152	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1153	. "ave(X)", zstophy,zout)
1154	c
1155	CALL histdef(nid_mth, "snowf", "Snow fall", "kg/(s*m2)",
1156	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1157	. "ave(X)", zstophy,zout)
1158	c
1159	c CALL histdef(nid_mth, "snow_mass", "Snow Mass", "kg/m2",
1160	c . iim,jj_nb,nhori, 1,1,1, -99, 32,
1161	c . "ave(X)", zstophy,zout)
1162	c
1163	CALL histdef(nid_mth, "evap",
1164	. "Evaporation plus sublimation rate", "kg/(s*m2)",
1165	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1166	. "ave(X)", zstophy,zout)
1167	c
1168	CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2",
1169	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1170	. "ave(X)", zstophy,zout)
1171	c
1172	CALL histdef(nid_mth, "tops0", "CS Solar rad. at TOA", "W/m2",
1173	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1174	. "ave(X)", zstophy,zout)
1175	c
1176	CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2",
1177	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1178	. "ave(X)", zstophy,zout)
1179	c
1180	CALL histdef(nid_mth, "topl0", "IR rad. at TOA", "W/m2",
1181	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1182	. "ave(X)", zstophy,zout)
1183	c
1184	CALL histdef(nid_mth, "SWupTOA", "SWup at TOA","W/m2",
1185	. iim,jj_nb,nhori, 1,1,1,-99,
1186	. 32, "ave(X)", zstophy,zout)
1187	c
1188	CALL histdef(nid_mth, "SWupTOAclr",
1189	. "SWup clear sky at TOA","W/m2",
1190	. iim,jj_nb,nhori, 1,1,1,-99,
1191	. 32, "ave(X)", zstophy,zout)
1192	c
1193	CALL histdef(nid_mth, "SWdnTOA", "SWdn at TOA","W/m2",
1194	. iim,jj_nb,nhori, 1,1,1,-99,
1195	. 32, "ave(X)", zstophy,zout)
1196	c
1197	CALL histdef(nid_mth, "SWdnTOAclr",
1198	. "SWdn clear sky at TOA","W/m2",
1199	. iim,jj_nb,nhori, 1,1,1,-99,
1200	. 32, "ave(X)", zstophy,zout)
1201	c
1202	CALL histdef(nid_mth, "SWup200", "SWup at 200mb","W/m2",
1203	. iim,jj_nb,nhori, 1,1,1,-99,
1204	. 32, "ave(X)", zstophy,zout)
1205	c
1206	CALL histdef(nid_mth, "SWup200clr",
1207	. "SWup clear sky at 200mb","W/m2",
1208	. iim,jj_nb,nhori, 1,1,1,-99,
1209	. 32, "ave(X)", zstophy,zout)
1210	c
1211	CALL histdef(nid_mth, "SWdn200", "SWdn at 200mb","W/m2",
1212	. iim,jj_nb,nhori, 1,1,1,-99,
1213	. 32, "ave(X)", zstophy,zout)
1214	c
1215	CALL histdef(nid_mth, "SWdn200clr",
1216	. "SWdn clear sky at 200mb","W/m2",
1217	. iim,jj_nb,nhori, 1,1,1,-99,
1218	. 32, "ave(X)", zstophy,zout)
1219	c
1220	CALL histdef(nid_mth, "LWup200", "LWup at 200mb","W/m2",
1221	. iim,jj_nb,nhori, 1,1,1,-99,
1222	. 32, "ave(X)", zstophy,zout)
1223	c
1224	CALL histdef(nid_mth, "LWup200clr",
1225	. "LWup clear sky at 200mb","W/m2",
1226	. iim,jj_nb,nhori, 1,1,1,-99,
1227	. 32, "ave(X)", zstophy,zout)
1228	c
1229	CALL histdef(nid_mth, "LWdn200", "LWdn at 200mb","W/m2",
1230	. iim,jj_nb,nhori, 1,1,1,-99,
1231	. 32, "ave(X)", zstophy,zout)
1232	c
1233	CALL histdef(nid_mth, "LWdn200clr",
1234	. "LWdn clear sky at 200mb","W/m2",
1235	. iim,jj_nb,nhori, 1,1,1,-99,
1236	. 32, "ave(X)", zstophy,zout)
1237	c
1238	CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2",
1239	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1240	. "ave(X)", zstophy,zout)
1241	c
1242	CALL histdef(nid_mth, "sols0", "Solar rad. at surf.", "W/m2",
1243	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1244	. "ave(X)", zstophy,zout)
1245	c
1246	CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2",
1247	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1248	. "ave(X)", zstophy,zout)
1249	c
1250	CALL histdef(nid_mth, "radsol", "Rayonnement au sol", "W/m2",
1251	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1252	. "ave(X)", zstophy,zout)
1253	c
1254	CALL histdef(nid_mth, "soll0", "IR rad. at surface", "W/m2",
1255	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1256	. "ave(X)", zstophy,zout)
1257	c
1258	CALL histdef(nid_mth, "SWupSFC", "SWup at surface","W/m2",
1259	. iim,jj_nb,nhori, 1,1,1,-99,
1260	. 32, "ave(X)", zstophy,zout)
1261	c
1262	CALL histdef(nid_mth, "SWupSFCclr",
1263	. "SWup clear sky at surface","W/m2",
1264	. iim,jj_nb,nhori, 1,1,1,-99,
1265	. 32, "ave(X)", zstophy,zout)
1266	c
1267	CALL histdef(nid_mth, "SWdnSFC", "SWdn at surface","W/m2",
1268	. iim,jj_nb,nhori, 1,1,1,-99,
1269	. 32, "ave(X)", zstophy,zout)
1270	c
1271	CALL histdef(nid_mth, "SWdnSFCclr",
1272	. "SWdn clear sky at surface","W/m2",
1273	. iim,jj_nb,nhori, 1,1,1,-99,
1274	. 32, "ave(X)", zstophy,zout)
1275	c
1276	CALL histdef(nid_mth,"LWupSFC","Upwd. IR rad. at surface",
1277	. "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32,
1278	. "ave(X)", zstophy,zout)
1279	c
1280	CALL histdef(nid_mth,"LWdnSFC","Down. IR rad. at surface",
1281	. "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32,
1282	. "ave(X)", zstophy,zout)
1283	c
1284	CALL histdef(nid_mth,"LWupSFCclr",
1285	. "CS Upwd. IR rad. at surface",
1286	. "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32,
1287	. "ave(X)", zstophy,zout)
1288	c
1289	CALL histdef(nid_mth,"LWdnSFCclr",
1290	. "Down. CS IR rad. at surface",
1291	. "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32,
1292	. "ave(X)", zstophy,zout)
1293	c
1294	CALL histdef(nid_mth, "bils", "Surf. total heat flux", "W/m2",
1295	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1296	. "ave(X)", zstophy,zout)
1297	c
1298	CALL histdef(nid_mth, "sens", "Sensible heat flux", "W/m2",
1299	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1300	. "ave(X)", zstophy,zout)
1301	c
1302	CALL histdef(nid_mth, "fder", "Heat flux derivation", "W/m2",
1303	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1304	. "ave(X)", zstophy,zout)
1305	c
1306	c CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa",
1307	c . iim,jj_nb,nhori, 1,1,1, -99, 32,
1308	c . "ave(X)", zstophy,zout)
1309	c
1310	c CALL histdef(nid_mth, "frtv", "Meridional wind stress", "Pa",
1311	c . iim,jj_nb,nhori, 1,1,1, -99, 32,
1312	c . "ave(X)", zstophy,zout)
1313	c
1314	CALL histdef(nid_mth, "ffonte","Thermal flux for snow melting",
1315	. "W/m2",iim,jj_nb,nhori, 1,1,1, -99, 32,
1316	. "ave(X)", zstophy,zout)
1317	c
1318	CALL histdef(nid_mth, "fqcalving","Ice Calving",
1319	. "kg/m2/s",iim,jj_nb,nhori, 1,1,1, -99, 32,
1320	. "ave(X)", zstophy,zout)
1321	c
1322	DO nsrf = 1, nbsrf
1323	call histdef(nid_mth, "taux_"//clnsurf(nsrf),
1324	$ "Zonal wind stress"//clnsurf(nsrf), "Pa",
1325	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
1326	$ "ave(X)", zstophy,zout)
1327
1328	call histdef(nid_mth, "tauy_"//clnsurf(nsrf),
1329	$ "Meridional wind stress "//clnsurf(nsrf), "Pa",
1330	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
1331	$ "ave(X)", zstophy,zout)
1332	ENDDO
1333	c
1334	call histdef(nid_mth, "pourc_"//clnsurf(is_oce),
1335	$ "% "//clnsurf(is_oce), "%",
1336	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
1337	$ "ave(X)", zstophy,zout)
1338	c
1339	call histdef(nid_mth, "fract_"//clnsurf(is_oce),
1340	$ "Fraction "//clnsurf(is_oce), "1",
1341	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
1342	$ "ave(X)", zstophy,zout)
1343	C
1344	DO nsrf = 1, nbsrf
1345	C
1346	call histdef(nid_mth, "tsol_"//clnsurf(nsrf),
1347	$ "Temperature "//clnsurf(nsrf), "K",
1348	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
1349	$ "ave(X)", zstophy,zout)
1350	C
1351	call histdef(nid_mth, "sens_"//clnsurf(nsrf),
1352	$ "Sensible heat flux "//clnsurf(nsrf), "W/m2",
1353	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
1354	$ "ave(X)", zstophy,zout)
1355	c
1356	call histdef(nid_mth, "lat_"//clnsurf(nsrf),
1357	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
1358	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
1359	$ "ave(X)", zstophy,zout)
1360	C
1361	call histdef(nid_mth, "flw_"//clnsurf(nsrf),
1362	$ "LW "//clnsurf(nsrf), "W/m2",
1363	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
1364	$ "ave(X)", zstophy,zout)
1365	c
1366	call histdef(nid_mth, "fsw_"//clnsurf(nsrf),
1367	$ "SW "//clnsurf(nsrf), "W/m2",
1368	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
1369	$ "ave(X)", zstophy,zout)
1370	C
1371	call histdef(nid_mth, "wbils_"//clnsurf(nsrf),
1372	$ "Bilan sol "//clnsurf(nsrf), "W/m2",
1373	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
1374	$ "ave(X)", zstophy,zout)
1375	C
1376	call histdef(nid_mth, "wbilo_"//clnsurf(nsrf),
1377	$ "Bilan eau "//clnsurf(nsrf), "kg/(m2*s)",
1378	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
1379	$ "ave(X)", zstophy,zout)
1380	C
1381	cIM 250105 BEG
1382	IF(1.EQ.0) THEN
1383	C
1384	call histdef(nid_mth, "evapour_"//clnsurf(nsrf),
1385	$ "Evappourc "//clnsurf(nsrf), "kg/(sm2)",
1386	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
1387	$ "ave(X)", zstophy,zout)
1388	C
1389	call histdef(nid_mth, "prepour_"//clnsurf(nsrf),
1390	$ "Precippourc "//clnsurf(nsrf), "kg/(sm2)",
1391	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
1392	$ "ave(X)", zstophy,zout)
1393	C
1394	ENDIF !(1.EQ.0) THEN
1395	cIM 250105 END
1396	C
1397	END DO
1398	c
1399	CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
1400	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1401	. "ave(X)", zstophy,zout)
1402	c
1403	CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
1404	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1405	. "ave(X)", zstophy,zout)
1406	c
1407	CALL histdef(nid_mth, "cldl", "Low-level cloudiness", "-",
1408	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1409	. "ave(X)", zstophy,zout)
1410	c
1411	CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", "-",
1412	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1413	. "ave(X)", zstophy,zout)
1414	c
1415	CALL histdef(nid_mth, "cldh", "High-level cloudiness", "-",
1416	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1417	. "ave(X)", zstophy,zout)
1418	c
1419	CALL histdef(nid_mth, "cldt", "Total cloudiness", "%",
1420	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1421	. "ave(X)", zstophy,zout)
1422	c
1423	CALL histdef(nid_mth,"cldq","Cloud liquid water path","kg/m2",
1424	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1425	. "ave(X)", zstophy,zout)
1426	c
1427	CALL histdef(nid_mth,"lwp","Cloud water path","kg/m2",
1428	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1429	. "ave(X)", zstophy,zout)
1430	c
1431	CALL histdef(nid_mth,"iwp","Cloud ice water path","kg/m2",
1432	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1433	. "ave(X)", zstophy,zout)
1434	c
1435	CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
1436	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1437	. "ave(X)", zstophy,zout)
1438	c
1439	CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
1440	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1441	. "ave(X)", zstophy,zout)
1442	c
1443	CALL histdef(nid_mth, "uq", "Zonal humidity transport", "-",
1444	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1445	. "ave(X)", zstophy,zout)
1446	c
1447	CALL histdef(nid_mth, "vq", "Merid humidity transport", "-",
1448	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1449	. "ave(X)", zstophy,zout)
1450	cKE43
1451	IF(iflag_con.GE.3) THEN ! sb
1452	c
1453	CALL histdef(nid_mth, "cape", "Conv avlbl pot ener", "J/kg",
1454	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1455	. "ave(X)", zstophy,zout)
1456	c
1457	CALL histdef(nid_mth, "pbase", "Cld base pressure", "mb",
1458	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1459	. "ave(X)", zstophy,zout)
1460	c
1461	CALL histdef(nid_mth, "ptop", "Cld top pressure", "mb",
1462	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1463	. "ave(X)", zstophy,zout)
1464	c
1465	CALL histdef(nid_mth, "fbase", "Cld base mass flux", "kg/m2/s",
1466	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1467	. "ave(X)", zstophy,zout)
1468	c
1469	CALL histdef(nid_mth, "prw", "Precipitable water", "kg/m2",
1470	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1471	. "ave(X)", zstophy,zout)
1472	c
1473	CALL histdef(nid_mth, "upwd", "saturated updraft", "kg/m2/s",
1474	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1475	. "ave(X)", zstophy,zout)
1476	c
1477	ENDIF !iflag_con .GE. 3
1478	c34EK
1479	cIM cf. AM 081204 BEG
1480	c34EK
1481	c HBTM2
1482	CALL histdef(nid_mth, "s_pblh", "Boundary Layer Height", "m",
1483	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1484	. "ave(X)", zstophy,zout)
1485	CALL histdef(nid_mth, "s_pblt", "T at Boundary Layer Height",
1486	. "K",
1487	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1488	. "ave(X)", zstophy,zout)
1489	CALL histdef(nid_mth, "s_lcl", "Condensation level", "m",
1490	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1491	. "ave(X)", zstophy,zout)
1492	CALL histdef(nid_mth,"s_capCL","Conv avlbl pot ener for ABL",
1493	. "J/m2",
1494	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1495	. "ave(X)", zstophy,zout)
1496	CALL histdef(nid_mth, "s_oliqCL", "Liq Water in BL", "kg/m2",
1497	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1498	. "ave(X)", zstophy,zout)
1499	CALL histdef(nid_mth,"s_cteiCL","Instability criteria (ABL)",
1500	. "K",
1501	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1502	. "ave(X)", zstophy,zout)
1503	CALL histdef(nid_mth, "s_therm", "Exces du thermique", "K",
1504	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1505	. "ave(X)", zstophy,zout)
1506	CALL histdef(nid_mth, "s_trmb1", "deep_cape(HBTM2)", "J/m2",
1507	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1508	. "ave(X)", zstophy,zout)
1509	CALL histdef(nid_mth, "s_trmb2", "inhibition (HBTM2)", "J/m2",
1510	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1511	. "ave(X)", zstophy,zout)
1512	CALL histdef(nid_mth, "s_trmb3", "Point Omega (HBTM2)", "m",
1513	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1514	. "ave(X)", zstophy,zout)
1515	c
1516	cIM cf. AM 081204 END
1517	c
1518	c Champs interpolles sur des niveaux de pression
1519	DO k=1, nlevSTD
1520	c
1521	c bb=clevSTD(k)
1522	c
1523	c IF(k.GE.2) THEN
1524	IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
1525	IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
1526	c aa=clevSTD(k)
1527	c bb=aa(1:lnblnk1(aa))
1528	c ENDIF
1529	c
1530	IF(bb2.EQ."850".OR.bb2.EQ."700"
1531	$ .OR.bb2.EQ."500".OR.bb2.EQ."200") THEN
1532	c
1533	CALL histdef(nid_mth, "u"//bb2,
1534	. "Zonal wind "//bb2//"mb","m/s",
1535	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1536	. "inst(X)", zout,zout)
1537	c
1538	CALL histdef(nid_mth, "v"//bb2,
1539	. "Meridional wind "//bb2//"mb","m/s",
1540	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1541	. "inst(X)", zout,zout)
1542	c
1543	CALL histdef(nid_mth, "w"//bb2,
1544	. "Vertical wind "//bb2//"mb","m/s",
1545	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1546	. "inst(X)", zout,zout)
1547	c
1548	CALL histdef(nid_mth, "phi"//bb2,
1549	. "Geopotential "//bb2//"mb","m",
1550	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1551	. "inst(X)", zout,zout)
1552	c
1553	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700"
1554	c
1555	ENDDO !nlevSTD
1556	c
1557	cIM diagnostiques flux ocean-atm ou ocean-glace de mer
1558	cIM pour utilisation dans un modele de "slab" ocean
1559	c
1560	CALL histdef(nid_mth, "fluxo",
1561	. "Flux turbulents ocean-atmosphere", "W/m2",
1562	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1563	. "ave(X)", zstophy,zout)
1564	c
1565	CALL histdef(nid_mth, "fluxg",
1566	. "Flux turbulents ocean-glace de mer","W/m2",
1567	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1568	. "ave(X)", zstophy,zout)
1569	c
1570	CALL histdef(nid_mth, "t_oce_sic",
1571	. "Temp. mixte oce-sic","K",
1572	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1573	. "ave(X)", zstophy,zout)
1574	c
1575	IF (OCEAN.EQ.'force ') THEN
1576	c
1577	CALL histdef(nid_mth, "lmt_bils",
1578	. "Bilan au sol atmosphere forcee", "W/m2",
1579	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1580	. "ave(X)", zstophy,zout)
1581	c
1582	ELSE IF (OCEAN.EQ.'slab ') THEN
1583	c
1584	CALL histdef(nid_mth, "slab_bils",
1585	. "Bilan au sol Slab", "W/m2",
1586	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1587	. "ave(X)", zstophy,zout)
1588	c
1589	CALL histdef(nid_mth, "tslab", "Slab SST ", "K",
1590	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1591	. "ave(X)", zstophy,zout)
1592	c
1593	CALL histdef(nid_mth, "seaice", "Slab seaice", "kg/m2",
1594	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1595	. "ave(X)", zstophy,zout)
1596	c
1597	CALL histdef(nid_mth, "siceh", "Slab seaice height", "m",
1598	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1599	. "ave(X)", zstophy,zout)
1600	c
1601	ENDIF !OCEAN.EQ.slab/force
1602	c
1603	ENDIF !lev_histmth.GE.1
1604	c
1605	IF(lev_histmth.GE.2) THEN
1606	c
1607	c Champs 3D:
1608	c
1609	c
1610	CALL histdef(nid_mth,"lwcon","Cloud liquid water content",
1611	. "kg/kg",
1612	. iim,jj_nb,nhori, klev,1,klev, nvert, 32,
1613	. "ave(X)", zstophy,zout)
1614	c
1615	CALL histdef(nid_mth,"iwcon","Cloud ice water content","kg/kg",
1616	. iim,jj_nb,nhori, klev,1,klev, nvert, 32,
1617	. "ave(X)", zstophy,zout)
1618	c
1619	CALL histdef(nid_mth, "temp", "Air temperature", "K",
1620	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1621	. "ave(X)", zstophy,zout)
1622	c
1623	CALL histdef(nid_mth, "ovap", "Specific humidity", "kg/kg",
1624	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1625	. "ave(X)", zstophy,zout)
1626	c
1627	c CALL histdef(nid_mth,"wvap","Water vapor mixing ratio","kg/kg",
1628	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1629	c . "ave(X)", zstophy,zout)
1630	c
1631	CALL histdef(nid_mth, "geop", "Geopotential height", "m2/s2",
1632	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1633	. "ave(X)", zstophy,zout)
1634	c
1635	CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s",
1636	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1637	. "ave(X)", zstophy,zout)
1638	c
1639	CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s",
1640	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1641	. "ave(X)", zstophy,zout)
1642	c
1643	CALL histdef(nid_mth, "vitw", "Vertical wind", "Pa/s",
1644	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1645	. "ave(X)", zstophy,zout)
1646	c
1647	CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
1648	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1649	. "ave(X)", zstophy,zout)
1650	c
1651	CALL histdef(nid_mth, "rneb", "Cloud fraction", "%",
1652	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1653	. "ave(X)", zstophy,zout)
1654	c
1655	CALL histdef(nid_mth, "rnebcon", "Convective Cloud Fraction"
1656	. , "-",
1657	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1658	. "ave(X)", zstophy,zout)
1659	c
1660	CALL histdef(nid_mth, "rhum", "Relative humidity", "-",
1661	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1662	. "ave(X)", zstophy,zout)
1663	c
1664	CALL histdef(nid_mth, "ozone", "Ozone concentration", "ppmv",
1665	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1666	. "ave(X)", zstophy,zout)
1667	c
1668	CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
1669	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1670	. "ave(X)", zstophy,zout)
1671	c
1672	CALL histdef(nid_mth, "dqphy","Physics dQ","(kg/kg)/s",
1673	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1674	. "ave(X)", zstophy,zout)
1675	c
1676	ENDIF !lev_histmth.GE.2
1677	c
1678	IF(lev_histmth.GE.3) THEN
1679	c
1680	DO nsrf=1, nbsrf
1681	c
1682	call histdef(nid_mth, "albe_"//clnsurf(nsrf),
1683	$ "Albedo surf. "//clnsurf(nsrf), "-",
1684	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
1685	$ "ave(X)", zstophy,zout)
1686	c
1687	call histdef(nid_mth, "rugs_"//clnsurf(nsrf),
1688	$ "Rugosite "//clnsurf(nsrf), "m",
1689	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
1690	$ "ave(X)", zstophy,zout)
1691	c
1692	CALL histdef(nid_mth, "ages_"//clnsurf(nsrf), "Snow age","day",
1693	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1694	. "ave(X)", zstophy,zout)
1695	c
1696	ENDDO !nsrf=1, nbsrf
1697	c
1698	CALL histdef(nid_mth, "albs", "Surface albedo", "-",
1699	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1700	. "ave(X)", zstophy,zout)
1701	CALL histdef(nid_mth, "albslw", "Surface albedo LW", "-",
1702	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1703	. "ave(X)", zstophy,zout)
1704	c
1705	ENDIF !lev_histmth.GE.3
1706	c
1707	IF(lev_histmth.GE.4) THEN
1708	c
1709	!FH Sorties pour la couche limite
1710	if (iflag_pbl>1) then
1711	CALL histdef(nid_mth, "tke","TKE","m2/s2",
1712	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1713	. "ave(X)", zstophy,zout)
1714	c
1715	CALL histdef(nid_mth, "tke_max","TKE max","m2/s2",
1716	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1717	. "t_max(X)", zstophy,zout)
1718	endif
1719	c
1720	CALL histdef(nid_mth, "kz","Kz melange","m2/s",
1721	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1722	. "ave(X)", zstophy,zout)
1723	c
1724	CALL histdef(nid_mth, "kz_max","Kz melange max","m2/s",
1725	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1726	. "t_max(X)", zstophy,zout)
1727
1728
1729
1730
1731	CALL histdef(nid_mth, "clwcon",
1732	. "Convective Cloud Liquid water content"
1733	. , "kg/kg",
1734	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1735	. "ave(X)", zstophy,zout)
1736	c
1737	CALL histdef(nid_mth,"Ma","undilute adiab updraft","kg/m2/s",
1738	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1739	. "ave(X)", zstophy,zout)
1740	c
1741	CALL histdef(nid_mth, "dnwd", "saturated downdraft","kg/m2/s",
1742	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1743	. "ave(X)", zstophy,zout)
1744	c
1745	CALL histdef(nid_mth, "dnwd0", "unsat. downdraft", "kg/m2/s",
1746	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1747	. "ave(X)", zstophy,zout)
1748	c
1749	CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
1750	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1751	. "ave(X)", zstophy,zout)
1752	c
1753	CALL histdef(nid_mth, "dqdyn","Dynamics dQ","(kg/kg)/s",
1754	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1755	. "ave(X)", zstophy,zout)
1756	c
1757	CALL histdef(nid_mth, "dqcon","Convection dQ","(kg/kg)/s",
1758	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1759	. "ave(X)", zstophy,zout)
1760	c
1761	CALL histdef(nid_mth, "dtlschr",
1762	$ "Large-scale condensational heating rate", "K/s",
1763	$ iim,jj_nb
1764	$ ,nhori, klev,1,klev,nvert, 32,"ave(X)", zstophy,zout)
1765	c
1766	CALL histdef(nid_mth, "dqlsc","Condensation dQ","(kg/kg)/s",
1767	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1768	. "ave(X)", zstophy,zout)
1769	c
1770	CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s",
1771	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1772	. "ave(X)", zstophy,zout)
1773	c
1774	CALL histdef(nid_mth, "dqvdf","Boundary-layer dQ","(kg/kg)/s",
1775	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1776	. "ave(X)", zstophy,zout)
1777	c
1778	CALL histdef(nid_mth, "dteva", "Reevaporation dT", "K/s",
1779	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1780	. "ave(X)", zstophy,zout)
1781	c
1782	CALL histdef(nid_mth,"dqeva","Reevaporation dQ","(kg/kg)/s",
1783	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1784	. "ave(X)", zstophy,zout)
1785
1786	CALL histdef(nid_mth, "ptconv", "POINTS CONVECTIFS"," ",
1787	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1788	. "ave(X)", zstophy,zout)
1789
1790	CALL histdef(nid_mth, "ratqs", "RATQS"," ",
1791	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1792	. "ave(X)", zstophy,zout)
1793
1794	c
1795	CALL histdef(nid_mth, "dqajs","Dry adjust. dQ","(kg/kg)/s",
1796	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1797	. "ave(X)", zstophy,zout)
1798	c
1799	CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s",
1800	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1801	. "ave(X)", zstophy,zout)
1802	c
1803	CALL histdef(nid_mth, "dtsw0", "CS SW radiation dT", "K/s",
1804	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1805	. "ave(X)", zstophy,zout)
1806	c
1807	CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
1808	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1809	. "ave(X)", zstophy,zout)
1810	c
1811	CALL histdef(nid_mth,"dtlw0","CS LW radiation dT","K/s",
1812	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1813	. "ave(X)", zstophy,zout)
1814	c
1815	CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
1816	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1817	. "ave(X)", zstophy,zout)
1818	c
1819	CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
1820	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1821	. "ave(X)", zstophy,zout)
1822	c
1823	CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
1824	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1825	. "ave(X)", zstophy,zout)
1826	c
1827	IF (ok_orodr) THEN
1828	CALL histdef(nid_mth, "duoro", "Orography dU", "m/s2",
1829	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1830	. "ave(X)", zstophy,zout)
1831	c
1832	CALL histdef(nid_mth, "dvoro", "Orography dV", "m/s2",
1833	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1834	. "ave(X)", zstophy,zout)
1835	c
1836	ENDIF
1837	C
1838	IF (ok_orolf) THEN
1839	CALL histdef(nid_mth, "dulif", "Orography dU", "m/s2",
1840	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1841	. "ave(X)", zstophy,zout)
1842	c
1843	CALL histdef(nid_mth, "dvlif", "Orography dV", "m/s2",
1844	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1845	. "ave(X)", zstophy,zout)
1846	ENDIF
1847	C
1848	#if defined(INCA_AER) && defined(CPP_COUPLE)
1849	CALL histdef(nid_mth, "tau", "cloud opt thickness", "-",
1850	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1851	. "ave(X)", zstophy,zout)
1852	CALL histdef(nid_mth, "CCM1", "CCM1", "",
1853	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1854	. "ave(X)", zstophy,zout)
1855	CALL histdef(nid_mth, "CCM2", "CCM2", "",
1856	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1857	. "ave(X)", zstophy,zout)
1858	CALL histdef(nid_mth, "taupi", "cloud opt thickn. (pi)", "-",
1859	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1860	. "ave(X)", zstophy,zout)
1861	CALL histdef(nid_mth, "cldemi", "emissivite IR", "-",
1862	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1863	. "ave(X)", zstophy,zout)
1864	CALL histdef(nid_mth, "SWTOAAS_ZERO", "?", "?",
1865	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1866	. "ave(X)", zstophy,zout)
1867	CALL histdef(nid_mth, "SWTOAAS_AER", "?", "?",
1868	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1869	. "ave(X)", zstophy,zout)
1870	CALL histdef(nid_mth, "SWTOAAS_NAT", "?", "?",
1871	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1872	. "ave(X)", zstophy,zout)
1873	CALL histdef(nid_mth, "SWTOAAS_AD", "?", "?",
1874	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1875	. "ave(X)", zstophy,zout)
1876	CALL histdef(nid_mth, "SWTOAAS_AI", "?", "?",
1877	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1878	. "ave(X)", zstophy,zout)
1879	CALL histdef(nid_mth, "SWTOACS_ZERO", "?", "?",
1880	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1881	. "ave(X)", zstophy,zout)
1882	CALL histdef(nid_mth, "SWTOACS_AER", "?", "?",
1883	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1884	. "ave(X)", zstophy,zout)
1885	CALL histdef(nid_mth, "SWTOACS_AD", "?", "?",
1886	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1887	. "ave(X)", zstophy,zout)
1888	CALL histdef(nid_mth, "SWTOACS_NAT", "?", "?",
1889	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1890	. "ave(X)", zstophy,zout)
1891	#endif
1892	#ifndef INCA
1893	if (nqmax.GE.3) THEN
1894	DO iq=3,nqmax
1895	iiq=niadv(iq)
1896	CALL histdef(nid_mth, tnom(iq), ttext(iiq), "-",
1897	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1898	. "ave(X)", zstophy,zout)
1899	ENDDO
1900	ENDIF
1901	#endif
1902
1903	ENDIF !lev_histmth.GE.4
1904	c
1905	ENDIF !type_run.EQ."CLIM".OR.type_run.EQ."ENSP")
1906	c
1907	CALL histend(nid_mth)
1908	c
1909	ndex2d = 0
1910	ndex3d = 0
1911	c
1912	c$OMP END MASTER
1913	ENDIF ! fin de test sur ok_mensuel

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

Download in other formats: