Context Navigation

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

Last change on this file since 556 was 511, checked in by lmdzadmin, 21 years ago
Enleve region "1d" des fichiers ini*.h IM
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 6.3 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, "SWdownOR",
114	. "Sfce incident SW radiation OR", "W/m^2",
115	. iim,jjmp1,nhori, 1,1,1, -99, 32,
116	. "ave(X)", zsto1,zout)
117	c
118	CALL histdef(nid_hf, "LWdownOR",
119	. "Sfce incident LW radiation OR", "W/m^2",
120	. iim,jjmp1,nhori, 1,1,1, -99, 32,
121	. "ave(X)", zsto1,zout)
122	ENDIF !lev_histhf.GE.2
123	c
124	cIM ENSEMBLES END
125	c
126	IF(lev_histhf.GE.3) THEN
127	c
128	DO k=1, nlevENS
129	c
130	CALL histdef(nid_hf, "t"//clev(k),
131	. "Temperature"//clev(k)//"mb","K",
132	. iim,jjmp1,nhori, 1,1,1, -99, 32,
133	. "inst(X)", zsto,zout)
134	c
135	IF(clev(k).NE."500") THEN
136	CALL histdef(nid_hf, "phi"//clev(k),
137	. "Geopotential"//clev(k)//"mb", "m2/s2",
138	. iim,jjmp1,nhori, 1,1,1, -99, 32,
139	. "inst(X)", zsto,zout)
140	ENDIF
141	c
142	CALL histdef(nid_hf, "q"//clev(k),
143	. "Specific humidity"//clev(k)//"mb","kg/kg",
144	. iim,jjmp1,nhori, 1,1,1, -99, 32,
145	. "inst(X)", zsto,zout)
146	c
147	IF(1.EQ.0) THEN
148	CALL histdef(nid_hf, "rh"//clev(k),
149	. "Relative humidity"//clev(k)//"mb", "%",
150	. iim,jjmp1,nhori, 1,1,1, -99, 32,
151	. "inst(X)", zsto,zout)
152	ENDIF
153	c
154	CALL histdef(nid_hf, "u"//clev(k),
155	. "Zonal wind"//clev(k)//"mb","m/s",
156	. iim,jjmp1,nhori, 1,1,1, -99, 32,
157	. "inst(X)", zsto,zout)
158	c
159	CALL histdef(nid_hf, "v"//clev(k),
160	. "Meridional wind"//clev(k)//"mb","m/s",
161	. iim,jjmp1,nhori, 1,1,1, -99, 32,
162	. "inst(X)", zsto,zout)
163	c
164	ENDDO !nlevENS
165	IF(1.EQ.0) THEN
166	CALL histdef(nid_hf, "cdrm", " Momentum drag coef.", "-",
167	. iim,jjmp1,nhori, 1,1,1, -99, 32,
168	. "inst(X)", zsto,zout)
169	c
170	CALL histdef(nid_hf, "cdrh", "Heat drag coef.", "-",
171	. iim,jjmp1,nhori, 1,1,1, -99, 32,
172	. "inst(X)", zsto,zout)
173	ENDIF !(1.EQ.0) THEN
174	c
175	ENDIF !lev_histhf.GE.3
176	c
177	IF(lev_histhf.GE.4) THEN
178	#define histhf3d
179	#ifdef histhf3d
180	#include "ini_histhf3d.h"
181	#endif
182	ENDIF !lev_histhf.GE.4
183	c
184	c#define histhf3d
185	c#ifdef histhf3d
186	c#include "ini_histhf3d.h"
187	c#endif
188	c
189	CALL histend(nid_hf)
190	c
191	endif ! ok_hf

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

Download in other formats: