Context Navigation

source: LMDZ4/trunk/libf/phylmd/ini_histmth.h @ 548

Last change on this file since 548 was 548, checked in by lmdzadmin, 20 years ago
This commit was generated by cvs2svn to compensate for changes in r547, which included commits to RCS files with non-trunk default branches.
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 29.1 KB

Line
1	!
2	! $Header$
3	!
4	IF (ok_mensuel) THEN
5	c
6	zsto = dtime
7	zout = dtime * ecrit_mth
8	c
9	c zsto1: pour des valeurs "instantannees" mensuelles
10	zsto1 = dtime * ecrit_mth
11	c zsto2: pour des flux radiatifs calcules tous les 2 heures
12	zsto2 = dtime * radpas
13	PRINT*,' zsto,zsto1,zsto2,zout=',zsto, zsto1, zsto2,zout
14	c
15	idayref = day_ref
16	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
17	c
18	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
19	DO i = 1, iim
20	zx_lon(i,1) = rlon(i+1)
21	zx_lon(i,jjmp1) = rlon(i+1)
22	ENDDO
23	DO ll=1,klev
24	znivsig(ll)=float(ll)
25	ENDDO
26	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
27	CALL histbeg("histmth.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
28	. 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
29	. nhori, nid_mth)
30	write(,)'Mensuel ', itau_phy, zjulian
31	CALL histvert(nid_mth, "presnivs", "Vertical levels", "mb",
32	. klev, presnivs, nvert)
33	c call histvert(nid_mth, 'sig_s', 'Niveaux sigma','-',
34	c . klev, znivsig, nvert)
35	c
36	c
37	IF(lev_histmth.GE.1) THEN
38	CALL histdef(nid_mth, "phis", "Surface geop. height", "-",
39	. iim,jjmp1,nhori, 1,1,1, -99, 32,
40	. "once", zsto,zout)
41	c
42	CALL histdef(nid_mth, "aire", "Grid area", "-",
43	. iim,jjmp1,nhori, 1,1,1, -99, 32,
44	. "once", zsto,zout)
45	c
46	call histdef(nid_mth, "pourc_"//clnsurf(is_ter),
47	$ "Fraction "//clnsurf(is_ter), "W/m2",
48	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
49	$ "once", zsto,zout)
50	c
51	c Champs 2D:
52	c
53	CALL histdef(nid_mth, "slp", "Sea Level Pressure", "Pa",
54	. iim,jjmp1,nhori, 1,1,1, -99, 32,
55	. "ave(X)", zsto,zout)
56	c
57	CALL histdef(nid_mth, "tsol", "Surface Temperature", "K",
58	. iim,jjmp1,nhori, 1,1,1, -99, 32,
59	. "ave(X)", zsto,zout)
60	c
61	CALL histdef(nid_mth, "t2m", "Temperature 2m", "K",
62	. iim,jjmp1,nhori, 1,1,1, -99, 32,
63	. "ave(X)", zsto,zout)
64	c ENSEMBLES BEG
65	CALL histdef(nid_mth, "t2m_min", "Temp. 2m min.",
66	. "K",
67	. iim,jjmp1,nhori, 1,1,1, -99, 32,
68	. t2mincels, zsto,zout)
69	c
70	CALL histdef(nid_mth, "t2m_max", "Temp. 2m max.",
71	. "K",
72	. iim,jjmp1,nhori, 1,1,1, -99, 32,
73	. t2maxcels, zsto,zout)
74	c
75	c CALL histdef(nid_mth, "tsoil", "Sfce soil Temperature",
76	c . "K",
77	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
78	c . "ave(X)", zsto,zout)
79	c
80	CALL histdef(nid_mth, "wind10m","10-m wind speed","m/s",
81	. iim,jjmp1,nhori, 1,1,1, -99, 32,
82	. "ave(X)", zsto,zout)
83	c
84	CALL histdef(nid_mth, "sicf", "Sea-ice fraction", "-",
85	. iim,jjmp1,nhori, 1,1,1, -99, 32,
86	. "ave(X)", zsto,zout)
87
88	c
89	c ENSEMBLES END
90	CALL histdef(nid_mth, "q2m", "Specific humidity 2m", "kg/kg",
91	. iim,jjmp1,nhori, 1,1,1, -99, 32,
92	. "ave(X)", zsto,zout)
93	c
94	CALL histdef(nid_mth, "u10m", "Vent zonal 10m", "m/s",
95	. iim,jjmp1,nhori, 1,1,1, -99, 32,
96	. "ave(X)", zsto,zout)
97	c
98	CALL histdef(nid_mth, "v10m", "Vent meridien 10m", "m/s",
99	. iim,jjmp1,nhori, 1,1,1, -99, 32,
100	. "ave(X)", zsto,zout)
101	c
102	CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa",
103	. iim,jjmp1,nhori, 1,1,1, -99, 32,
104	. "ave(X)", zsto,zout)
105	c
106	CALL histdef(nid_mth, "qsurf", "Surface Air humidity", "kg/kg",
107	. iim,jjmp1,nhori, 1,1,1, -99, 32,
108	. "ave(X)", zsto,zout)
109	c
110	if (.not. ok_veget) then
111	CALL histdef(nid_mth, "qsol", "Soil watter content", "mm",
112	. iim,jjmp1,nhori, 1,1,1, -99, 32,
113	. "ave(X)", zsto,zout)
114	endif
115	c
116	CALL histdef(nid_mth, "ndayrain",
117	. "Number of day with rain (liq+sol)", "-",
118	. iim,jjmp1,nhori, 1,1,1, -99, 32,
119	. "inst(X)", zsto1,zout)
120	c
121	CALL histdef(nid_mth, "precip", "Precipitation Totale liq+sol",
122	. "kg/(s*m2)",
123	. iim,jjmp1,nhori, 1,1,1, -99, 32,
124	. "ave(X)", zsto,zout)
125	c
126	CALL histdef(nid_mth, "plul", "Large-scale Precip.",
127	. "kg/(s*m2)",
128	. iim,jjmp1,nhori, 1,1,1, -99, 32,
129	. "ave(X)", zsto,zout)
130	c
131	CALL histdef(nid_mth, "pluc", "Convective Precip.",
132	. "kg/(s*m2)",
133	. iim,jjmp1,nhori, 1,1,1, -99, 32,
134	. "ave(X)", zsto,zout)
135	c
136	CALL histdef(nid_mth, "snow", "Snow fall", "kg/(s*m2)",
137	. iim,jjmp1,nhori, 1,1,1, -99, 32,
138	. "ave(X)", zsto,zout)
139	c
140	CALL histdef(nid_mth, "snow_mass", "Snow Mass", "kg/m2",
141	. iim,jjmp1,nhori, 1,1,1, -99, 32,
142	. "ave(X)", zsto,zout)
143	c
144	CALL histdef(nid_mth, "evap", "Evaporation", "kg/(s*m2)",
145	. iim,jjmp1,nhori, 1,1,1, -99, 32,
146	. "ave(X)", zsto,zout)
147	c
148	CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2",
149	. iim,jjmp1,nhori, 1,1,1, -99, 32,
150	. "ave(X)", zsto,zout)
151	c
152	CALL histdef(nid_mth, "tops0", "Solar rad. at TOA", "W/m2",
153	. iim,jjmp1,nhori, 1,1,1, -99, 32,
154	. "ave(X)", zsto,zout)
155	c
156	CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2",
157	. iim,jjmp1,nhori, 1,1,1, -99, 32,
158	. "ave(X)", zsto,zout)
159	c
160	CALL histdef(nid_mth, "topl0", "IR rad. at TOA", "W/m2",
161	. iim,jjmp1,nhori, 1,1,1, -99, 32,
162	. "ave(X)", zsto,zout)
163	c
164	CALL histdef(nid_mth, "SWupTOA", "SWup at TOA","W/m2",
165	. iim,jjmp1,nhori, 1,1,1,-99,
166	. 32, "ave(X)", zsto2,zout)
167	c
168	CALL histdef(nid_mth, "SWupTOAclr",
169	. "SWup clear sky at TOA","W/m2",
170	. iim,jjmp1,nhori, 1,1,1,-99,
171	. 32, "ave(X)", zsto2,zout)
172	c
173	CALL histdef(nid_mth, "SWdnTOA", "SWdn at TOA","W/m2",
174	. iim,jjmp1,nhori, 1,1,1,-99,
175	. 32, "ave(X)", zsto2,zout)
176	c
177	CALL histdef(nid_mth, "SWdnTOAclr",
178	. "SWdn clear sky at TOA","W/m2",
179	. iim,jjmp1,nhori, 1,1,1,-99,
180	. 32, "ave(X)", zsto2,zout)
181	c
182	CALL histdef(nid_mth, "SWup200", "SWup at 200hPa","W/m2",
183	. iim,jjmp1,nhori, 1,1,1,-99,
184	. 32, "ave(X)", zsto2,zout)
185	c
186	CALL histdef(nid_mth, "SWup200clr",
187	. "SWup clear sky at 200hPa","W/m2",
188	. iim,jjmp1,nhori, 1,1,1,-99,
189	. 32, "ave(X)", zsto2,zout)
190	c
191	CALL histdef(nid_mth, "SWdn200", "SWdn at 200hPa","W/m2",
192	. iim,jjmp1,nhori, 1,1,1,-99,
193	. 32, "ave(X)", zsto2,zout)
194	c
195	CALL histdef(nid_mth, "SWdn200clr",
196	. "SWdn clear sky at 200hPa","W/m2",
197	. iim,jjmp1,nhori, 1,1,1,-99,
198	. 32, "ave(X)", zsto2,zout)
199	c
200	CALL histdef(nid_mth, "LWup200", "LWup at 200hPa","W/m2",
201	. iim,jjmp1,nhori, 1,1,1,-99,
202	. 32, "ave(X)", zsto2,zout)
203	c
204	CALL histdef(nid_mth, "LWup200clr",
205	. "LWup clear sky at 200hPa","W/m2",
206	. iim,jjmp1,nhori, 1,1,1,-99,
207	. 32, "ave(X)", zsto2,zout)
208	c
209	CALL histdef(nid_mth, "LWdn200", "LWdn at 200hPa","W/m2",
210	. iim,jjmp1,nhori, 1,1,1,-99,
211	. 32, "ave(X)", zsto2,zout)
212	c
213	CALL histdef(nid_mth, "LWdn200clr",
214	. "LWdn clear sky at 200hPa","W/m2",
215	. iim,jjmp1,nhori, 1,1,1,-99,
216	. 32, "ave(X)", zsto2,zout)
217	c
218	CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2",
219	. iim,jjmp1,nhori, 1,1,1, -99, 32,
220	. "ave(X)", zsto,zout)
221	c
222	CALL histdef(nid_mth, "sols0", "Solar rad. at surf.", "W/m2",
223	. iim,jjmp1,nhori, 1,1,1, -99, 32,
224	. "ave(X)", zsto,zout)
225	c
226	CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2",
227	. iim,jjmp1,nhori, 1,1,1, -99, 32,
228	. "ave(X)", zsto,zout)
229	c
230	CALL histdef(nid_mth, "soll0", "IR rad. at surface", "W/m2",
231	. iim,jjmp1,nhori, 1,1,1, -99, 32,
232	. "ave(X)", zsto,zout)
233	c
234	CALL histdef(nid_mth, "SWupSFC", "SWup at surface","W/m2",
235	. iim,jjmp1,nhori, 1,1,1,-99,
236	. 32, "ave(X)", zsto2,zout)
237	c
238	CALL histdef(nid_mth, "SWupSFCclr",
239	. "SWup clear sky at surface","W/m2",
240	. iim,jjmp1,nhori, 1,1,1,-99,
241	. 32, "ave(X)", zsto2,zout)
242	c
243	CALL histdef(nid_mth, "SWdnSFC", "SWdn at surface","W/m2",
244	. iim,jjmp1,nhori, 1,1,1,-99,
245	. 32, "ave(X)", zsto2,zout)
246	c
247	CALL histdef(nid_mth, "SWdnSFCclr",
248	. "SWdn clear sky at surface","W/m2",
249	. iim,jjmp1,nhori, 1,1,1,-99,
250	. 32, "ave(X)", zsto2,zout)
251	c
252	CALL histdef(nid_mth,"LWupSFC","Upwd. IR rad. at surface",
253	. "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
254	. "ave(X)", zsto2,zout)
255	c
256	CALL histdef(nid_mth,"LWdnSFC","Down. IR rad. at surface",
257	. "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
258	. "ave(X)", zsto2,zout)
259	c
260	CALL histdef(nid_mth,"LWupSFCclr",
261	. "CS Upwd. IR rad. at surface",
262	. "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
263	. "ave(X)", zsto2,zout)
264	c
265	CALL histdef(nid_mth,"LWdnSFCclr",
266	. "Down. CS IR rad. at surface",
267	. "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
268	. "ave(X)", zsto2,zout)
269	c
270	CALL histdef(nid_mth, "bils", "Surf. total heat flux", "W/m2",
271	. iim,jjmp1,nhori, 1,1,1, -99, 32,
272	. "ave(X)", zsto,zout)
273	c
274	CALL histdef(nid_mth, "sens", "Sensible heat flux", "W/m2",
275	. iim,jjmp1,nhori, 1,1,1, -99, 32,
276	. "ave(X)", zsto,zout)
277	c
278	CALL histdef(nid_mth, "fder", "Heat flux derivation", "W/m2",
279	. iim,jjmp1,nhori, 1,1,1, -99, 32,
280	. "ave(X)", zsto,zout)
281	c
282	c Effets des aerosols
283	c
284	c IF (ok_ade.OR.ok_aie) THEN
285	CALL histdef(nid_mth, "topsad", "ADE at TOA", "W/m2",
286	. iim,jjmp1,nhori, 1,1,1, -99, 32,
287	. "ave(X)", zsto,zout)
288	c
289	CALL histdef(nid_mth, "solsad", "ADE at sfc", "W/m2",
290	. iim,jjmp1,nhori, 1,1,1, -99, 32,
291	. "ave(X)", zsto,zout)
292	c
293	CALL histdef(nid_mth, "topsai", "AIE at TOA", "W/m2",
294	. iim,jjmp1,nhori, 1,1,1, -99, 32,
295	. "ave(X)", zsto,zout)
296	c
297	CALL histdef(nid_mth, "solsai", "AIE at sfc", "W/m2",
298	. iim,jjmp1,nhori, 1,1,1, -99, 32,
299	. "ave(X)", zsto,zout)
300	c endif
301	c
302
303	c
304	c CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa",
305	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
306	c . "ave(X)", zsto,zout)
307	c
308	c CALL histdef(nid_mth, "frtv", "Meridional wind stress", "Pa",
309	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
310	c . "ave(X)", zsto,zout)
311	c
312	CALL histdef(nid_mth, "ffonte","Thermal flux for snow melting",
313	. "W/m2",iim,jjmp1,nhori, 1,1,1, -99, 32,
314	. "ave(X)", zsto,zout)
315	c
316	CALL histdef(nid_mth, "fqcalving","Ice Calving",
317	. "kg/m2/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
318	. "ave(X)", zsto,zout)
319	cIM: 171003
320	DO nsrf = 1, nbsrf
321	call histdef(nid_mth, "taux_"//clnsurf(nsrf),
322	$ "Zonal wind stress"//clnsurf(nsrf), "Pa",
323	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
324	$ "ave(X)", zsto,zout)
325
326	call histdef(nid_mth, "tauy_"//clnsurf(nsrf),
327	$ "Meridional xind stress "//clnsurf(nsrf), "Pa",
328	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
329	$ "ave(X)", zsto,zout)
330	ENDDO
331	cIM: 171003
332	c
333	DO nsrf = 1, nbsrf
334	C
335	IF(nsrf.GT.1) THEN
336	call histdef(nid_mth, "pourc_"//clnsurf(nsrf),
337	$ "Fraction "//clnsurf(nsrf), "W/m2",
338	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
339	$ "ave(X)", zsto,zout)
340	ENDIF !nsrf.GT.1
341	C
342	call histdef(nid_mth, "tsol_"//clnsurf(nsrf),
343	$ "Fraction "//clnsurf(nsrf), "W/m2",
344	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
345	$ "ave(X)", zsto,zout)
346	C
347	call histdef(nid_mth, "sens_"//clnsurf(nsrf),
348	$ "Sensible heat flux "//clnsurf(nsrf), "W/m2",
349	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
350	$ "ave(X)", zsto,zout)
351	c
352	call histdef(nid_mth, "lat_"//clnsurf(nsrf),
353	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
354	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
355	$ "ave(X)", zsto,zout)
356	C
357	call histdef(nid_mth, "flw_"//clnsurf(nsrf),
358	$ "LW "//clnsurf(nsrf), "W/m2",
359	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
360	$ "ave(X)", zsto2,zout)
361	c
362	call histdef(nid_mth, "fsw_"//clnsurf(nsrf),
363	$ "SW "//clnsurf(nsrf), "W/m2",
364	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
365	$ "ave(X)", zsto2,zout)
366	C
367	call histdef(nid_mth, "wbils_"//clnsurf(nsrf),
368	$ "Bilan sol "//clnsurf(nsrf), "W/m2",
369	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
370	$ "ave(X)", zsto,zout)
371	C
372	END DO
373	c
374	CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
375	. iim,jjmp1,nhori, 1,1,1, -99, 32,
376	. "ave(X)", zsto,zout)
377	c
378	CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
379	. iim,jjmp1,nhori, 1,1,1, -99, 32,
380	. "ave(X)", zsto,zout)
381	c
382	CALL histdef(nid_mth, "cldl", "Low-level cloudiness", "-",
383	. iim,jjmp1,nhori, 1,1,1, -99, 32,
384	. "ave(X)", zsto,zout)
385	c
386	CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", "-",
387	. iim,jjmp1,nhori, 1,1,1, -99, 32,
388	. "ave(X)", zsto,zout)
389	c
390	CALL histdef(nid_mth, "cldh", "High-level cloudiness", "-",
391	. iim,jjmp1,nhori, 1,1,1, -99, 32,
392	. "ave(X)", zsto,zout)
393	c
394	CALL histdef(nid_mth, "cldt", "Total cloudiness", "-",
395	. iim,jjmp1,nhori, 1,1,1, -99, 32,
396	. "ave(X)", zsto,zout)
397	c
398	CALL histdef(nid_mth,"cldq","Cloud liquid water path","kg/m2",
399	. iim,jjmp1,nhori, 1,1,1, -99, 32,
400	. "ave(X)", zsto,zout)
401	cIM: 071003
402	CALL histdef(nid_mth,"lwp","Cloud water path","kg/m2",
403	. iim,jjmp1,nhori, 1,1,1, -99, 32,
404	. "ave(X)", zsto,zout)
405	c
406	CALL histdef(nid_mth,"iwp","Cloud ice water path","kg/m2",
407	. iim,jjmp1,nhori, 1,1,1, -99, 32,
408	. "ave(X)", zsto,zout)
409	cIM: 071003
410	c
411	CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
412	. iim,jjmp1,nhori, 1,1,1, -99, 32,
413	. "ave(X)", zsto,zout)
414	c
415	CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
416	. iim,jjmp1,nhori, 1,1,1, -99, 32,
417	. "ave(X)", zsto,zout)
418	c
419	CALL histdef(nid_mth, "uq", "Zonal humidity transport", "-",
420	. iim,jjmp1,nhori, 1,1,1, -99, 32,
421	. "ave(X)", zsto,zout)
422	c
423	CALL histdef(nid_mth, "vq", "Merid humidity transport", "-",
424	. iim,jjmp1,nhori, 1,1,1, -99, 32,
425	. "ave(X)", zsto,zout)
426	cKE43
427	IF(iflag_con.GE.3) THEN ! sb
428	c
429	CALL histdef(nid_mth, "cape", "Conv avlbl pot ener", "J/kg",
430	. iim,jjmp1,nhori, 1,1,1, -99, 32,
431	. "ave(X)", zsto,zout)
432	c
433	CALL histdef(nid_mth, "pbase", "Cld base pressure", "hPa",
434	. iim,jjmp1,nhori, 1,1,1, -99, 32,
435	. "ave(X)", zsto,zout)
436	c
437	CALL histdef(nid_mth, "ptop", "Cld top pressure", "hPa",
438	. iim,jjmp1,nhori, 1,1,1, -99, 32,
439	. "ave(X)", zsto,zout)
440	c
441	CALL histdef(nid_mth, "fbase", "Cld base mass flux", "kg/m2/s",
442	. iim,jjmp1,nhori, 1,1,1, -99, 32,
443	. "ave(X)", zsto,zout)
444	c
445	CALL histdef(nid_mth, "prw", "Precipitable water", "kg/m2",
446	. iim,jjmp1,nhori, 1,1,1, -99, 32,
447	. "ave(X)", zsto,zout)
448	c
449	ENDIF !iflag_con .GE. 3
450	c34EK
451	c
452	c Champs interpolles sur des niveaux de pression
453	DO k=1, nlevENS
454	CALL histdef(nid_mth, "u"//clev(k),
455	. "Zonal wind"//clev(k)//"mb","m/s",
456	. iim,jjmp1,nhori, 1,1,1, -99, 32,
457	. "ave(X)", zsto,zout)
458	c
459	CALL histdef(nid_mth, "v"//clev(k),
460	. "Meridional wind"//clev(k)//"mb","m/s",
461	. iim,jjmp1,nhori, 1,1,1, -99, 32,
462	. "ave(X)", zsto,zout)
463	c
464	CALL histdef(nid_mth, "w"//clev(k),
465	. "Vertical wind"//clev(k)//"mb","m/s",
466	. iim,jjmp1,nhori, 1,1,1, -99, 32,
467	. "ave(X)", zsto,zout)
468	c
469	CALL histdef(nid_mth, "phi"//clev(k),
470	. "Geopotential"//clev(k)//"mb","m2/s2",
471	. iim,jjmp1,nhori, 1,1,1, -99, 32,
472	. "ave(X)", zsto,zout)
473	c
474	ENDDO
475	ENDIF !lev_histmth.GE.1
476	c
477	IF(lev_histmth.GE.2) THEN
478	c
479	c Champs 3D:
480	c
481	c
482	CALL histdef(nid_mth,"lwcon","Cloud water content","kg/kg",
483	. iim,jjmp1,nhori, klev,1,klev, nvert, 32,
484	. "ave(X)", zsto,zout)
485	c
486	CALL histdef(nid_mth,"iwcon","Cloud ice water content","kg/kg",
487	. iim,jjmp1,nhori, klev,1,klev, nvert, 32,
488	. "ave(X)", zsto,zout)
489	c
490	CALL histdef(nid_mth, "temp", "Air temperature", "K",
491	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
492	. "ave(X)", zsto,zout)
493	c
494	CALL histdef(nid_mth, "ovap", "Specific humidity", "kg/kg",
495	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
496	. "ave(X)", zsto,zout)
497	c
498	c CALL histdef(nid_mth,"wvap","Water vapor mixing ratio","kg/kg",
499	c . iim,jjmp1,nhori, klev,1,klev,nvert, 32,
500	c . "ave(X)", zsto,zout)
501	c
502	CALL histdef(nid_mth, "geop", "Geopotential height", "m",
503	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
504	. "ave(X)", zsto,zout)
505	c
506	CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s",
507	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
508	. "ave(X)", zsto,zout)
509	c
510	CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s",
511	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
512	. "ave(X)", zsto,zout)
513	c
514	CALL histdef(nid_mth, "vitw", "Vertical wind", "m/s",
515	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
516	. "ave(X)", zsto,zout)
517	c
518	CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
519	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
520	. "ave(X)", zsto,zout)
521	c
522	CALL histdef(nid_mth, "rneb", "Cloud fraction", "-",
523	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
524	. "ave(X)", zsto,zout)
525	c
526	CALL histdef(nid_mth, "rnebcon", "Convective Cloud Fraction"
527	. , "-",
528	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
529	. "ave(X)", zsto,zout)
530	c
531	CALL histdef(nid_mth, "rhum", "Relative humidity", "-",
532	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
533	. "ave(X)", zsto,zout)
534	c
535	CALL histdef(nid_mth, "ozone", "Ozone concentration", "-",
536	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
537	. "ave(X)", zsto,zout)
538	c
539	CALL histdef(nid_mth, "upwd", "saturated updraft", "kg/m2/s",
540	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
541	. "ave(X)", zsto,zout)
542	c
543	CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
544	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
545	. "ave(X)", zsto,zout)
546	c
547	CALL histdef(nid_mth, "dqphy", "Physics dQ", "kg/kg/s",
548	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
549	. "ave(X)", zsto,zout)
550	c
551	c#define histmthNMC
552	c#ifdef histmthNMC
553	c#include "ini_histmthNMC.h"
554	c#endif
555	c
556	ENDIF !lev_histmth.GE.2
557	c
558	IF(lev_histmth.GE.3) THEN
559	c
560	DO nsrf=1, nbsrf
561	c
562	call histdef(nid_mth, "albe_"//clnsurf(nsrf),
563	$ "Albedo surf. "//clnsurf(nsrf), "-",
564	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
565	$ "ave(X)", zsto,zout)
566	c
567	call histdef(nid_mth, "rugs_"//clnsurf(nsrf),
568	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
569	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
570	$ "ave(X)", zsto,zout)
571	c
572	CALL histdef(nid_mth, "ages_"//clnsurf(nsrf), "Snow age","day",
573	. iim,jjmp1,nhori, 1,1,1, -99, 32,
574	. "ave(X)", zsto,zout)
575	c
576	ENDDO !nsrf=1, nbsrf
577	c
578	CALL histdef(nid_mth, "albs", "Surface albedo", "-",
579	. iim,jjmp1,nhori, 1,1,1, -99, 32,
580	. "ave(X)", zsto,zout)
581	CALL histdef(nid_mth, "albslw", "Surface albedo LW", "-",
582	. iim,jjmp1,nhori, 1,1,1, -99, 32,
583	. "ave(X)", zsto,zout)
584	c
585	ENDIF !lev_histmth.GE.3
586	c
587	IF(lev_histmth.GE.4) THEN
588	c
589	CALL histdef(nid_mth, "clwcon",
590	. "Convective Cloud Liquid water content"
591	. , "kg/kg",
592	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
593	. "ave(X)", zsto,zout)
594	c
595	CALL histdef(nid_mth,"Ma","undilute adiab updraft","kg/m2/s",
596	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
597	. "ave(X)", zsto,zout)
598	c
599	CALL histdef(nid_mth, "dnwd", "saturated downdraft","kg/m2/s",
600	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
601	. "ave(X)", zsto,zout)
602	c
603	CALL histdef(nid_mth, "dnwd0", "unsat. downdraft", "kg/m2/s",
604	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
605	. "ave(X)", zsto,zout)
606	c
607	CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
608	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
609	. "ave(X)", zsto,zout)
610	c
611	CALL histdef(nid_mth, "dqdyn", "Dynamics dQ", "kg/kg/s",
612	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
613	. "ave(X)", zsto,zout)
614	c
615	CALL histdef(nid_mth, "dtcon", "Convection dT", "K/s",
616	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
617	. "ave(X)", zsto,zout)
618	c
619	c CALL histdef(nid_mth, "ducon", "Convection du", "m/s2",
620	c . iim,jjmp1,nhori, klev,1,klev,nvert, 32,
621	c . "ave(X)", zsto,zout)
622	c
623	CALL histdef(nid_mth, "dqcon", "Convection dQ", "kg/kg/s",
624	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
625	. "ave(X)", zsto,zout)
626	c
627	CALL histdef(nid_mth, "dtlsc", "Condensation dT", "K/s",
628	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
629	. "ave(X)", zsto,zout)
630	cIM: 071003
631	CALL histdef(nid_mth, "dtlschr",
632	$ "Large-scale condensational heating rate", "K/s",iim,jjmp1
633	$ ,nhori, klev,1,klev,nvert, 32,"ave(X)", zsto,zout)
634	c
635	CALL histdef(nid_mth, "dqlsc", "Condensation dQ", "kg/kg/s",
636	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
637	. "ave(X)", zsto,zout)
638	c
639	CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s",
640	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
641	. "ave(X)", zsto,zout)
642	c
643	CALL histdef(nid_mth, "dqvdf", "Boundary-layer dQ", "kg/kg/s",
644	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
645	. "ave(X)", zsto,zout)
646	c
647	CALL histdef(nid_mth, "dteva", "Reevaporation dT", "K/s",
648	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
649	. "ave(X)", zsto,zout)
650	c
651	CALL histdef(nid_mth, "dqeva", "Reevaporation dQ", "kg/kg/s",
652	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
653	. "ave(X)", zsto,zout)
654
655	CALL histdef(nid_mth, "ptconv", "POINTS CONVECTIFS"," ",
656	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
657	. "ave(X)", zsto,zout)
658
659	CALL histdef(nid_mth, "ratqs", "RATQS"," ",
660	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
661	. "ave(X)", zsto,zout)
662
663	c
664	CALL histdef(nid_mth, "dtajs", "Dry adjust. dT", "K/s",
665	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
666	. "ave(X)", zsto,zout)
667
668	CALL histdef(nid_mth, "dqajs", "Dry adjust. dQ", "kg/kg/s",
669	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
670	. "ave(X)", zsto,zout)
671	c
672	CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s",
673	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
674	. "ave(X)", zsto,zout)
675	c
676	CALL histdef(nid_mth, "dtsw0", "CS SW radiation dT", "K/s",
677	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
678	. "ave(X)", zsto,zout)
679	c
680	CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
681	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
682	. "ave(X)", zsto,zout)
683	c
684	CALL histdef(nid_mth,"dtlw0","CS LW radiation dT","K/s",
685	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
686	. "ave(X)", zsto,zout)
687	c
688	CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
689	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
690	. "ave(X)", zsto,zout)
691	c
692	CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
693	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
694	. "ave(X)", zsto,zout)
695	c
696	CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
697	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
698	. "ave(X)", zsto,zout)
699	c
700	IF (ok_orodr) THEN
701	CALL histdef(nid_mth, "duoro", "Orography dU", "m/s2",
702	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
703	. "ave(X)", zsto,zout)
704	c
705	CALL histdef(nid_mth, "dvoro", "Orography dV", "m/s2",
706	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
707	. "ave(X)", zsto,zout)
708	c
709	ENDIF
710	C
711	IF (ok_orolf) THEN
712	CALL histdef(nid_mth, "dulif", "Orography dU", "m/s2",
713	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
714	. "ave(X)", zsto,zout)
715	c
716	CALL histdef(nid_mth, "dvlif", "Orography dV", "m/s2",
717	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
718	. "ave(X)", zsto,zout)
719	ENDIF
720	c
721	c Effets des aerosols
722	c
723	c IF (ok_ade.OR.ok_aie) THEN
724	CALL histdef(nid_mth, "re", "CDR", "um",
725	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
726	. "ave(X)", zsto,zout)
727	c
728	CALL histdef(nid_mth, "redenom", "CDR denominator", "-",
729	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
730	. "ave(X)", zsto,zout)
731	c
732	CALL histdef(nid_mth, "tau", "cloud opt thickness", "-",
733	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
734	. "ave(X)", zsto,zout)
735	c
736	CALL histdef(nid_mth, "taupi", "cloud opt thickn. (pi)", "-",
737	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
738	. "ave(X)", zsto,zout)
739	c endif
740	c
741	c CALL histdef(nid_mth, "ozone", "Ozone concentration", "-",
742	c . iim,jjmp1,nhori, klev,1,klev,nvert, 32,
743	c . "ave(X)", zsto,zout)
744	c
745	if (nqmax.GE.3) THEN
746	DO iq=1,nqmax-2
747	IF (iq.LE.99) THEN
748	WRITE(str2,'(i2.2)') iq
749	CALL histdef(nid_mth, "trac"//str2, "Tracer No."//str2, "-",
750	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
751	. "ave(X)", zsto,zout)
752	ELSE
753	PRINT*, "Trop de traceurs"
754	CALL abort
755	ENDIF
756	ENDDO
757	ENDIF
758	c
759	ENDIF !lev_histmth.GE.4
760	c
761	CALL histend(nid_mth)
762	c
763	ndex2d = 0
764	ndex3d = 0
765	c
766	ENDIF ! fin de test sur ok_mensuel

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

Download in other formats: