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

Last change on this file since 860 was 860, checked in by lmdzadmin, 17 years ago

Unite de vitw etait fausse dans histday.nc et histmth.nc
IM
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 68.4 KB
Line 
1c
2c $Header$
3c
4      IF (ok_mensuel) THEN
5c$OMP MASTER
6c
7       zstophy = dtime
8       zstomth = ecrit_mth                       !valeurs mensuelles "instantannees"
9       zstorad = dtime * FLOAT(radpas)    !variables stockees toutes les 2h
10       zstohf = ecrit_hf                         !variables stockees toutes les 6h
11       zout = ecrit_mth                          !frequence des sorties : 1 val. par mois
12c
13         idayref = day_ref
14         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
15c
16cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
17cym         DO i = 1, iim
18cym            zx_lon(i,1) = rlon(i+1)
19cym            zx_lon(i,jjmp1) = rlon(i+1)
20cym         ENDDO
21         DO ll=1,klev
22            znivsig(ll)=float(ll)
23         ENDDO
24cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
25cym         CALL histbeg("histmth.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
26cym     .                 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 
27cym     .                 nhori, nid_mth)
28
29           CALL histbeg_phy("histmth",itau_phy, zjulian, dtime, 
30     .                 nhori, nid_mth)
31     
32       write(*,*)'Mensuel ', itau_phy, zjulian
33         CALL histvert(nid_mth, "presnivs", "Vertical levels", "mb",
34     .                 klev, presnivs/100., nvert)
35c
36      IF(type_run.EQ."CLIM".OR.type_run.EQ."ENSP") THEN
37c
38      IF(lev_histmth.GE.1) THEN
39c
40         CALL histdef(nid_mth, "phis", 
41     .                "Surface geop. height", "m2/s2",
42     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
43     .                "once",  zstophy,zout)
44c
45         CALL histdef(nid_mth, "aire", "Grid area", "-",
46     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
47     .                "once",  zstophy,zout)
48c
49         CALL histdef(nid_mth, "pourc_"//clnsurf(is_ter),
50     .         "% "//clnsurf(is_ter), "%", 
51     .         iim,jj_nb,nhori, 1,1,1, -99, 32,
52     .         "once", zstophy,zout)
53c
54         CALL histdef(nid_mth, "fract_"//clnsurf(is_ter),
55     .         "Fraction "//clnsurf(is_ter), "1", 
56     .         iim,jj_nb,nhori, 1,1,1, -99, 32,
57     .         "once", zstophy,zout)
58c
59c Champs 2D:
60c
61         CALL histdef(nid_mth, "slp", "Sea Level Pressure", "Pa",
62     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
63     .                "ave(X)", zstophy,zout)
64c
65         CALL histdef(nid_mth, "tsol", "Surface Temperature", "K",
66     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
67     .                "ave(X)", zstophy,zout)
68c
69         CALL histdef(nid_mth, "t2m", "Temperature 2m", "K",
70     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
71     .                "ave(X)", zstophy,zout)
72c
73c ENSEMBLES BEG
74         CALL histdef(nid_mth, "t2m_min", "Temp. 2m min.",
75     .                "K",
76     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
77     .                t2mincels, zstophy,zout)
78c
79         CALL histdef(nid_mth, "t2m_max", "Temp. 2m max.",
80     .                "K",
81     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
82     .                t2maxcels, zstophy,zout)
83c
84c        CALL histdef(nid_mth, "tsoil", "Sfce soil Temperature",
85c    .                "K",
86c    .                iim,jj_nb,nhori, 1,1,1, -99, 32,
87c    .                "ave(X)", zstophy,zout)
88c
89         CALL histdef(nid_mth, "wind10m","10-m wind speed","m/s",
90     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
91     .                "ave(X)", zstophy,zout)
92c
93         CALL histdef(nid_mth, "sicf", "Sea-ice fraction", "-",
94     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
95     .                "ave(X)", zstophy,zout)
96c
97c ENSEMBLES END
98         CALL histdef(nid_mth, "q2m", "Specific humidity 2m", "kg/kg",
99     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
100     .                "ave(X)", zstophy,zout)
101c
102         CALL histdef(nid_mth, "u10m", "Vent zonal 10m", "m/s",
103     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
104     .                "ave(X)", zstophy,zout)
105c
106         CALL histdef(nid_mth, "v10m", "Vent meridien 10m", "m/s",
107     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
108     .                "ave(X)", zstophy,zout)
109c
110         CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa",
111     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
112     .                "ave(X)", zstophy,zout)
113c
114         CALL histdef(nid_mth, "qsurf", "Surface Air humidity", "kg/kg",
115     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
116     .                "ave(X)", zstophy,zout)
117c
118      if (.not. ok_veget) then
119         CALL histdef(nid_mth, "qsol", "Soil watter content", "mm",
120     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
121     .                "ave(X)", zstophy,zout)
122      endif
123c
124         CALL histdef(nid_mth, "ndayrain", 
125     .                "Number of day with rain (liq+sol)", "-",
126     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
127     .                "inst(X)", zstomth,zout)
128c
129         CALL histdef(nid_mth, "precip", "Precipitation Totale liq+sol", 
130     .                "kg/(s*m2)",
131     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
132     .                "ave(X)", zstophy,zout)
133c
134         CALL histdef(nid_mth, "plul", "Large-scale Precip.", 
135     .   "kg/(s*m2)",
136     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
137     .                "ave(X)", zstophy,zout)
138c
139         CALL histdef(nid_mth, "pluc", "Convective Precip.", 
140     .   "kg/(s*m2)",
141     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
142     .                "ave(X)", zstophy,zout)
143c
144         CALL histdef(nid_mth, "snow", "Snow fall", "kg/(s*m2)",
145     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
146     .                "ave(X)", zstophy,zout)
147c
148c        CALL histdef(nid_mth, "snow_mass", "Snow Mass", "kg/m2",
149c    .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
150c    .                "ave(X)", zstophy,zout)
151c
152         CALL histdef(nid_mth, "evap", "Evaporation", "kg/(s*m2)",
153     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
154     .                "ave(X)", zstophy,zout)
155c
156         CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2",
157     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
158     .                "ave(X)", zstophy,zout)
159c
160         CALL histdef(nid_mth, "tops0", "CS Solar rad. at TOA", "W/m2",
161     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
162     .                "ave(X)", zstophy,zout)
163c
164         CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2",
165     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
166     .                "ave(X)", zstophy,zout)
167c
168         CALL histdef(nid_mth, "topl0", "IR rad. at TOA", "W/m2",
169     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
170     .                "ave(X)", zstophy,zout)
171c
172         CALL histdef(nid_mth, "SWupTOA", "SWup at TOA","W/m2",
173     .                iim,jj_nb,nhori, 1,1,1,-99,
174     .                32, "ave(X)", zstophy,zout)
175c
176         CALL histdef(nid_mth, "SWupTOAclr", 
177     .                "SWup clear sky at TOA","W/m2",
178     .                iim,jj_nb,nhori, 1,1,1,-99,
179     .                32, "ave(X)", zstophy,zout)
180c
181         CALL histdef(nid_mth, "SWdnTOA", "SWdn at TOA","W/m2",
182     .                iim,jj_nb,nhori, 1,1,1,-99,
183     .                32, "ave(X)", zstophy,zout)
184c
185         CALL histdef(nid_mth, "SWdnTOAclr", 
186     .                "SWdn clear sky at TOA","W/m2",
187     .                iim,jj_nb,nhori, 1,1,1,-99,
188     .                32, "ave(X)", zstophy,zout)
189c
190         CALL histdef(nid_mth, "SWup200", "SWup at 200mb","W/m2",
191     .                iim,jj_nb,nhori, 1,1,1,-99,
192     .                32, "ave(X)", zstophy,zout)
193c
194         CALL histdef(nid_mth, "SWup200clr", 
195     .                "SWup clear sky at 200mb","W/m2",
196     .                iim,jj_nb,nhori, 1,1,1,-99,
197     .                32, "ave(X)", zstophy,zout)
198c
199         CALL histdef(nid_mth, "SWdn200", "SWdn at 200mb","W/m2",
200     .                iim,jj_nb,nhori, 1,1,1,-99,
201     .                32, "ave(X)", zstophy,zout)
202c
203         CALL histdef(nid_mth, "SWdn200clr", 
204     .                "SWdn clear sky at 200mb","W/m2",
205     .                iim,jj_nb,nhori, 1,1,1,-99,
206     .                32, "ave(X)", zstophy,zout)
207c
208         CALL histdef(nid_mth, "LWup200", "LWup at 200mb","W/m2",
209     .                iim,jj_nb,nhori, 1,1,1,-99,
210     .                32, "ave(X)", zstophy,zout)
211c
212         CALL histdef(nid_mth, "LWup200clr", 
213     .                "LWup clear sky at 200mb","W/m2",
214     .                iim,jj_nb,nhori, 1,1,1,-99,
215     .                32, "ave(X)", zstophy,zout)
216c
217         CALL histdef(nid_mth, "LWdn200", "LWdn at 200mb","W/m2",
218     .                iim,jj_nb,nhori, 1,1,1,-99,
219     .                32, "ave(X)", zstophy,zout)
220c
221         CALL histdef(nid_mth, "LWdn200clr", 
222     .                "LWdn clear sky at 200mb","W/m2",
223     .                iim,jj_nb,nhori, 1,1,1,-99,
224     .                32, "ave(X)", zstophy,zout)
225c
226         CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2",
227     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
228     .                "ave(X)", zstophy,zout)
229c
230         CALL histdef(nid_mth, "sols0", "Solar rad. at surf.", "W/m2",
231     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
232     .                "ave(X)", zstophy,zout)
233c
234         CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2",
235     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
236     .                "ave(X)", zstophy,zout)
237c
238         CALL histdef(nid_mth, "radsol", "Rayonnement au sol", "W/m2",
239     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
240     .                "ave(X)", zstophy,zout)
241c
242         CALL histdef(nid_mth, "soll0", "IR rad. at surface", "W/m2",
243     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
244     .                "ave(X)", zstophy,zout)
245c
246         CALL histdef(nid_mth, "SWupSFC", "SWup at surface","W/m2",
247     .                iim,jj_nb,nhori, 1,1,1,-99,
248     .                32, "ave(X)", zstophy,zout)
249c
250         CALL histdef(nid_mth, "SWupSFCclr", 
251     .                "SWup clear sky at surface","W/m2",
252     .                iim,jj_nb,nhori, 1,1,1,-99,
253     .                32, "ave(X)", zstophy,zout)
254c
255         CALL histdef(nid_mth, "SWdnSFC", "SWdn at surface","W/m2",
256     .                iim,jj_nb,nhori, 1,1,1,-99,
257     .                32, "ave(X)", zstophy,zout)
258c
259         CALL histdef(nid_mth, "SWdnSFCclr", 
260     .                "SWdn clear sky at surface","W/m2",
261     .                iim,jj_nb,nhori, 1,1,1,-99,
262     .                32, "ave(X)", zstophy,zout)
263c
264         CALL histdef(nid_mth,"LWupSFC","Upwd. IR rad. at surface", 
265     .                "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32, 
266     .                "ave(X)", zstophy,zout)
267c
268         CALL histdef(nid_mth,"LWdnSFC","Down. IR rad. at surface", 
269     .                "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32, 
270     .                "ave(X)", zstophy,zout)
271c
272         CALL histdef(nid_mth,"LWupSFCclr",
273     .                "CS Upwd. IR rad. at surface", 
274     .                "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32, 
275     .                "ave(X)", zstophy,zout)
276c
277         CALL histdef(nid_mth,"LWdnSFCclr",
278     .                "Down. CS IR rad. at surface", 
279     .                "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32, 
280     .                "ave(X)", zstophy,zout)
281c
282         CALL histdef(nid_mth, "bils", "Surf. total heat flux", "W/m2",
283     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
284     .                "ave(X)", zstophy,zout)
285c
286         CALL histdef(nid_mth, "sens", "Sensible heat flux", "W/m2",
287     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
288     .                "ave(X)", zstophy,zout)
289c
290         CALL histdef(nid_mth, "fder", "Heat flux derivation", "W/m2",
291     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
292     .                "ave(X)", zstophy,zout)
293c
294c          CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa",
295c    .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
296c    .                "ave(X)", zstophy,zout)
297c
298c        CALL histdef(nid_mth, "frtv", "Meridional wind stress", "Pa",
299c    .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
300c    .                "ave(X)", zstophy,zout)
301c
302         CALL histdef(nid_mth, "ffonte","Thermal flux for snow melting",
303     .                "W/m2",iim,jj_nb,nhori, 1,1,1, -99, 32,
304     .                "ave(X)", zstophy,zout)
305c
306         CALL histdef(nid_mth, "fqcalving","Ice Calving",
307     .                "kg/m2/s",iim,jj_nb,nhori, 1,1,1, -99, 32,
308     .                "ave(X)", zstophy,zout)
309c
310         CALL histdef(nid_mth, "fqfonte","Land ice melt",
311     .                "kg/m2/s",iim,jj_nb,nhori, 1,1,1, -99, 32,
312     .                "ave(X)", zstophy,zout)
313
314         DO nsrf = 1, nbsrf
315           call histdef(nid_mth, "taux_"//clnsurf(nsrf),
316     $         "Zonal wind stress"//clnsurf(nsrf), "Pa", 
317     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
318     $         "ave(X)", zstophy,zout)
319
320           call histdef(nid_mth, "tauy_"//clnsurf(nsrf),
321     $         "Meridional wind stress "//clnsurf(nsrf), "Pa", 
322     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
323     $         "ave(X)", zstophy,zout)
324         ENDDO
325c
326         DO nsrf = 1, nbsrf
327C
328          IF(nsrf.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", "Pa/s",
633     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
634     .                "ave(X)", zstophy,zout)
635c
636         CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
637     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
638     .                "ave(X)", zstophy,zout)
639c
640        CALL histdef(nid_mth, "rneb", "Cloud fraction", "-",
641     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
642     .                "ave(X)", zstophy,zout)
643c
644         CALL histdef(nid_mth, "rnebcon", "Convective Cloud Fraction"
645     .                , "-",
646     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
647     .                "ave(X)", zstophy,zout)
648c
649         CALL histdef(nid_mth, "rhum", "Relative humidity", "-",
650     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
651     .                "ave(X)", zstophy,zout)
652c
653         CALL histdef(nid_mth, "ozone", "Ozone concentration", "ppmv",
654     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
655     .                "ave(X)", zstophy,zout)
656c
657         CALL histdef(nid_mth, "upwd", "saturated updraft", "kg/m2/s",
658     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
659     .                "ave(X)", zstophy,zout)
660c
661         CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
662     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
663     .                "ave(X)", zstophy,zout)
664c
665         CALL histdef(nid_mth, "dqphy", "Physics dQ", "(kg/kg)/s",
666     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
667     .                "ave(X)", zstophy,zout)
668c
669      ENDIF !lev_histmth.GE.2
670c
671      IF(lev_histmth.GE.3) THEN
672c
673        DO nsrf=1, nbsrf
674c
675           call histdef(nid_mth, "albe_"//clnsurf(nsrf),
676     $         "Albedo surf. "//clnsurf(nsrf), "-", 
677     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
678     $         "ave(X)", zstophy,zout)
679c
680           call histdef(nid_mth, "rugs_"//clnsurf(nsrf),
681     $         "Latent heat flux "//clnsurf(nsrf), "W/m2", 
682     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
683     $         "ave(X)", zstophy,zout)
684c
685         CALL histdef(nid_mth, "ages_"//clnsurf(nsrf), "Snow age","day",
686     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
687     .                "ave(X)", zstophy,zout)
688c
689         ENDDO !nsrf=1, nbsrf
690c
691         CALL histdef(nid_mth, "albs", "Surface albedo", "-",
692     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
693     .                "ave(X)", zstophy,zout)
694         CALL histdef(nid_mth, "albslw", "Surface albedo LW", "-",
695     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
696     .                "ave(X)", zstophy,zout)
697c
698      ENDIF !lev_histmth.GE.3
699c
700      IF(lev_histmth.GE.4) THEN
701c
702         CALL histdef(nid_mth, "clwcon", 
703     .                "Convective Cloud Liquid water content"
704     .                , "kg/kg",
705     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
706     .                "ave(X)", zstophy,zout)
707c
708         CALL histdef(nid_mth,"Ma","undilute adiab updraft","kg/m2/s",
709     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
710     .                "ave(X)", zstophy,zout)
711c
712         CALL histdef(nid_mth, "dnwd", "saturated downdraft","kg/m2/s",
713     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
714     .                "ave(X)", zstophy,zout)
715c
716         CALL histdef(nid_mth, "dnwd0", "unsat. downdraft", "kg/m2/s",
717     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
718     .                "ave(X)", zstophy,zout)
719c
720         CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
721     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
722     .                "ave(X)", zstophy,zout)
723c
724         CALL histdef(nid_mth, "dqdyn", "Dynamics dQ", "(kg/kg)/s",
725     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
726     .                "ave(X)", zstophy,zout)
727c
728         CALL histdef(nid_mth, "dtcon", "Convection dT", "K/s",
729     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
730     .                "ave(X)", zstophy,zout)
731c
732c        CALL histdef(nid_mth, "ducon", "Convection du", "m/s2",
733c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
734c    .                "ave(X)", zstophy,zout)
735c
736         CALL histdef(nid_mth, "dqcon", "Convection dQ", "(kg/kg)/s",
737     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
738     .                "ave(X)", zstophy,zout)
739c
740         CALL histdef(nid_mth, "dtlsc", "Condensation dT", "K/s",
741     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
742     .                "ave(X)", zstophy,zout)
743c
744         CALL histdef(nid_mth, "dtlschr",
745     $     "Large-scale condensational heating rate", "K/s",iim,jj_nb
746     $     ,nhori, klev,1,klev,nvert, 32,"ave(X)", zstophy,zout)
747c
748         CALL histdef(nid_mth, "dqlsc", "Condensation dQ", "(kg/kg)/s",
749     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
750     .                "ave(X)", zstophy,zout)
751c
752         CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s",
753     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
754     .                "ave(X)", zstophy,zout)
755c
756         CALL histdef(nid_mth, "dqvdf","Boundary-layer dQ","(kg/kg)/s",
757     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
758     .                "ave(X)", zstophy,zout)
759c
760         CALL histdef(nid_mth, "dteva", "Reevaporation dT", "K/s",
761     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
762     .                "ave(X)", zstophy,zout)
763c
764         CALL histdef(nid_mth, "dqeva","Reevaporation dQ","(kg/kg)/s",
765     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
766     .                "ave(X)", zstophy,zout)
767
768         CALL histdef(nid_mth, "ptconv", "POINTS CONVECTIFS"," ",
769     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
770     .                "ave(X)", zstophy,zout)
771
772         CALL histdef(nid_mth, "ratqs", "RATQS"," ",
773     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
774     .                "ave(X)", zstophy,zout)
775
776c
777         CALL histdef(nid_mth, "dtajs", "Dry adjust. dT", "K/s",
778     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
779     .                "ave(X)", zstophy,zout)
780
781         CALL histdef(nid_mth,"dqajs","Dry adjust. dQ","(kg/kg)/s",
782     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
783     .                "ave(X)", zstophy,zout)
784c
785         CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s",
786     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
787     .                "ave(X)", zstophy,zout)
788c
789         CALL histdef(nid_mth, "dtsw0", "CS SW radiation dT", "K/s",
790     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
791     .                "ave(X)", zstophy,zout)
792c
793         CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
794     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
795     .                "ave(X)", zstophy,zout)
796c
797         CALL histdef(nid_mth,"dtlw0","CS LW radiation dT","K/s",
798     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
799     .                "ave(X)", zstophy,zout)
800c
801         CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
802     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
803     .                "ave(X)", zstophy,zout)
804c
805         CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
806     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
807     .                "ave(X)", zstophy,zout)
808c
809         CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
810     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
811     .                "ave(X)", zstophy,zout)
812c
813         IF (ok_orodr) THEN
814         CALL histdef(nid_mth, "duoro", "Orography dU", "m/s2",
815     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
816     .                "ave(X)", zstophy,zout)
817c
818         CALL histdef(nid_mth, "dvoro", "Orography dV", "m/s2",
819     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
820     .                "ave(X)", zstophy,zout)
821c
822         ENDIF
823C
824         IF (ok_orolf) THEN
825         CALL histdef(nid_mth, "dulif", "Orography dU", "m/s2",
826     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
827     .                "ave(X)", zstophy,zout)
828c
829         CALL histdef(nid_mth, "dvlif", "Orography dV", "m/s2",
830     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
831     .                "ave(X)", zstophy,zout)
832         ENDIF
833C
834         if (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", "Pa/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.