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

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

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