source: LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histhf.h @ 1546

Last change on this file since 1546 was 598, checked in by Laurent Fairhead, 20 years ago

SWnetOR non defini mais utilise dans write_histhf.h AC
LF

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 6.5 KB
Line 
1
2      IF (ok_hf) THEN
3c
4        zout = dtime * REAL(NINT(86400./dtime*ecrit_hf))
5        zsto = zout
6        zsto1 = dtime
7        PRINT*, 'La frequence de sortie instant. est de ', ecrit_hf
8c
9         idayref = day_ref
10         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
11
12c
13         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
14         DO i = 1, iim
15            zx_lon(i,1) = rlon(i+1)
16            zx_lon(i,jjmp1) = rlon(i+1)
17         ENDDO
18
19         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
20
21cccIM      CALL histbeg("histhf", iim,zx_lon, jjmp1,zx_lat,
22         CALL histbeg("histhf", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
23     .                 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 
24     .                 nhori, nid_hf)
25
26         CALL histvert(nid_hf, "presnivs", "Vertical levels", "mb",
27     .                 klev, presnivs, nvert)
28c
29
30      IF(lev_histhf.GE.1) THEN
31c
32c        CALL histdef(nid_hf, "phis", "Surface geop. height", "-",
33c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
34c    .                "once", zsto,zout)
35c
36c        CALL histdef(nid_hf, "aire", "Grid area", "-",
37c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
38c    .                "once", zsto,zout)
39c
40         CALL histdef(nid_hf, "aireTER","Grid area CONT","-",
41     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
42     .                "once", zsto,zout)
43c
44         CALL histdef(nid_hf, "contfracATM","% sfce ter+lic ","-",
45     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
46     .                "once", zsto,zout)
47c
48         CALL histdef(nid_hf, "contfracOR","% sfce terre OR", "-",
49     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
50     .                "once", zsto,zout)
51c
52c Champs 2D:
53c
54         CALL histdef(nid_hf, "t2m", "Temperature 2m", "K",
55     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
56     .                "inst(X)", zsto,zout)
57c
58         CALL histdef(nid_hf, "q2m", "Specific humidity", "kg/kg",
59     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
60     .                "inst(X)", zsto,zout)
61c
62         CALL histdef(nid_hf, "psol", "Surface Pressure", "Pa",
63     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
64     .                "inst(X)", zsto,zout)
65c
66         print*,'ATTENTION METTRE AVE(X) POUR LES PRECIPS'
67
68         CALL histdef(nid_hf, "rain", "Precipitation", "kg/m^2s",
69     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
70     .                "ave(X)", zsto1,zout)
71c
72cIM ENSEMBLES BEG
73c
74         CALL histdef(nid_hf, "tsol", "Surface Temperature", "K",
75     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
76     .                "inst(X)", zsto,zout)
77c
78         CALL histdef(nid_hf, "slp", "Sea Level Pressure", "Pa",
79     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
80     .                "inst(X)", zsto,zout)
81c
82c
83         CALL histdef(nid_hf, "u10m", "Vent zonal 10m", "m/s",
84     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
85     .                "inst(X)", zsto,zout)
86c
87         CALL histdef(nid_hf, "v10m", "Vent meridien 10m", "m/s",
88     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
89     .                "inst(X)", zsto,zout)
90c
91         CALL histdef(nid_hf, "wind10m","10-m wind speed","m/s",
92     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
93     .                "inst(X)", zsto,zout)
94c
95         DO k=1, nlevENS
96          IF(clev(k).EQ."500") THEN
97          CALL histdef(nid_hf, "phi"//clev(k),
98     .                "Geopotential"//clev(k)//"mb", "m2/s2",
99     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
100     .                "inst(X)", zsto,zout)
101          ENDIF !clev(k).EQ."500"
102         ENDDO !k=1, nlevENS
103c
104      ENDIF !lev_histhf.GE.1
105c
106      IF(lev_histhf.GE.2) THEN
107c
108         CALL histdef(nid_hf, "cldt", "Total cloudiness", "-",
109     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
110     .                "ave(X)", zsto1,zout)
111c
112c
113         CALL histdef(nid_hf, "SWnetOR",
114     .                "Sfce net SW radiation OR", "W/m2",
115     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
116     .                "ave(X)", zsto1,zout)
117c
118c
119         CALL histdef(nid_hf, "SWdownOR",
120     .                "Sfce incident SW radiation OR", "W/m^2",
121     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
122     .                "ave(X)", zsto1,zout)
123c
124         CALL histdef(nid_hf, "LWdownOR",
125     .                "Sfce incident LW radiation OR", "W/m^2",
126     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
127     .                "ave(X)", zsto1,zout)
128      ENDIF !lev_histhf.GE.2
129c
130cIM ENSEMBLES END
131c
132      IF(lev_histhf.GE.3) THEN
133c
134         DO k=1, nlevENS
135c
136          CALL histdef(nid_hf, "t"//clev(k),
137     .                 "Temperature"//clev(k)//"mb","K",
138     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
139     .                "inst(X)", zsto,zout)
140c
141          IF(clev(k).NE."500") THEN
142          CALL histdef(nid_hf, "phi"//clev(k),
143     .                "Geopotential"//clev(k)//"mb", "m2/s2",
144     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
145     .                "inst(X)", zsto,zout)
146          ENDIF
147c
148          CALL histdef(nid_hf, "q"//clev(k),
149     .                 "Specific humidity"//clev(k)//"mb","kg/kg",
150     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
151     .                "inst(X)", zsto,zout)
152c
153         IF(1.EQ.0) THEN
154         CALL histdef(nid_hf, "rh"//clev(k),
155     .                 "Relative humidity"//clev(k)//"mb", "%",
156     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
157     .                "inst(X)", zsto,zout)
158         ENDIF
159c
160          CALL histdef(nid_hf, "u"//clev(k),
161     .                 "Zonal wind"//clev(k)//"mb","m/s",
162     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
163     .                "inst(X)", zsto,zout)
164c
165          CALL histdef(nid_hf, "v"//clev(k),
166     .                 "Meridional wind"//clev(k)//"mb","m/s",
167     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
168     .                "inst(X)", zsto,zout)
169c
170         ENDDO !nlevENS
171         IF(1.EQ.0) THEN
172         CALL histdef(nid_hf, "cdrm", " Momentum drag coef.", "-",
173     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
174     .                "inst(X)", zsto,zout)
175c
176         CALL histdef(nid_hf, "cdrh", "Heat drag coef.", "-",
177     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
178     .                "inst(X)", zsto,zout)
179         ENDIF !(1.EQ.0) THEN
180c
181      ENDIF !lev_histhf.GE.3
182c
183      IF(lev_histhf.GE.4) THEN
184#define histhf3d
185#ifdef histhf3d
186#include "ini_histhf3d.h"
187#endif
188      ENDIF !lev_histhf.GE.4
189c
190c#define histhf3d
191c#ifdef histhf3d
192c#include "ini_histhf3d.h"
193c#endif
194c
195         CALL histend(nid_hf)
196c
197      endif ! ok_hf
Note: See TracBrowser for help on using the repository browser.