Context Navigation

ini_histhf.h @ 2622

Last change on this file since 2622 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
3	c
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
8	c
9	idayref = day_ref
10	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
11
12	c
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
21	cccIM 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)
28	c
29
30	IF(lev_histhf.GE.1) THEN
31	c
32	c CALL histdef(nid_hf, "phis", "Surface geop. height", "-",
33	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
34	c . "once", zsto,zout)
35	c
36	c CALL histdef(nid_hf, "aire", "Grid area", "-",
37	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
38	c . "once", zsto,zout)
39	c
40	CALL histdef(nid_hf, "aireTER","Grid area CONT","-",
41	. iim,jjmp1,nhori, 1,1,1, -99, 32,
42	. "once", zsto,zout)
43	c
44	CALL histdef(nid_hf, "contfracATM","% sfce ter+lic ","-",
45	. iim,jjmp1,nhori, 1,1,1, -99, 32,
46	. "once", zsto,zout)
47	c
48	CALL histdef(nid_hf, "contfracOR","% sfce terre OR", "-",
49	. iim,jjmp1,nhori, 1,1,1, -99, 32,
50	. "once", zsto,zout)
51	c
52	c Champs 2D:
53	c
54	CALL histdef(nid_hf, "t2m", "Temperature 2m", "K",
55	. iim,jjmp1,nhori, 1,1,1, -99, 32,
56	. "inst(X)", zsto,zout)
57	c
58	CALL histdef(nid_hf, "q2m", "Specific humidity", "kg/kg",
59	. iim,jjmp1,nhori, 1,1,1, -99, 32,
60	. "inst(X)", zsto,zout)
61	c
62	CALL histdef(nid_hf, "psol", "Surface Pressure", "Pa",
63	. iim,jjmp1,nhori, 1,1,1, -99, 32,
64	. "inst(X)", zsto,zout)
65	c
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)
71	c
72	cIM ENSEMBLES BEG
73	c
74	CALL histdef(nid_hf, "tsol", "Surface Temperature", "K",
75	. iim,jjmp1,nhori, 1,1,1, -99, 32,
76	. "inst(X)", zsto,zout)
77	c
78	CALL histdef(nid_hf, "slp", "Sea Level Pressure", "Pa",
79	. iim,jjmp1,nhori, 1,1,1, -99, 32,
80	. "inst(X)", zsto,zout)
81	c
82	c
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)
86	c
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)
90	c
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)
94	c
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
103	c
104	ENDIF !lev_histhf.GE.1
105	c
106	IF(lev_histhf.GE.2) THEN
107	c
108	CALL histdef(nid_hf, "cldt", "Total cloudiness", "-",
109	. iim,jjmp1,nhori, 1,1,1, -99, 32,
110	. "ave(X)", zsto1,zout)
111	c
112	c
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)
117	c
118	c
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)
123	c
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
129	c
130	cIM ENSEMBLES END
131	c
132	IF(lev_histhf.GE.3) THEN
133	c
134	DO k=1, nlevENS
135	c
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)
140	c
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
147	c
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)
152	c
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
159	c
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)
164	c
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)
169	c
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)
175	c
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
180	c
181	ENDIF !lev_histhf.GE.3
182	c
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
189	c
190	c#define histhf3d
191	c#ifdef histhf3d
192	c#include "ini_histhf3d.h"
193	c#endif
194	c
195	CALL histend(nid_hf)
196	c
197	endif ! ok_hf

Note: See TracBrowser for help on using the repository browser.

Context Navigation

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

Download in other formats: