source: LMDZ4/trunk/libf/phytherm/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: 70.5 KB
Line 
1c
2c $Header$
3c
4      IF (ok_mensuel) THEN
5c$OMP MASTER
6c
7       zstophy = dtime
8       zstomth = ecrit_mth                       !valeurs mensuelles "instantannees"
9       zstorad = dtime * FLOAT(radpas)    !variables stockees toutes les 2h
10       zstohf = ecrit_hf                         !variables stockees toutes les 6h
11       zout = ecrit_mth                          !frequence des sorties : 1 val. par mois
12c
13         idayref = day_ref
14         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
15c
16cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
17cym         DO i = 1, iim
18cym            zx_lon(i,1) = rlon(i+1)
19cym            zx_lon(i,jjmp1) = rlon(i+1)
20cym         ENDDO
21         DO ll=1,klev
22            znivsig(ll)=float(ll)
23         ENDDO
24cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
25cym         CALL histbeg("histmth.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
26cym     .                 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 
27cym     .                 nhori, nid_mth)
28
29           CALL histbeg_phy("histmth",itau_phy, zjulian, dtime, 
30     .                 nhori, nid_mth)
31     
32       write(*,*)'Mensuel ', itau_phy, zjulian
33         CALL histvert(nid_mth, "presnivs", "Vertical levels", "mb",
34     .                 klev, presnivs/100., nvert)
35c
36      IF(type_run.EQ."CLIM".OR.type_run.EQ."ENSP") THEN
37c
38      IF(lev_histmth.GE.1) THEN
39c
40         CALL histdef(nid_mth, "phis", 
41     .                "Surface geop. height", "m2/s2",
42     .                iim,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!FH Sorties pour la couche limite
703      if (iflag_pbl>1) then
704         CALL histdef(nid_mth, "tke","TKE","m2/s2",
705     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
706     .                "ave(X)", zstophy,zout)
707c
708         CALL histdef(nid_mth, "tke_max","TKE max","m2/s2",
709     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
710     .                "t_max(X)", zstophy,zout)
711      endif
712c
713         CALL histdef(nid_mth, "kz","Kz melange","m2/s",
714     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
715     .                "ave(X)", zstophy,zout)
716c
717         CALL histdef(nid_mth, "kz_max","Kz melange max","m2/s",
718     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
719     .                "t_max(X)", zstophy,zout)
720
721
722         CALL histdef(nid_mth, "clwcon", 
723     .                "Convective Cloud Liquid water content"
724     .                , "kg/kg",
725     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
726     .                "ave(X)", zstophy,zout)
727c
728         CALL histdef(nid_mth,"Ma","undilute adiab updraft","kg/m2/s",
729     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
730     .                "ave(X)", zstophy,zout)
731c
732         CALL histdef(nid_mth, "dnwd", "saturated downdraft","kg/m2/s",
733     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
734     .                "ave(X)", zstophy,zout)
735c
736         CALL histdef(nid_mth, "dnwd0", "unsat. downdraft", "kg/m2/s",
737     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
738     .                "ave(X)", zstophy,zout)
739c
740         CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
741     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
742     .                "ave(X)", zstophy,zout)
743c
744         CALL histdef(nid_mth, "dqdyn", "Dynamics dQ", "(kg/kg)/s",
745     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
746     .                "ave(X)", zstophy,zout)
747c
748         CALL histdef(nid_mth, "dtcon", "Convection dT", "K/s",
749     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
750     .                "ave(X)", zstophy,zout)
751c
752
753
754c
755c        CALL histdef(nid_mth, "ducon", "Convection du", "m/s2",
756c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
757c    .                "ave(X)", zstophy,zout)
758c
759         CALL histdef(nid_mth, "dqcon", "Convection dQ", "(kg/kg)/s",
760     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
761     .                "ave(X)", zstophy,zout)
762c
763         CALL histdef(nid_mth, "dtlsc", "Condensation dT", "K/s",
764     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
765     .                "ave(X)", zstophy,zout)
766c
767         CALL histdef(nid_mth, "dtlschr",
768     $     "Large-scale condensational heating rate", "K/s",iim,jj_nb
769     $     ,nhori, klev,1,klev,nvert, 32,"ave(X)", zstophy,zout)
770c
771         CALL histdef(nid_mth, "dqlsc", "Condensation dQ", "(kg/kg)/s",
772     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
773     .                "ave(X)", zstophy,zout)
774c
775         CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s",
776     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
777     .                "ave(X)", zstophy,zout)
778c
779         CALL histdef(nid_mth, "dqvdf","Boundary-layer dQ","(kg/kg)/s",
780     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
781     .                "ave(X)", zstophy,zout)
782c
783         CALL histdef(nid_mth, "dteva", "Reevaporation dT", "K/s",
784     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
785     .                "ave(X)", zstophy,zout)
786c
787         CALL histdef(nid_mth, "dqeva","Reevaporation dQ","(kg/kg)/s",
788     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
789     .                "ave(X)", zstophy,zout)
790
791         CALL histdef(nid_mth, "ptconv", "POINTS CONVECTIFS"," ",
792     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
793     .                "ave(X)", zstophy,zout)
794
795         CALL histdef(nid_mth, "ratqs", "RATQS"," ",
796     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
797     .                "ave(X)", zstophy,zout)
798
799c
800         CALL histdef(nid_mth, "dtthe", "Dry adjust. dT", "K/s",
801     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
802     .                "ave(X)", zstophy,zout)
803
804         CALL histdef(nid_mth,"dqthe","Dry adjust. dQ","(kg/kg)/s",
805     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
806     .                "ave(X)", zstophy,zout)
807c
808         CALL histdef(nid_mth, "dtajs", "Dry adjust. dT", "K/s",
809     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
810     .                "ave(X)", zstophy,zout)
811
812         CALL histdef(nid_mth,"dqajs","Dry adjust. dQ","(kg/kg)/s",
813     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
814     .                "ave(X)", zstophy,zout)
815c
816         CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s",
817     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
818     .                "ave(X)", zstophy,zout)
819c
820         CALL histdef(nid_mth, "dtsw0", "CS SW radiation dT", "K/s",
821     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
822     .                "ave(X)", zstophy,zout)
823c
824         CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
825     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
826     .                "ave(X)", zstophy,zout)
827c
828         CALL histdef(nid_mth,"dtlw0","CS LW radiation dT","K/s",
829     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
830     .                "ave(X)", zstophy,zout)
831c
832         CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
833     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
834     .                "ave(X)", zstophy,zout)
835c
836         CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
837     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
838     .                "ave(X)", zstophy,zout)
839c
840         CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
841     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
842     .                "ave(X)", zstophy,zout)
843c
844         IF (ok_orodr) THEN
845         CALL histdef(nid_mth, "duoro", "Orography dU", "m/s2",
846     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
847     .                "ave(X)", zstophy,zout)
848c
849         CALL histdef(nid_mth, "dvoro", "Orography dV", "m/s2",
850     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
851     .                "ave(X)", zstophy,zout)
852c
853         ENDIF
854C
855         IF (ok_orolf) THEN
856         CALL histdef(nid_mth, "dulif", "Orography dU", "m/s2",
857     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
858     .                "ave(X)", zstophy,zout)
859c
860         CALL histdef(nid_mth, "dvlif", "Orography dV", "m/s2",
861     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
862     .                "ave(X)", zstophy,zout)
863         ENDIF
864C
865         if (nqmax.GE.3) THEN
866           DO iq=3,nqmax
867             iiq=niadv(iq)
868             CALL histdef(nid_mth, tnom(iq), ttext(iiq), "-",
869     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
870     .                "ave(X)", zstophy,zout)
871             ENDDO
872           ENDIF
873c
874      ENDIF !lev_histmth.GE.4
875c
876      ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
877c
878      IF(lev_histmth.GE.1) THEN
879c
880         CALL histdef(nid_mth, "phis", "Surface geop. height", "m",
881     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
882     .                "once",  zstophy,zout)
883c
884         CALL histdef(nid_mth, "aire", "Grid area", "-",
885     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
886     .                "once",  zstophy,zout)
887c
888         CALL histdef(nid_mth, "pourc_"//clnsurf(is_ter),
889     .         "% "//clnsurf(is_ter), "%", 
890     .         iim,jj_nb,nhori, 1,1,1, -99, 32,
891     .         "once", zstophy,zout)
892c
893         CALL histdef(nid_mth, "fract_"//clnsurf(is_ter),
894     .         "Fraction "//clnsurf(is_ter), "1", 
895     .         iim,jj_nb,nhori, 1,1,1, -99, 32,
896     .         "once", zstophy,zout)
897c
898         CALL histdef(nid_mth, "pourc_"//clnsurf(is_lic),
899     .         "% "//clnsurf(is_lic), "%", 
900     .         iim,jj_nb,nhori, 1,1,1, -99, 32,
901     .         "once", zstophy,zout)
902c
903         CALL histdef(nid_mth, "fract_"//clnsurf(is_lic),
904     .         "Fraction "//clnsurf(is_lic), "1", 
905     .         iim,jj_nb,nhori, 1,1,1, -99, 32,
906     .         "once", zstophy,zout)
907c
908         CALL histdef(nid_mth, "pourc_"//clnsurf(is_sic),
909     .         "% "//clnsurf(is_sic), "%",
910     .         iim,jj_nb,nhori, 1,1,1, -99, 32,
911     .                "ave(X)", zstophy,zout)
912c
913         CALL histdef(nid_mth, "fract_"//clnsurf(is_sic),
914     .         "Fraction "//clnsurf(is_sic), "1",
915     .         iim,jj_nb,nhori, 1,1,1, -99, 32,
916     .                "ave(X)", zstophy,zout)
917c
918         CALL histdef(nid_mth, "flat", "Latent heat flux", "W/m2",
919     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
920     .                "ave(X)", zstophy,zout)
921c
922         CALL histdef(nid_mth, "bilTOA",
923     .                "Net radiation at TOA", "W/m2",
924     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
925     .                "ave(X)", zstophy,zout)
926c
927c
928c        CALL histdef(nid_mth, "tops01",
929c    .                "TOA reflected clear-sky shortwave
930c    .                radiation (method 1)", "W/m2",
931c    .                iim,jj_nb,nhori, 1,1,1, -99, 32,
932c    .                "ave(X)", zstophy,zout)
933c
934          CALL histdef(nid_mth, "dtdia","dtdia","K/s",
935     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
936     .                "ave(X)", zstophy,zout)
937c
938          CALL histdef(nid_mth, "dtsw","dtsw","K/s",
939     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
940     .                "ave(X)", zstophy,zout)
941c
942          CALL histdef(nid_mth, "dtlw","dtlw","K/s",
943     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
944     .                "ave(X)", zstophy,zout)
945c
946          CALL histdef(nid_mth, "dtcon","dtcon","K/s",
947     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
948     .                "ave(X)", zstophy,zout)
949c
950         CALL histdef(nid_mth, "dtthe", "Dry adjust. dT", "K/s",
951     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
952     .                "ave(X)", zstophy,zout)
953
954         CALL histdef(nid_mth,"dqthe","Dry adjust. dQ","(kg/kg)/s",
955     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
956     .                "ave(X)", zstophy,zout)
957c
958          CALL histdef(nid_mth, "dtajs","dtajs","K/s",
959     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
960     .                "ave(X)", zstophy,zout)
961c
962          CALL histdef(nid_mth, "dtlsc","dtlsc","K/s",
963     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
964     .                "ave(X)", zstophy,zout)
965c
966          CALL histdef(nid_mth, "dqdia","dqdia","(kg/kg)/s",
967     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
968     .                "ave(X)", zstophy,zout)
969c
970          CALL histdef(nid_mth, "duoli",
971     .         "duoli","m/s2",
972     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
973     .                "ave(X)", zstophy,zout)
974c
975          CALL histdef(nid_mth, "dvoli",
976     .         "dvoli","m/s2",
977     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
978     .                "ave(X)", zstophy,zout)
979c
980          CALL histdef(nid_mth, "ducon",
981     .         "ducon","m/s2",
982     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
983     .                "ave(X)", zstophy,zout)
984c
985          CALL histdef(nid_mth, "dvcon",
986     .         "dvcon","m/s2",
987     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
988     .                "ave(X)", zstophy,zout)
989c
990          CALL histdef(nid_mth, "dutot",
991     .         "dutot","m/s2",
992     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
993     .                "ave(X)", zstophy,zout)
994c
995          CALL histdef(nid_mth, "dvtot",
996     .         "dvtot","m/s2",
997     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
998     .                "ave(X)", zstophy,zout)
999c
1000c         CALL histdef(nid_mth, "cldtrue",
1001c    .         "TRUE Cloud fraction ","%",
1002c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1003c    .         "inst(X)", zout,zout)
1004c
1005c         CALL histdef(nid_mth, "cldTOA",
1006c    .         "Cloud amount (satellite view)","%",
1007c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1008c    .         "ave(X)", zstophy,zout)
1009c
1010c         CALL histdef(nid_mth, "cldSFC",
1011c    .         "Cloud amount (surface view)","%",
1012c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1013c    .         "ave(X)", zstophy,zout)
1014c
1015          CALL histdef(nid_mth, "GWDtaux",
1016     .         "GWD induced Ewd wind stress","Pa",
1017     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1018     .         "ave(X)", zstophy,zout)
1019c
1020          CALL histdef(nid_mth, "GWDtauy",
1021     .         "GWD induced Nwd wind stress","Pa",
1022     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1023     .         "ave(X)", zstophy,zout)
1024c
1025c Champs 2D:
1026c
1027         CALL histdef(nid_mth, "slp", "Sea Level Pressure", "Pa",
1028     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1029     .                "ave(X)", zstorad,zout)
1030c
1031         CALL histdef(nid_mth, "tsol", "Surface Temperature", "K",
1032     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1033     .                "ave(X)", zstophy,zout)
1034c
1035         CALL histdef(nid_mth, "t2m", "Temperature 2m", "K",
1036     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1037     .                "ave(X)", zstorad,zout)
1038c
1039         CALL histdef(nid_mth, "t2m_min", "Temp. 2m min.",
1040     .                "K",
1041     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1042     .                t2mincels, zstophy,zout)
1043c
1044         CALL histdef(nid_mth, "t2m_max", "Temp. 2m max.",
1045     .                "K",
1046     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1047     .                t2maxcels, zstophy,zout)
1048c
1049c        CALL histdef(nid_mth, "tsoil", "Sfce soil Temperature",
1050c    .                "K",
1051c    .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1052c    .                "ave(X)", zstophy,zout)
1053c
1054         CALL histdef(nid_mth, "wind10m","10-m wind speed","m/s",
1055     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1056     .                "ave(X)", zstorad,zout)
1057c
1058         CALL histdef(nid_mth, "sicf", "Sea-ice fraction", "-",
1059     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1060     .                "ave(X)", zstophy,zout)
1061
1062c
1063         CALL histdef(nid_mth, "q2m","Specific humidity 2m","kg/kg",
1064     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1065     .                "ave(X)", zstorad,zout)
1066c
1067         CALL histdef(nid_mth, "u10m", "Vent zonal 10m", "m/s",
1068     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1069     .                "ave(X)", zstorad,zout)
1070c
1071         CALL histdef(nid_mth, "v10m", "Vent meridien 10m", "m/s",
1072     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1073     .                "ave(X)", zstorad,zout)
1074c
1075         CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa",
1076     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1077     .                "ave(X)", zstorad,zout)
1078c
1079         CALL histdef(nid_mth, "qsurf","Surface Air humidity","kg/kg",
1080     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1081     .                "ave(X)", zstophy,zout)
1082c
1083      if (.not. ok_veget) then
1084         CALL histdef(nid_mth, "qsol", "Soil watter content", "mm",
1085     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1086     .                "ave(X)", zstophy,zout)
1087      endif
1088c
1089         CALL histdef(nid_mth, "ndayrain", 
1090     .                "Number of day with rain (liq+sol)", "-",
1091     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1092     .                "inst(X)", zstomth,zout)
1093c
1094         CALL histdef(nid_mth, "precip", "Precipitation Totale liq+sol", 
1095     .                "kg/(s*m2)",
1096     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1097     .                "ave(X)", zstophy,zout)
1098c
1099         CALL histdef(nid_mth, "plul", "Large-scale Precip.", 
1100     .   "kg/(s*m2)",
1101     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1102     .                "ave(X)", zstophy,zout)
1103c
1104         CALL histdef(nid_mth, "pluc", "Convective Precip.", 
1105     .   "kg/(s*m2)",
1106     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1107     .                "ave(X)", zstophy,zout)
1108c
1109         CALL histdef(nid_mth, "snowf", "Snow fall", "kg/(s*m2)",
1110     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1111     .                "ave(X)", zstophy,zout)
1112c
1113c        CALL histdef(nid_mth, "snow_mass", "Snow Mass", "kg/m2",
1114c    .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1115c    .                "ave(X)", zstophy,zout)
1116c
1117         CALL histdef(nid_mth, "evap", 
1118     .                "Evaporation plus sublimation rate", "kg/(s*m2)",
1119     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1120     .                "ave(X)", zstophy,zout)
1121c
1122         CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2",
1123     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1124     .                "ave(X)", zstophy,zout)
1125c
1126         CALL histdef(nid_mth, "tops0", "CS Solar rad. at TOA", "W/m2",
1127     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1128     .                "ave(X)", zstophy,zout)
1129c
1130         CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2",
1131     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1132     .                "ave(X)", zstophy,zout)
1133c
1134         CALL histdef(nid_mth, "topl0", "IR rad. at TOA", "W/m2",
1135     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1136     .                "ave(X)", zstophy,zout)
1137c
1138         CALL histdef(nid_mth, "SWupTOA", "SWup at TOA","W/m2",
1139     .                iim,jj_nb,nhori, 1,1,1,-99,
1140     .                32, "ave(X)", zstophy,zout)
1141c
1142         CALL histdef(nid_mth, "SWupTOAclr", 
1143     .                "SWup clear sky at TOA","W/m2",
1144     .                iim,jj_nb,nhori, 1,1,1,-99,
1145     .                32, "ave(X)", zstophy,zout)
1146c
1147         CALL histdef(nid_mth, "SWdnTOA", "SWdn at TOA","W/m2",
1148     .                iim,jj_nb,nhori, 1,1,1,-99,
1149     .                32, "ave(X)", zstophy,zout)
1150c
1151         CALL histdef(nid_mth, "SWdnTOAclr", 
1152     .                "SWdn clear sky at TOA","W/m2",
1153     .                iim,jj_nb,nhori, 1,1,1,-99,
1154     .                32, "ave(X)", zstophy,zout)
1155c
1156         CALL histdef(nid_mth, "SWup200", "SWup at 200mb","W/m2",
1157     .                iim,jj_nb,nhori, 1,1,1,-99,
1158     .                32, "ave(X)", zstophy,zout)
1159c
1160         CALL histdef(nid_mth, "SWup200clr", 
1161     .                "SWup clear sky at 200mb","W/m2",
1162     .                iim,jj_nb,nhori, 1,1,1,-99,
1163     .                32, "ave(X)", zstophy,zout)
1164c
1165         CALL histdef(nid_mth, "SWdn200", "SWdn at 200mb","W/m2",
1166     .                iim,jj_nb,nhori, 1,1,1,-99,
1167     .                32, "ave(X)", zstophy,zout)
1168c
1169         CALL histdef(nid_mth, "SWdn200clr", 
1170     .                "SWdn clear sky at 200mb","W/m2",
1171     .                iim,jj_nb,nhori, 1,1,1,-99,
1172     .                32, "ave(X)", zstophy,zout)
1173c
1174         CALL histdef(nid_mth, "LWup200", "LWup at 200mb","W/m2",
1175     .                iim,jj_nb,nhori, 1,1,1,-99,
1176     .                32, "ave(X)", zstophy,zout)
1177c
1178         CALL histdef(nid_mth, "LWup200clr", 
1179     .                "LWup clear sky at 200mb","W/m2",
1180     .                iim,jj_nb,nhori, 1,1,1,-99,
1181     .                32, "ave(X)", zstophy,zout)
1182c
1183         CALL histdef(nid_mth, "LWdn200", "LWdn at 200mb","W/m2",
1184     .                iim,jj_nb,nhori, 1,1,1,-99,
1185     .                32, "ave(X)", zstophy,zout)
1186c
1187         CALL histdef(nid_mth, "LWdn200clr", 
1188     .                "LWdn clear sky at 200mb","W/m2",
1189     .                iim,jj_nb,nhori, 1,1,1,-99,
1190     .                32, "ave(X)", zstophy,zout)
1191c
1192         CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2",
1193     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1194     .                "ave(X)", zstophy,zout)
1195c
1196         CALL histdef(nid_mth, "sols0", "Solar rad. at surf.", "W/m2",
1197     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1198     .                "ave(X)", zstophy,zout)
1199c
1200         CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2",
1201     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1202     .                "ave(X)", zstophy,zout)
1203c
1204         CALL histdef(nid_mth, "radsol", "Rayonnement au sol", "W/m2",
1205     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1206     .                "ave(X)", zstophy,zout)
1207c
1208         CALL histdef(nid_mth, "soll0", "IR rad. at surface", "W/m2",
1209     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1210     .                "ave(X)", zstophy,zout)
1211c
1212         CALL histdef(nid_mth, "SWupSFC", "SWup at surface","W/m2",
1213     .                iim,jj_nb,nhori, 1,1,1,-99,
1214     .                32, "ave(X)", zstophy,zout)
1215c
1216         CALL histdef(nid_mth, "SWupSFCclr", 
1217     .                "SWup clear sky at surface","W/m2",
1218     .                iim,jj_nb,nhori, 1,1,1,-99,
1219     .                32, "ave(X)", zstophy,zout)
1220c
1221         CALL histdef(nid_mth, "SWdnSFC", "SWdn at surface","W/m2",
1222     .                iim,jj_nb,nhori, 1,1,1,-99,
1223     .                32, "ave(X)", zstophy,zout)
1224c
1225         CALL histdef(nid_mth, "SWdnSFCclr", 
1226     .                "SWdn clear sky at surface","W/m2",
1227     .                iim,jj_nb,nhori, 1,1,1,-99,
1228     .                32, "ave(X)", zstophy,zout)
1229c
1230         CALL histdef(nid_mth,"LWupSFC","Upwd. IR rad. at surface", 
1231     .                "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32, 
1232     .                "ave(X)", zstophy,zout)
1233c
1234         CALL histdef(nid_mth,"LWdnSFC","Down. IR rad. at surface", 
1235     .                "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32, 
1236     .                "ave(X)", zstophy,zout)
1237c
1238         CALL histdef(nid_mth,"LWupSFCclr",
1239     .                "CS Upwd. IR rad. at surface", 
1240     .                "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32, 
1241     .                "ave(X)", zstophy,zout)
1242c
1243         CALL histdef(nid_mth,"LWdnSFCclr",
1244     .                "Down. CS IR rad. at surface", 
1245     .                "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32, 
1246     .                "ave(X)", zstophy,zout)
1247c
1248         CALL histdef(nid_mth, "bils", "Surf. total heat flux", "W/m2",
1249     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1250     .                "ave(X)", zstophy,zout)
1251c
1252         CALL histdef(nid_mth, "sens", "Sensible heat flux", "W/m2",
1253     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1254     .                "ave(X)", zstophy,zout)
1255c
1256         CALL histdef(nid_mth, "fder", "Heat flux derivation", "W/m2",
1257     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1258     .                "ave(X)", zstophy,zout)
1259c
1260c          CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa",
1261c    .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1262c    .                "ave(X)", zstophy,zout)
1263c
1264c        CALL histdef(nid_mth, "frtv", "Meridional wind stress", "Pa",
1265c    .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1266c    .                "ave(X)", zstophy,zout)
1267c
1268         CALL histdef(nid_mth, "ffonte","Thermal flux for snow melting",
1269     .                "W/m2",iim,jj_nb,nhori, 1,1,1, -99, 32,
1270     .                "ave(X)", zstophy,zout)
1271c
1272         CALL histdef(nid_mth, "fqcalving","Ice Calving",
1273     .                "kg/m2/s",iim,jj_nb,nhori, 1,1,1, -99, 32,
1274     .                "ave(X)", zstophy,zout)
1275c
1276         DO nsrf = 1, nbsrf
1277           call histdef(nid_mth, "taux_"//clnsurf(nsrf),
1278     $         "Zonal wind stress"//clnsurf(nsrf), "Pa", 
1279     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
1280     $         "ave(X)", zstophy,zout)
1281
1282           call histdef(nid_mth, "tauy_"//clnsurf(nsrf),
1283     $         "Meridional wind stress "//clnsurf(nsrf), "Pa", 
1284     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
1285     $         "ave(X)", zstophy,zout)
1286         ENDDO
1287c
1288           call histdef(nid_mth, "pourc_"//clnsurf(is_oce),
1289     $         "% "//clnsurf(is_oce), "%", 
1290     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
1291     $         "ave(X)", zstophy,zout)
1292c
1293           call histdef(nid_mth, "fract_"//clnsurf(is_oce),
1294     $         "Fraction "//clnsurf(is_oce), "1", 
1295     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
1296     $         "ave(X)", zstophy,zout)
1297C
1298         DO nsrf = 1, nbsrf
1299C
1300           call histdef(nid_mth, "tsol_"//clnsurf(nsrf),
1301     $         "Temperature "//clnsurf(nsrf), "K", 
1302     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
1303     $         "ave(X)", zstophy,zout)
1304C
1305           call histdef(nid_mth, "sens_"//clnsurf(nsrf),
1306     $         "Sensible heat flux "//clnsurf(nsrf), "W/m2", 
1307     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
1308     $         "ave(X)", zstophy,zout)
1309c
1310           call histdef(nid_mth, "lat_"//clnsurf(nsrf),
1311     $         "Latent heat flux "//clnsurf(nsrf), "W/m2", 
1312     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
1313     $         "ave(X)", zstophy,zout)
1314C
1315           call histdef(nid_mth, "flw_"//clnsurf(nsrf),
1316     $         "LW "//clnsurf(nsrf), "W/m2",
1317     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
1318     $         "ave(X)", zstophy,zout)
1319c
1320           call histdef(nid_mth, "fsw_"//clnsurf(nsrf),
1321     $         "SW "//clnsurf(nsrf), "W/m2",
1322     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
1323     $         "ave(X)", zstophy,zout)
1324C
1325           call histdef(nid_mth, "wbils_"//clnsurf(nsrf),
1326     $         "Bilan sol "//clnsurf(nsrf), "W/m2",
1327     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
1328     $         "ave(X)", zstophy,zout)
1329C
1330           call histdef(nid_mth, "wbilo_"//clnsurf(nsrf),
1331     $         "Bilan eau "//clnsurf(nsrf), "kg/(m2*s)",
1332     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
1333     $         "ave(X)", zstophy,zout)
1334C
1335cIM 250105 BEG
1336          IF(1.EQ.0) THEN
1337C
1338           call histdef(nid_mth, "evapour_"//clnsurf(nsrf),
1339     $         "Evap*pourc "//clnsurf(nsrf), "kg/(s*m2)",
1340     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
1341     $         "ave(X)", zstophy,zout)
1342C
1343           call histdef(nid_mth, "prepour_"//clnsurf(nsrf),
1344     $         "Precip*pourc "//clnsurf(nsrf), "kg/(s*m2)",
1345     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
1346     $         "ave(X)", zstophy,zout)
1347C
1348          ENDIF !(1.EQ.0) THEN
1349cIM 250105 END
1350C
1351         END DO
1352c
1353         CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
1354     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1355     .                "ave(X)", zstophy,zout)
1356c
1357         CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
1358     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1359     .                "ave(X)", zstophy,zout)
1360c
1361         CALL histdef(nid_mth, "cldl", "Low-level cloudiness", "-",
1362     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1363     .                "ave(X)", zstophy,zout)
1364c
1365         CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", "-",
1366     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1367     .                "ave(X)", zstophy,zout)
1368c
1369         CALL histdef(nid_mth, "cldh", "High-level cloudiness", "-",
1370     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1371     .                "ave(X)", zstophy,zout)
1372c
1373         CALL histdef(nid_mth, "cldt", "Total cloudiness", "%",
1374     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1375     .                "ave(X)", zstophy,zout)
1376c
1377         CALL histdef(nid_mth,"cldq","Cloud liquid water path","kg/m2",
1378     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1379     .                "ave(X)", zstophy,zout)
1380c
1381         CALL histdef(nid_mth,"lwp","Cloud water path","kg/m2",
1382     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1383     .                "ave(X)", zstophy,zout)
1384c
1385         CALL histdef(nid_mth,"iwp","Cloud ice water path","kg/m2",
1386     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1387     .                "ave(X)", zstophy,zout)
1388c
1389         CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
1390     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1391     .                "ave(X)", zstophy,zout)
1392c
1393         CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
1394     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1395     .                "ave(X)", zstophy,zout)
1396c
1397         CALL histdef(nid_mth, "uq", "Zonal humidity transport", "-",
1398     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1399     .                "ave(X)", zstophy,zout)
1400c
1401         CALL histdef(nid_mth, "vq", "Merid humidity transport", "-",
1402     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1403     .                "ave(X)", zstophy,zout)
1404cKE43
1405      IF(iflag_con.GE.3) THEN ! sb
1406c
1407         CALL histdef(nid_mth, "cape", "Conv avlbl pot ener", "J/kg",
1408     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1409     .                "ave(X)", zstophy,zout)
1410c
1411         CALL histdef(nid_mth, "pbase", "Cld base pressure", "mb",
1412     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1413     .                "ave(X)", zstophy,zout)
1414c
1415         CALL histdef(nid_mth, "ptop", "Cld top pressure", "mb",
1416     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1417     .                "ave(X)", zstophy,zout)
1418c
1419         CALL histdef(nid_mth, "fbase", "Cld base mass flux", "kg/m2/s",
1420     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1421     .                "ave(X)", zstophy,zout)
1422c
1423         CALL histdef(nid_mth, "prw", "Precipitable water", "kg/m2",
1424     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1425     .                "ave(X)", zstophy,zout)
1426c
1427         CALL histdef(nid_mth, "upwd", "saturated updraft", "kg/m2/s",
1428     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1429     .                "ave(X)", zstophy,zout)
1430c
1431      ENDIF !iflag_con .GE. 3
1432c34EK
1433cIM cf. AM 081204 BEG
1434c34EK
1435c HBTM2
1436         CALL histdef(nid_mth, "s_pblh", "Boundary Layer Height", "m",
1437     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1438     .                "ave(X)", zstophy,zout)
1439         CALL histdef(nid_mth, "s_pblt", "T at Boundary Layer Height",
1440     .                "K",
1441     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1442     .                "ave(X)", zstophy,zout)
1443         CALL histdef(nid_mth, "s_lcl", "Condensation level", "m",
1444     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1445     .                "ave(X)", zstophy,zout)
1446         CALL histdef(nid_mth,"s_capCL","Conv avlbl pot ener for ABL",
1447     .                "J/m2",
1448     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1449     .                "ave(X)", zstophy,zout)
1450         CALL histdef(nid_mth, "s_oliqCL", "Liq Water in BL", "kg/m2",
1451     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1452     .                "ave(X)", zstophy,zout)
1453         CALL histdef(nid_mth,"s_cteiCL","Instability criteria (ABL)",
1454     .                "K",
1455     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1456     .                "ave(X)", zstophy,zout)
1457         CALL histdef(nid_mth, "s_therm", "Exces du thermique", "K",
1458     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1459     .                "ave(X)", zstophy,zout)
1460         CALL histdef(nid_mth, "s_trmb1", "deep_cape(HBTM2)", "J/m2",
1461     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1462     .                "ave(X)", zstophy,zout)
1463         CALL histdef(nid_mth, "s_trmb2", "inhibition (HBTM2)", "J/m2",
1464     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1465     .                "ave(X)", zstophy,zout)
1466         CALL histdef(nid_mth, "s_trmb3", "Point Omega (HBTM2)", "m",
1467     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1468     .                "ave(X)", zstophy,zout)
1469c
1470cIM cf. AM 081204 END
1471c
1472c Champs interpolles sur des niveaux de pression
1473         DO k=1, nlevSTD
1474c
1475c        bb=clevSTD(k)
1476c
1477c        IF(k.GE.2) THEN
1478         IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
1479         IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
1480c         aa=clevSTD(k)
1481c         bb=aa(1:lnblnk1(aa))
1482c        ENDIF
1483c
1484       IF(bb2.EQ."850".OR.bb2.EQ."700"
1485     $   .OR.bb2.EQ."500".OR.bb2.EQ."200") THEN
1486c
1487          CALL histdef(nid_mth, "u"//bb2,
1488     .                 "Zonal wind "//bb2//"mb","m/s",
1489     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1490     .                "inst(X)", zout,zout)
1491c
1492          CALL histdef(nid_mth, "v"//bb2,
1493     .                 "Meridional wind "//bb2//"mb","m/s",
1494     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1495     .                "inst(X)", zout,zout)
1496c
1497          CALL histdef(nid_mth, "w"//bb2,
1498     .                 "Vertical wind "//bb2//"mb","m/s",
1499     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1500     .                "inst(X)", zout,zout)
1501c
1502          CALL histdef(nid_mth, "phi"//bb2,
1503     .                 "Geopotential "//bb2//"mb","m",
1504     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
1505     .                "inst(X)", zout,zout)
1506c
1507         ENDIF !(bb2.EQ."850".OR.bb2.EQ."700"
1508c
1509       ENDDO !nlevSTD
1510c
1511cIM diagnostiques flux ocean-atm ou ocean-glace de mer
1512cIM pour utilisation dans un modele de "slab" ocean
1513c
1514       CALL histdef(nid_mth, "fluxo",
1515     .              "Flux turbulents ocean-atmosphere", "W/m2",
1516     .              iim,jj_nb,nhori, 1,1,1, -99, 32,
1517     .              "ave(X)", zstophy,zout)
1518c
1519       CALL histdef(nid_mth, "fluxg",
1520     .              "Flux turbulents ocean-glace de mer","W/m2",
1521     .              iim,jj_nb,nhori, 1,1,1, -99, 32,
1522     .              "ave(X)", zstophy,zout)
1523c
1524       CALL histdef(nid_mth, "t_oce_sic",
1525     .              "Temp. mixte oce-sic","K",
1526     .              iim,jj_nb,nhori, 1,1,1, -99, 32,
1527     .              "ave(X)", zstophy,zout)
1528c
1529       IF (OCEAN.EQ.'force ') THEN
1530c
1531        CALL histdef(nid_mth, "lmt_bils",
1532     .       "Bilan au sol atmosphere forcee", "W/m2",
1533     .       iim,jj_nb,nhori, 1,1,1, -99, 32,
1534     .       "ave(X)", zstophy,zout)
1535c
1536       ELSE IF (OCEAN.EQ.'slab  ') THEN
1537c
1538        CALL histdef(nid_mth, "slab_bils",
1539     .              "Bilan au sol Slab", "W/m2",
1540     .              iim,jj_nb,nhori, 1,1,1, -99, 32,
1541     .              "ave(X)", zstophy,zout)
1542c
1543        CALL histdef(nid_mth, "tslab", "Slab SST ", "K",
1544     .              iim,jj_nb,nhori, 1,1,1, -99, 32,
1545     .              "ave(X)", zstophy,zout)
1546c
1547        CALL histdef(nid_mth, "seaice", "Slab seaice", "kg/m2",
1548     .              iim,jj_nb,nhori, 1,1,1, -99, 32,
1549     .              "ave(X)", zstophy,zout)
1550c
1551        CALL histdef(nid_mth, "siceh", "Slab seaice height", "m",
1552     .              iim,jj_nb,nhori, 1,1,1, -99, 32,
1553     .              "ave(X)", zstophy,zout)
1554c
1555       ENDIF !OCEAN.EQ.slab/force
1556c
1557      ENDIF !lev_histmth.GE.1
1558c
1559      IF(lev_histmth.GE.2) THEN
1560c
1561c Champs 3D:
1562c
1563c
1564         CALL histdef(nid_mth,"lwcon","Cloud liquid water content",
1565     .                "kg/kg",
1566     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32, 
1567     .                "ave(X)", zstophy,zout)
1568c
1569         CALL histdef(nid_mth,"iwcon","Cloud ice water content","kg/kg",
1570     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32, 
1571     .                "ave(X)", zstophy,zout)
1572c
1573         CALL histdef(nid_mth, "temp", "Air temperature", "K",
1574     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1575     .                "ave(X)", zstophy,zout)
1576c
1577         CALL histdef(nid_mth, "ovap", "Specific humidity", "kg/kg",
1578     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1579     .                "ave(X)", zstophy,zout)
1580c
1581c        CALL histdef(nid_mth,"wvap","Water vapor mixing ratio","kg/kg",
1582c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1583c    .                "ave(X)", zstophy,zout)
1584c
1585         CALL histdef(nid_mth, "geop", "Geopotential height", "m2/s2",
1586     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1587     .                "ave(X)", zstophy,zout)
1588c
1589         CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s",
1590     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1591     .                "ave(X)", zstophy,zout)
1592c
1593         CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s",
1594     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1595     .                "ave(X)", zstophy,zout)
1596c
1597         CALL histdef(nid_mth, "vitw", "Vertical wind", "Pa/s",
1598     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1599     .                "ave(X)", zstophy,zout)
1600c
1601         CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
1602     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1603     .                "ave(X)", zstophy,zout)
1604c
1605        CALL histdef(nid_mth, "rneb", "Cloud fraction", "%",
1606     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1607     .                "ave(X)", zstophy,zout)
1608c
1609         CALL histdef(nid_mth, "rnebcon", "Convective Cloud Fraction"
1610     .                , "-",
1611     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1612     .                "ave(X)", zstophy,zout)
1613c
1614         CALL histdef(nid_mth, "rhum", "Relative humidity", "-",
1615     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1616     .                "ave(X)", zstophy,zout)
1617c
1618         CALL histdef(nid_mth, "ozone", "Ozone concentration", "ppmv",
1619     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1620     .                "ave(X)", zstophy,zout)
1621c
1622         CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
1623     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1624     .                "ave(X)", zstophy,zout)
1625c
1626         CALL histdef(nid_mth, "dqphy","Physics dQ","(kg/kg)/s",
1627     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1628     .                "ave(X)", zstophy,zout)
1629c
1630      ENDIF !lev_histmth.GE.2
1631c
1632      IF(lev_histmth.GE.3) THEN
1633c
1634        DO nsrf=1, nbsrf
1635c
1636           call histdef(nid_mth, "albe_"//clnsurf(nsrf),
1637     $         "Albedo surf. "//clnsurf(nsrf), "-", 
1638     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
1639     $         "ave(X)", zstophy,zout)
1640c
1641           call histdef(nid_mth, "rugs_"//clnsurf(nsrf),
1642     $         "Rugosite "//clnsurf(nsrf), "m", 
1643     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
1644     $         "ave(X)", zstophy,zout)
1645c
1646         CALL histdef(nid_mth, "ages_"//clnsurf(nsrf), "Snow age","day",
1647     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1648     .                "ave(X)", zstophy,zout)
1649c
1650         ENDDO !nsrf=1, nbsrf
1651c
1652         CALL histdef(nid_mth, "albs", "Surface albedo", "-",
1653     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1654     .                "ave(X)", zstophy,zout)
1655         CALL histdef(nid_mth, "albslw", "Surface albedo LW", "-",
1656     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
1657     .                "ave(X)", zstophy,zout)
1658c
1659      ENDIF !lev_histmth.GE.3
1660c
1661      IF(lev_histmth.GE.4) THEN
1662c
1663!FH Sorties pour la couche limite
1664      if (iflag_pbl>1) then
1665         CALL histdef(nid_mth, "tke","TKE","m2/s2",
1666     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1667     .                "ave(X)", zstophy,zout)
1668c
1669         CALL histdef(nid_mth, "tke_max","TKE max","m2/s2",
1670     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1671     .                "t_max(X)", zstophy,zout)
1672      endif
1673c
1674         CALL histdef(nid_mth, "kz","Kz melange","m2/s",
1675     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1676     .                "ave(X)", zstophy,zout)
1677c
1678         CALL histdef(nid_mth, "kz_max","Kz melange max","m2/s",
1679     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1680     .                "t_max(X)", zstophy,zout)
1681
1682
1683
1684
1685         CALL histdef(nid_mth, "clwcon", 
1686     .                "Convective Cloud Liquid water content"
1687     .                , "kg/kg",
1688     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1689     .                "ave(X)", zstophy,zout)
1690c
1691         CALL histdef(nid_mth,"Ma","undilute adiab updraft","kg/m2/s",
1692     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1693     .                "ave(X)", zstophy,zout)
1694c
1695         CALL histdef(nid_mth, "dnwd", "saturated downdraft","kg/m2/s",
1696     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1697     .                "ave(X)", zstophy,zout)
1698c
1699         CALL histdef(nid_mth, "dnwd0", "unsat. downdraft", "kg/m2/s",
1700     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1701     .                "ave(X)", zstophy,zout)
1702c
1703         CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
1704     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1705     .                "ave(X)", zstophy,zout)
1706c
1707         CALL histdef(nid_mth, "dqdyn","Dynamics dQ","(kg/kg)/s",
1708     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1709     .                "ave(X)", zstophy,zout)
1710c
1711         CALL histdef(nid_mth, "dqcon","Convection dQ","(kg/kg)/s",
1712     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1713     .                "ave(X)", zstophy,zout)
1714c
1715         CALL histdef(nid_mth, "dtlschr",
1716     $       "Large-scale condensational heating rate", "K/s",
1717     $       iim,jj_nb
1718     $       ,nhori, klev,1,klev,nvert, 32,"ave(X)", zstophy,zout)
1719c
1720         CALL histdef(nid_mth, "dqlsc","Condensation dQ","(kg/kg)/s",
1721     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1722     .                "ave(X)", zstophy,zout)
1723c
1724         CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s",
1725     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1726     .                "ave(X)", zstophy,zout)
1727c
1728         CALL histdef(nid_mth, "dqvdf","Boundary-layer dQ","(kg/kg)/s",
1729     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1730     .                "ave(X)", zstophy,zout)
1731c
1732         CALL histdef(nid_mth, "dteva", "Reevaporation dT", "K/s",
1733     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1734     .                "ave(X)", zstophy,zout)
1735c
1736         CALL histdef(nid_mth,"dqeva","Reevaporation dQ","(kg/kg)/s",
1737     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1738     .                "ave(X)", zstophy,zout)
1739
1740         CALL histdef(nid_mth, "ptconv", "POINTS CONVECTIFS"," ",
1741     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1742     .                "ave(X)", zstophy,zout)
1743
1744         CALL histdef(nid_mth, "ratqs", "RATQS"," ",
1745     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1746     .                "ave(X)", zstophy,zout)
1747
1748c
1749         CALL histdef(nid_mth, "dqajs","Dry adjust. dQ","(kg/kg)/s",
1750     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1751     .                "ave(X)", zstophy,zout)
1752c
1753         CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s",
1754     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1755     .                "ave(X)", zstophy,zout)
1756c
1757         CALL histdef(nid_mth, "dtsw0", "CS SW radiation dT", "K/s",
1758     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1759     .                "ave(X)", zstophy,zout)
1760c
1761         CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
1762     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1763     .                "ave(X)", zstophy,zout)
1764c
1765         CALL histdef(nid_mth,"dtlw0","CS LW radiation dT","K/s",
1766     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1767     .                "ave(X)", zstophy,zout)
1768c
1769         CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
1770     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1771     .                "ave(X)", zstophy,zout)
1772c
1773         CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
1774     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1775     .                "ave(X)", zstophy,zout)
1776c
1777         CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
1778     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1779     .                "ave(X)", zstophy,zout)
1780c
1781         IF (ok_orodr) THEN
1782         CALL histdef(nid_mth, "duoro", "Orography dU", "m/s2",
1783     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1784     .                "ave(X)", zstophy,zout)
1785c
1786         CALL histdef(nid_mth, "dvoro", "Orography dV", "m/s2",
1787     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1788     .                "ave(X)", zstophy,zout)
1789c
1790         ENDIF
1791C
1792         IF (ok_orolf) THEN
1793         CALL histdef(nid_mth, "dulif", "Orography dU", "m/s2",
1794     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1795     .                "ave(X)", zstophy,zout)
1796c
1797         CALL histdef(nid_mth, "dvlif", "Orography dV", "m/s2",
1798     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1799     .                "ave(X)", zstophy,zout)
1800         ENDIF
1801C
1802         if (nqmax.GE.3) THEN
1803         DO iq=3,nqmax
1804         iiq=niadv(iq)
1805         CALL histdef(nid_mth, tnom(iq), ttext(iiq), "-",
1806     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
1807     .                "ave(X)", zstophy,zout)
1808         ENDDO
1809         ENDIF
1810
1811      ENDIF !lev_histmth.GE.4
1812c
1813      ENDIF !type_run.EQ."CLIM".OR.type_run.EQ."ENSP")
1814c
1815         CALL histend(nid_mth)
1816c
1817         ndex2d = 0
1818         ndex3d = 0
1819c
1820c$OMP END MASTER
1821      ENDIF ! fin de test sur ok_mensuel
Note: See TracBrowser for help on using the repository browser.