Context Navigation

source: LMDZ4/trunk/libf/phytherm/ini_histmth.h @ 1072

Last change on this file since 1072 was 864, checked in by lsce, 17 years ago
ACo : ajout clef CPP (COUPLE + AER) pour le modele couple avec les aerosol (IPSL_ESM_V1)
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	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 CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa",
295	c . iim,jj_nb,nhori, 1,1,1, -99, 32,
296	c . "ave(X)", zstophy,zout)
297	c
298	c CALL histdef(nid_mth, "frtv", "Meridional wind stress", "Pa",
299	c . iim,jj_nb,nhori, 1,1,1, -99, 32,
300	c . "ave(X)", zstophy,zout)
301	c
302	CALL histdef(nid_mth, "ffonte","Thermal flux for snow melting",
303	. "W/m2",iim,jj_nb,nhori, 1,1,1, -99, 32,
304	. "ave(X)", zstophy,zout)
305	c
306	CALL histdef(nid_mth, "fqcalving","Ice Calving",
307	. "kg/m2/s",iim,jj_nb,nhori, 1,1,1, -99, 32,
308	. "ave(X)", zstophy,zout)
309	c
310	CALL histdef(nid_mth, "fqfonte","Land ice melt",
311	. "kg/m2/s",iim,jj_nb,nhori, 1,1,1, -99, 32,
312	. "ave(X)", zstophy,zout)
313
314	DO nsrf = 1, nbsrf
315	call histdef(nid_mth, "taux_"//clnsurf(nsrf),
316	$ "Zonal wind stress"//clnsurf(nsrf), "Pa",
317	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
318	$ "ave(X)", zstophy,zout)
319
320	call histdef(nid_mth, "tauy_"//clnsurf(nsrf),
321	$ "Meridional wind stress "//clnsurf(nsrf), "Pa",
322	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
323	$ "ave(X)", zstophy,zout)
324	ENDDO
325	c
326	DO nsrf = 1, nbsrf
327	C
328	IF(nsrf.GT.2) THEN
329	call histdef(nid_mth, "pourc_"//clnsurf(nsrf),
330	$ "% "//clnsurf(nsrf), "%",
331	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
332	$ "ave(X)", zstophy,zout)
333	c
334	call histdef(nid_mth, "fract_"//clnsurf(nsrf),
335	$ "Fraction "//clnsurf(nsrf), "1",
336	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
337	$ "ave(X)", zstophy,zout)
338	ENDIF !nsrf.GT.2
339	C
340	call histdef(nid_mth, "tsol_"//clnsurf(nsrf),
341	$ "Temperature "//clnsurf(nsrf), "K",
342	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
343	$ "ave(X)", zstophy,zout)
344	C
345	call histdef(nid_mth, "sens_"//clnsurf(nsrf),
346	$ "Sensible heat flux "//clnsurf(nsrf), "W/m2",
347	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
348	$ "ave(X)", zstophy,zout)
349	c
350	call histdef(nid_mth, "lat_"//clnsurf(nsrf),
351	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
352	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
353	$ "ave(X)", zstophy,zout)
354	C
355	call histdef(nid_mth, "flw_"//clnsurf(nsrf),
356	$ "LW "//clnsurf(nsrf), "W/m2",
357	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
358	$ "ave(X)", zstophy,zout)
359	c
360	call histdef(nid_mth, "fsw_"//clnsurf(nsrf),
361	$ "SW "//clnsurf(nsrf), "W/m2",
362	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
363	$ "ave(X)", zstophy,zout)
364	C
365	call histdef(nid_mth, "wbils_"//clnsurf(nsrf),
366	$ "Bilan sol "//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, "wbilo_"//clnsurf(nsrf),
371	$ "Bilan eau "//clnsurf(nsrf), "kg/(m2*s)",
372	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
373	$ "ave(X)", zstophy,zout)
374	C
375	cIM 250105 BEG
376	IF(1.EQ.0) THEN
377	C
378	call histdef(nid_mth, "evapour_"//clnsurf(nsrf),
379	$ "Evappourc "//clnsurf(nsrf), "kg/(sm2)",
380	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
381	$ "ave(X)", zstophy,zout)
382	C
383	call histdef(nid_mth, "prepour_"//clnsurf(nsrf),
384	$ "Precippourc "//clnsurf(nsrf), "kg/(sm2)",
385	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
386	$ "ave(X)", zstophy,zout)
387	C
388	ENDIF !(1.EQ.0) THEN
389	cIM 250105 END
390	C
391	END DO
392	c
393	CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
394	. iim,jj_nb,nhori, 1,1,1, -99, 32,
395	. "ave(X)", zstophy,zout)
396	c
397	CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
398	. iim,jj_nb,nhori, 1,1,1, -99, 32,
399	. "ave(X)", zstophy,zout)
400	c
401	CALL histdef(nid_mth, "cldl", "Low-level cloudiness", "-",
402	. iim,jj_nb,nhori, 1,1,1, -99, 32,
403	. "ave(X)", zstophy,zout)
404	c
405	CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", "-",
406	. iim,jj_nb,nhori, 1,1,1, -99, 32,
407	. "ave(X)", zstophy,zout)
408	c
409	CALL histdef(nid_mth, "cldh", "High-level cloudiness", "-",
410	. iim,jj_nb,nhori, 1,1,1, -99, 32,
411	. "ave(X)", zstophy,zout)
412	c
413	CALL histdef(nid_mth, "cldt", "Total cloudiness", "-",
414	. iim,jj_nb,nhori, 1,1,1, -99, 32,
415	. "ave(X)", zstophy,zout)
416	c
417	CALL histdef(nid_mth,"cldq","Cloud liquid water path","kg/m2",
418	. iim,jj_nb,nhori, 1,1,1, -99, 32,
419	. "ave(X)", zstophy,zout)
420	c
421	CALL histdef(nid_mth,"lwp","Cloud water path","kg/m2",
422	. iim,jj_nb,nhori, 1,1,1, -99, 32,
423	. "ave(X)", zstophy,zout)
424	c
425	CALL histdef(nid_mth,"iwp","Cloud ice 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, "ue", "Zonal energy transport", "-",
430	. iim,jj_nb,nhori, 1,1,1, -99, 32,
431	. "ave(X)", zstophy,zout)
432	c
433	CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
434	. iim,jj_nb,nhori, 1,1,1, -99, 32,
435	. "ave(X)", zstophy,zout)
436	c
437	CALL histdef(nid_mth, "uq", "Zonal humidity transport", "-",
438	. iim,jj_nb,nhori, 1,1,1, -99, 32,
439	. "ave(X)", zstophy,zout)
440	c
441	CALL histdef(nid_mth, "vq", "Merid humidity transport", "-",
442	. iim,jj_nb,nhori, 1,1,1, -99, 32,
443	. "ave(X)", zstophy,zout)
444	cKE43
445	IF(iflag_con.GE.3) THEN ! sb
446	c
447	CALL histdef(nid_mth, "cape", "Conv avlbl pot ener", "J/kg",
448	. iim,jj_nb,nhori, 1,1,1, -99, 32,
449	. "ave(X)", zstophy,zout)
450	c
451	CALL histdef(nid_mth, "pbase", "Cld base pressure", "mb",
452	. iim,jj_nb,nhori, 1,1,1, -99, 32,
453	. "ave(X)", zstophy,zout)
454	c
455	CALL histdef(nid_mth, "ptop", "Cld top pressure", "mb",
456	. iim,jj_nb,nhori, 1,1,1, -99, 32,
457	. "ave(X)", zstophy,zout)
458	c
459	CALL histdef(nid_mth, "fbase", "Cld base mass flux", "kg/m2/s",
460	. iim,jj_nb,nhori, 1,1,1, -99, 32,
461	. "ave(X)", zstophy,zout)
462	c
463	CALL histdef(nid_mth, "prw", "Precipitable water", "kg/m2",
464	. iim,jj_nb,nhori, 1,1,1, -99, 32,
465	. "ave(X)", zstophy,zout)
466	c
467	ENDIF !iflag_con .GE. 3
468	c34EK
469	cIM cf. AM 081204 BEG
470	c34EK
471	c HBTM2
472	CALL histdef(nid_mth, "s_pblh", "Boundary Layer Height", "m",
473	. iim,jj_nb,nhori, 1,1,1, -99, 32,
474	. "ave(X)", zstophy,zout)
475	CALL histdef(nid_mth, "s_pblt", "t at Boundary Layer Height",
476	. "K",
477	. iim,jj_nb,nhori, 1,1,1, -99, 32,
478	. "ave(X)", zstophy,zout)
479	CALL histdef(nid_mth, "s_lcl", "Condensation level", "m",
480	. iim,jj_nb,nhori, 1,1,1, -99, 32,
481	. "ave(X)", zstophy,zout)
482	CALL histdef(nid_mth, "s_capCL", "Conv avlbl pot ener
483	. for ABL", "J/m2",
484	. iim,jj_nb,nhori, 1,1,1, -99, 32,
485	. "ave(X)", zstophy,zout)
486	CALL histdef(nid_mth, "s_oliqCL", "Liq Water in BL", "kg/m2",
487	. iim,jj_nb,nhori, 1,1,1, -99, 32,
488	. "ave(X)", zstophy,zout)
489	CALL histdef(nid_mth, "s_cteiCL", "Instability criteria
490	. (ABL)", "K",
491	. iim,jj_nb,nhori, 1,1,1, -99, 32,
492	. "ave(X)", zstophy,zout)
493	CALL histdef(nid_mth, "s_therm", "Exces du thermique", "K",
494	. iim,jj_nb,nhori, 1,1,1, -99, 32,
495	. "ave(X)", zstophy,zout)
496	CALL histdef(nid_mth, "s_trmb1", "deep_cape(HBTM2)", "J/m2",
497	. iim,jj_nb,nhori, 1,1,1, -99, 32,
498	. "ave(X)", zstophy,zout)
499	CALL histdef(nid_mth, "s_trmb2", "inhibition (HBTM2)", "J/m2",
500	. iim,jj_nb,nhori, 1,1,1, -99, 32,
501	. "ave(X)", zstophy,zout)
502	CALL histdef(nid_mth, "s_trmb3", "Point Omega (HBTM2)", "m",
503	. iim,jj_nb,nhori, 1,1,1, -99, 32,
504	. "ave(X)", zstophy,zout)
505	c
506	cIM cf. AM 081204 END
507	c
508	c Champs interpolles sur des niveaux de pression
509	DO k=1, nlevSTD
510	c
511	c bb=clevSTD(k)
512	c
513	c IF(k.GE.2) THEN
514	IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
515	IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
516	c aa=clevSTD(k)
517	c bb=aa(1:lnblnk1(aa))
518	c ENDIF
519	c
520	IF(bb2.EQ."850".OR.bb2.EQ."700".OR.
521	$ bb2.EQ."500".OR.bb2.EQ."200") THEN
522	c
523	CALL histdef(nid_mth, "u"//bb2,
524	. "Zonal wind "//bb2//"mb","m/s",
525	. iim,jj_nb,nhori, 1,1,1, -99, 32,
526	. "inst(X)", zout,zout)
527	c
528	CALL histdef(nid_mth, "v"//bb2,
529	. "Meridional wind "//bb2//"mb","m/s",
530	. iim,jj_nb,nhori, 1,1,1, -99, 32,
531	. "inst(X)", zout,zout)
532	c
533	CALL histdef(nid_mth, "w"//bb2,
534	. "Vertical wind "//bb2//"mb","m/s",
535	. iim,jj_nb,nhori, 1,1,1, -99, 32,
536	. "inst(X)", zout,zout)
537	c
538	CALL histdef(nid_mth, "phi"//bb2,
539	. "Geopotential "//bb2//"mb","m",
540	. iim,jj_nb,nhori, 1,1,1, -99, 32,
541	. "inst(X)", zout,zout)
542	c
543	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.
544	ENDDO
545	c
546	cIM diagnostiques flux ocean-atm ou ocean-glace de mer
547	cIM pour utilisation dans un modele de "slab" ocean
548	c
549	CALL histdef(nid_mth, "fluxo",
550	. "Flux turbulents ocean-atmosphere", "W/m2",
551	. iim,jj_nb,nhori, 1,1,1, -99, 32,
552	. "ave(X)", zstophy,zout)
553	c
554	CALL histdef(nid_mth, "fluxg",
555	. "Flux turbulents ocean-glace de mer","W/m2",
556	. iim,jj_nb,nhori, 1,1,1, -99, 32,
557	. "ave(X)", zstophy,zout)
558	c
559	CALL histdef(nid_mth, "t_oce_sic",
560	. "Temp. mixte oce-sic","K",
561	. iim,jj_nb,nhori, 1,1,1, -99, 32,
562	. "ave(X)", zstophy,zout)
563	c
564	IF (OCEAN.EQ.'force ') THEN
565	c
566	CALL histdef(nid_mth, "lmt_bils",
567	. "Bilan au sol atmosphere forcee", "W/m2",
568	. iim,jj_nb,nhori, 1,1,1, -99, 32,
569	. "ave(X)", zstophy,zout)
570	c
571	ELSE IF (OCEAN.EQ.'slab ') THEN
572	c
573	CALL histdef(nid_mth, "slab_bils",
574	. "Bilan au sol Slab", "W/m2",
575	. iim,jj_nb,nhori, 1,1,1, -99, 32,
576	. "ave(X)", zstophy,zout)
577	c
578	CALL histdef(nid_mth, "tslab", "Slab SST ", "K",
579	. iim,jj_nb,nhori, 1,1,1, -99, 32,
580	. "ave(X)", zstophy,zout)
581	c
582	CALL histdef(nid_mth, "seaice", "Slab seaice", "kg/m2",
583	. iim,jj_nb,nhori, 1,1,1, -99, 32,
584	. "ave(X)", zstophy,zout)
585	c
586	CALL histdef(nid_mth, "siceh", "Slab seaice height", "m",
587	. iim,jj_nb,nhori, 1,1,1, -99, 32,
588	. "ave(X)", zstophy,zout)
589	c
590	ENDIF
591	c
592	ENDIF !lev_histmth.GE.1
593	c
594	IF(lev_histmth.GE.2) THEN
595	c
596	c Champs 3D:
597	c
598	c
599	CALL histdef(nid_mth,"lwcon","Cloud liquid water content",
600	. "kg/kg",
601	. iim,jj_nb,nhori, klev,1,klev, nvert, 32,
602	. "ave(X)", zstophy,zout)
603	c
604	CALL histdef(nid_mth,"iwcon","Cloud ice water content","kg/kg",
605	. iim,jj_nb,nhori, klev,1,klev, nvert, 32,
606	. "ave(X)", zstophy,zout)
607	c
608	CALL histdef(nid_mth, "temp", "Air temperature", "K",
609	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
610	. "ave(X)", zstophy,zout)
611	c
612	CALL histdef(nid_mth, "ovap", "Specific humidity", "kg/kg",
613	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
614	. "ave(X)", zstophy,zout)
615	c
616	c CALL histdef(nid_mth,"wvap","Water vapor mixing ratio","kg/kg",
617	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
618	c . "ave(X)", zstophy,zout)
619	c
620	CALL histdef(nid_mth, "geop", "Geopotential height", "m2/s2",
621	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
622	. "ave(X)", zstophy,zout)
623	c
624	CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s",
625	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
626	. "ave(X)", zstophy,zout)
627	c
628	CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s",
629	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
630	. "ave(X)", zstophy,zout)
631	c
632	CALL histdef(nid_mth, "vitw", "Vertical wind", "Pa/s",
633	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
634	. "ave(X)", zstophy,zout)
635	c
636	CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
637	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
638	. "ave(X)", zstophy,zout)
639	c
640	CALL histdef(nid_mth, "rneb", "Cloud fraction", "-",
641	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
642	. "ave(X)", zstophy,zout)
643	c
644	CALL histdef(nid_mth, "rnebcon", "Convective Cloud Fraction"
645	. , "-",
646	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
647	. "ave(X)", zstophy,zout)
648	c
649	CALL histdef(nid_mth, "rhum", "Relative humidity", "-",
650	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
651	. "ave(X)", zstophy,zout)
652	c
653	CALL histdef(nid_mth, "ozone", "Ozone concentration", "ppmv",
654	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
655	. "ave(X)", zstophy,zout)
656	c
657	CALL histdef(nid_mth, "upwd", "saturated updraft", "kg/m2/s",
658	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
659	. "ave(X)", zstophy,zout)
660	c
661	CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
662	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
663	. "ave(X)", zstophy,zout)
664	c
665	CALL histdef(nid_mth, "dqphy", "Physics dQ", "(kg/kg)/s",
666	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
667	. "ave(X)", zstophy,zout)
668	c
669	ENDIF !lev_histmth.GE.2
670	c
671	IF(lev_histmth.GE.3) THEN
672	c
673	DO nsrf=1, nbsrf
674	c
675	call histdef(nid_mth, "albe_"//clnsurf(nsrf),
676	$ "Albedo surf. "//clnsurf(nsrf), "-",
677	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
678	$ "ave(X)", zstophy,zout)
679	c
680	call histdef(nid_mth, "rugs_"//clnsurf(nsrf),
681	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
682	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
683	$ "ave(X)", zstophy,zout)
684	c
685	CALL histdef(nid_mth, "ages_"//clnsurf(nsrf), "Snow age","day",
686	. iim,jj_nb,nhori, 1,1,1, -99, 32,
687	. "ave(X)", zstophy,zout)
688	c
689	ENDDO !nsrf=1, nbsrf
690	c
691	CALL histdef(nid_mth, "albs", "Surface albedo", "-",
692	. iim,jj_nb,nhori, 1,1,1, -99, 32,
693	. "ave(X)", zstophy,zout)
694	CALL histdef(nid_mth, "albslw", "Surface albedo LW", "-",
695	. iim,jj_nb,nhori, 1,1,1, -99, 32,
696	. "ave(X)", zstophy,zout)
697	c
698	ENDIF !lev_histmth.GE.3
699	c
700	IF(lev_histmth.GE.4) THEN
701	c
702	!FH Sorties pour la couche limite
703	if (iflag_pbl>1) then
704	CALL histdef(nid_mth, "tke","TKE","m2/s2",
705	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
706	. "ave(X)", zstophy,zout)
707	c
708	CALL histdef(nid_mth, "tke_max","TKE max","m2/s2",
709	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
710	. "t_max(X)", zstophy,zout)
711	endif
712	c
713	CALL histdef(nid_mth, "kz","Kz melange","m2/s",
714	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
715	. "ave(X)", zstophy,zout)
716	c
717	CALL histdef(nid_mth, "kz_max","Kz melange max","m2/s",
718	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
719	. "t_max(X)", zstophy,zout)
720
721
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
753
754	c
755	c CALL histdef(nid_mth, "ducon", "Convection du", "m/s2",
756	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
757	c . "ave(X)", zstophy,zout)
758	c
759	CALL histdef(nid_mth, "dqcon", "Convection dQ", "(kg/kg)/s",
760	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
761	. "ave(X)", zstophy,zout)
762	c
763	CALL histdef(nid_mth, "dtlsc", "Condensation dT", "K/s",
764	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
765	. "ave(X)", zstophy,zout)
766	c
767	CALL histdef(nid_mth, "dtlschr",
768	$ "Large-scale condensational heating rate", "K/s",iim,jj_nb
769	$ ,nhori, klev,1,klev,nvert, 32,"ave(X)", zstophy,zout)
770	c
771	CALL histdef(nid_mth, "dqlsc", "Condensation dQ", "(kg/kg)/s",
772	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
773	. "ave(X)", zstophy,zout)
774	c
775	CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s",
776	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
777	. "ave(X)", zstophy,zout)
778	c
779	CALL histdef(nid_mth, "dqvdf","Boundary-layer 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, "dteva", "Reevaporation 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, "dqeva","Reevaporation dQ","(kg/kg)/s",
788	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
789	. "ave(X)", zstophy,zout)
790
791	CALL histdef(nid_mth, "ptconv", "POINTS CONVECTIFS"," ",
792	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
793	. "ave(X)", zstophy,zout)
794
795	CALL histdef(nid_mth, "ratqs", "RATQS"," ",
796	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
797	. "ave(X)", zstophy,zout)
798
799	c
800	CALL histdef(nid_mth, "dtthe", "Dry adjust. dT", "K/s",
801	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
802	. "ave(X)", zstophy,zout)
803
804	CALL histdef(nid_mth,"dqthe","Dry adjust. dQ","(kg/kg)/s",
805	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
806	. "ave(X)", zstophy,zout)
807	c
808	CALL histdef(nid_mth, "dtajs", "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,"dqajs","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, "dtswr", "SW radiation dT", "K/s",
817	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
818	. "ave(X)", zstophy,zout)
819	c
820	CALL histdef(nid_mth, "dtsw0", "CS SW radiation dT", "K/s",
821	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
822	. "ave(X)", zstophy,zout)
823	c
824	CALL histdef(nid_mth, "dtlwr", "LW 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,"dtlw0","CS LW 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, "dtec", "Cinetic dissip 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, "duvdf", "Boundary-layer dU", "m/s2",
837	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
838	. "ave(X)", zstophy,zout)
839	c
840	CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
841	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
842	. "ave(X)", zstophy,zout)
843	c
844	IF (ok_orodr) THEN
845	CALL histdef(nid_mth, "duoro", "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, "dvoro", "Orography dV", "m/s2",
850	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
851	. "ave(X)", zstophy,zout)
852	c
853	ENDIF
854	C
855	IF (ok_orolf) THEN
856	CALL histdef(nid_mth, "dulif", "Orography dU", "m/s2",
857	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
858	. "ave(X)", zstophy,zout)
859	c
860	CALL histdef(nid_mth, "dvlif", "Orography dV", "m/s2",
861	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
862	. "ave(X)", zstophy,zout)
863	ENDIF
864	C
865	#if defined(INCA_AER) && defined(CPP_COUPLE)
866	CALL histdef(nid_mth, "tau", "cloud opt thickness", "-",
867	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
868	. "ave(X)", zstophy,zout)
869	CALL histdef(nid_mth, "CCM1", "CCM1", "",
870	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
871	. "ave(X)", zstophy,zout)
872	CALL histdef(nid_mth, "CCM2", "CCM2", "",
873	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
874	. "ave(X)", zstophy,zout)
875	CALL histdef(nid_mth, "taupi", "cloud opt thickn. (pi)", "-",
876	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
877	. "ave(X)", zstophy,zout)
878	CALL histdef(nid_mth, "cldemi", "emissivite IR", "-",
879	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
880	. "ave(X)", zstophy,zout)
881	CALL histdef(nid_mth, "SWTOAAS_ZERO", "?", "?",
882	. iim,jj_nb,nhori, 1,1,1, -99, 32,
883	. "ave(X)", zstophy,zout)
884	CALL histdef(nid_mth, "SWTOAAS_AER", "?", "?",
885	. iim,jj_nb,nhori, 1,1,1, -99, 32,
886	. "ave(X)", zstophy,zout)
887	CALL histdef(nid_mth, "SWTOAAS_NAT", "?", "?",
888	. iim,jj_nb,nhori, 1,1,1, -99, 32,
889	. "ave(X)", zstophy,zout)
890	CALL histdef(nid_mth, "SWTOAAS_AD", "?", "?",
891	. iim,jj_nb,nhori, 1,1,1, -99, 32,
892	. "ave(X)", zstophy,zout)
893	CALL histdef(nid_mth, "SWTOAAS_AI", "?", "?",
894	. iim,jj_nb,nhori, 1,1,1, -99, 32,
895	. "ave(X)", zstophy,zout)
896	CALL histdef(nid_mth, "SWTOACS_ZERO", "?", "?",
897	. iim,jj_nb,nhori, 1,1,1, -99, 32,
898	. "ave(X)", zstophy,zout)
899	CALL histdef(nid_mth, "SWTOACS_AER", "?", "?",
900	. iim,jj_nb,nhori, 1,1,1, -99, 32,
901	. "ave(X)", zstophy,zout)
902	CALL histdef(nid_mth, "SWTOACS_AD", "?", "?",
903	. iim,jj_nb,nhori, 1,1,1, -99, 32,
904	. "ave(X)", zstophy,zout)
905	CALL histdef(nid_mth, "SWTOACS_NAT", "?", "?",
906	. iim,jj_nb,nhori, 1,1,1, -99, 32,
907	. "ave(X)", zstophy,zout)
908	#endif
909	#ifndef INCA
910	if (nqmax.GE.3) THEN
911	DO iq=3,nqmax
912	iiq=niadv(iq)
913	CALL histdef(nid_mth, tnom(iq), ttext(iiq), "-",
914	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
915	. "ave(X)", zstophy,zout)
916	ENDDO
917	ENDIF
918	#endif
919	c
920	ENDIF !lev_histmth.GE.4
921	c
922	ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
923	c
924	IF(lev_histmth.GE.1) THEN
925	c
926	CALL histdef(nid_mth, "phis", "Surface geop. height", "m",
927	. iim,jj_nb,nhori, 1,1,1, -99, 32,
928	. "once", zstophy,zout)
929	c
930	CALL histdef(nid_mth, "aire", "Grid area", "-",
931	. iim,jj_nb,nhori, 1,1,1, -99, 32,
932	. "once", zstophy,zout)
933	c
934	CALL histdef(nid_mth, "pourc_"//clnsurf(is_ter),
935	. "% "//clnsurf(is_ter), "%",
936	. iim,jj_nb,nhori, 1,1,1, -99, 32,
937	. "once", zstophy,zout)
938	c
939	CALL histdef(nid_mth, "fract_"//clnsurf(is_ter),
940	. "Fraction "//clnsurf(is_ter), "1",
941	. iim,jj_nb,nhori, 1,1,1, -99, 32,
942	. "once", zstophy,zout)
943	c
944	CALL histdef(nid_mth, "pourc_"//clnsurf(is_lic),
945	. "% "//clnsurf(is_lic), "%",
946	. iim,jj_nb,nhori, 1,1,1, -99, 32,
947	. "once", zstophy,zout)
948	c
949	CALL histdef(nid_mth, "fract_"//clnsurf(is_lic),
950	. "Fraction "//clnsurf(is_lic), "1",
951	. iim,jj_nb,nhori, 1,1,1, -99, 32,
952	. "once", zstophy,zout)
953	c
954	CALL histdef(nid_mth, "pourc_"//clnsurf(is_sic),
955	. "% "//clnsurf(is_sic), "%",
956	. iim,jj_nb,nhori, 1,1,1, -99, 32,
957	. "ave(X)", zstophy,zout)
958	c
959	CALL histdef(nid_mth, "fract_"//clnsurf(is_sic),
960	. "Fraction "//clnsurf(is_sic), "1",
961	. iim,jj_nb,nhori, 1,1,1, -99, 32,
962	. "ave(X)", zstophy,zout)
963	c
964	CALL histdef(nid_mth, "flat", "Latent heat flux", "W/m2",
965	. iim,jj_nb,nhori, 1,1,1, -99, 32,
966	. "ave(X)", zstophy,zout)
967	c
968	CALL histdef(nid_mth, "bilTOA",
969	. "Net radiation at TOA", "W/m2",
970	. iim,jj_nb,nhori, 1,1,1, -99, 32,
971	. "ave(X)", zstophy,zout)
972	c
973	c
974	c CALL histdef(nid_mth, "tops01",
975	c . "TOA reflected clear-sky shortwave
976	c . radiation (method 1)", "W/m2",
977	c . iim,jj_nb,nhori, 1,1,1, -99, 32,
978	c . "ave(X)", zstophy,zout)
979	c
980	CALL histdef(nid_mth, "dtdia","dtdia","K/s",
981	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
982	. "ave(X)", zstophy,zout)
983	c
984	CALL histdef(nid_mth, "dtsw","dtsw","K/s",
985	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
986	. "ave(X)", zstophy,zout)
987	c
988	CALL histdef(nid_mth, "dtlw","dtlw","K/s",
989	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
990	. "ave(X)", zstophy,zout)
991	c
992	CALL histdef(nid_mth, "dtcon","dtcon","K/s",
993	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
994	. "ave(X)", zstophy,zout)
995	c
996	CALL histdef(nid_mth, "dtthe", "Dry adjust. dT", "K/s",
997	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
998	. "ave(X)", zstophy,zout)
999
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	CALL histdef(nid_mth, "t2m_min", "Temp. 2m min.",
1086	. "K",
1087	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1088	. t2mincels, zstophy,zout)
1089	c
1090	CALL histdef(nid_mth, "t2m_max", "Temp. 2m max.",
1091	. "K",
1092	. iim,jj_nb,nhori, 1,1,1, -99, 32,
1093	. 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: