source: LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histday.h @ 5456

Last change on this file since 5456 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: 17.7 KB
Line 
1      IF (ok_journe) THEN
2c
3         zsto = dtime
4         zout = dtime * FLOAT(ecrit_day)
5         zsto1= dtime * FLOAT(ecrit_day)
6c         zout = dtime * REAL(ecrit_day)
7c         zsto1= dtime * REAL(ecrit_day)
8c
9         idayref = day_ref
10         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
11c
12         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
13         DO i = 1, iim
14            zx_lon(i,1) = rlon(i+1)
15            zx_lon(i,jjmp1) = rlon(i+1)
16         ENDDO
17         DO ll=1,klev
18            znivsig(ll)=float(ll)
19         ENDDO
20         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
21         CALL histbeg("histday", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
22     .                 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 
23     .                 nhori, nid_day)
24         write(*,*)'Journee ', itau_phy, zjulian
25         CALL histvert(nid_day, "presnivs", "Vertical levels", "mb",
26     .                 klev, presnivs, nvert)
27c
28       IF(lev_histday.GE.1) THEN
29c
30         CALL histdef(nid_day, "phis", "Surface geop. height", "-",
31     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
32     .                "once", zsto,zout)
33c
34         CALL histdef(nid_day, "aire", "Grid area", "-",
35     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
36     .                "once", zsto,zout)
37c
38         CALL histdef(nid_day, "contfracATM","% sfce ter+lic ","-",
39     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
40     .                "once", zsto,zout)
41c
42         CALL histdef(nid_day, "contfracOR","% sfce terre OR", "-",
43     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
44     .                "once", zsto,zout)
45c
46c Champs 2D:
47c
48         CALL histdef(nid_day, "tsol", "Surface Temperature", "K",
49     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
50     .                "ave(X)", zsto,zout)
51c
52         CALL histdef(nid_day, "t2m", "Temperature 2m", "K",
53     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
54     .                "ave(X)", zsto,zout)
55c
56c Champs retires momentannéement en attendant un hypothetique
57c debugage
58
59         CALL histdef(nid_day, "t2m_min", "Temp. 2m min.",
60     .                "K",
61     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
62     .                t2mincels, zsto,zout)
63c
64         CALL histdef(nid_day, "t2m_max", "Temp. 2m max.",
65     .                "K",
66     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
67     .                t2maxcels, zsto,zout)
68c
69         CALL histdef(nid_day, "plul", "Large-scale Precip.",
70     .   "kg/(s*m2)",
71     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
72     .                "ave(X)", zsto,zout)
73c
74         CALL histdef(nid_day, "pluc", "Convective Precip.",
75     .   "kg/(s*m2)",
76     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
77     .                "ave(X)", zsto,zout)
78c
79         CALL histdef(nid_day, "snowl", "Solid Large-scale Precip.",
80     .   "kg/(s*m2)",
81     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
82     .                "ave(X)", zsto,zout)
83c
84         IF(1.EQ.0) THEN !snowc=0.
85         CALL histdef(nid_day, "snowc", "Solid Convective Precip.",
86     .   "kg/(s*m2)",
87     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
88     .                "ave(X)", zsto,zout)
89         ENDIF !snowc=0.
90c
91         CALL histdef(nid_day, "flat", "Latent heat flux", "W/m2",
92     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
93     .                "ave(X)", zsto,zout)
94c
95         CALL histdef(nid_day, "sicf", "Sea-ice fraction", "-",
96     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
97     .                "inst(X)", zsto1,zout)
98c
99         CALL histdef(nid_day, "q2m", "Specific humidity", "kg/kg",
100     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
101     .                "ave(X)", zsto,zout)
102c
103         CALL histdef(nid_day, "u10m", "Vent zonal 10m", "m/s",
104     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
105     .                "ave(X)", zsto,zout)
106c
107         CALL histdef(nid_day, "v10m", "Vent meridien 10m", "m/s",
108     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
109     .                "ave(X)", zsto,zout)
110c
111         CALL histdef(nid_day, "wind10m","10-m wind speed","m/s",
112     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
113     .                "ave(X)", zsto,zout)
114c
115         CALL histdef(nid_day, "wind10max", "10-m wind speed max.",
116     .                "m/s",
117     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
118     .                t2maxcels, zsto,zout)
119c
120         CALL histdef(nid_day, "psol", "Surface Pressure", "Pa",
121     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
122     .                "ave(X)", zsto,zout)
123c
124         CALL histdef(nid_day, "precip","Precipitation Totale liq+sol"
125     .                , "kg/(s*m2)",
126     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
127     .                "ave(X)", zsto,zout)
128c
129         CALL histdef(nid_day, "snowf", "Snow fall", "kg/(s*m2)",
130     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
131     .                "ave(X)", zsto,zout)
132c
133         CALL histdef(nid_day, "snow_mass", "Snow Mass", "kg/m2",
134     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
135     .                "ave(X)", zsto,zout)
136c
137         CALL histdef(nid_day, "evap", "Evaporation", "kg/(s*m2)",
138     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
139     .                "ave(X)", zsto,zout)
140c
141         CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",
142     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
143     .                "ave(X)", zsto,zout)
144c
145         CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2",
146     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
147     .                "ave(X)", zsto,zout)
148c
149         CALL histdef(nid_day, "sols", "Net Solar rad. at surf.", 
150     .                "W/m2",
151     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
152     .                "ave(X)", zsto,zout)
153c
154         CALL histdef(nid_day, "soll", "Net IR rad. at surface", "W/m2",
155     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
156     .                "ave(X)", zsto,zout)
157c
158         CALL histdef(nid_day, "cldl", "Low-level cloudiness", "-",
159     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
160     .                "ave(X)", zsto,zout)
161c
162         CALL histdef(nid_day, "cldm", "Mid-level cloudiness", "-",
163     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
164     .                "ave(X)", zsto,zout)
165c
166         CALL histdef(nid_day, "cldh", "High-level cloudiness", "-",
167     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
168     .                "ave(X)", zsto,zout)
169c
170         CALL histdef(nid_day, "cldt", "Total cloudiness", "-",
171     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
172     .                "ave(X)", zsto,zout)
173c
174         CALL histdef(nid_day, "cldq", "Cloud liquid water path", 
175     .                "kg/m2",
176     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
177     .                "ave(X)", zsto,zout)
178c
179         CALL histdef(nid_day, "prw", "Precipitable water", "kg/m2",
180     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
181     .                "ave(X)", zsto,zout)
182c
183c  Champs dynamiques sur niveaux de pression
184         DO k=1, nlevENS
185c
186          CALL histdef(nid_day, "u"//clev(k),
187     .                 "Zonal wind"//clev(k)//"mb","m/s",
188     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
189     .                "ave(X)", zsto,zout)
190c
191          CALL histdef(nid_day, "v"//clev(k),
192     .                 "Meridional wind"//clev(k)//"mb","m/s",
193     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
194     .                "ave(X)", zsto,zout)
195c
196         ENDDO !nlevENS
197c
198         CALL histdef(nid_day, "w500", "Verical wind 500mb", "m/s",
199     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
200     .                "ave(X)", zsto,zout)
201c
202         CALL histdef(nid_day,"phi500", "Geopotentiel a 500mb","m2/s2",
203     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
204     .                "ave(X)", zsto,zout)
205c
206         CALL histdef(nid_day, "slp", "Sea Level Pressure", "Pa",
207     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
208     .                "ave(X)", zsto,zout)
209c
210         CALL histdef(nid_day, "cape_max", "CAPE max.",
211     .                "J/kg",
212     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
213     .                capemaxcels, zsto,zout)
214c
215         CALL histdef(nid_day, "solldown", "Down. IR rad. at surface", 
216     .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32, 
217     .                "ave(X)", zsto,zout)
218c
219         CALL histdef(nid_day, "sens", "Sensible heat flux", "W/m2",
220     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
221     .                "ave(X)", zsto,zout)
222c
223         CALL histdef(nid_day, "SWdnSFC", "SWdn at surface","W/m2",
224     .                iim,jjmp1,nhori, 1,1,1,-99,
225     .                32, "ave(X)", zsto,zout)
226c
227      ENDIF !lev_histday.GE.1
228c
229      IF (lev_histday.GE.2) THEN
230c
231         CALL histdef(nid_day, "bils", "Surf. total heat flux", "W/m2",
232     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
233     .                "ave(X)", zsto,zout)
234c
235         CALL histdef(nid_day, "fder", "Heat flux derivation", "W/m2/K",
236     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
237     .                "ave(X)", zsto,zout)
238c
239      ENDIF !lev_histday.GE.2
240c=================================================================
241      IF(lev_histday.GE.3) THEN
242c=================================================================
243c   INITIALISATION DES CHAMPS 3D
244c=================================================================
245c=================================================================
246c=================================================================
247c Champs 3D:
248c
249         CALL histdef(nid_day, "temp", "Air temperature", "K",
250     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
251     .                "ave(X)", zsto,zout)
252c
253         CALL histdef(nid_day, "ovap", "Specific humidity", "kg/kg",
254     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
255     .                "ave(X)", zsto,zout)
256c
257         CALL histdef(nid_day, "geop", "Geopotential height", "m",
258     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
259     .                "ave(X)", zsto,zout)
260c
261         CALL histdef(nid_day, "vitu", "Zonal wind", "m/s",
262     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
263     .                "ave(X)", zsto,zout)
264c
265         CALL histdef(nid_day, "vitv", "Meridional wind", "m/s",
266     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
267     .                "ave(X)", zsto,zout)
268c
269         CALL histdef(nid_day, "vitw", "Vertical wind", "Pa/s",
270     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
271     .                "ave(X)", zsto,zout)
272c
273         CALL histdef(nid_day, "pres", "Air pressure", "Pa",
274     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
275     .                "ave(X)", zsto,zout)
276c
277c=================================================================
278c   FIN INITIALISATION DES CHAMPS 3D
279c=================================================================
280      ENDIF !  lev_histday.GE.3
281c
282      IF (lev_histday.GE.4) THEN
283c=================================================================
284c
285c   INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
286c
287c=================================================================
288c
289         CALL histdef(nid_day, "SWupTOA", "SWup at TOA","W/m2",
290     .                iim,jjmp1,nhori, 1,1,1,-99,
291     .                32, "ave(X)", zsto,zout)
292c
293         CALL histdef(nid_day, "SWupSFC", "SWup at surface","W/m2",
294     .                iim,jjmp1,nhori, 1,1,1,-99,
295     .                32, "ave(X)", zsto,zout)
296c
297         CALL histdef(nid_day, "SWdnTOA", "SWdn at TOA","W/m2",
298     .                iim,jjmp1,nhori, 1,1,1,-99,
299     .                32, "ave(X)", zsto,zout)
300c
301         CALL histdef(nid_day, "SWupTOAclr",
302     .                "SWup clear sky at TOA","W/m2",
303     .                iim,jjmp1,nhori, 1,1,1,-99,
304     .                32, "ave(X)", zsto,zout)
305
306         CALL histdef(nid_day, "SWupSFCclr",
307     .                "SWup clear sky at surface","W/m2",
308     .                iim,jjmp1,nhori, 1,1,1,-99,
309     .                32, "ave(X)", zsto,zout)
310
311         CALL histdef(nid_day, "SWdnTOAclr",
312     .                "SWdn clear sky at TOA","W/m2",
313     .                iim,jjmp1,nhori, 1,1,1,-99,
314     .                32, "ave(X)", zsto,zout)
315
316         CALL histdef(nid_day, "SWdnSFCclr",
317     .                "SWdn clear sky at surface","W/m2",
318     .                iim,jjmp1,nhori, 1,1,1,-99,
319     .                32, "ave(X)", zsto,zout)
320c
321c=================================================================
322c=================================================================
323c=================================================================
324c   INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
325c=================================================================
326c
327         CALL histdef(nid_day, "tter", "Surface Temperature", "K",
328     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
329     .                "ave(X)", zsto,zout)
330c
331         CALL histdef(nid_day, "tlic", "Surface Temperature", "K",
332     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
333     .                "ave(X)", zsto,zout)
334c
335         CALL histdef(nid_day, "toce", "Surface Temperature", "K",
336     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
337     .                "ave(X)", zsto,zout)
338c
339         CALL histdef(nid_day, "tsic", "Surface Temperature", "K",
340     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
341     .                "ave(X)", zsto,zout)
342c
343         CALL histdef(nid_day, "t2mter", "Temp.terre 2m", "K",
344     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
345     .                "ave(X)", zsto,zout)
346c
347         CALL histdef(nid_day, "t2mlic", "Temp.lic 2m", "K",
348     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
349     .                "ave(X)", zsto,zout)
350c
351         CALL histdef(nid_day, "t2moce", "Temp.oce 2m", "K",
352     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
353     .                "ave(X)", zsto,zout)
354c
355         CALL histdef(nid_day, "t2msic", "Temp.sic 2m", "K",
356     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
357     .                "ave(X)", zsto,zout)
358c
359         CALL histdef(nid_day, "t2mter_min", "Temp.terre 2m min.",
360     .                "K",
361     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
362     .                t2mincels, zsto,zout)
363c
364         CALL histdef(nid_day, "t2mter_max", "Temp.terre 2m max.",
365     .                "K",
366     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
367     .                t2maxcels, zsto,zout)
368c
369         CALL histdef(nid_day, "u10mter", "Vent zonal ter 10m", "m/s",
370     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
371     .                "ave(X)", zsto,zout)
372c
373         CALL histdef(nid_day, "u10mlic", "Vent zonal lic 10m", "m/s",
374     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
375     .                "ave(X)", zsto,zout)
376c
377         CALL histdef(nid_day, "u10moce", "Vent zonal oce 10m", "m/s",
378     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
379     .                "ave(X)", zsto,zout)
380c
381         CALL histdef(nid_day, "u10msic", "Vent zonal sic 10m",
382     .                 "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
383     .                "ave(X)", zsto,zout)
384c
385         CALL histdef(nid_day, "v10mter", "Vent meridien ter 10m", 
386     .                "m/s", iim,jjmp1,nhori, 1,1,1, -99, 32,
387     .                "ave(X)", zsto,zout)
388c
389         CALL histdef(nid_day, "v10mlic", "Vent meridien lic 10m",
390     .                 "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
391     .                "ave(X)", zsto,zout)
392c
393         CALL histdef(nid_day, "v10moce", "Vent meridien oce 10m",
394     .                 "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
395     .                "ave(X)", zsto,zout)
396c
397         CALL histdef(nid_day, "v10msic", "Vent meridien sic 10m",
398     .                 "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
399     .                "ave(X)", zsto,zout)
400c
401         DO nsrf = 1, nbsrf
402C
403           call histdef(nid_day, "pourc_"//clnsurf(nsrf),
404     $         "Fraction"//clnsurf(nsrf), "W/m2", 
405     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
406     $         "ave(X)", zsto,zout)
407C
408           call histdef(nid_day, "tsol_"//clnsurf(nsrf),
409     $         "Fraction"//clnsurf(nsrf), "W/m2", 
410     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
411     $         "ave(X)", zsto,zout)
412C
413           call histdef(nid_day, "sens_"//clnsurf(nsrf),
414     $         "Sensible heat flux "//clnsurf(nsrf), "W/m2", 
415     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
416     $         "ave(X)", zsto,zout)
417c
418           call histdef(nid_day, "lat_"//clnsurf(nsrf),
419     $         "Latent heat flux "//clnsurf(nsrf), "W/m2", 
420     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
421     $         "ave(X)", zsto,zout)
422C
423           call histdef(nid_day, "taux_"//clnsurf(nsrf),
424     $         "Zonal wind stress"//clnsurf(nsrf),"Pa",
425     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
426     $         "ave(X)", zsto,zout)
427
428           call histdef(nid_day, "tauy_"//clnsurf(nsrf),
429     $         "Meridional xind stress "//clnsurf(nsrf), "Pa", 
430     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
431     $         "ave(X)", zsto,zout)
432C
433           call histdef(nid_day, "albe_"//clnsurf(nsrf),
434     $         "Albedo surf. "//clnsurf(nsrf), "W/m2", 
435     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
436     $         "ave(X)", zsto,zout)
437C
438           call histdef(nid_day, "rugs_"//clnsurf(nsrf),
439     $         "Latent heat flux "//clnsurf(nsrf), "W/m2", 
440     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
441     $         "ave(X)", zsto,zout)
442C
443         END DO
444C           
445c=================================================================
446c
447c   FIN INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
448c
449c=================================================================
450      ENDIF !lev_histday.GE.4
451c=================================================================
452c
453         CALL histend(nid_day)
454c
455         ndex2d = 0
456         ndex3d = 0
457c
458c=================================================================
459      ENDIF ! fin de test sur ok_journe
Note: See TracBrowser for help on using the repository browser.