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

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

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