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

Last change on this file since 463 was 463, checked in by lmdzadmin, 21 years ago

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