Context Navigation

source: LMDZ4/trunk/libf/phylmd/ini_bilKP_ins.h @ 758

Last change on this file since 758 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: 10.4 KB

Line
1	c
2	c $Header$
3	c
4	IF (ok_journe) THEN
5	c
6	zsto = dtime
7	zout = dtime
8	typeval=tinst
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	c
25	cIM 280405 BEG
26	c
27	cIM cf. AM 081204 BEG region
28	imin_ins=1
29	imax_ins=iim
30	jmin_ins=1
31	jmax_ins=jjmp1
32	do i=1,iim-1
33	if(zx_lon(i,1).lt.lonmin_ins) imin_ins=i
34	if(zx_lon(i,1).le.lonmax_ins) imax_ins=i+1
35	enddo
36	do j=1,jjmp1
37	if(zx_lat(1,j).ge.latmin_ins) jmax_ins=j
38	if(zx_lat(1,j).gt.latmax_ins) jmin_ins=j
39	enddo
40	c
41	print*,'On stoke le fichier bilKP instantanne sur ',
42	s imin_ins,imax_ins,jmin_ins,jmax_ins
43	print*,'On stoke le fichier bilKP instantanne sur ',
44	s zx_lon(imin_ins,1),zx_lon(imax_ins,1),
45	s zx_lat(1,jmin_ins),zx_lat(1,jmax_ins)
46	cIM cf. AM 081204 END region
47	c
48	cIM 280405 END
49	c
50	IF(1.EQ.0) THEN
51	CALL histbeg("histbilKP_ins", iim,zx_lon(:,1), jjmp1,
52	. zx_lat(1,:),
53	. 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
54	. nhori, nid_bilKPins)
55	ENDIF
56	c
57	cIM 280405 BEG
58	c
59	cIM cf. AM 081204 BEG region
60	CALL histbeg("histbilKP_ins", iim,zx_lon(:,1),
61	. jjmp1,zx_lat(1,:),
62	. imin_ins,imax_ins-imin_ins+1,
63	. jmin_ins,jmax_ins-jmin_ins+1,
64	. itau_phy, zjulian, dtime,
65	. nhori, nid_bilKPins)
66	cIM 081204 END
67	c
68	cIM 280405 END
69	c
70	write(,)'Journee ', itau_phy, zjulian
71	CALL histvert(nid_bilKPins, "presnivs",
72	. "Vertical levels","mb",
73	. klev, presnivs/100., nvert)
74	c
75	c Champs 3D:
76	c
77	CALL histdef(nid_bilKPins,"ue",
78	. "Zonal energy transport","-",
79	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
80	. typeval, zsto,zout)
81	c
82	CALL histdef(nid_bilKPins,"ve",
83	. "Merid energy transport","-",
84	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
85	. typeval, zsto,zout)
86	c
87	CALL histdef(nid_bilKPins,"uq",
88	. "Zonal humidity transport","-",
89	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
90	. typeval, zsto,zout)
91	c
92	CALL histdef(nid_bilKPins,"vq",
93	. "Merid humidity transport","-",
94	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
95	. typeval, zsto,zout)
96	c
97	c Champs 3D:
98	c
99	CALL histdef(nid_bilKPins, "temp",
100	. "Air temperature", "K",
101	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
102	. typeval, zsto,zout)
103	c
104	CALL histdef(nid_bilKPins,"ovap",
105	. "Specific humidity","Kg/Kg",
106	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
107	. typeval, zsto,zout)
108	c
109	CALL histdef(nid_bilKPins,"geop",
110	. "Geopotential height", "m",
111	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
112	. typeval, zsto,zout)
113	c
114	CALL histdef(nid_bilKPins,"vitu",
115	. "Zonal wind", "m/s",
116	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
117	. typeval, zsto,zout)
118	c
119	CALL histdef(nid_bilKPins,"vitv",
120	. "Meridional wind", "m/s",
121	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
122	. typeval, zsto,zout)
123	c
124	CALL histdef(nid_bilKPins, "vitw",
125	. "Vertical wind", "m/s",
126	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
127	. typeval, zsto,zout)
128	c
129	CALL histdef(nid_bilKPins, "pres",
130	. "Inter-Layer Air pressure",
131	. "Pa",
132	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
133	. typeval, zsto,zout)
134	c
135	CALL histdef(nid_bilKPins, "play",
136	. "Mean-Layer Air pressure",
137	. "Pa",
138	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
139	. typeval, zsto,zout)
140	c
141	CALL histdef(nid_bilKPins, "oliq",
142	. "Liquid water content",
143	. "kg/kg",
144	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
145	. typeval, zsto,zout)
146	c
147	CALL histdef(nid_bilKPins, "dtdyn",
148	. "Dynamics dT", "K/s",
149	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
150	. typeval, zsto,zout)
151	c
152	CALL histdef(nid_bilKPins, "dqdyn",
153	. "Dynamics dQ", "Kg/Kg/s",
154	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
155	. typeval, zsto,zout)
156	c
157	CALL histdef(nid_bilKPins, "dtcon",
158	. "Convection dT", "K/s",
159	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
160	. typeval, zsto,zout)
161	c
162	CALL histdef(nid_bilKPins, "ducon",
163	. "Convection du", "m/s2",
164	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
165	. typeval, zsto,zout)
166	c
167	CALL histdef(nid_bilKPins, "dvcon",
168	. "Convection dv", "m/s2",
169	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
170	. typeval, zsto,zout)
171	c
172	CALL histdef(nid_bilKPins,"dqcon",
173	. "Convection dQ","Kg/Kg/s",
174	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
175	. typeval, zsto,zout)
176	c
177	CALL histdef(nid_bilKPins, "dtlsc",
178	. "Condensation dT", "K/s",
179	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
180	. typeval, zsto,zout)
181	c
182	CALL histdef(nid_bilKPins,"dqlsc",
183	. "Condensation dQ","Kg/Kg/s",
184	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
185	. typeval, zsto,zout)
186	c
187	CALL histdef(nid_bilKPins,"dtvdf",
188	. "Boundary-layer dT","K/s",
189	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
190	. typeval, zsto,zout)
191	c
192	CALL histdef(nid_bilKPins, "dqvdf",
193	. "Boundary-layer dQ",
194	. "Kg/Kg/s",
195	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
196	. typeval, zsto,zout)
197	c
198	CALL histdef(nid_bilKPins,"dtajs",
199	. "Ajustement sec dT","K/s",
200	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
201	. typeval, zsto,zout)
202	c
203	CALL histdef(nid_bilKPins,"dqajs",
204	. "Ajustement sec dQ",
205	. "Kg/Kg/s",
206	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
207	. typeval, zsto,zout)
208	c
209	CALL histdef(nid_bilKPins,"dteva",
210	. "Reevaporation dT","K/s",
211	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
212	. typeval, zsto,zout)
213	c
214	CALL histdef(nid_bilKPins,"dqeva",
215	. "Reevaporation dQ",
216	. "Kg/Kg/s",
217	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
218	. typeval, zsto,zout)
219
220	c
221	CALL histdef(nid_bilKPins, "dtswr",
222	. "SW radiation dT", "K/s",
223	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
224	. typeval, zsto,zout)
225	c
226	CALL histdef(nid_bilKPins, "dtsw0",
227	. "SW radiation dT", "K/s",
228	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
229	. typeval, zsto,zout)
230	c
231	CALL histdef(nid_bilKPins, "dtlwr",
232	. "LW radiation dT", "K/s",
233	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
234	. typeval, zsto,zout)
235	c
236	CALL histdef(nid_bilKPins, "dtlw0",
237	. "LW radiation dT", "K/s",
238	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
239	. typeval, zsto,zout)
240	c
241	CALL histdef(nid_bilKPins,"duvdf",
242	. "Boundary-layer dU","m/s2",
243	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
244	. typeval, zsto,zout)
245	c
246	CALL histdef(nid_bilKPins,"dvvdf",
247	. "Boundary-layer dV","m/s2",
248	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
249	. typeval, zsto,zout)
250	c
251	IF (ok_orodr) THEN
252	IF (ok_orolf) THEN
253	CALL histdef(nid_bilKPins, "duoli",
254	. "Orography dU", "m/s2",
255	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
256	. typeval, zsto,zout)
257	c
258	CALL histdef(nid_bilKPins, "dvoli",
259	. "Orography dV", "m/s2",
260	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
261	. typeval, zsto,zout)
262	c
263	ENDIF
264	ENDIF
265	C
266	CALL histdef(nid_bilKPins, "duphy",
267	. "Physiq dU","-",
268	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
269	. typeval, zsto,zout)
270	C
271	CALL histdef(nid_bilKPins, "dvphy",
272	. "Physiq dV","-",
273	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
274	. typeval, zsto,zout)
275	C
276	CALL histdef(nid_bilKPins, "dtphy",
277	. "Physiq dT","-",
278	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
279	. typeval, zsto,zout)
280	C
281	CALL histdef(nid_bilKPins, "dqphy",
282	. "Physiq dQ","-",
283	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
284	. typeval, zsto,zout)
285	C
286	CALL histdef(nid_bilKPins, "dqlphy",
287	. "Physiq dQl","-",
288	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
289	. typeval, zsto,zout)
290	cIM 280405 BEG
291	c
292	c Champs 2D:
293	c
294	c u850, v850
295	c DO k=1, nlevSTD
296	DO k=1, 12
297	c
298	IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
299	IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
300	c
301	IF(bb2.EQ."850") THEN
302	c
303	CALL histdef(nid_bilKPins, "u"//bb2,
304	. "Zonal wind "//bb2//"mb","m/s",
305	. iim,jjmp1,nhori, 1,1,1, -99, 32,
306	. typeval, zsto,zout)
307	c
308	CALL histdef(nid_bilKPins, "v"//bb2,
309	. "Meridional wind "//bb2//"mb","m/s",
310	. iim,jjmp1,nhori, 1,1,1, -99, 32,
311	. typeval, zsto,zout)
312	c
313	ENDIF !(bb2.EQ."850")
314	c
315	ENDDO !k=1, 12
316	c
317	cIM 280405 END
318	c
319	CALL histend(nid_bilKPins)
320	c
321	ndex2d = 0
322	ndex3d = 0
323	c
324	ENDIF ! fin de test sur ok_journe

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

Download in other formats: