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

Last change on this file since 720 was 704, checked in by Laurent Fairhead, 19 years ago

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