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

Last change on this file since 892 was 892, checked in by Laurent Fairhead, 16 years ago

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