source: LMDZ4/trunk/libf/phylmd/ini_histday.h @ 895

Last change on this file since 895 was 895, checked in by Laurent Fairhead, 16 years ago

Un peu de menage sur les fichiers de sortie apres les tests aquaplanete FH
LF
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 28.1 KB
Line 
1c
2c $Header$
3c
4      IF (ok_journe) THEN
5
6c$OMP MASTER
7c
8       zstophy = dtime
9       zstoday= ecrit_day
10       zout = ecrit_day
11c
12         idayref = day_ref
13         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
14c
15cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
16cym         DO i = 1, iim
17cym            zx_lon(i,1) = rlon(i+1)
18cym            zx_lon(i,jjmp1) = rlon(i+1)
19cym         ENDDO
20         DO ll=1,klev
21            znivsig(ll)=float(ll)
22         ENDDO
23cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
24cym         CALL histbeg("histday", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
25cym     .                 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 
26cym     .                 nhori, nid_day)
27
28         CALL histbeg_phy("histday",itau_phy, zjulian, dtime, 
29     .                    nhori, nid_day)
30
31         write(*,*)'Journee ', itau_phy, zjulian
32         CALL histvert(nid_day, "presnivs", "Vertical levels", "mb",
33     .                 klev, presnivs/100., nvert)
34c
35       IF(lev_histday.GE.1) THEN
36c
37         CALL histdef(nid_day, "phis", "Surface geop. height", "m2/s2",
38     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
39     .                "once", zstophy,zout)
40c
41         CALL histdef(nid_day, "aire", "Grid area", "-",
42     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
43     .                "once", zstophy,zout)
44c
45         CALL histdef(nid_day, "contfracATM","% sfce ter+lic ","-",
46     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
47     .                "once", zstophy,zout)
48c
49         CALL histdef(nid_day, "contfracOR","% sfce terre OR", "-",
50     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
51     .                "once", zstophy,zout)
52c
53c Champs 2D:
54c
55         CALL histdef(nid_day, "weakinv", "Weak inversion", "",
56     s           iim,jjmp1,nhori, 1,1,1, -99, 32,
57     s           "ave(X)", zstophy,zout)
58
59         CALL histdef(nid_day, "dthmin", "dTheta mini", "K/m",
60     s           iim,jjmp1,nhori, 1,1,1, -99, 32,
61     s           "ave(X)", zstophy,zout)
62
63
64         CALL histdef(nid_day, "tsol", "Surface Temperature", "K",
65     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
66     .                "ave(X)", zstophy,zout)
67c
68         CALL histdef(nid_day, "t2m", "Temperature 2m", "K",
69     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
70     .                "ave(X)", zstophy,zout)
71c
72c Champs retires momentannement en attendant un hypothetique
73c debugage
74
75         CALL histdef(nid_day, "t2m_min", "Temp. 2m min.",
76     .                "K",
77     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
78     .                t2mincels, zstophy,zout)
79c
80         CALL histdef(nid_day, "t2m_max", "Temp. 2m max.",
81     .                "K",
82     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
83     .                t2maxcels, zstophy,zout)
84c
85!        CALL histdef(nid_day, "plul", "Large-scale Precip.",
86!    .   "kg/(s*m2)",
87!    .                iim,jj_nb,nhori, 1,1,1, -99, 32,
88!    .                "ave(X)", zstophy,zout)
89c
90         CALL histdef(nid_day, "pluc", "Convective Precip.",
91     .   "kg/(s*m2)",
92     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
93     .                "ave(X)", zstophy,zout)
94c
95         CALL histdef(nid_day, "snowl", "Solid Large-scale Precip.",
96     .   "kg/(m2*s)",
97     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
98     .                "ave(X)", zstophy,zout)
99c
100         IF(1.EQ.0) THEN !snowc=0.
101         CALL histdef(nid_day, "snowc", "Solid Convective Precip.",
102     .   "kg/(m2*s)",
103     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
104     .                "ave(X)", zstophy,zout)
105         ENDIF !snowc=0.
106c
107         CALL histdef(nid_day, "flat", "Latent heat flux", "W/m2",
108     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
109     .                "ave(X)", zstophy,zout)
110c
111         CALL histdef(nid_day, "sicf", "Sea-ice fraction", "-",
112     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
113     .                "inst(X)", zstoday,zout)
114c
115         CALL histdef(nid_day, "q2m", "Specific humidity", "kg/kg",
116     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
117     .                "ave(X)", zstophy,zout)
118c
119         CALL histdef(nid_day, "u10m", "Vent zonal 10m", "m/s",
120     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
121     .                "ave(X)", zstophy,zout)
122c
123         CALL histdef(nid_day, "v10m", "Vent meridien 10m", "m/s",
124     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
125     .                "ave(X)", zstophy,zout)
126c
127         CALL histdef(nid_day, "wind10m","10-m wind speed","m/s",
128     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
129     .                "ave(X)", zstophy,zout)
130c
131         CALL histdef(nid_day, "wind10max", "10-m wind speed max.",
132     .                "m/s",
133     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
134     .                t2maxcels, zstophy,zout)
135c
136         CALL histdef(nid_day, "psol", "Surface Pressure", "Pa",
137     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
138     .                "ave(X)", zstophy,zout)
139c
140         CALL histdef(nid_day, "precip","Precipitation Totale liq+sol"
141     .                , "kg/(s*m2)",
142     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
143     .                "ave(X)", zstophy,zout)
144c
145         CALL histdef(nid_day, "plul","Precipitation ls liq+sol"
146     .                , "kg/(s*m2)",
147     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
148     .                "ave(X)", zstophy,zout)
149c
150         CALL histdef(nid_day, "snowf", "Snow fall", "kg/(m2*s)",
151     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
152     .                "ave(X)", zstophy,zout)
153c
154c        CALL histdef(nid_day, "snow_mass", "Snow Mass", "kg/m2",
155c    .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
156c    .                "ave(X)", zstophy,zout)
157c
158         CALL histdef(nid_day, "evap", "Evaporation", "kg/(m2*s)",
159     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
160     .                "ave(X)", zstophy,zout)
161c
162         CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",
163     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
164     .                "ave(X)", zstophy,zout)
165c
166         CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2",
167     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
168     .                "ave(X)", zstophy,zout)
169c
170         CALL histdef(nid_day, "sols", "Net Solar rad. at surf.", 
171     .                "W/m2",
172     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
173     .                "ave(X)", zstophy,zout)
174c
175         CALL histdef(nid_day, "soll", "Net IR rad. at surface", "W/m2",
176     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
177     .                "ave(X)", zstophy,zout)
178c
179         CALL histdef(nid_day, "radsol", "Rayonnement au sol", "W/m2",
180     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
181     .                "ave(X)", zstophy,zout)
182c
183         CALL histdef(nid_day, "cldl", "Low-level cloudiness", "-",
184     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
185     .                "ave(X)", zstophy,zout)
186c
187         CALL histdef(nid_day, "cldm", "Mid-level cloudiness", "-",
188     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
189     .                "ave(X)", zstophy,zout)
190c
191         CALL histdef(nid_day, "cldh", "High-level cloudiness", "-",
192     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
193     .                "ave(X)", zstophy,zout)
194c
195         CALL histdef(nid_day, "cldt", "Total cloudiness", "-",
196     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
197     .                "ave(X)", zstophy,zout)
198c
199         CALL histdef(nid_day, "cldq", "Cloud liquid water path", 
200     .                "kg/m2",
201     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
202     .                "ave(X)", zstophy,zout)
203c
204         CALL histdef(nid_day, "prw", "Precipitable water", "kg/m2",
205     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
206     .                "ave(X)", zstophy,zout)
207c
208c  Champs dynamiques sur niveaux de pression
209c
210c        DO k=1, nlevENS
211         DO k=1, nlevSTD
212c
213c        bb=clevSTD(k)
214c
215         IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
216         IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
217c         aa=clevSTD(k)
218c         bb2=clevSTD(k)
219c        ENDIF
220c
221         IF(bb2.EQ."850".OR.bb2.EQ."700"
222     .   .OR.bb2.EQ."500".OR.bb2.EQ."200") THEN
223c
224          CALL histdef(nid_day, "u"//bb2,
225     .                 "Zonal wind "//bb2//"mb","m/s",
226     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
227     .                "inst(X)", zout,zout)
228c
229          CALL histdef(nid_day, "v"//bb2,
230     .                 "Meridional wind "//bb2//"mb","m/s",
231     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
232     .                "inst(X)", zout,zout)
233c
234         ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.bb2.EQ."500".OR.bb2.EQ."200") THEN
235c
236         IF(bb2.EQ."500") THEN
237c
238         CALL histdef(nid_day, "w"//bb2, "Vertical wind "//bb2//"mb",
239     .                "m/s",
240     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
241     .                "inst(X)", zout,zout)
242c
243         CALL histdef(nid_day,"phi"//bb2, "Geopotential "//bb2//"mb",
244     .                "m",
245     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
246     .                "inst(X)", zout,zout)
247c
248         ENDIF !(bb2.EQ."500")
249c
250         ENDDO !nlevSTD
251c
252         CALL histdef(nid_day, "slp", "Sea Level Pressure", "Pa",
253     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
254     .                "ave(X)", zstophy,zout)
255c
256         CALL histdef(nid_day, "cape_max", "CAPE max.",
257     .                "J/kg",
258     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
259     .                capemaxcels, zstophy,zout)
260c
261         CALL histdef(nid_day, "solldown", "Down. IR rad. at surface", 
262     .                "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32, 
263     .                "ave(X)", zstophy,zout)
264c
265         CALL histdef(nid_day, "sens", "Sensible heat flux", "W/m2",
266     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
267     .                "ave(X)", zstophy,zout)
268c
269         CALL histdef(nid_day, "SWdnSFC", "SWdn at surface","W/m2",
270     .                iim,jj_nb,nhori, 1,1,1,-99,
271     .                32, "ave(X)", zstophy,zout)
272c
273cIM diagnostiques flux ocean-atm ou ocean-glace de mer
274cIM pour utilisation dans un modele de "slab" ocean
275c
276         IF (OCEAN.EQ.'force ') THEN
277c
278          CALL histdef(nid_day, "lmt_bils", 
279     .         "Bilan au sol atmosphere forcee", "W/m2",
280     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
281     .                "ave(X)", zstophy,zout)
282c
283         ELSE IF (OCEAN.EQ.'slab  ') THEN
284c
285          CALL histdef(nid_day, "slab_bils",
286     .                "Bilan au sol slab", "W/m2",
287     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
288     .                "ave(X)", zstophy,zout)
289c
290          CALL histdef(nid_day, "tslab", "Slab SST ", "K",
291     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
292     .                "ave(X)", zstophy,zout)
293c
294          CALL histdef(nid_day, "seaice", "Slab seaice", "kg/m2",
295     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
296     .                "ave(X)", zstophy,zout)
297c
298          CALL histdef(nid_day, "siceh", "Slab seaice height", "m",
299     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
300     .                "ave(X)", zstophy,zout)
301c
302         ENDIF !OCEAN.EQ.force/slab
303c
304         CALL histdef(nid_day, "fluxo",
305     .                "Flux turbulents ocean-atmosphere", "W/m2",
306     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
307     .                "ave(X)", zstophy,zout)
308c
309         CALL histdef(nid_day, "fluxg",
310     .                "Flux ocean-glace de mer", "W/m2",
311     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
312     .                "ave(X)", zstophy,zout)
313c
314      ENDIF !lev_histday.GE.1
315c
316      IF (lev_histday.GE.2) THEN
317c
318         CALL histdef(nid_day, "bils", "Surf. total heat flux", "W/m2",
319     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
320     .                "ave(X)", zstophy,zout)
321c
322         CALL histdef(nid_day, "fder", "Heat flux derivation", "W/m2/K",
323     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
324     .                "ave(X)", zstophy,zout)
325c
326      ENDIF !lev_histday.GE.2
327c=================================================================
328      IF(lev_histday.GE.3) THEN
329c=================================================================
330c   INITIALISATION DES CHAMPS 3D
331c=================================================================
332c=================================================================
333c=================================================================
334c Champs 3D:
335c
336         CALL histdef(nid_day, "temp", "Air temperature", "K",
337     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
338     .                "ave(X)", zstophy,zout)
339c
340         CALL histdef(nid_day, "rhum", "relative humidity", " ",
341     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
342     .                "ave(X)", zstophy,zout)
343c
344         CALL histdef(nid_day, "ovap", "Specific humidity", "kg/kg",
345     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
346     .                "ave(X)", zstophy,zout)
347c
348         CALL histdef(nid_day, "geop", "Geopotential height", "m2/s2",
349     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
350     .                "ave(X)", zstophy,zout)
351c
352         CALL histdef(nid_day, "vitu", "Zonal wind", "m/s",
353     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
354     .                "ave(X)", zstophy,zout)
355c
356         CALL histdef(nid_day, "vitv", "Meridional wind", "m/s",
357     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
358     .                "ave(X)", zstophy,zout)
359c
360         CALL histdef(nid_day, "vitw", "Vertical wind", "Pa/s",
361     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
362     .                "ave(X)", zstophy,zout)
363c
364         CALL histdef(nid_day, "pres", "Air pressure", "Pa",
365     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
366     .                "ave(X)", zstophy,zout)
367c
368c=================================================================
369c   FIN INITIALISATION DES CHAMPS 3D
370c=================================================================
371      ENDIF !  lev_histday.GE.3
372c
373      IF (lev_histday.GE.4) THEN
374c=================================================================
375c
376c   INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
377c
378c=================================================================
379c
380         CALL histdef(nid_day, "SWupTOA", "SWup at TOA","W/m2",
381     .                iim,jj_nb,nhori, 1,1,1,-99,
382     .                32, "ave(X)", zstophy,zout)
383c
384         CALL histdef(nid_day, "SWupSFC", "SWup at surface","W/m2",
385     .                iim,jj_nb,nhori, 1,1,1,-99,
386     .                32, "ave(X)", zstophy,zout)
387c
388         CALL histdef(nid_day, "SWdnTOA", "SWdn at TOA","W/m2",
389     .                iim,jj_nb,nhori, 1,1,1,-99,
390     .                32, "ave(X)", zstophy,zout)
391c
392         CALL histdef(nid_day, "SWupTOAclr",
393     .                "SWup clear sky at TOA","W/m2",
394     .                iim,jj_nb,nhori, 1,1,1,-99,
395     .                32, "ave(X)", zstophy,zout)
396
397         CALL histdef(nid_day, "SWupSFCclr",
398     .                "SWup clear sky at surface","W/m2",
399     .                iim,jj_nb,nhori, 1,1,1,-99,
400     .                32, "ave(X)", zstophy,zout)
401
402         CALL histdef(nid_day, "SWdnTOAclr",
403     .                "SWdn clear sky at TOA","W/m2",
404     .                iim,jj_nb,nhori, 1,1,1,-99,
405     .                32, "ave(X)", zstophy,zout)
406
407         CALL histdef(nid_day, "SWdnSFCclr",
408     .                "SWdn clear sky at surface","W/m2",
409     .                iim,jj_nb,nhori, 1,1,1,-99,
410     .                32, "ave(X)", zstophy,zout)
411c
412cIM AMMA-MIP
413         CALL histdef(nid_day, "LWdnSFC",
414     .                "LW down at surface","W/m2",
415     .                iim,jj_nb,nhori, 1,1,1,-99,
416     .                32, "ave(X)", zstophy,zout)
417c
418         CALL histdef(nid_day, "LWupSFC",
419     .                "LW down at surface","W/m2",
420     .                iim,jj_nb,nhori, 1,1,1,-99,
421     .                32, "ave(X)", zstophy,zout)
422c
423cIM     IF (iflag_con.GE.3) THEN
424          if (ok_ISCCP) then
425            CALL histdef(nid_day, "ptop",
426c    .                "Cld top pressure","mb",
427     .                "Cld top pressure (cf ISCCP simulator)","Pa",
428     .                iim,jj_nb,nhori, 1,1,1,-99,
429     .                32, "ave(X)", zstophy,zout)
430          ENDIF
431cIM     ENDIF !(iflag_con.GE.3) THEN
432c
433c=================================================================
434c=================================================================
435c=================================================================
436c   INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
437c=================================================================
438c
439         CALL histdef(nid_day, "tter", "Surface Temperature ter", "K",
440     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
441     .                "ave(X)", zstophy,zout)
442c
443         CALL histdef(nid_day, "tlic", "Surface Temperature lic", "K",
444     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
445     .                "ave(X)", zstophy,zout)
446c
447         CALL histdef(nid_day, "toce", "Surface Temperature oce", "K",
448     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
449     .                "ave(X)", zstophy,zout)
450c
451         CALL histdef(nid_day, "tsic", "Surface Temperature sic", "K",
452     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
453     .                "ave(X)", zstophy,zout)
454c
455         CALL histdef(nid_day, "t2mter", "Temp.tere 2m", "K",
456     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
457     .                "ave(X)", zstophy,zout)
458c
459         CALL histdef(nid_day, "t2mlic", "Temp.lic 2m", "K",
460     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
461     .                "ave(X)", zstophy,zout)
462c
463         CALL histdef(nid_day, "t2moce", "Temp.oce 2m", "K",
464     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
465     .                "ave(X)", zstophy,zout)
466c
467         CALL histdef(nid_day, "t2msic", "Temp.sic 2m", "K",
468     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
469     .                "ave(X)", zstophy,zout)
470c
471         CALL histdef(nid_day, "t2mter_min", "Temp.terre 2m min.",
472     .                "K",
473     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
474     .                t2mincels, zstophy,zout)
475c
476         CALL histdef(nid_day, "t2mter_max", "Temp.terre 2m max.",
477     .                "K",
478     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
479     .                t2maxcels, zstophy,zout)
480c
481         CALL histdef(nid_day, "u10mter", "Vent zonal ter 10m", "m/s",
482     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
483     .                "ave(X)", zstophy,zout)
484c
485         CALL histdef(nid_day, "u10mlic", "Vent zonal lic 10m", "m/s",
486     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
487     .                "ave(X)", zstophy,zout)
488c
489         CALL histdef(nid_day, "u10moce", "Vent zonal oce 10m", "m/s",
490     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
491     .                "ave(X)", zstophy,zout)
492c
493         CALL histdef(nid_day, "u10msic", "Vent zonal sic 10m",
494     .                 "m/s",iim,jj_nb,nhori, 1,1,1, -99, 32,
495     .                "ave(X)", zstophy,zout)
496c
497         CALL histdef(nid_day, "v10mter", "Vent meridien ter 10m", 
498     .                "m/s", iim,jj_nb,nhori, 1,1,1, -99, 32,
499     .                "ave(X)", zstophy,zout)
500c
501         CALL histdef(nid_day, "v10mlic", "Vent meridien lic 10m",
502     .                 "m/s",iim,jj_nb,nhori, 1,1,1, -99, 32,
503     .                "ave(X)", zstophy,zout)
504c
505         CALL histdef(nid_day, "v10moce", "Vent meridien oce 10m",
506     .                 "m/s",iim,jj_nb,nhori, 1,1,1, -99, 32,
507     .                "ave(X)", zstophy,zout)
508c
509         CALL histdef(nid_day, "v10msic", "Vent meridien sic 10m",
510     .                 "m/s",iim,jj_nb,nhori, 1,1,1, -99, 32,
511     .                "ave(X)", zstophy,zout)
512c
513         DO nsrf = 1, nbsrf
514C
515           call histdef(nid_day, "pourc_"//clnsurf(nsrf),
516     $         "% "//clnsurf(nsrf), "%", 
517     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
518     $         "ave(X)", zstophy,zout)
519C
520           call histdef(nid_day, "fract_"//clnsurf(nsrf),
521     $         "Fraction "//clnsurf(nsrf), "1", 
522     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
523     $         "ave(X)", zstophy,zout)
524C
525           call histdef(nid_day, "tsol_"//clnsurf(nsrf),
526     $         "Temperature "//clnsurf(nsrf), "K", 
527     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
528     $         "ave(X)", zstophy,zout)
529C
530           call histdef(nid_day, "sens_"//clnsurf(nsrf),
531     $         "Sensible heat flux "//clnsurf(nsrf), "W/m2", 
532     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
533     $         "ave(X)", zstophy,zout)
534c
535           call histdef(nid_day, "lat_"//clnsurf(nsrf),
536     $         "Latent heat flux "//clnsurf(nsrf), "W/m2", 
537     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
538     $         "ave(X)", zstophy,zout)
539C
540           call histdef(nid_day, "taux_"//clnsurf(nsrf),
541     $         "Zonal wind stress "//clnsurf(nsrf),"Pa",
542     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
543     $         "ave(X)", zstophy,zout)
544
545           call histdef(nid_day, "tauy_"//clnsurf(nsrf),
546     $         "Meridional wind stress "//clnsurf(nsrf), "Pa", 
547     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
548     $         "ave(X)", zstophy,zout)
549C
550           call histdef(nid_day, "albe_"//clnsurf(nsrf),
551     $         "Albedo surf. "//clnsurf(nsrf), "-", 
552     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
553     $         "ave(X)", zstophy,zout)
554C
555           call histdef(nid_day, "rugs_"//clnsurf(nsrf),
556     $         "Rugosite "//clnsurf(nsrf), "m", 
557     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
558     $         "ave(X)", zstophy,zout)
559C
560
561
562! FH Sorties specifiques pour Mellor et Yamada
563      if (iflag_pbl>1 .and. lev_histday.gt.10 ) then
564           call histdef(nid_day, "tke_"//clnsurf(nsrf),
565     $         "Max Turb. Kinetic Energy "//clnsurf(nsrf), "-", 
566     $         iim,jj_nb,nhori, klev,1,klev,nvert, 32,
567     $         "ave(X)", zstophy,zout)
568
569           call histdef(nid_day, "tke_max_"//clnsurf(nsrf),
570     $         "Max Turb. Kinetic Energy "//clnsurf(nsrf), "-", 
571     $         iim,jj_nb,nhori, klev,1,klev,nvert, 32,
572     $         "t_max(X)", zstophy,zout)
573      endif     
574
575C
576         END DO
577C           
578c=================================================================
579c
580c   FIN INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
581c
582c=================================================================
583      ENDIF !lev_histday.GE.4
584c
585      IF(lev_histday.GE.5) THEN !lev_histday.GE.5
586c
587c rajout sorties F. Aires
588         CALL histdef(nid_day, "tops0", "CS Solar rad. at TOA", 
589     .                "W/m2",
590     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
591     .                "ave(X)", zstophy,zout)
592c
593         CALL histdef(nid_day, "topl0", "CS IR rad. at TOA", 
594     .                "W/m2",
595     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
596     .                "ave(X)", zstophy,zout)
597c
598         CALL histdef(nid_day, "sols0", "CS Net Solar rad. at surf.", 
599     .                "W/m2",
600     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
601     .                "ave(X)", zstophy,zout)
602c
603         CALL histdef(nid_day, "soll0", "CS Net IR rad. at surface", 
604     .                "W/m2",
605     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
606     .                "ave(X)", zstophy,zout)
607c
608         CALL histdef(nid_day, "rneb", "Cloud fraction", "-",
609     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
610     .                "ave(X)", zstophy,zout)
611c
612         CALL histdef(nid_day, "dqpbl", "PBL dT", "K/s",
613     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
614     .                "ave(X)", zstophy,zout)
615c
616         CALL histdef(nid_day, "rnebcon", "Convective Cloud Fraction"
617     .                , "-",
618     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
619     .                "ave(X)", zstophy,zout)
620c
621         CALL histdef(nid_day, "dtwak", "Wake dT", "K/s",
622     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
623     .                "ave(X)", zstophy,zout)
624c
625         CALL histdef(nid_day,"lwcon","Cloud liquid water content",
626     .                "kg/kg",
627     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32,
628     .                "ave(X)", zstophy,zout)
629c
630        CALL histdef(nid_day,"iwcon","Cloud ice water content","kg/kg",
631     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32,
632     .                "ave(X)", zstophy,zout)
633c
634         CALL histdef(nid_day,"lwp","Cloud water path","kg/m2",
635     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
636     .                "ave(X)", zstophy,zout)
637c
638         CALL histdef(nid_day,"iwp","Cloud ice water path","kg/m2",
639     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
640     .                "ave(X)", zstophy,zout)
641c
642
643         CALL histdef(nid_day, "cape", "Conv avlbl pot ener", "J/kg",
644     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
645     .                "ave(X)", zstophy,zout)
646
647         CALL histdef(nid_day, "wh", "Wake height", "m",
648     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
649     .                "ave(X)", zstophy,zout)
650c
651         CALL histdef(nid_day, "ws", "Wake surface", "m2",
652     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
653     .                "ave(X)", zstophy,zout)
654c
655        CALL histdef(nid_day,"meantaucld",
656     .                "ISCCP mean cloud optical thickness","1",
657     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
658     .                "ave(X)", zstophy,zout)
659c
660        CALL histdef(nid_day,"cldtau","Cloud optical thickness","1",
661     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32,
662     .                "ave(X)", zstophy,zout)
663c
664        CALL histdef(nid_day,"cldemi","Cloud optical emissivity","1",
665     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32,
666     .                "ave(X)", zstophy,zout)
667c
668c
669      CALL histdef(nid_day, "cin", "Conv Inhibition", "J/kg",
670     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
671     .                "ave(X)", zstophy,zout)     
672c
673         CALL histdef(nid_day, "wale", "Wake Available Energy", "J/kg",
674     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
675     .                "ave(X)", zstophy,zout)
676c
677         CALL histdef(nid_day, "walp", 
678     .                "Available Lifting Energy due to wake", "W/m2",
679     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
680     .                "ave(X)", zstophy,zout)
681c
682         CALL histdef(nid_day, "blale", "PBL Available Energy", "J/kg",
683     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
684     .                "ave(X)", zstophy,zout)
685c
686         CALL histdef(nid_day, "blalp", 
687     .                "Available Lifting Energy due to PBL", "W/m2",
688     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
689     .                "ave(X)", zstophy,zout)
690c
691         CALL histdef(nid_day, "wdt1", 
692     .                "Temp diff wake layer 1", "K",
693     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
694     .                "ave(X)", zstophy,zout)
695c
696         CALL histdef(nid_day, "wdq1", 
697     .                "Temp diff wake layer 1", "g/kg",
698     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
699     .                "ave(X)", zstophy,zout)
700cIM rajout AMMA-MIP
701c
702         CALL histdef(nid_day, "dtcon", "Convection dT", "K/s",
703     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
704     .                "ave(X)", zstophy,zout)
705c
706         CALL histdef(nid_day, "dqcon", "Convection dQ", "(kg/kg)/s",
707     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
708     .                "ave(X)", zstophy,zout)
709c
710         CALL histdef(nid_day, "rh2m", "Relative humidity at 2m", "%",
711     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
712     .                "ave(X)", zstophy,zout)
713c
714         CALL histdef(nid_day, "qsat2m", "Saturant humidity at 2m", "%",
715     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
716     .                "ave(X)", zstophy,zout)
717c
718         CALL histdef(nid_day, "tpot", 
719     .                "Surface air potential temperature", "K",
720     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
721     .                "ave(X)", zstophy,zout)
722c
723         CALL histdef(nid_day, "tpote", 
724     .        "Surface air equivalent potential temperature", "K",
725     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
726     .                "ave(X)", zstophy,zout)
727c
728      ENDIF !lev_histday.GE.5
729c=================================================================
730c
731         CALL histend(nid_day)
732c
733         ndex2d = 0
734         ndex3d = 0
735c
736c=================================================================
737c$OMP END MASTER
738      ENDIF ! fin de test sur ok_journe
Note: See TracBrowser for help on using the repository browser.