source: LMDZ4/branches/V3_test/libf/phylmd/ini_histmth.h @ 4462

Last change on this file since 4462 was 735, checked in by Laurent Fairhead, 18 years ago

Modifications version parallele
YM/LF
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 69.5 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,jjphy_nb,nhori, 1,1,1, -99, 32, 
43     .                "once",  zstophy,zout)
44c
45         CALL histdef(nid_mth, "aire", "Grid area", "-",
46     .                iim,jjphy_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,jjphy_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,jjphy_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,jjphy_nb,nhori, 1,1,1, -99, 32,
63     .                "ave(X)", zstophy,zout)
64c
65         CALL histdef(nid_mth, "tsol", "Surface Temperature", "K",
66     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
67     .                "ave(X)", zstophy,zout)
68c
69         CALL histdef(nid_mth, "t2m", "Temperature 2m", "K",
70     .                iim,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_nb,nhori, 1,1,1, -99, 32,
91     .                "ave(X)", zstophy,zout)
92c
93         CALL histdef(nid_mth, "sicf", "Sea-ice fraction", "-",
94     .                iim,jjphy_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,jjphy_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,jjphy_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,jjphy_nb,nhori, 1,1,1, -99, 32,
108     .                "ave(X)", zstophy,zout)
109c
110         CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa",
111     .                iim,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjmp1,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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_nb,nhori, 1,1,1, -99, 32,
312     .                "ave(X)", zsto,zout)
313
314         DO nsrf = 1, nbsrf
315           call histdef(nid_mth, "taux_"//clnsurf(nsrf),
316     $         "Zonal wind stress"//clnsurf(nsrf), "Pa", 
317     $         iim,jjphy_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,jjphy_nb,nhori, 1,1,1, -99, 32,
323     $         "ave(X)", zstophy,zout)
324         ENDDO
325c
326         DO nsrf = 1, nbsrf
327C
328          IF(nsrf.GT.2) THEN
329           call histdef(nid_mth, "pourc_"//clnsurf(nsrf),
330     $         "% "//clnsurf(nsrf), "%", 
331     $         iim,jjphy_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,jjphy_nb,nhori, 1,1,1, -99, 32,
337     $         "ave(X)", zstophy,zout)
338          ENDIF !nsrf.GT.2
339C
340           call histdef(nid_mth, "tsol_"//clnsurf(nsrf),
341     $         "Temperature "//clnsurf(nsrf), "K", 
342     $         iim,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_nb,nhori, 1,1,1, -99, 32, 
395     .                "ave(X)", zstophy,zout)
396c
397         CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
398     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
399     .                "ave(X)", zstophy,zout)
400c
401         CALL histdef(nid_mth, "cldl", "Low-level cloudiness", "-",
402     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
403     .                "ave(X)", zstophy,zout)
404c
405         CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", "-",
406     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
407     .                "ave(X)", zstophy,zout)
408c
409         CALL histdef(nid_mth, "cldh", "High-level cloudiness", "-",
410     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
411     .                "ave(X)", zstophy,zout)
412c
413         CALL histdef(nid_mth, "cldt", "Total cloudiness", "-",
414     .                iim,jjphy_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,jjphy_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,jjphy_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,jjphy_nb,nhori, 1,1,1, -99, 32, 
427     .                "ave(X)", zstophy,zout)
428c
429         CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
430     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
431     .                "ave(X)", zstophy,zout)
432c
433         CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
434     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
435     .                "ave(X)", zstophy,zout)
436c
437         CALL histdef(nid_mth, "uq", "Zonal humidity transport", "-",
438     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
439     .                "ave(X)", zstophy,zout)
440c
441         CALL histdef(nid_mth, "vq", "Merid humidity transport", "-",
442     .                iim,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_nb,nhori, 1,1,1, -99, 32,
478     .                "ave(X)", zstophy,zout)
479         CALL histdef(nid_mth, "s_lcl", "Condensation level", "m",
480     .                iim,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_nb,nhori, 1,1,1, -99, 32,
576     .              "ave(X)", zstophy,zout)
577c
578        CALL histdef(nid_mth, "tslab", "Slab SST ", "K",
579     .              iim,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_nb,nhori, klev,1,klev, nvert, 32, 
606     .                "ave(X)", zstophy,zout)
607c
608         CALL histdef(nid_mth, "temp", "Air temperature", "K",
609     .                iim,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_nb,nhori, klev,1,klev,nvert, 32,
630     .                "ave(X)", zstophy,zout)
631c
632         CALL histdef(nid_mth, "vitw", "Vertical wind", "m/s",
633     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
634     .                "ave(X)", zstophy,zout)
635c
636         CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
637     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
638     .                "ave(X)", zstophy,zout)
639c
640        CALL histdef(nid_mth, "rneb", "Cloud fraction", "-",
641     .                iim,jjphy_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,jjphy_nb,nhori, klev,1,klev,nvert, 32,
647     .                "ave(X)", zstophy,zout)
648c
649         CALL histdef(nid_mth, "rhum", "Relative humidity", "-",
650     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
651     .                "ave(X)", zstophy,zout)
652c
653         CALL histdef(nid_mth, "ozone", "Ozone concentration", "ppmv",
654     .                iim,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_nb,nhori, 1,1,1, -99, 32, 
693     .                "ave(X)", zstophy,zout)
694         CALL histdef(nid_mth, "albslw", "Surface albedo LW", "-",
695     .                iim,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_nb,nhori, klev,1,klev,nvert, 32,
766     .                "ave(X)", zstophy,zout)
767
768         CALL histdef(nid_mth, "ptconv", "POINTS CONVECTIFS"," ",
769     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
770     .                "ave(X)", zstophy,zout)
771
772         CALL histdef(nid_mth, "ratqs", "RATQS"," ",
773     .                iim,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_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,jjphy_nb,nhori, klev,1,klev,nvert, 32,
831     .                "ave(X)", zstophy,zout)
832         ENDIF
833C
834         if (nqmax.GE.3) THEN
835         DO iq=1,nqmax-2
836         IF (iq.LE.99) THEN
837         WRITE(str2,'(i2.2)') iq
838         CALL histdef(nid_mth, "trac"//str2, "Tracer No."//str2, "-",
839     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
840     .                "ave(X)", zstophy,zout)
841         ELSE
842         PRINT*, "Trop de traceurs"
843         CALL abort
844         ENDIF
845         ENDDO
846         ENDIF
847c
848      ENDIF !lev_histmth.GE.4
849c
850      ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
851c
852      IF(lev_histmth.GE.1) THEN
853c
854         CALL histdef(nid_mth, "phis", "Surface geop. height", "m",
855     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
856     .                "once",  zstophy,zout)
857c
858         CALL histdef(nid_mth, "aire", "Grid area", "-",
859     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
860     .                "once",  zstophy,zout)
861c
862         CALL histdef(nid_mth, "pourc_"//clnsurf(is_ter),
863     .         "% "//clnsurf(is_ter), "%", 
864     .         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
865     .         "once", zstophy,zout)
866c
867         CALL histdef(nid_mth, "fract_"//clnsurf(is_ter),
868     .         "Fraction "//clnsurf(is_ter), "1", 
869     .         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
870     .         "once", zstophy,zout)
871c
872         CALL histdef(nid_mth, "pourc_"//clnsurf(is_lic),
873     .         "% "//clnsurf(is_lic), "%", 
874     .         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
875     .         "once", zstophy,zout)
876c
877         CALL histdef(nid_mth, "fract_"//clnsurf(is_lic),
878     .         "Fraction "//clnsurf(is_lic), "1", 
879     .         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
880     .         "once", zstophy,zout)
881c
882         CALL histdef(nid_mth, "pourc_"//clnsurf(is_sic),
883     .         "% "//clnsurf(is_sic), "%",
884     .         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
885     .                "ave(X)", zstophy,zout)
886c
887         CALL histdef(nid_mth, "fract_"//clnsurf(is_sic),
888     .         "Fraction "//clnsurf(is_sic), "1",
889     .         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
890     .                "ave(X)", zstophy,zout)
891c
892         CALL histdef(nid_mth, "flat", "Latent heat flux", "W/m2",
893     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
894     .                "ave(X)", zstophy,zout)
895c
896         CALL histdef(nid_mth, "bilTOA",
897     .                "Net radiation at TOA", "W/m2",
898     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
899     .                "ave(X)", zstophy,zout)
900c
901c
902c        CALL histdef(nid_mth, "tops01",
903c    .                "TOA reflected clear-sky shortwave
904c    .                radiation (method 1)", "W/m2",
905c    .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
906c    .                "ave(X)", zstophy,zout)
907c
908          CALL histdef(nid_mth, "dtdia","dtdia","K/s",
909     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
910     .                "ave(X)", zstophy,zout)
911c
912          CALL histdef(nid_mth, "dtsw","dtsw","K/s",
913     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
914     .                "ave(X)", zstophy,zout)
915c
916          CALL histdef(nid_mth, "dtlw","dtlw","K/s",
917     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
918     .                "ave(X)", zstophy,zout)
919c
920          CALL histdef(nid_mth, "dtcon","dtcon","K/s",
921     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
922     .                "ave(X)", zstophy,zout)
923c
924          CALL histdef(nid_mth, "dtajs","dtajs","K/s",
925     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
926     .                "ave(X)", zstophy,zout)
927c
928          CALL histdef(nid_mth, "dtlsc","dtlsc","K/s",
929     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
930     .                "ave(X)", zstophy,zout)
931c
932          CALL histdef(nid_mth, "dqdia","dqdia","(kg/kg)/s",
933     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
934     .                "ave(X)", zstophy,zout)
935c
936          CALL histdef(nid_mth, "duoli",
937     .         "duoli","m/s2",
938     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
939     .                "ave(X)", zstophy,zout)
940c
941          CALL histdef(nid_mth, "dvoli",
942     .         "dvoli","m/s2",
943     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
944     .                "ave(X)", zstophy,zout)
945c
946          CALL histdef(nid_mth, "ducon",
947     .         "ducon","m/s2",
948     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
949     .                "ave(X)", zstophy,zout)
950c
951          CALL histdef(nid_mth, "dvcon",
952     .         "dvcon","m/s2",
953     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
954     .                "ave(X)", zstophy,zout)
955c
956          CALL histdef(nid_mth, "dutot",
957     .         "dutot","m/s2",
958     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
959     .                "ave(X)", zstophy,zout)
960c
961          CALL histdef(nid_mth, "dvtot",
962     .         "dvtot","m/s2",
963     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
964     .                "ave(X)", zstophy,zout)
965c
966c         CALL histdef(nid_mth, "cldtrue",
967c    .         "TRUE Cloud fraction ","%",
968c    .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
969c    .         "inst(X)", zout,zout)
970c
971c         CALL histdef(nid_mth, "cldTOA",
972c    .         "Cloud amount (satellite view)","%",
973c    .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
974c    .         "ave(X)", zstophy,zout)
975c
976c         CALL histdef(nid_mth, "cldSFC",
977c    .         "Cloud amount (surface view)","%",
978c    .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
979c    .         "ave(X)", zstophy,zout)
980c
981          CALL histdef(nid_mth, "GWDtaux",
982     .         "GWD induced Ewd wind stress","Pa",
983     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
984     .         "ave(X)", zstophy,zout)
985c
986          CALL histdef(nid_mth, "GWDtauy",
987     .         "GWD induced Nwd wind stress","Pa",
988     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
989     .         "ave(X)", zstophy,zout)
990c
991c Champs 2D:
992c
993         CALL histdef(nid_mth, "slp", "Sea Level Pressure", "Pa",
994     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
995     .                "ave(X)", zstorad,zout)
996c
997         CALL histdef(nid_mth, "tsol", "Surface Temperature", "K",
998     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
999     .                "ave(X)", zstophy,zout)
1000c
1001         CALL histdef(nid_mth, "t2m", "Temperature 2m", "K",
1002     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1003     .                "ave(X)", zstorad,zout)
1004c
1005         CALL histdef(nid_mth, "t2m_min", "Temp. 2m min.",
1006     .                "K",
1007     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1008     .                t2mincels, zstophy,zout)
1009c
1010         CALL histdef(nid_mth, "t2m_max", "Temp. 2m max.",
1011     .                "K",
1012     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1013     .                t2maxcels, zstophy,zout)
1014c
1015c        CALL histdef(nid_mth, "tsoil", "Sfce soil Temperature",
1016c    .                "K",
1017c    .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1018c    .                "ave(X)", zstophy,zout)
1019c
1020         CALL histdef(nid_mth, "wind10m","10-m wind speed","m/s",
1021     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1022     .                "ave(X)", zstorad,zout)
1023c
1024         CALL histdef(nid_mth, "sicf", "Sea-ice fraction", "-",
1025     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1026     .                "ave(X)", zstophy,zout)
1027
1028c
1029         CALL histdef(nid_mth, "q2m","Specific humidity 2m","kg/kg",
1030     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1031     .                "ave(X)", zstorad,zout)
1032c
1033         CALL histdef(nid_mth, "u10m", "Vent zonal 10m", "m/s",
1034     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1035     .                "ave(X)", zstorad,zout)
1036c
1037         CALL histdef(nid_mth, "v10m", "Vent meridien 10m", "m/s",
1038     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1039     .                "ave(X)", zstorad,zout)
1040c
1041         CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa",
1042     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1043     .                "ave(X)", zstorad,zout)
1044c
1045         CALL histdef(nid_mth, "qsurf","Surface Air humidity","kg/kg",
1046     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1047     .                "ave(X)", zstophy,zout)
1048c
1049      if (.not. ok_veget) then
1050         CALL histdef(nid_mth, "qsol", "Soil watter content", "mm",
1051     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1052     .                "ave(X)", zstophy,zout)
1053      endif
1054c
1055         CALL histdef(nid_mth, "ndayrain", 
1056     .                "Number of day with rain (liq+sol)", "-",
1057     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1058     .                "inst(X)", zstomth,zout)
1059c
1060         CALL histdef(nid_mth, "precip", "Precipitation Totale liq+sol", 
1061     .                "kg/(s*m2)",
1062     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1063     .                "ave(X)", zstophy,zout)
1064c
1065         CALL histdef(nid_mth, "plul", "Large-scale Precip.", 
1066     .   "kg/(s*m2)",
1067     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1068     .                "ave(X)", zstophy,zout)
1069c
1070         CALL histdef(nid_mth, "pluc", "Convective Precip.", 
1071     .   "kg/(s*m2)",
1072     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1073     .                "ave(X)", zstophy,zout)
1074c
1075         CALL histdef(nid_mth, "snowf", "Snow fall", "kg/(s*m2)",
1076     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1077     .                "ave(X)", zstophy,zout)
1078c
1079c        CALL histdef(nid_mth, "snow_mass", "Snow Mass", "kg/m2",
1080c    .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1081c    .                "ave(X)", zstophy,zout)
1082c
1083         CALL histdef(nid_mth, "evap", 
1084     .                "Evaporation plus sublimation rate", "kg/(s*m2)",
1085     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1086     .                "ave(X)", zstophy,zout)
1087c
1088         CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2",
1089     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1090     .                "ave(X)", zstophy,zout)
1091c
1092         CALL histdef(nid_mth, "tops0", "CS Solar rad. at TOA", "W/m2",
1093     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1094     .                "ave(X)", zstophy,zout)
1095c
1096         CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2",
1097     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1098     .                "ave(X)", zstophy,zout)
1099c
1100         CALL histdef(nid_mth, "topl0", "IR rad. at TOA", "W/m2",
1101     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1102     .                "ave(X)", zstophy,zout)
1103c
1104         CALL histdef(nid_mth, "SWupTOA", "SWup at TOA","W/m2",
1105     .                iim,jjphy_nb,nhori, 1,1,1,-99,
1106     .                32, "ave(X)", zstophy,zout)
1107c
1108         CALL histdef(nid_mth, "SWupTOAclr", 
1109     .                "SWup clear sky at TOA","W/m2",
1110     .                iim,jjphy_nb,nhori, 1,1,1,-99,
1111     .                32, "ave(X)", zstophy,zout)
1112c
1113         CALL histdef(nid_mth, "SWdnTOA", "SWdn at TOA","W/m2",
1114     .                iim,jjphy_nb,nhori, 1,1,1,-99,
1115     .                32, "ave(X)", zstophy,zout)
1116c
1117         CALL histdef(nid_mth, "SWdnTOAclr", 
1118     .                "SWdn clear sky at TOA","W/m2",
1119     .                iim,jjphy_nb,nhori, 1,1,1,-99,
1120     .                32, "ave(X)", zstophy,zout)
1121c
1122         CALL histdef(nid_mth, "SWup200", "SWup at 200mb","W/m2",
1123     .                iim,jjphy_nb,nhori, 1,1,1,-99,
1124     .                32, "ave(X)", zstophy,zout)
1125c
1126         CALL histdef(nid_mth, "SWup200clr", 
1127     .                "SWup clear sky at 200mb","W/m2",
1128     .                iim,jjphy_nb,nhori, 1,1,1,-99,
1129     .                32, "ave(X)", zstophy,zout)
1130c
1131         CALL histdef(nid_mth, "SWdn200", "SWdn at 200mb","W/m2",
1132     .                iim,jjphy_nb,nhori, 1,1,1,-99,
1133     .                32, "ave(X)", zstophy,zout)
1134c
1135         CALL histdef(nid_mth, "SWdn200clr", 
1136     .                "SWdn clear sky at 200mb","W/m2",
1137     .                iim,jjphy_nb,nhori, 1,1,1,-99,
1138     .                32, "ave(X)", zstophy,zout)
1139c
1140         CALL histdef(nid_mth, "LWup200", "LWup at 200mb","W/m2",
1141     .                iim,jjphy_nb,nhori, 1,1,1,-99,
1142     .                32, "ave(X)", zstophy,zout)
1143c
1144         CALL histdef(nid_mth, "LWup200clr", 
1145     .                "LWup clear sky at 200mb","W/m2",
1146     .                iim,jjphy_nb,nhori, 1,1,1,-99,
1147     .                32, "ave(X)", zstophy,zout)
1148c
1149         CALL histdef(nid_mth, "LWdn200", "LWdn at 200mb","W/m2",
1150     .                iim,jjphy_nb,nhori, 1,1,1,-99,
1151     .                32, "ave(X)", zstophy,zout)
1152c
1153         CALL histdef(nid_mth, "LWdn200clr", 
1154     .                "LWdn clear sky at 200mb","W/m2",
1155     .                iim,jjphy_nb,nhori, 1,1,1,-99,
1156     .                32, "ave(X)", zstophy,zout)
1157c
1158         CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2",
1159     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1160     .                "ave(X)", zstophy,zout)
1161c
1162         CALL histdef(nid_mth, "sols0", "Solar rad. at surf.", "W/m2",
1163     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1164     .                "ave(X)", zstophy,zout)
1165c
1166         CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2",
1167     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1168     .                "ave(X)", zstophy,zout)
1169c
1170         CALL histdef(nid_mth, "radsol", "Rayonnement au sol", "W/m2",
1171     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1172     .                "ave(X)", zstophy,zout)
1173c
1174         CALL histdef(nid_mth, "soll0", "IR rad. at surface", "W/m2",
1175     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1176     .                "ave(X)", zstophy,zout)
1177c
1178         CALL histdef(nid_mth, "SWupSFC", "SWup at surface","W/m2",
1179     .                iim,jjphy_nb,nhori, 1,1,1,-99,
1180     .                32, "ave(X)", zstophy,zout)
1181c
1182         CALL histdef(nid_mth, "SWupSFCclr", 
1183     .                "SWup clear sky at surface","W/m2",
1184     .                iim,jjphy_nb,nhori, 1,1,1,-99,
1185     .                32, "ave(X)", zstophy,zout)
1186c
1187         CALL histdef(nid_mth, "SWdnSFC", "SWdn at surface","W/m2",
1188     .                iim,jjphy_nb,nhori, 1,1,1,-99,
1189     .                32, "ave(X)", zstophy,zout)
1190c
1191         CALL histdef(nid_mth, "SWdnSFCclr", 
1192     .                "SWdn clear sky at surface","W/m2",
1193     .                iim,jjphy_nb,nhori, 1,1,1,-99,
1194     .                32, "ave(X)", zstophy,zout)
1195c
1196         CALL histdef(nid_mth,"LWupSFC","Upwd. IR rad. at surface", 
1197     .                "W/m2", iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1198     .                "ave(X)", zstophy,zout)
1199c
1200         CALL histdef(nid_mth,"LWdnSFC","Down. IR rad. at surface", 
1201     .                "W/m2", iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1202     .                "ave(X)", zstophy,zout)
1203c
1204         CALL histdef(nid_mth,"LWupSFCclr",
1205     .                "CS Upwd. IR rad. at surface", 
1206     .                "W/m2", iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1207     .                "ave(X)", zstophy,zout)
1208c
1209         CALL histdef(nid_mth,"LWdnSFCclr",
1210     .                "Down. CS IR rad. at surface", 
1211     .                "W/m2", iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1212     .                "ave(X)", zstophy,zout)
1213c
1214         CALL histdef(nid_mth, "bils", "Surf. total heat flux", "W/m2",
1215     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1216     .                "ave(X)", zstophy,zout)
1217c
1218         CALL histdef(nid_mth, "sens", "Sensible heat flux", "W/m2",
1219     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1220     .                "ave(X)", zstophy,zout)
1221c
1222         CALL histdef(nid_mth, "fder", "Heat flux derivation", "W/m2",
1223     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1224     .                "ave(X)", zstophy,zout)
1225c
1226c          CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa",
1227c    .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1228c    .                "ave(X)", zstophy,zout)
1229c
1230c        CALL histdef(nid_mth, "frtv", "Meridional wind stress", "Pa",
1231c    .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1232c    .                "ave(X)", zstophy,zout)
1233c
1234         CALL histdef(nid_mth, "ffonte","Thermal flux for snow melting",
1235     .                "W/m2",iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1236     .                "ave(X)", zstophy,zout)
1237c
1238         CALL histdef(nid_mth, "fqcalving","Ice Calving",
1239     .                "kg/m2/s",iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1240     .                "ave(X)", zstophy,zout)
1241c
1242         DO nsrf = 1, nbsrf
1243           call histdef(nid_mth, "taux_"//clnsurf(nsrf),
1244     $         "Zonal wind stress"//clnsurf(nsrf), "Pa", 
1245     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1246     $         "ave(X)", zstophy,zout)
1247
1248           call histdef(nid_mth, "tauy_"//clnsurf(nsrf),
1249     $         "Meridional wind stress "//clnsurf(nsrf), "Pa", 
1250     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1251     $         "ave(X)", zstophy,zout)
1252         ENDDO
1253c
1254           call histdef(nid_mth, "pourc_"//clnsurf(is_oce),
1255     $         "% "//clnsurf(is_oce), "%", 
1256     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1257     $         "ave(X)", zstophy,zout)
1258c
1259           call histdef(nid_mth, "fract_"//clnsurf(is_oce),
1260     $         "Fraction "//clnsurf(is_oce), "1", 
1261     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1262     $         "ave(X)", zstophy,zout)
1263C
1264         DO nsrf = 1, nbsrf
1265C
1266           call histdef(nid_mth, "tsol_"//clnsurf(nsrf),
1267     $         "Temperature "//clnsurf(nsrf), "K", 
1268     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1269     $         "ave(X)", zstophy,zout)
1270C
1271           call histdef(nid_mth, "sens_"//clnsurf(nsrf),
1272     $         "Sensible heat flux "//clnsurf(nsrf), "W/m2", 
1273     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1274     $         "ave(X)", zstophy,zout)
1275c
1276           call histdef(nid_mth, "lat_"//clnsurf(nsrf),
1277     $         "Latent heat flux "//clnsurf(nsrf), "W/m2", 
1278     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1279     $         "ave(X)", zstophy,zout)
1280C
1281           call histdef(nid_mth, "flw_"//clnsurf(nsrf),
1282     $         "LW "//clnsurf(nsrf), "W/m2",
1283     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1284     $         "ave(X)", zstophy,zout)
1285c
1286           call histdef(nid_mth, "fsw_"//clnsurf(nsrf),
1287     $         "SW "//clnsurf(nsrf), "W/m2",
1288     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1289     $         "ave(X)", zstophy,zout)
1290C
1291           call histdef(nid_mth, "wbils_"//clnsurf(nsrf),
1292     $         "Bilan sol "//clnsurf(nsrf), "W/m2",
1293     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1294     $         "ave(X)", zstophy,zout)
1295C
1296           call histdef(nid_mth, "wbilo_"//clnsurf(nsrf),
1297     $         "Bilan eau "//clnsurf(nsrf), "kg/(m2*s)",
1298     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1299     $         "ave(X)", zstophy,zout)
1300C
1301cIM 250105 BEG
1302          IF(1.EQ.0) THEN
1303C
1304           call histdef(nid_mth, "evapour_"//clnsurf(nsrf),
1305     $         "Evap*pourc "//clnsurf(nsrf), "kg/(s*m2)",
1306     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1307     $         "ave(X)", zstophy,zout)
1308C
1309           call histdef(nid_mth, "prepour_"//clnsurf(nsrf),
1310     $         "Precip*pourc "//clnsurf(nsrf), "kg/(s*m2)",
1311     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1312     $         "ave(X)", zstophy,zout)
1313C
1314          ENDIF !(1.EQ.0) THEN
1315cIM 250105 END
1316C
1317         END DO
1318c
1319         CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
1320     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1321     .                "ave(X)", zstophy,zout)
1322c
1323         CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
1324     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1325     .                "ave(X)", zstophy,zout)
1326c
1327         CALL histdef(nid_mth, "cldl", "Low-level cloudiness", "-",
1328     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1329     .                "ave(X)", zstophy,zout)
1330c
1331         CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", "-",
1332     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1333     .                "ave(X)", zstophy,zout)
1334c
1335         CALL histdef(nid_mth, "cldh", "High-level cloudiness", "-",
1336     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1337     .                "ave(X)", zstophy,zout)
1338c
1339         CALL histdef(nid_mth, "cldt", "Total cloudiness", "%",
1340     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1341     .                "ave(X)", zstophy,zout)
1342c
1343         CALL histdef(nid_mth,"cldq","Cloud liquid water path","kg/m2",
1344     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1345     .                "ave(X)", zstophy,zout)
1346c
1347         CALL histdef(nid_mth,"lwp","Cloud water path","kg/m2",
1348     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1349     .                "ave(X)", zstophy,zout)
1350c
1351         CALL histdef(nid_mth,"iwp","Cloud ice water path","kg/m2",
1352     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1353     .                "ave(X)", zstophy,zout)
1354c
1355         CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
1356     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1357     .                "ave(X)", zstophy,zout)
1358c
1359         CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
1360     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1361     .                "ave(X)", zstophy,zout)
1362c
1363         CALL histdef(nid_mth, "uq", "Zonal humidity transport", "-",
1364     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1365     .                "ave(X)", zstophy,zout)
1366c
1367         CALL histdef(nid_mth, "vq", "Merid humidity transport", "-",
1368     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1369     .                "ave(X)", zstophy,zout)
1370cKE43
1371      IF(iflag_con.GE.3) THEN ! sb
1372c
1373         CALL histdef(nid_mth, "cape", "Conv avlbl pot ener", "J/kg",
1374     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1375     .                "ave(X)", zstophy,zout)
1376c
1377         CALL histdef(nid_mth, "pbase", "Cld base pressure", "mb",
1378     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1379     .                "ave(X)", zstophy,zout)
1380c
1381         CALL histdef(nid_mth, "ptop", "Cld top pressure", "mb",
1382     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1383     .                "ave(X)", zstophy,zout)
1384c
1385         CALL histdef(nid_mth, "fbase", "Cld base mass flux", "kg/m2/s",
1386     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1387     .                "ave(X)", zstophy,zout)
1388c
1389         CALL histdef(nid_mth, "prw", "Precipitable water", "kg/m2",
1390     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1391     .                "ave(X)", zstophy,zout)
1392c
1393         CALL histdef(nid_mth, "upwd", "saturated updraft", "kg/m2/s",
1394     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1395     .                "ave(X)", zstophy,zout)
1396c
1397      ENDIF !iflag_con .GE. 3
1398c34EK
1399cIM cf. AM 081204 BEG
1400c34EK
1401c HBTM2
1402         CALL histdef(nid_mth, "s_pblh", "Boundary Layer Height", "m",
1403     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1404     .                "ave(X)", zstophy,zout)
1405         CALL histdef(nid_mth, "s_pblt", "T at Boundary Layer Height",
1406     .                "K",
1407     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1408     .                "ave(X)", zstophy,zout)
1409         CALL histdef(nid_mth, "s_lcl", "Condensation level", "m",
1410     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1411     .                "ave(X)", zstophy,zout)
1412         CALL histdef(nid_mth,"s_capCL","Conv avlbl pot ener for ABL",
1413     .                "J/m2",
1414     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1415     .                "ave(X)", zstophy,zout)
1416         CALL histdef(nid_mth, "s_oliqCL", "Liq Water in BL", "kg/m2",
1417     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1418     .                "ave(X)", zstophy,zout)
1419         CALL histdef(nid_mth,"s_cteiCL","Instability criteria (ABL)",
1420     .                "K",
1421     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1422     .                "ave(X)", zstophy,zout)
1423         CALL histdef(nid_mth, "s_therm", "Exces du thermique", "K",
1424     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1425     .                "ave(X)", zstophy,zout)
1426         CALL histdef(nid_mth, "s_trmb1", "deep_cape(HBTM2)", "J/m2",
1427     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1428     .                "ave(X)", zstophy,zout)
1429         CALL histdef(nid_mth, "s_trmb2", "inhibition (HBTM2)", "J/m2",
1430     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1431     .                "ave(X)", zstophy,zout)
1432         CALL histdef(nid_mth, "s_trmb3", "Point Omega (HBTM2)", "m",
1433     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1434     .                "ave(X)", zstophy,zout)
1435c
1436cIM cf. AM 081204 END
1437c
1438c Champs interpolles sur des niveaux de pression
1439         DO k=1, nlevSTD
1440c
1441c        bb=clevSTD(k)
1442c
1443c        IF(k.GE.2) THEN
1444         IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
1445         IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
1446c         aa=clevSTD(k)
1447c         bb=aa(1:lnblnk1(aa))
1448c        ENDIF
1449c
1450       IF(bb2.EQ."850".OR.bb2.EQ."700"
1451     $   .OR.bb2.EQ."500".OR.bb2.EQ."200") THEN
1452c
1453          CALL histdef(nid_mth, "u"//bb2,
1454     .                 "Zonal wind "//bb2//"mb","m/s",
1455     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1456     .                "inst(X)", zout,zout)
1457c
1458          CALL histdef(nid_mth, "v"//bb2,
1459     .                 "Meridional wind "//bb2//"mb","m/s",
1460     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1461     .                "inst(X)", zout,zout)
1462c
1463          CALL histdef(nid_mth, "w"//bb2,
1464     .                 "Vertical wind "//bb2//"mb","m/s",
1465     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1466     .                "inst(X)", zout,zout)
1467c
1468          CALL histdef(nid_mth, "phi"//bb2,
1469     .                 "Geopotential "//bb2//"mb","m",
1470     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1471     .                "inst(X)", zout,zout)
1472c
1473         ENDIF !(bb2.EQ."850".OR.bb2.EQ."700"
1474c
1475       ENDDO !nlevSTD
1476c
1477cIM diagnostiques flux ocean-atm ou ocean-glace de mer
1478cIM pour utilisation dans un modele de "slab" ocean
1479c
1480       CALL histdef(nid_mth, "fluxo",
1481     .              "Flux turbulents ocean-atmosphere", "W/m2",
1482     .              iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1483     .              "ave(X)", zstophy,zout)
1484c
1485       CALL histdef(nid_mth, "fluxg",
1486     .              "Flux turbulents ocean-glace de mer","W/m2",
1487     .              iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1488     .              "ave(X)", zstophy,zout)
1489c
1490       CALL histdef(nid_mth, "t_oce_sic",
1491     .              "Temp. mixte oce-sic","K",
1492     .              iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1493     .              "ave(X)", zstophy,zout)
1494c
1495       IF (OCEAN.EQ.'force ') THEN
1496c
1497        CALL histdef(nid_mth, "lmt_bils",
1498     .       "Bilan au sol atmosphere forcee", "W/m2",
1499     .       iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1500     .       "ave(X)", zstophy,zout)
1501c
1502       ELSE IF (OCEAN.EQ.'slab  ') THEN
1503c
1504        CALL histdef(nid_mth, "slab_bils",
1505     .              "Bilan au sol Slab", "W/m2",
1506     .              iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1507     .              "ave(X)", zstophy,zout)
1508c
1509        CALL histdef(nid_mth, "tslab", "Slab SST ", "K",
1510     .              iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1511     .              "ave(X)", zstophy,zout)
1512c
1513        CALL histdef(nid_mth, "seaice", "Slab seaice", "kg/m2",
1514     .              iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1515     .              "ave(X)", zstophy,zout)
1516c
1517        CALL histdef(nid_mth, "siceh", "Slab seaice height", "m",
1518     .              iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1519     .              "ave(X)", zstophy,zout)
1520c
1521       ENDIF !OCEAN.EQ.slab/force
1522c
1523      ENDIF !lev_histmth.GE.1
1524c
1525      IF(lev_histmth.GE.2) THEN
1526c
1527c Champs 3D:
1528c
1529c
1530         CALL histdef(nid_mth,"lwcon","Cloud liquid water content",
1531     .                "kg/kg",
1532     .                iim,jjphy_nb,nhori, klev,1,klev, nvert, 32, 
1533     .                "ave(X)", zstophy,zout)
1534c
1535         CALL histdef(nid_mth,"iwcon","Cloud ice water content","kg/kg",
1536     .                iim,jjphy_nb,nhori, klev,1,klev, nvert, 32, 
1537     .                "ave(X)", zstophy,zout)
1538c
1539         CALL histdef(nid_mth, "temp", "Air temperature", "K",
1540     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1541     .                "ave(X)", zstophy,zout)
1542c
1543         CALL histdef(nid_mth, "ovap", "Specific humidity", "kg/kg",
1544     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1545     .                "ave(X)", zstophy,zout)
1546c
1547c        CALL histdef(nid_mth,"wvap","Water vapor mixing ratio","kg/kg",
1548c    .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1549c    .                "ave(X)", zstophy,zout)
1550c
1551         CALL histdef(nid_mth, "geop", "Geopotential height", "m2/s2",
1552     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1553     .                "ave(X)", zstophy,zout)
1554c
1555         CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s",
1556     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1557     .                "ave(X)", zstophy,zout)
1558c
1559         CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s",
1560     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1561     .                "ave(X)", zstophy,zout)
1562c
1563         CALL histdef(nid_mth, "vitw", "Vertical wind", "m/s",
1564     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1565     .                "ave(X)", zstophy,zout)
1566c
1567         CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
1568     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1569     .                "ave(X)", zstophy,zout)
1570c
1571        CALL histdef(nid_mth, "rneb", "Cloud fraction", "%",
1572     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1573     .                "ave(X)", zstophy,zout)
1574c
1575         CALL histdef(nid_mth, "rnebcon", "Convective Cloud Fraction"
1576     .                , "-",
1577     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1578     .                "ave(X)", zstophy,zout)
1579c
1580         CALL histdef(nid_mth, "rhum", "Relative humidity", "-",
1581     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1582     .                "ave(X)", zstophy,zout)
1583c
1584         CALL histdef(nid_mth, "ozone", "Ozone concentration", "ppmv",
1585     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1586     .                "ave(X)", zstophy,zout)
1587c
1588         CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
1589     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1590     .                "ave(X)", zstophy,zout)
1591c
1592         CALL histdef(nid_mth, "dqphy","Physics dQ","(kg/kg)/s",
1593     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1594     .                "ave(X)", zstophy,zout)
1595c
1596      ENDIF !lev_histmth.GE.2
1597c
1598      IF(lev_histmth.GE.3) THEN
1599c
1600        DO nsrf=1, nbsrf
1601c
1602           call histdef(nid_mth, "albe_"//clnsurf(nsrf),
1603     $         "Albedo surf. "//clnsurf(nsrf), "-", 
1604     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1605     $         "ave(X)", zstophy,zout)
1606c
1607           call histdef(nid_mth, "rugs_"//clnsurf(nsrf),
1608     $         "Rugosite "//clnsurf(nsrf), "m", 
1609     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
1610     $         "ave(X)", zstophy,zout)
1611c
1612         CALL histdef(nid_mth, "ages_"//clnsurf(nsrf), "Snow age","day",
1613     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1614     .                "ave(X)", zstophy,zout)
1615c
1616         ENDDO !nsrf=1, nbsrf
1617c
1618         CALL histdef(nid_mth, "albs", "Surface albedo", "-",
1619     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1620     .                "ave(X)", zstophy,zout)
1621         CALL histdef(nid_mth, "albslw", "Surface albedo LW", "-",
1622     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32, 
1623     .                "ave(X)", zstophy,zout)
1624c
1625      ENDIF !lev_histmth.GE.3
1626c
1627      IF(lev_histmth.GE.4) THEN
1628c
1629         CALL histdef(nid_mth, "clwcon", 
1630     .                "Convective Cloud Liquid water content"
1631     .                , "kg/kg",
1632     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1633     .                "ave(X)", zstophy,zout)
1634c
1635         CALL histdef(nid_mth,"Ma","undilute adiab updraft","kg/m2/s",
1636     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1637     .                "ave(X)", zstophy,zout)
1638c
1639         CALL histdef(nid_mth, "dnwd", "saturated downdraft","kg/m2/s",
1640     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1641     .                "ave(X)", zstophy,zout)
1642c
1643         CALL histdef(nid_mth, "dnwd0", "unsat. downdraft", "kg/m2/s",
1644     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1645     .                "ave(X)", zstophy,zout)
1646c
1647         CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
1648     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1649     .                "ave(X)", zstophy,zout)
1650c
1651         CALL histdef(nid_mth, "dqdyn","Dynamics dQ","(kg/kg)/s",
1652     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1653     .                "ave(X)", zstophy,zout)
1654c
1655         CALL histdef(nid_mth, "dqcon","Convection dQ","(kg/kg)/s",
1656     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1657     .                "ave(X)", zstophy,zout)
1658c
1659         CALL histdef(nid_mth, "dtlschr",
1660     $       "Large-scale condensational heating rate", "K/s",
1661     $       iim,jjphy_nb
1662     $       ,nhori, klev,1,klev,nvert, 32,"ave(X)", zstophy,zout)
1663c
1664         CALL histdef(nid_mth, "dqlsc","Condensation dQ","(kg/kg)/s",
1665     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1666     .                "ave(X)", zstophy,zout)
1667c
1668         CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s",
1669     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1670     .                "ave(X)", zstophy,zout)
1671c
1672         CALL histdef(nid_mth, "dqvdf","Boundary-layer dQ","(kg/kg)/s",
1673     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1674     .                "ave(X)", zstophy,zout)
1675c
1676         CALL histdef(nid_mth, "dteva", "Reevaporation dT", "K/s",
1677     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1678     .                "ave(X)", zstophy,zout)
1679c
1680         CALL histdef(nid_mth,"dqeva","Reevaporation dQ","(kg/kg)/s",
1681     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1682     .                "ave(X)", zstophy,zout)
1683
1684         CALL histdef(nid_mth, "ptconv", "POINTS CONVECTIFS"," ",
1685     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1686     .                "ave(X)", zstophy,zout)
1687
1688         CALL histdef(nid_mth, "ratqs", "RATQS"," ",
1689     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1690     .                "ave(X)", zstophy,zout)
1691
1692c
1693         CALL histdef(nid_mth, "dqajs","Dry adjust. dQ","(kg/kg)/s",
1694     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1695     .                "ave(X)", zstophy,zout)
1696c
1697         CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s",
1698     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1699     .                "ave(X)", zstophy,zout)
1700c
1701         CALL histdef(nid_mth, "dtsw0", "CS SW radiation dT", "K/s",
1702     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1703     .                "ave(X)", zstophy,zout)
1704c
1705         CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
1706     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1707     .                "ave(X)", zstophy,zout)
1708c
1709         CALL histdef(nid_mth,"dtlw0","CS LW radiation dT","K/s",
1710     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1711     .                "ave(X)", zstophy,zout)
1712c
1713         CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
1714     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1715     .                "ave(X)", zstophy,zout)
1716c
1717         CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
1718     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1719     .                "ave(X)", zstophy,zout)
1720c
1721         CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
1722     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1723     .                "ave(X)", zstophy,zout)
1724c
1725         IF (ok_orodr) THEN
1726         CALL histdef(nid_mth, "duoro", "Orography dU", "m/s2",
1727     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1728     .                "ave(X)", zstophy,zout)
1729c
1730         CALL histdef(nid_mth, "dvoro", "Orography dV", "m/s2",
1731     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1732     .                "ave(X)", zstophy,zout)
1733c
1734         ENDIF
1735C
1736         IF (ok_orolf) THEN
1737         CALL histdef(nid_mth, "dulif", "Orography dU", "m/s2",
1738     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1739     .                "ave(X)", zstophy,zout)
1740c
1741         CALL histdef(nid_mth, "dvlif", "Orography dV", "m/s2",
1742     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1743     .                "ave(X)", zstophy,zout)
1744         ENDIF
1745C
1746         if (nqmax.GE.3) THEN
1747         DO iq=3,nqmax
1748         iiq=niadv(iq)
1749         CALL histdef(nid_mth, tnom(iq), ttext(iiq), "-",
1750     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
1751     .                "ave(X)", zstophy,zout)
1752         ENDDO
1753         ENDIF
1754
1755      ENDIF !lev_histmth.GE.4
1756c
1757      ENDIF !type_run.EQ."CLIM".OR.type_run.EQ."ENSP")
1758c
1759         CALL histend(nid_mth)
1760c
1761         ndex2d = 0
1762         ndex3d = 0
1763c
1764c$OMP END MASTER
1765      ENDIF ! fin de test sur ok_mensuel
Note: See TracBrowser for help on using the repository browser.