Context Navigation

source: LMDZ4/branches/pre_V3/libf/phylmd/ini_bilKP_ave.h @ 5451

Last change on this file since 5451 was 684, checked in by lmdzadmin, 19 years ago
Ajout nouveaux flags (ip_ebil_phy, ok_slab_sicOBS) ds physiq.def + Lecture freq. ecriture en nombre de jours IM
Property svn:eol-style set to `native` Property svn:executable set to ``* Property svn:keywords set to `Author Date Id Revision`
File size: 8.6 KB

Line
1	c
2	c $Header$
3	c
4	IF (ok_journe) THEN
5	c
6	zsto = dtime
7	zout = ecrit_day
8	typeval=tave
9	c
10	idayref = day_ref
11	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
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	DO ll=1,klev
19	znivsig(ll)=float(ll)
20	ENDDO
21	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
22	write(,)'zx_lon = ',zx_lon(:,1)
23	write(,)'zx_lat = ',zx_lat(1,:)
24	CALL histbeg("histbilKP_ave", iim,zx_lon(:,1), jjmp1,
25	. zx_lat(1,:),
26	. 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
27	. nhori, nid_bilKPave)
28	write(,)'Journee ', itau_phy, zjulian
29	CALL histvert(nid_bilKPave, "presnivs",
30	. "Vertical levels","mb",
31	. klev, presnivs/100., nvert)
32	c
33	c
34	c Champs 3D:
35	c
36	CALL histdef(nid_bilKPave,"ue",
37	. "Zonal energy transport","-",
38	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
39	. typeval, zsto,zout)
40	c
41	CALL histdef(nid_bilKPave,"ve",
42	. "Merid energy transport","-",
43	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
44	. typeval, zsto,zout)
45	c
46	CALL histdef(nid_bilKPave,"uq",
47	. "Zonal humidity transport","-",
48	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
49	. typeval, zsto,zout)
50	c
51	CALL histdef(nid_bilKPave,"vq",
52	. "Merid humidity transport","-",
53	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
54	. typeval, zsto,zout)
55	c
56	c Champs 3D:
57	c
58	CALL histdef(nid_bilKPave,"temp",
59	. "Air temperature","K",
60	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
61	. typeval, zsto,zout)
62	c
63	CALL histdef(nid_bilKPave,"ovap",
64	. "Specific humidity","Kg/Kg",
65	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
66	. typeval, zsto,zout)
67	c
68	CALL histdef(nid_bilKPave,"geop",
69	. "Geopotential height","m",
70	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
71	. typeval, zsto,zout)
72	c
73	CALL histdef(nid_bilKPave,"vitu",
74	. "Zonal wind","m/s",
75	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
76	. typeval, zsto,zout)
77	c
78	CALL histdef(nid_bilKPave,"vitv",
79	. "Meridional wind","m/s",
80	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
81	. typeval, zsto,zout)
82	c
83	CALL histdef(nid_bilKPave, "vitw",
84	. "Vertical wind", "m/s",
85	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
86	. typeval, zsto,zout)
87	c
88	CALL histdef(nid_bilKPave, "pres",
89	. "Inter-Layer Air pressure",
90	. "Pa",
91	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
92	. typeval, zsto,zout)
93	c
94	CALL histdef(nid_bilKPave, "play",
95	. "Mean-Layer Air pressure",
96	. "Pa",
97	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
98	. typeval, zsto,zout)
99	c
100	CALL histdef(nid_bilKPave, "oliq",
101	. "Liquid water content",
102	. "kg/kg",
103	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
104	. typeval, zsto,zout)
105	c
106	CALL histdef(nid_bilKPave, "dtdyn",
107	. "Dynamics dT", "K/s",
108	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
109	. typeval, zsto,zout)
110	c
111	CALL histdef(nid_bilKPave, "dqdyn",
112	. "Dynamics dQ", "Kg/Kg/s",
113	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
114	. typeval, zsto,zout)
115	c
116	CALL histdef(nid_bilKPave, "dtcon",
117	. "Convection dT", "K/s",
118	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
119	. typeval, zsto,zout)
120	c
121	CALL histdef(nid_bilKPave, "ducon",
122	. "Convection du", "m/s2",
123	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
124	. typeval, zsto,zout)
125	c
126	CALL histdef(nid_bilKPave, "dvcon",
127	. "Convection dv", "m/s2",
128	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
129	. typeval, zsto,zout)
130	c
131	CALL histdef(nid_bilKPave,"dqcon",
132	. "Convection dQ","Kg/Kg/s",
133	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
134	. typeval, zsto,zout)
135	c
136	CALL histdef(nid_bilKPave, "dtlsc",
137	. "Condensation dT", "K/s",
138	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
139	. typeval, zsto,zout)
140	c
141	CALL histdef(nid_bilKPave,"dqlsc",
142	. "Condensation dQ","Kg/Kg/s",
143	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
144	. typeval, zsto,zout)
145	c
146	CALL histdef(nid_bilKPave,"dtvdf",
147	. "Boundary-layer dT","K/s",
148	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
149	. typeval, zsto,zout)
150	c
151	CALL histdef(nid_bilKPave, "dqvdf",
152	. "Boundary-layer dQ",
153	. "Kg/Kg/s",
154	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
155	. typeval, zsto,zout)
156	c
157	CALL histdef(nid_bilKPave,"dtajs",
158	. "Ajustement sec dT","K/s",
159	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
160	. typeval, zsto,zout)
161	c
162	CALL histdef(nid_bilKPave, "dqajs",
163	. "Ajustement sec dQ",
164	. "Kg/Kg/s",
165	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
166	. typeval, zsto,zout)
167	c
168	CALL histdef(nid_bilKPave, "dteva",
169	. "Reevaporation dT", "K/s",
170	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
171	. typeval, zsto,zout)
172	c
173	CALL histdef(nid_bilKPave,"dqeva",
174	. "Reevaporation dQ",
175	. "Kg/Kg/s",
176	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
177	. typeval, zsto,zout)
178
179	c
180	CALL histdef(nid_bilKPave, "dtswr",
181	. "SW radiation dT", "K/s",
182	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
183	. typeval, zsto,zout)
184	c
185	CALL histdef(nid_bilKPave, "dtsw0",
186	. "SW radiation dT", "K/s",
187	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
188	. typeval, zsto,zout)
189	c
190	CALL histdef(nid_bilKPave, "dtlwr",
191	. "LW radiation dT", "K/s",
192	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
193	. typeval, zsto,zout)
194	c
195	CALL histdef(nid_bilKPave, "dtlw0",
196	. "LW radiation dT", "K/s",
197	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
198	. typeval, zsto,zout)
199	c
200	CALL histdef(nid_bilKPave,"duvdf",
201	. "Boundary-layer dU","m/s2",
202	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
203	. typeval, zsto,zout)
204	c
205	CALL histdef(nid_bilKPave,"dvvdf",
206	. "Boundary-layer dV","m/s2",
207	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
208	. typeval, zsto,zout)
209	c
210	IF (ok_orodr) THEN
211	IF (ok_orolf) THEN
212	CALL histdef(nid_bilKPave, "duoli",
213	. "Orography dU", "m/s2",
214	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
215	. typeval, zsto,zout)
216	c
217	CALL histdef(nid_bilKPave, "dvoli",
218	. "Orography dV", "m/s2",
219	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
220	. typeval, zsto,zout)
221	c
222	ENDIF
223	ENDIF
224	C
225	CALL histdef(nid_bilKPave, "duphy",
226	. "Physiq dU","-",
227	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
228	. typeval, zsto,zout)
229	C
230	CALL histdef(nid_bilKPave, "dvphy",
231	. "Physiq dV","-",
232	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
233	. typeval, zsto,zout)
234	C
235	CALL histdef(nid_bilKPave, "dtphy",
236	. "Physiq dT","-",
237	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
238	. typeval, zsto,zout)
239	C
240	CALL histdef(nid_bilKPave, "dqphy",
241	. "Physiq dQ","-",
242	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
243	. typeval, zsto,zout)
244	C
245	CALL histdef(nid_bilKPave, "dqlphy",
246	. "Physiq dQl","-",
247	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
248	. typeval, zsto,zout)
249	C
250	C
251	CALL histend(nid_bilKPave)
252	c
253	ndex2d = 0
254	ndex3d = 0
255	c
256	ENDIF ! fin de test sur ok_journe

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

Download in other formats: