source: LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histmth.h @ 506

Last change on this file since 506 was 506, checked in by lmdzadmin, 20 years ago

IM: ajout lev_hist(hf,day,mth) niveaux output definis dans physiq.def et

variables ENSEMBLES + rearrangement
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 27.9 KB
Line 
1      IF (ok_mensuel) THEN
2c
3        IF(region.EQ.'3d') THEN
4         zsto = dtime
5         zout = dtime * ecrit_mth
6c
7c zsto1: pour des valeurs "instantannees" mensuelles
8         zsto1 = dtime * ecrit_mth
9c zsto2: pour des flux radiatifs calcules tous les 2 heures
10         zsto2 = dtime * radpas
11         PRINT*,' zsto,zsto1,zsto2,zout=',zsto, zsto1, zsto2,zout
12        ELSEIF(region.EQ.'1d') THEN
13          zsto = dtime
14c         zout = dtime
15          zout = dtime *4.
16        ENDIF
17c
18         idayref = day_ref
19         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
20c
21         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
22         DO i = 1, iim
23            zx_lon(i,1) = rlon(i+1)
24            zx_lon(i,jjmp1) = rlon(i+1)
25         ENDDO
26         DO ll=1,klev
27            znivsig(ll)=float(ll)
28         ENDDO
29         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
30         CALL histbeg("histmth.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
31c071103
32c    .                 1,iim,1,jjmp1, 
33     .                 imin_debut,nbpti,
34     .                 jmin_debut,nbptj,
35c071103
36     .                 itau_phy, zjulian, dtime, 
37     .                 nhori, nid_mth)
38         write(*,*)'Mensuel ', itau_phy, zjulian
39         CALL histvert(nid_mth, "presnivs", "Vertical levels", "mb",
40     .                 klev, presnivs, nvert)
41c        call histvert(nid_mth, 'sig_s', 'Niveaux sigma','-',
42c    .              klev, znivsig, nvert)
43c
44c
45      IF(lev_histmth.GE.1) THEN
46         CALL histdef(nid_mth, "phis", "Surface geop. height", "-",
47     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
48     .                "once",  zsto,zout)
49c
50         CALL histdef(nid_mth, "aire", "Grid area", "-",
51     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
52     .                "once",  zsto,zout)
53c
54           call histdef(nid_mth, "pourc_"//clnsurf(is_ter),
55     $         "Fraction "//clnsurf(is_ter), "W/m2", 
56     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
57     $         "once", zsto,zout)
58c
59c Champs 2D:
60c
61         CALL histdef(nid_mth, "slp", "Sea Level Pressure", "Pa",
62     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
63     .                "ave(X)", zsto,zout)
64c
65         CALL histdef(nid_mth, "tsol", "Surface Temperature", "K",
66     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
67     .                "ave(X)", zsto,zout)
68c
69         CALL histdef(nid_mth, "t2m", "Temperature 2m", "K",
70     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
71     .                "ave(X)", zsto,zout)
72c ENSEMBLES BEG
73         CALL histdef(nid_mth, "t2m_min", "Temp. 2m min.",
74     .                "K",
75     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
76     .                t2mincels, zsto,zout)
77c
78         CALL histdef(nid_mth, "t2m_max", "Temp. 2m max.",
79     .                "K",
80     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
81     .                t2maxcels, zsto,zout)
82c
83c        CALL histdef(nid_mth, "tsoil", "Sfce soil Temperature",
84c    .                "K",
85c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
86c    .                "ave(X)", zsto,zout)
87c
88         CALL histdef(nid_mth, "wind10m","10-m wind speed","m/s",
89     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
90     .                "ave(X)", zsto,zout)
91c
92         CALL histdef(nid_mth, "sicf", "Sea-ice fraction", "-",
93     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
94     .                "ave(X)", zsto,zout)
95
96c
97c ENSEMBLES END
98         CALL histdef(nid_mth, "q2m", "Specific humidity 2m", "Kg/Kg",
99     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
100     .                "ave(X)", zsto,zout)
101c
102         CALL histdef(nid_mth, "u10m", "Vent zonal 10m", "m/s",
103     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
104     .                "ave(X)", zsto,zout)
105c
106         CALL histdef(nid_mth, "v10m", "Vent meridien 10m", "m/s",
107     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
108     .                "ave(X)", zsto,zout)
109c
110         CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa",
111     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
112     .                "ave(X)", zsto,zout)
113c
114         CALL histdef(nid_mth, "qsurf", "Surface Air humidity", "Kg/Kg",
115     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
116     .                "ave(X)", zsto,zout)
117c
118      if (.not. ok_veget) then
119         CALL histdef(nid_mth, "qsol", "Soil watter content", "mm",
120     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
121     .                "ave(X)", zsto,zout)
122      endif
123c
124         CALL histdef(nid_mth, "ndayrain", 
125     .                "Number of day with rain (liq+sol)", "-",
126     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
127     .                "inst(X)", zsto1,zout)
128c
129         CALL histdef(nid_mth, "precip", "Precipitation Totale liq+sol", 
130     .                "Kg/(s*m2)",
131     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
132     .                "ave(X)", zsto,zout)
133c
134         CALL histdef(nid_mth, "plul", "Large-scale Precip.", 
135     .   "Kg/(s*m2)",
136     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
137     .                "ave(X)", zsto,zout)
138c
139         CALL histdef(nid_mth, "pluc", "Convective Precip.", 
140     .   "Kg/(s*m2)",
141     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
142     .                "ave(X)", zsto,zout)
143c
144         CALL histdef(nid_mth, "snow", "Snow fall", "Kg/(s*m2)",
145     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
146     .                "ave(X)", zsto,zout)
147c
148         CALL histdef(nid_mth, "snow_mass", "Snow Mass", "Kg/m2",
149     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
150     .                "ave(X)", zsto,zout)
151c
152cIM: 071003 evap(Kg/m2/s) ==> evap(mm/day)
153        CALL histdef(nid_mth, "evap", "Evaporation", "kg/(s*m2)",
154c       CALL histdef(nid_mth, "evap", "Evaporation", "mm/day", 
155     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
156     .                "ave(X)", zsto,zout)
157c
158         CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2",
159     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
160     .                "ave(X)", zsto,zout)
161c
162         CALL histdef(nid_mth, "tops0", "Solar rad. at TOA", "W/m2",
163     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
164     .                "ave(X)", zsto,zout)
165c
166         CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2",
167     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
168     .                "ave(X)", zsto,zout)
169c
170         CALL histdef(nid_mth, "topl0", "IR rad. at TOA", "W/m2",
171     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
172     .                "ave(X)", zsto,zout)
173c
174         CALL histdef(nid_mth, "SWupTOA", "SWup at TOA","W/m2",
175     .                iim,jjmp1,nhori, 1,1,1,-99,
176     .                32, "ave(X)", zsto2,zout)
177c
178         CALL histdef(nid_mth, "SWupTOAclr", 
179     .                "SWup clear sky at TOA","W/m2",
180     .                iim,jjmp1,nhori, 1,1,1,-99,
181     .                32, "ave(X)", zsto2,zout)
182c
183         CALL histdef(nid_mth, "SWdnTOA", "SWdn at TOA","W/m2",
184     .                iim,jjmp1,nhori, 1,1,1,-99,
185     .                32, "ave(X)", zsto2,zout)
186c
187         CALL histdef(nid_mth, "SWdnTOAclr", 
188     .                "SWdn clear sky at TOA","W/m2",
189     .                iim,jjmp1,nhori, 1,1,1,-99,
190     .                32, "ave(X)", zsto2,zout)
191c
192         CALL histdef(nid_mth, "SWup200", "SWup at 200hPa","W/m2",
193     .                iim,jjmp1,nhori, 1,1,1,-99,
194     .                32, "ave(X)", zsto2,zout)
195c
196         CALL histdef(nid_mth, "SWup200clr", 
197     .                "SWup clear sky at 200hPa","W/m2",
198     .                iim,jjmp1,nhori, 1,1,1,-99,
199     .                32, "ave(X)", zsto2,zout)
200c
201         CALL histdef(nid_mth, "SWdn200", "SWdn at 200hPa","W/m2",
202     .                iim,jjmp1,nhori, 1,1,1,-99,
203     .                32, "ave(X)", zsto2,zout)
204c
205         CALL histdef(nid_mth, "SWdn200clr", 
206     .                "SWdn clear sky at 200hPa","W/m2",
207     .                iim,jjmp1,nhori, 1,1,1,-99,
208     .                32, "ave(X)", zsto2,zout)
209c
210         CALL histdef(nid_mth, "LWup200", "LWup at 200hPa","W/m2",
211     .                iim,jjmp1,nhori, 1,1,1,-99,
212     .                32, "ave(X)", zsto2,zout)
213c
214         CALL histdef(nid_mth, "LWup200clr", 
215     .                "LWup clear sky at 200hPa","W/m2",
216     .                iim,jjmp1,nhori, 1,1,1,-99,
217     .                32, "ave(X)", zsto2,zout)
218c
219         CALL histdef(nid_mth, "LWdn200", "LWdn at 200hPa","W/m2",
220     .                iim,jjmp1,nhori, 1,1,1,-99,
221     .                32, "ave(X)", zsto2,zout)
222c
223         CALL histdef(nid_mth, "LWdn200clr", 
224     .                "LWdn clear sky at 200hPa","W/m2",
225     .                iim,jjmp1,nhori, 1,1,1,-99,
226     .                32, "ave(X)", zsto2,zout)
227c
228         CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2",
229     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
230     .                "ave(X)", zsto,zout)
231c
232         CALL histdef(nid_mth, "sols0", "Solar rad. at surf.", "W/m2",
233     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
234     .                "ave(X)", zsto,zout)
235c
236         CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2",
237     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
238     .                "ave(X)", zsto,zout)
239c
240         CALL histdef(nid_mth, "soll0", "IR rad. at surface", "W/m2",
241     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
242     .                "ave(X)", zsto,zout)
243c
244         CALL histdef(nid_mth, "SWupSFC", "SWup at surface","W/m2",
245     .                iim,jjmp1,nhori, 1,1,1,-99,
246     .                32, "ave(X)", zsto2,zout)
247c
248         CALL histdef(nid_mth, "SWupSFCclr", 
249     .                "SWup clear sky at surface","W/m2",
250     .                iim,jjmp1,nhori, 1,1,1,-99,
251     .                32, "ave(X)", zsto2,zout)
252c
253         CALL histdef(nid_mth, "SWdnSFC", "SWdn at surface","W/m2",
254     .                iim,jjmp1,nhori, 1,1,1,-99,
255     .                32, "ave(X)", zsto2,zout)
256c
257         CALL histdef(nid_mth, "SWdnSFCclr", 
258     .                "SWdn clear sky at surface","W/m2",
259     .                iim,jjmp1,nhori, 1,1,1,-99,
260     .                32, "ave(X)", zsto2,zout)
261c
262         CALL histdef(nid_mth,"LWupSFC","Upwd. IR rad. at surface", 
263     .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32, 
264     .                "ave(X)", zsto2,zout)
265c
266         CALL histdef(nid_mth,"LWdnSFC","Down. IR rad. at surface", 
267     .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32, 
268     .                "ave(X)", zsto2,zout)
269c
270         CALL histdef(nid_mth,"LWupSFCclr",
271     .                "CS Upwd. IR rad. at surface", 
272     .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32, 
273     .                "ave(X)", zsto2,zout)
274c
275         CALL histdef(nid_mth,"LWdnSFCclr",
276     .                "Down. CS IR rad. at surface", 
277     .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32, 
278     .                "ave(X)", zsto2,zout)
279c
280         CALL histdef(nid_mth, "bils", "Surf. total heat flux", "W/m2",
281     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
282     .                "ave(X)", zsto,zout)
283c
284         CALL histdef(nid_mth, "sens", "Sensible heat flux", "W/m2",
285     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
286     .                "ave(X)", zsto,zout)
287c
288         CALL histdef(nid_mth, "fder", "Heat flux derivation", "W/m2",
289     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
290     .                "ave(X)", zsto,zout)
291c
292c          CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa",
293c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
294c    .                "ave(X)", zsto,zout)
295c
296c        CALL histdef(nid_mth, "frtv", "Meridional wind stress", "Pa",
297c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
298c    .                "ave(X)", zsto,zout)
299c
300         CALL histdef(nid_mth, "ffonte","Thermal flux for snow melting",
301     .                "W/m2",iim,jjmp1,nhori, 1,1,1, -99, 32,
302     .                "ave(X)", zsto,zout)
303c
304         CALL histdef(nid_mth, "fqcalving","Ice Calving",
305     .                "Kg/m2/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
306     .                "ave(X)", zsto,zout)
307cIM: 171003
308         DO nsrf = 1, nbsrf
309           call histdef(nid_mth, "taux_"//clnsurf(nsrf),
310     $         "Zonal wind stress"//clnsurf(nsrf), "Pa", 
311     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
312     $         "ave(X)", zsto,zout)
313
314           call histdef(nid_mth, "tauy_"//clnsurf(nsrf),
315     $         "Meridional xind stress "//clnsurf(nsrf), "Pa", 
316     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
317     $         "ave(X)", zsto,zout)
318         ENDDO
319cIM: 171003
320c
321cIM       if ( lev_histday.gt.1 ) then
322         DO nsrf = 1, nbsrf
323C
324          IF(nsrf.GT.1) THEN
325           call histdef(nid_mth, "pourc_"//clnsurf(nsrf),
326     $         "Fraction "//clnsurf(nsrf), "W/m2", 
327     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
328     $         "ave(X)", zsto,zout)
329          ENDIF !nsrf.GT.1
330C
331           call histdef(nid_mth, "tsol_"//clnsurf(nsrf),
332     $         "Fraction "//clnsurf(nsrf), "W/m2", 
333     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
334     $         "ave(X)", zsto,zout)
335C
336           call histdef(nid_mth, "sens_"//clnsurf(nsrf),
337     $         "Sensible heat flux "//clnsurf(nsrf), "W/m2", 
338     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
339     $         "ave(X)", zsto,zout)
340c
341           call histdef(nid_mth, "lat_"//clnsurf(nsrf),
342     $         "Latent heat flux "//clnsurf(nsrf), "W/m2", 
343     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
344     $         "ave(X)", zsto,zout)
345C
346           call histdef(nid_mth, "flw_"//clnsurf(nsrf),
347     $         "LW "//clnsurf(nsrf), "W/m2",
348     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
349     $         "ave(X)", zsto2,zout)
350c
351           call histdef(nid_mth, "fsw_"//clnsurf(nsrf),
352     $         "SW "//clnsurf(nsrf), "W/m2",
353     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
354     $         "ave(X)", zsto2,zout)
355C
356           call histdef(nid_mth, "wbils_"//clnsurf(nsrf),
357     $         "Bilan sol "//clnsurf(nsrf), "W/m2",
358     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
359     $         "ave(X)", zsto,zout)
360C
361         END DO
362cIM         endif !lev_histday
363c
364         CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
365     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
366     .                "ave(X)", zsto,zout)
367c
368         CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
369     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
370     .                "ave(X)", zsto,zout)
371c
372         CALL histdef(nid_mth, "cldl", "Low-level cloudiness", "-",
373     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
374     .                "ave(X)", zsto,zout)
375c
376         CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", "-",
377     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
378     .                "ave(X)", zsto,zout)
379c
380         CALL histdef(nid_mth, "cldh", "High-level cloudiness", "-",
381     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
382     .                "ave(X)", zsto,zout)
383c
384         CALL histdef(nid_mth, "cldt", "Total cloudiness", "-",
385     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
386     .                "ave(X)", zsto,zout)
387c
388         CALL histdef(nid_mth,"cldq","Cloud liquid water path","Kg/m2",
389     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
390     .                "ave(X)", zsto,zout)
391cIM: 071003
392         CALL histdef(nid_mth,"lwp","Cloud water path","Kg/m2",
393     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
394     .                "ave(X)", zsto,zout)
395c
396         CALL histdef(nid_mth,"iwp","Cloud ice water path","Kg/m2",
397     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
398     .                "ave(X)", zsto,zout)
399cIM: 071003
400c
401         CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
402     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
403     .                "ave(X)", zsto,zout)
404c
405         CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
406     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
407     .                "ave(X)", zsto,zout)
408c
409         CALL histdef(nid_mth, "uq", "Zonal humidity transport", "-",
410     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
411     .                "ave(X)", zsto,zout)
412c
413         CALL histdef(nid_mth, "vq", "Merid humidity transport", "-",
414     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
415     .                "ave(X)", zsto,zout)
416cKE43
417      IF(iflag_con.GE.3) THEN ! sb
418c
419         CALL histdef(nid_mth, "cape", "Conv avlbl pot ener", "J/Kg",
420     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
421     .                "ave(X)", zsto,zout)
422c
423         CALL histdef(nid_mth, "pbase", "Cld base pressure", "hPa",
424     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
425     .                "ave(X)", zsto,zout)
426c
427         CALL histdef(nid_mth, "ptop", "Cld top pressure", "hPa",
428     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
429     .                "ave(X)", zsto,zout)
430c
431         CALL histdef(nid_mth, "fbase", "Cld base mass flux", "Kg/m2/s",
432     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
433     .                "ave(X)", zsto,zout)
434c
435         CALL histdef(nid_mth, "prw", "Precipitable water", "Kg/m2",
436     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
437     .                "ave(X)", zsto,zout)
438c
439      ENDIF !iflag_con .GE. 3
440c34EK
441c
442c Champs interpolles sur des niveaux de pression
443         DO k=1, nlevENS
444          CALL histdef(nid_mth, "u"//clev(k),
445     .                 "Zonal wind"//clev(k)//"mb","m/s",
446     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
447     .                "ave(X)", zsto,zout)
448c
449          CALL histdef(nid_mth, "v"//clev(k),
450     .                 "Meridional wind"//clev(k)//"mb","m/s",
451     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
452     .                "ave(X)", zsto,zout)
453c
454          CALL histdef(nid_mth, "w"//clev(k),
455     .                 "Vertical wind"//clev(k)//"mb","m/s",
456     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
457     .                "ave(X)", zsto,zout)
458c
459          CALL histdef(nid_mth, "phi"//clev(k),
460     .                 "Geopotential"//clev(k)//"mb","m2/s2",
461     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
462     .                "ave(X)", zsto,zout)
463c
464         ENDDO
465      ENDIF !lev_histmth.GE.1
466c
467      IF(lev_histmth.GE.2) THEN
468c
469c Champs 3D:
470c
471c
472         CALL histdef(nid_mth,"lwcon","Cloud water content","Kg/Kg",
473     .                iim,jjmp1,nhori, klev,1,klev, nvert, 32, 
474     .                "ave(X)", zsto,zout)
475c
476         CALL histdef(nid_mth,"iwcon","Cloud ice water content","Kg/Kg",
477     .                iim,jjmp1,nhori, klev,1,klev, nvert, 32, 
478     .                "ave(X)", zsto,zout)
479c
480         CALL histdef(nid_mth, "temp", "Air temperature", "K",
481     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
482     .                "ave(X)", zsto,zout)
483c
484         CALL histdef(nid_mth, "ovap", "Specific humidity", "Kg/Kg",
485     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
486     .                "ave(X)", zsto,zout)
487c
488c        CALL histdef(nid_mth,"wvap","Water vapor mixing ratio","Kg/Kg",
489c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
490c    .                "ave(X)", zsto,zout)
491c
492         CALL histdef(nid_mth, "geop", "Geopotential height", "m",
493     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
494     .                "ave(X)", zsto,zout)
495c
496         CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s",
497     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
498     .                "ave(X)", zsto,zout)
499c
500         CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s",
501     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
502     .                "ave(X)", zsto,zout)
503c
504         CALL histdef(nid_mth, "vitw", "Vertical wind", "m/s",
505     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
506     .                "ave(X)", zsto,zout)
507c
508         CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
509     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
510     .                "ave(X)", zsto,zout)
511c
512        CALL histdef(nid_mth, "rneb", "Cloud fraction", "-",
513     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
514     .                "ave(X)", zsto,zout)
515c
516         CALL histdef(nid_mth, "rnebcon", "Convective Cloud Fraction"
517     .                , "-",
518     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
519     .                "ave(X)", zsto,zout)
520c
521         CALL histdef(nid_mth, "rhum", "Relative humidity", "-",
522     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
523     .                "ave(X)", zsto,zout)
524c
525         CALL histdef(nid_mth, "ozone", "Ozone concentration", "-",
526     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
527     .                "ave(X)", zsto,zout)
528c
529         CALL histdef(nid_mth, "upwd", "saturated updraft", "Kg/m2/s",
530     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
531     .                "ave(X)", zsto,zout)
532c
533         CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
534     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
535     .                "ave(X)", zsto,zout)
536c
537         CALL histdef(nid_mth, "dqphy", "Physics dQ", "Kg/Kg/s",
538     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
539     .                "ave(X)", zsto,zout)
540c
541c#define histmthNMC
542c#ifdef histmthNMC
543c#include "ini_histmthNMC.h"
544c#endif
545c
546      ENDIF !lev_histmth.GE.2
547c
548      IF(lev_histmth.GE.3) THEN
549c
550        DO nsrf=1, nbsrf
551c
552           call histdef(nid_mth, "albe_"//clnsurf(nsrf),
553     $         "Albedo surf. "//clnsurf(nsrf), "-", 
554     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
555     $         "ave(X)", zsto,zout)
556c
557           call histdef(nid_mth, "rugs_"//clnsurf(nsrf),
558     $         "Latent heat flux "//clnsurf(nsrf), "W/m2", 
559     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
560     $         "ave(X)", zsto,zout)
561c
562         CALL histdef(nid_mth, "ages_"//clnsurf(nsrf), "Snow age","day",
563     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
564     .                "ave(X)", zsto,zout)
565c
566         ENDDO !nsrf=1, nbsrf
567c
568         CALL histdef(nid_mth, "albs", "Surface albedo", "-",
569     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
570     .                "ave(X)", zsto,zout)
571         CALL histdef(nid_mth, "albslw", "Surface albedo LW", "-",
572     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
573     .                "ave(X)", zsto,zout)
574c
575      ENDIF !lev_histmth.GE.3
576c
577      IF(lev_histmth.GE.4) THEN
578c
579         CALL histdef(nid_mth, "clwcon", 
580     .                "Convective Cloud Liquid water content"
581     .                , "Kg/Kg",
582     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
583     .                "ave(X)", zsto,zout)
584c
585         CALL histdef(nid_mth,"Ma","undilute adiab updraft","Kg/m2/s",
586     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
587     .                "ave(X)", zsto,zout)
588c
589         CALL histdef(nid_mth, "dnwd", "saturated downdraft","Kg/m2/s",
590     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
591     .                "ave(X)", zsto,zout)
592c
593         CALL histdef(nid_mth, "dnwd0", "unsat. downdraft", "Kg/m2/s",
594     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
595     .                "ave(X)", zsto,zout)
596c
597         CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
598     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
599     .                "ave(X)", zsto,zout)
600c
601         CALL histdef(nid_mth, "dqdyn", "Dynamics dQ", "Kg/Kg/s",
602     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
603     .                "ave(X)", zsto,zout)
604c
605         CALL histdef(nid_mth, "dtcon", "Convection dT", "K/s",
606     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
607     .                "ave(X)", zsto,zout)
608c
609c        CALL histdef(nid_mth, "ducon", "Convection du", "m/s2",
610c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
611c    .                "ave(X)", zsto,zout)
612c
613         CALL histdef(nid_mth, "dqcon", "Convection dQ", "Kg/Kg/s",
614     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
615     .                "ave(X)", zsto,zout)
616c
617         CALL histdef(nid_mth, "dtlsc", "Condensation dT", "K/s",
618     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
619     .                "ave(X)", zsto,zout)
620cIM: 071003
621         CALL histdef(nid_mth, "dtlschr",
622     $       "Large-scale condensational heating rate", "K/s",iim,jjmp1
623     $       ,nhori, klev,1,klev,nvert, 32,"ave(X)", zsto,zout)
624c
625         CALL histdef(nid_mth, "dqlsc", "Condensation dQ", "Kg/Kg/s",
626     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
627     .                "ave(X)", zsto,zout)
628c
629         CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s",
630     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
631     .                "ave(X)", zsto,zout)
632c
633         CALL histdef(nid_mth, "dqvdf", "Boundary-layer dQ", "Kg/Kg/s",
634     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
635     .                "ave(X)", zsto,zout)
636c
637         CALL histdef(nid_mth, "dteva", "Reevaporation dT", "K/s",
638     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
639     .                "ave(X)", zsto,zout)
640c
641         CALL histdef(nid_mth, "dqeva", "Reevaporation dQ", "Kg/Kg/s",
642     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
643     .                "ave(X)", zsto,zout)
644
645         CALL histdef(nid_mth, "ptconv", "POINTS CONVECTIFS"," ",
646     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
647     .                "ave(X)", zsto,zout)
648
649         CALL histdef(nid_mth, "ratqs", "RATQS"," ",
650     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
651     .                "ave(X)", zsto,zout)
652
653c
654         CALL histdef(nid_mth, "dtajs", "Dry adjust. dT", "K/s",
655     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
656     .                "ave(X)", zsto,zout)
657
658         CALL histdef(nid_mth, "dqajs", "Dry adjust. dQ", "Kg/Kg/s",
659     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
660     .                "ave(X)", zsto,zout)
661c
662         CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s",
663     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
664     .                "ave(X)", zsto,zout)
665c
666         CALL histdef(nid_mth, "dtsw0", "CS SW radiation dT", "K/s",
667     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
668     .                "ave(X)", zsto,zout)
669c
670         CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
671     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
672     .                "ave(X)", zsto,zout)
673c
674         CALL histdef(nid_mth,"dtlw0","CS LW radiation dT","K/s",
675     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
676     .                "ave(X)", zsto,zout)
677c
678         CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
679     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
680     .                "ave(X)", zsto,zout)
681c
682         CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
683     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
684     .                "ave(X)", zsto,zout)
685c
686         CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
687     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
688     .                "ave(X)", zsto,zout)
689c
690         IF (ok_orodr) THEN
691         CALL histdef(nid_mth, "duoro", "Orography dU", "m/s2",
692     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
693     .                "ave(X)", zsto,zout)
694c
695         CALL histdef(nid_mth, "dvoro", "Orography dV", "m/s2",
696     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
697     .                "ave(X)", zsto,zout)
698c
699         ENDIF
700C
701         IF (ok_orolf) THEN
702         CALL histdef(nid_mth, "dulif", "Orography dU", "m/s2",
703     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
704     .                "ave(X)", zsto,zout)
705c
706         CALL histdef(nid_mth, "dvlif", "Orography dV", "m/s2",
707     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
708     .                "ave(X)", zsto,zout)
709         ENDIF
710C
711         if (nqmax.GE.3) THEN
712         DO iq=1,nqmax-2
713         IF (iq.LE.99) THEN
714         WRITE(str2,'(i2.2)') iq
715         CALL histdef(nid_mth, "trac"//str2, "Tracer No."//str2, "-",
716     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
717     .                "ave(X)", zsto,zout)
718         ELSE
719         PRINT*, "Trop de traceurs"
720         CALL abort
721         ENDIF
722         ENDDO
723         ENDIF
724c
725      ENDIF !lev_histmth.GE.4
726c
727         CALL histend(nid_mth)
728c
729         ndex2d = 0
730         ndex3d = 0
731c
732      ENDIF ! fin de test sur ok_mensuel
Note: See TracBrowser for help on using the repository browser.