Context Navigation

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

Last change on this file since 892 was 890, checked in by Laurent Fairhead, 17 years ago
Mise sous flag ok_ISCCP d'une sortie d'une variable dependant du simulateur ISCCP Correction coquille sur des sorties sous-surfaces SD/LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 74.8 KB

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

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

Download in other formats: