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

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

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