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

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

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