Context Navigation

source: LMDZ4/trunk/libf/phylmd/ini_bilKP_ave.h @ 644

Last change on this file since 644 was 644, checked in by Laurent Fairhead, 19 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: 8.6 KB

Line
1	c
2	c $Header$
3	c
4	IF (ok_journe) THEN
5	c
6	zsto = dtime
7	zout = dtime * FLOAT(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: