Context Navigation

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

Last change on this file since 673 was 644, checked in by Laurent Fairhead, 20 years ago
Synchronisation avec tous les diagnostiques de Ionela IM Inclusion du slab ocean IM LF
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-1
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: