Context Navigation

ini_histmth.h @ 759

Last change on this file since 759 was 740, checked in by lmdzadmin, 18 years ago
Correction bogues: les ecrit_ sont des REALs lus dans conf_phys.F90 en nombre de jours sauf pour ecrit_ins et ecrit_tra en secondes! Les ecrit_ sont initialises dans conf_phys.F90 et peuvent etre modifies dans physiq.def. IM
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 68.2 KB

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

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

Context Navigation

source: LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histmth.h @ 759

Download in other formats: