Context Navigation

source: LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/ini_histmth.h @ 897

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

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

Download in other formats: