source: LMDZ4/trunk/libf/phylmd/ini_histmth.h @ 644

Last change on this file since 644 was 644, checked in by Laurent Fairhead, 20 years ago

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