source: LMDZ4/branches/V3_test/libf/phylmd/ini_histday.h @ 715

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

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