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

Last change on this file since 889 was 879, checked in by Laurent Fairhead, 18 years ago

Suite de la bascule vers une physique avec thermiques, nouvelle convection, poche froide ...
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         CALL histdef(nid_day, "ptop",
425c    .                "Cld top pressure","mb",
426     .                "Cld top pressure (cf ISCCP simulator)","Pa",
427     .                iim,jj_nb,nhori, 1,1,1,-99,
428     .                32, "ave(X)", zstophy,zout)
429cIM     ENDIF !(iflag_con.GE.3) THEN
430c
431c=================================================================
432c=================================================================
433c=================================================================
434c   INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
435c=================================================================
436c
437         CALL histdef(nid_day, "tter", "Surface Temperature ter", "K",
438     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
439     .                "ave(X)", zstophy,zout)
440c
441         CALL histdef(nid_day, "tlic", "Surface Temperature lic", "K",
442     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
443     .                "ave(X)", zstophy,zout)
444c
445         CALL histdef(nid_day, "toce", "Surface Temperature oce", "K",
446     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
447     .                "ave(X)", zstophy,zout)
448c
449         CALL histdef(nid_day, "tsic", "Surface Temperature sic", "K",
450     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
451     .                "ave(X)", zstophy,zout)
452c
453         CALL histdef(nid_day, "t2mter", "Temp.tere 2m", "K",
454     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
455     .                "ave(X)", zstophy,zout)
456c
457         CALL histdef(nid_day, "t2mlic", "Temp.lic 2m", "K",
458     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
459     .                "ave(X)", zstophy,zout)
460c
461         CALL histdef(nid_day, "t2moce", "Temp.oce 2m", "K",
462     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
463     .                "ave(X)", zstophy,zout)
464c
465         CALL histdef(nid_day, "t2msic", "Temp.sic 2m", "K",
466     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
467     .                "ave(X)", zstophy,zout)
468c
469         CALL histdef(nid_day, "t2mter_min", "Temp.terre 2m min.",
470     .                "K",
471     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
472     .                t2mincels, zstophy,zout)
473c
474         CALL histdef(nid_day, "t2mter_max", "Temp.terre 2m max.",
475     .                "K",
476     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
477     .                t2maxcels, zstophy,zout)
478c
479         CALL histdef(nid_day, "u10mter", "Vent zonal ter 10m", "m/s",
480     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
481     .                "ave(X)", zstophy,zout)
482c
483         CALL histdef(nid_day, "u10mlic", "Vent zonal lic 10m", "m/s",
484     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
485     .                "ave(X)", zstophy,zout)
486c
487         CALL histdef(nid_day, "u10moce", "Vent zonal oce 10m", "m/s",
488     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
489     .                "ave(X)", zstophy,zout)
490c
491         CALL histdef(nid_day, "u10msic", "Vent zonal sic 10m",
492     .                 "m/s",iim,jj_nb,nhori, 1,1,1, -99, 32,
493     .                "ave(X)", zstophy,zout)
494c
495         CALL histdef(nid_day, "v10mter", "Vent meridien ter 10m", 
496     .                "m/s", iim,jj_nb,nhori, 1,1,1, -99, 32,
497     .                "ave(X)", zstophy,zout)
498c
499         CALL histdef(nid_day, "v10mlic", "Vent meridien lic 10m",
500     .                 "m/s",iim,jj_nb,nhori, 1,1,1, -99, 32,
501     .                "ave(X)", zstophy,zout)
502c
503         CALL histdef(nid_day, "v10moce", "Vent meridien oce 10m",
504     .                 "m/s",iim,jj_nb,nhori, 1,1,1, -99, 32,
505     .                "ave(X)", zstophy,zout)
506c
507         CALL histdef(nid_day, "v10msic", "Vent meridien sic 10m",
508     .                 "m/s",iim,jj_nb,nhori, 1,1,1, -99, 32,
509     .                "ave(X)", zstophy,zout)
510c
511         DO nsrf = 1, nbsrf
512C
513           call histdef(nid_day, "pourc_"//clnsurf(nsrf),
514     $         "% "//clnsurf(nsrf), "%", 
515     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
516     $         "ave(X)", zstophy,zout)
517C
518           call histdef(nid_day, "fract_"//clnsurf(nsrf),
519     $         "Fraction "//clnsurf(nsrf), "1", 
520     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
521     $         "ave(X)", zstophy,zout)
522C
523           call histdef(nid_day, "tsol_"//clnsurf(nsrf),
524     $         "Temperature "//clnsurf(nsrf), "K", 
525     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
526     $         "ave(X)", zstophy,zout)
527C
528           call histdef(nid_day, "sens_"//clnsurf(nsrf),
529     $         "Sensible heat flux "//clnsurf(nsrf), "W/m2", 
530     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
531     $         "ave(X)", zstophy,zout)
532c
533           call histdef(nid_day, "lat_"//clnsurf(nsrf),
534     $         "Latent heat flux "//clnsurf(nsrf), "W/m2", 
535     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
536     $         "ave(X)", zstophy,zout)
537C
538           call histdef(nid_day, "taux_"//clnsurf(nsrf),
539     $         "Zonal wind stress "//clnsurf(nsrf),"Pa",
540     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
541     $         "ave(X)", zstophy,zout)
542
543           call histdef(nid_day, "tauy_"//clnsurf(nsrf),
544     $         "Meridional wind stress "//clnsurf(nsrf), "Pa", 
545     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
546     $         "ave(X)", zstophy,zout)
547C
548           call histdef(nid_day, "albe_"//clnsurf(nsrf),
549     $         "Albedo surf. "//clnsurf(nsrf), "-", 
550     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
551     $         "ave(X)", zstophy,zout)
552C
553           call histdef(nid_day, "rugs_"//clnsurf(nsrf),
554     $         "Rugosite "//clnsurf(nsrf), "m", 
555     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
556     $         "ave(X)", zstophy,zout)
557C
558
559
560! FH Sorties specifiques pour Mellor et Yamada
561      if (iflag_pbl>1 .and. lev_histday.gt.10 ) then
562           call histdef(nid_day, "tke_"//clnsurf(nsrf),
563     $         "Max Turb. Kinetic Energy "//clnsurf(nsrf), "-", 
564     $         iim,jj_nb,nhori, klev,1,klev,nvert, 32,
565     $         "ave(X)", zstophy,zout)
566
567           call histdef(nid_day, "tke_max_"//clnsurf(nsrf),
568     $         "Max Turb. Kinetic Energy "//clnsurf(nsrf), "-", 
569     $         iim,jj_nb,nhori, klev,1,klev,nvert, 32,
570     $         "t_max(X)", zstophy,zout)
571      endif     
572
573C
574         END DO
575C           
576c=================================================================
577c
578c   FIN INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
579c
580c=================================================================
581      ENDIF !lev_histday.GE.4
582c
583      IF(lev_histday.GE.5) THEN !lev_histday.GE.5
584c
585c rajout sorties F. Aires
586         CALL histdef(nid_day, "tops0", "CS Solar rad. at TOA", 
587     .                "W/m2",
588     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
589     .                "ave(X)", zstophy,zout)
590c
591         CALL histdef(nid_day, "topl0", "CS IR rad. at TOA", 
592     .                "W/m2",
593     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
594     .                "ave(X)", zstophy,zout)
595c
596         CALL histdef(nid_day, "sols0", "CS Net Solar rad. at surf.", 
597     .                "W/m2",
598     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
599     .                "ave(X)", zstophy,zout)
600c
601         CALL histdef(nid_day, "soll0", "CS Net IR rad. at surface", 
602     .                "W/m2",
603     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
604     .                "ave(X)", zstophy,zout)
605c
606         CALL histdef(nid_day, "rneb", "Cloud fraction", "-",
607     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
608     .                "ave(X)", zstophy,zout)
609c
610         CALL histdef(nid_day, "dqpbl", "PBL dT", "K/s",
611     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
612     .                "ave(X)", zstophy,zout)
613c
614         CALL histdef(nid_day, "rnebcon", "Convective Cloud Fraction"
615     .                , "-",
616     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
617     .                "ave(X)", zstophy,zout)
618c
619         CALL histdef(nid_day, "dtwak", "Wake dT", "K/s",
620     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
621     .                "ave(X)", zstophy,zout)
622c
623         CALL histdef(nid_day,"lwcon","Cloud liquid water content",
624     .                "kg/kg",
625     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32,
626     .                "ave(X)", zstophy,zout)
627c
628        CALL histdef(nid_day,"iwcon","Cloud ice water content","kg/kg",
629     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32,
630     .                "ave(X)", zstophy,zout)
631c
632         CALL histdef(nid_day,"lwp","Cloud water path","kg/m2",
633     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
634     .                "ave(X)", zstophy,zout)
635c
636         CALL histdef(nid_day,"iwp","Cloud ice water path","kg/m2",
637     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
638     .                "ave(X)", zstophy,zout)
639c
640
641         CALL histdef(nid_day, "cape", "Conv avlbl pot ener", "J/kg",
642     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
643     .                "ave(X)", zstophy,zout)
644
645         CALL histdef(nid_day, "wh", "Wake height", "m",
646     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
647     .                "ave(X)", zstophy,zout)
648c
649         CALL histdef(nid_day, "ws", "Wake surface", "m2",
650     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
651     .                "ave(X)", zstophy,zout)
652c
653        CALL histdef(nid_day,"meantaucld",
654     .                "ISCCP mean cloud optical thickness","1",
655     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
656     .                "ave(X)", zstophy,zout)
657c
658        CALL histdef(nid_day,"cldtau","Cloud optical thickness","1",
659     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32,
660     .                "ave(X)", zstophy,zout)
661c
662        CALL histdef(nid_day,"cldemi","Cloud optical emissivity","1",
663     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32,
664     .                "ave(X)", zstophy,zout)
665c
666c
667      CALL histdef(nid_day, "cin", "Conv Inhibition", "J/kg",
668     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
669     .                "ave(X)", zstophy,zout)     
670c
671         CALL histdef(nid_day, "wale", "Wake Available Energy", "J/kg",
672     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
673     .                "ave(X)", zstophy,zout)
674c
675         CALL histdef(nid_day, "walp", 
676     .                "Available Lifting Energy due to wake", "W/m2",
677     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
678     .                "ave(X)", zstophy,zout)
679c
680         CALL histdef(nid_day, "blale", "PBL Available Energy", "J/kg",
681     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
682     .                "ave(X)", zstophy,zout)
683c
684         CALL histdef(nid_day, "blalp", 
685     .                "Available Lifting Energy due to PBL", "W/m2",
686     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
687     .                "ave(X)", zstophy,zout)
688c
689         CALL histdef(nid_day, "wdt1", 
690     .                "Temp diff wake layer 1", "K",
691     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
692     .                "ave(X)", zstophy,zout)
693c
694         CALL histdef(nid_day, "wdq1", 
695     .                "Temp diff wake layer 1", "g/kg",
696     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
697     .                "ave(X)", zstophy,zout)
698cIM rajout AMMA-MIP
699c
700         CALL histdef(nid_day, "dtcon", "Convection dT", "K/s",
701     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
702     .                "ave(X)", zstophy,zout)
703c
704         CALL histdef(nid_day, "dqcon", "Convection dQ", "(kg/kg)/s",
705     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
706     .                "ave(X)", zstophy,zout)
707c
708         CALL histdef(nid_day, "rh2m", "Relative humidity at 2m", "%",
709     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
710     .                "ave(X)", zstophy,zout)
711c
712         CALL histdef(nid_day, "qsat2m", "Saturant humidity at 2m", "%",
713     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
714     .                "ave(X)", zstophy,zout)
715c
716         CALL histdef(nid_day, "tpot", 
717     .                "Surface air potential temperature", "K",
718     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
719     .                "ave(X)", zstophy,zout)
720c
721         CALL histdef(nid_day, "tpote", 
722     .        "Surface air equivalent potential temperature", "K",
723     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
724     .                "ave(X)", zstophy,zout)
725c
726      ENDIF !lev_histday.GE.5
727c=================================================================
728c
729         CALL histend(nid_day)
730c
731         ndex2d = 0
732         ndex3d = 0
733c
734c=================================================================
735c$OMP END MASTER
736      ENDIF ! fin de test sur ok_journe
Note: See TracBrowser for help on using the repository browser.