Context Navigation

source: LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histmth.h @ 517

Last change on this file since 517 was 517, checked in by lmdzadmin, 20 years ago
Inclusion des modifications de O. Boucher et de J. Quaas pour le calcul des premiers effets directs et indirects dus aux aerosols LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 29.1 KB

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

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

Download in other formats: