source: LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/ini_histhf.h @ 5306

Last change on this file since 5306 was 918, checked in by Laurent Fairhead, 17 years ago
  • correction du bug ISCCP (n'ecrire ptop que quand ISCCP est appelé)
  • petite inversion de boucle dans isccp_cloud_types.F pour aller + vite
  • "CFisation" d'un certain nombre d'unités pour les hist*
  • les suggestions de JL pour rugoro

SD
LF
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 12.8 KB
Line 
1!
2! $Header$
3!
4
5      IF (ok_hf) THEN
6c$OMP MASTER
7c
8       zstohf = ecrit_hf
9       zstophy = dtime
10       zstoday = ecrit_day
11       zout = ecrit_hf
12c
13         idayref = day_ref
14         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
15
16c
17cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
18cym         DO i = 1, iim
19cym            zx_lon(i,1) = rlon(i+1)
20cym            zx_lon(i,jjmp1) = rlon(i+1)
21cym         ENDDO
22
23cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
24
25cccIM      CALL histbeg("histhf", iim,zx_lon, jjmp1,zx_lat,
26cym         CALL histbeg("histhf", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
27cym     .                 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 
28cym     .                 nhori, nid_hf)
29
30         CALL histbeg_phy("histhf", itau_phy, zjulian, dtime, 
31     .                 nhori, nid_hf)
32
33         CALL histvert(nid_hf, "presnivs", "Vertical levels", "mb",
34     .                 klev, presnivs/100., nvert)
35c
36
37      IF(type_run.EQ."CLIM".OR.type_run.EQ."ENSP") THEN
38c
39      IF(lev_histhf.GE.1) THEN
40c
41         CALL histdef(nid_hf, "aireTER","Grid area CONT","m2",
42     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
43     .                "once", zstohf,zout)
44c
45         CALL histdef(nid_hf, "contfracATM","% sfce ter+lic ","1",
46     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
47     .                "once", zstohf,zout)
48c
49         CALL histdef(nid_hf, "contfracOR","% sfce terre OR", "1",
50     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
51     .                "once", zstohf,zout)
52c
53c Champs 2D:
54c
55         CALL histdef(nid_hf, "t2m", "Temperature 2m", "K",
56     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
57     .                "inst(X)", zstohf,zout)
58c
59         CALL histdef(nid_hf, "q2m", "Specific humidity", "kg/kg",
60     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
61     .                "inst(X)", zstohf,zout)
62c
63         CALL histdef(nid_hf, "psol", "Surface Pressure", "Pa",
64     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
65     .                "inst(X)", zstohf,zout)
66c
67         CALL histdef(nid_hf, "precip", "Precipitation", "kg/(s*m2)",
68     .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
69     .                "ave(X)", zstophy,zout)
70c
71cIM ENSEMBLES BEG
72         CALL histdef(nid_hf, "tsol", "Surface Temperature", "K",
73     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
74     .                "inst(X)", zstohf,zout)
75c
76         CALL histdef(nid_hf, "slp", "Sea Level Pressure", "Pa",
77     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
78     .                "inst(X)", zstohf,zout)
79c
80         CALL histdef(nid_hf, "u10m", "Vent zonal 10m", "m/s",
81     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
82     .                "inst(X)", zstohf,zout)
83c
84         CALL histdef(nid_hf, "v10m", "Vent meridien 10m", "m/s",
85     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
86     .                "inst(X)", zstohf,zout)
87c
88         CALL histdef(nid_hf, "wind10m","10-m wind speed","m/s",
89     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
90     .                "inst(X)", zstohf,zout)
91c
92         DO k=1, nlevSTD
93c
94         bb2=clevSTD(k)
95
96         IF(bb2.EQ."500") THEN
97         CALL histdef(nid_hf, "phi"//bb2,
98     .                "Geopotential "//bb2//"mb", "m",
99     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
100     .                "inst(X)", zstohf,zout)
101         ENDIF !bb2.EQ."500"
102         ENDDO !k=1, nlevSTD
103c
104      ENDIF !lev_histhf.GE.1
105c
106      IF(lev_histhf.GE.2) THEN
107c
108         CALL histdef(nid_hf, "cldt", "Total cloudiness", "1",
109     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
110     .                "ave(X)", zstophy,zout)
111c
112         CALL histdef(nid_hf, "SWnetOR",
113     .                "Sfce net SW radiation OR", "W/m2",
114     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
115     .                "ave(X)", zstophy,zout)
116c
117         CALL histdef(nid_hf, "SWdownOR",
118     .                "Sfce incident SW radiation OR", "W/m2",
119     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
120     .                "ave(X)", zstophy,zout)
121c
122         CALL histdef(nid_hf, "LWdownOR",
123     .                "Sfce incident LW radiation OR", "W/m2",
124     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
125     .                "ave(X)", zstophy,zout)
126      ENDIF !lev_histhf.GE.2
127c
128cIM ENSEMBLES END
129c
130      IF(lev_histhf.GE.3) THEN
131c
132         DO k=1, nlevSTD
133         bb2=clevSTD(k)
134         IF(bb2.EQ."1000".OR.bb2.EQ."925".OR.
135     $      bb2.EQ."850".OR.bb2.EQ."700".OR.
136     $      bb2.EQ."500".OR.bb2.EQ."300".OR.
137     $      bb2.EQ."250".OR.bb2.EQ."200".OR.
138     $      bb2.EQ."150".OR.bb2.EQ."100".OR.
139     $      bb2.EQ."50".OR.bb2.EQ."30".OR.
140     $      bb2.EQ."10")THEN
141c
142          CALL histdef(nid_hf, "t"//bb2,
143     .                 "Temperature "//bb2//"mb","K",
144     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
145     .                "inst(X)", zstohf,zout)
146c
147          IF(bb2.NE."500") THEN
148          CALL histdef(nid_hf, "phi"//bb2,
149     .                "Geopotential "//bb2//"mb", "m",
150     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
151     .                "inst(X)", zstohf,zout)
152          ENDIF
153c
154          IF(bb2.GE."200") THEN
155          CALL histdef(nid_hf, "q"//bb2,
156     .                 "Specific humidity "//bb2//"mb","kg/kg",
157     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
158     .                "inst(X)", zstohf,zout)
159          ENDIF
160c
161          CALL histdef(nid_hf, "u"//bb2,
162     .                 "Zonal wind "//bb2//"mb","m/s",
163     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
164     .                "inst(X)", zstohf,zout)
165c
166          CALL histdef(nid_hf, "v"//bb2,
167     .                 "Meridional wind "//bb2//"mb","m/s",
168     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
169     .                "inst(X)", zstohf,zout)
170c
171          CALL histdef(nid_hf, "w"//bb2,
172     .                 "Vertical motion "//bb2//"mb","Pa/s",
173     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
174     .                "inst(X)", zstohf,zstohf)
175         ENDIF
176c
177         ENDDO !nlevSTD
178c
179      ENDIF !lev_histhf.GE.3
180c
181      IF(lev_histhf.GE.4) THEN
182c
183cIM sorties hf 3D
184c
185#undef histhf3d
186#define histhf3d
187#ifdef histhf3d
188#include "ini_histhf3d.h"
189#endif
190      ENDIF !lev_histhf.GE.4
191c
192      ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
193c
194      IF(lev_histhf.GE.1) THEN
195c
196        DO k=1, nlevSTD
197c
198c
199c        bb=clevSTD(k)
200c
201c        IF(k.GE.2) THEN
202         IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
203         IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
204c         aa=clevSTD(k)
205c         bb=aa(1:lnblnk1(aa))
206c        ENDIF
207c
208         IF(bb2.EQ."850".or.bb2.EQ."200") THEN
209c
210          CALL histdef(nid_hf, "u"//bb2,
211     .                 "Zonal wind "//bb2//"mb","m/s",
212     .         iim,jj_nb,nhori, 1,1,1, -99, 32,
213     .         "inst(X)", zstohf,zstohf)
214c
215          CALL histdef(nid_hf, "v"//bb2,
216     .                 "Meridional wind "//bb2//"mb","m/s",
217     .         iim,jj_nb,nhori, 1,1,1, -99, 32,
218     .         "inst(X)", zstohf,zstohf)
219c
220         ENDIF !bb2.EQ."850".or.bb2.EQ."200"
221        ENDDO !k=1, nlevSTD
222c
223        CALL histdef(nid_hf, "topl", 
224     .       "Outgoing longwave radiation", "W/m2",
225     .       iim,jj_nb,nhori, 1,1,1, -99, 32,
226     .       "ave(X)", zstophy,zout)
227c
228        CALL histdef(nid_hf, "precip","Total precipitation rate",
229     .       "kg/(s*m2)",iim,jj_nb,nhori, 1,1,1, -99, 32, 
230     .       "ave(X)", zstophy,zout)
231c
232        CALL histdef(nid_hf, "slp", 
233     .       "Mean sea-level pressure", "Pa",
234     .       iim,jj_nb,nhori, 1,1,1, -99, 32,
235     .       "inst(X)", zstohf,zout)
236c
237cIM 280504 BEG
238      ENDIF !(lev_histhf.GE.1) THEN
239c
240      IF(lev_histhf.GE.2) THEN
241c
242        DO k=1, nlevSTD
243c
244c        bb=clevSTD(k)
245c
246c        IF(k.GE.2) THEN
247         IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
248         IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
249c         aa=clevSTD(k)
250c         bb=aa(1:lnblnk1(aa))
251c        ENDIF
252c
253         IF(bb2.EQ."850".OR.bb2.EQ."500") THEN
254c
255          CALL histdef(nid_hf, "t"//bb2,
256     .                 "Temperature "//bb2//"mb","K",
257     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
258     .                "inst(X)", zstohf,zstohf)
259c
260          CALL histdef(nid_hf, "q"//bb2,
261     .                 "Specific humidity "//bb2//"mb","kg/kg",
262     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
263     .                "inst(X)", zstohf,zstohf)
264c
265         ENDIF !bb2
266c
267         IF(bb2.EQ."500") THEN
268c
269          CALL histdef(nid_hf, "phi"//bb2,
270     .                "Geopotential "//bb2//"mb", "m",
271     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
272     .                "inst(X)", zstohf,zstohf)
273c
274          CALL histdef(nid_hf, "w"//bb2,
275     .                 "Vertical motion "//bb2//"mb","Pa/s",
276     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
277     .                "inst(X)", zstohf,zstohf)
278         ENDIF !bb2.EQ."500"
279c
280         IF(bb3.EQ."50") THEN
281c
282          CALL histdef(nid_hf, "t"//bb3,
283     .                 "Temperature "//bb3//"mb","K",
284     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
285     .                "inst(X)", zstohf,zstohf)
286c
287          CALL histdef(nid_hf, "u"//bb3,
288     .                 "Zonal wind "//bb3//"mb","m/s",
289     .         iim,jj_nb,nhori, 1,1,1, -99, 32,
290     .         "inst(X)", zstohf,zstohf)
291c
292          CALL histdef(nid_hf, "v"//bb3,
293     .                 "Meridional wind "//bb3//"mb","m/s",
294     .         iim,jj_nb,nhori, 1,1,1, -99, 32,
295     .         "inst(X)", zstohf,zstohf)
296c
297         ENDIF !bb3
298c
299        ENDDO !k=1, nlevSTD
300c
301        CALL histdef(nid_hf, "u10m", "Vent zonal 10m", "m/s",
302     .               iim,jj_nb,nhori, 1,1,1, -99, 32,
303     .               "inst(X)", zstohf,zout)
304c
305        CALL histdef(nid_hf, "v10m", "Vent meridien 10m", "m/s",
306     .               iim,jj_nb,nhori, 1,1,1, -99, 32,
307     .               "inst(X)", zstohf,zout)
308c
309        CALL histdef(nid_hf, "t2m", "Temperature 2m", "K",
310     .               iim,jj_nb,nhori, 1,1,1, -99, 32,
311     .               "inst(X)", zstohf,zout)
312c
313        CALL histdef(nid_hf, "q2m", "Specific humidity", "kg/kg",
314     .               iim,jj_nb,nhori, 1,1,1, -99, 32,
315     .               "inst(X)", zstohf,zout)
316c
317        CALL histdef(nid_hf, "psol", "Surface Pressure", "Pa",
318     .               iim,jj_nb,nhori, 1,1,1, -99, 32, 
319     .               "inst(X)", zstohf,zout)
320c
321        CALL histdef(nid_hf, "cldt", "Total cloudiness", "1",
322     .               iim,jj_nb,nhori, 1,1,1, -99, 32,
323     .               "ave(X)", zstohf,zout)
324c
325        CALL histdef(nid_hf, "taux",
326     $         "Zonal wind stress","Pa",
327     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
328     $         "ave(X)", zstohf,zout)
329c
330        CALL histdef(nid_hf, "tauy",
331     $         "Meridional wind stress ", "Pa",
332     $         iim,jj_nb,nhori, 1,1,1, -99, 32,
333     $         "ave(X)", zstohf,zout)
334c
335        CALL histdef(nid_hf, "prw", "Precipitable water", "kg/m2",
336     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
337cIM 090904  .                "ave(X)", zstohf,zout)
338     .                "ave(X)", zstophy,zout)
339c
340        CALL histdef(nid_hf, "sens", "Sensible heat flux", "W/m2",
341     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
342cIM 090904  .                "ave(X)", zstohf,zout)
343     .                "ave(X)", zstophy,zout)
344c
345        CALL histdef(nid_hf, "flat", "Latent heat flux", "W/m2",
346     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
347     .                "ave(X)", zstohf,zout)
348c
349c       CALL histdef(nid_hf, "snow_mass", "Snow Mass", "kg/m2",
350c    .                iim,jj_nb,nhori, 1,1,1, -99, 32,
351c    .                "ave(X)", zstohf,zout)
352c
353        CALL histdef(nid_hf, "SWdnSFC", "SWdn at surface","W/m2",
354     .                iim,jj_nb,nhori, 1,1,1,-99,
355     .                32, "ave(X)", zstohf,zout)
356c
357        CALL histdef(nid_hf, "SWupSFC", "SWup at surface","W/m2",
358     .                iim,jj_nb,nhori, 1,1,1,-99,
359     .                32, "ave(X)", zstohf,zout)
360c
361        CALL histdef(nid_hf,"LWdnSFC","Down. IR rad. at surface",
362     .                "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32,
363     .                "ave(X)", zstohf,zout)
364c
365        CALL histdef(nid_hf,"LWupSFC","Upwd. IR rad. at surface",
366     .                "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32,
367     .                "ave(X)", zstohf,zout)
368c
369        CALL histdef(nid_hf, "SWdnTOA", "SWdn at TOA","W/m2",
370     .                iim,jj_nb,nhori, 1,1,1,-99,
371     .                32, "ave(X)", zstohf,zout)
372c
373        CALL histdef(nid_hf, "SWupTOA", "SWup at TOA","W/m2",
374     .                iim,jj_nb,nhori, 1,1,1,-99,
375     .                32, "ave(X)", zstohf,zout)
376c
377        DO k=1, nbteta
378         CALL histdef(nid_hf, "PV"//ctetaSTD(k),
379     $               "Vorticite potentielle "//ctetaSTD(k)//"K",
380     $               "1/(Pa*s)",iim,jj_nb,nhori, 1,1,1, -99,
381     $                 32, "inst(X)", zstohf,zout) 
382        ENDDO !k=1, nbteta
383c
384c planetary boundary layer height
385c
386c       ENDIF !1.EQ.0
387c
388      ENDIF !(lev_histhf.GE.2) THEN
389cIM 280504 END
390c
391      ENDIF !type_run
392c
393      CALL histend(nid_hf)
394c
395c$OMP END MASTER
396      ENDIF ! ok_hf
Note: See TracBrowser for help on using the repository browser.