Context Navigation

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

Last change on this file since 905 was 905, checked in by Laurent Fairhead, 16 years ago
Comme quoi faut faire gaffe quand on veut reintegrer des vieilles modifs MAFO LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 76.3 KB

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

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

Download in other formats: