Context Navigation

ini_histmth.h @ 4588

Last change on this file since 4588 was 918, checked in by Laurent Fairhead, 17 years ago

correction du bug ISCCP (n'ecrire ptop que quand ISCCP est appelé)
petite inversion de boucle dans isccp_cloud_types.F pour aller + vite
"CFisation" d'un certain nombre d'unités pour les hist*
les suggestions de JL pour rugoro

SD
LF

Property svn:eol-style set to native
Property svn:keywords set to Author Date Id Revision

File size: 74.2 KB

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

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

Context Navigation

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

Download in other formats: