Context Navigation

source: LMDZ4/trunk/libf/phylmd/ini_histhf.h @ 642

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

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

Download in other formats: