Context Navigation

ini_histmth.h @ 710

Last change on this file since 710 was 710, checked in by Laurent Fairhead, 18 years ago
Integration des modifications d'OM et JLD pour une meilleure prise en compte de la conservation du flux d'eau LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 29.3 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	c
320	CALL histdef(nid_mth, "fqfonte","Land ice melt",
321	. "kg/m2/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
322	. "ave(X)", zsto,zout)
323	cIM: 171003
324	DO nsrf = 1, nbsrf
325	call histdef(nid_mth, "taux_"//clnsurf(nsrf),
326	$ "Zonal wind stress"//clnsurf(nsrf), "Pa",
327	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
328	$ "ave(X)", zsto,zout)
329
330	call histdef(nid_mth, "tauy_"//clnsurf(nsrf),
331	$ "Meridional xind stress "//clnsurf(nsrf), "Pa",
332	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
333	$ "ave(X)", zsto,zout)
334	ENDDO
335	cIM: 171003
336	c
337	DO nsrf = 1, nbsrf
338	C
339	IF(nsrf.GT.1) THEN
340	call histdef(nid_mth, "pourc_"//clnsurf(nsrf),
341	$ "Fraction "//clnsurf(nsrf), "W/m2",
342	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
343	$ "ave(X)", zsto,zout)
344	ENDIF !nsrf.GT.1
345	C
346	call histdef(nid_mth, "tsol_"//clnsurf(nsrf),
347	$ "Fraction "//clnsurf(nsrf), "W/m2",
348	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
349	$ "ave(X)", zsto,zout)
350	C
351	call histdef(nid_mth, "sens_"//clnsurf(nsrf),
352	$ "Sensible heat flux "//clnsurf(nsrf), "W/m2",
353	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
354	$ "ave(X)", zsto,zout)
355	c
356	call histdef(nid_mth, "lat_"//clnsurf(nsrf),
357	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
358	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
359	$ "ave(X)", zsto,zout)
360	C
361	call histdef(nid_mth, "flw_"//clnsurf(nsrf),
362	$ "LW "//clnsurf(nsrf), "W/m2",
363	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
364	$ "ave(X)", zsto2,zout)
365	c
366	call histdef(nid_mth, "fsw_"//clnsurf(nsrf),
367	$ "SW "//clnsurf(nsrf), "W/m2",
368	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
369	$ "ave(X)", zsto2,zout)
370	C
371	call histdef(nid_mth, "wbils_"//clnsurf(nsrf),
372	$ "Bilan sol "//clnsurf(nsrf), "W/m2",
373	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
374	$ "ave(X)", zsto,zout)
375	C
376	END DO
377	c
378	CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
379	. iim,jjmp1,nhori, 1,1,1, -99, 32,
380	. "ave(X)", zsto,zout)
381	c
382	CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
383	. iim,jjmp1,nhori, 1,1,1, -99, 32,
384	. "ave(X)", zsto,zout)
385	c
386	CALL histdef(nid_mth, "cldl", "Low-level cloudiness", "-",
387	. iim,jjmp1,nhori, 1,1,1, -99, 32,
388	. "ave(X)", zsto,zout)
389	c
390	CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", "-",
391	. iim,jjmp1,nhori, 1,1,1, -99, 32,
392	. "ave(X)", zsto,zout)
393	c
394	CALL histdef(nid_mth, "cldh", "High-level cloudiness", "-",
395	. iim,jjmp1,nhori, 1,1,1, -99, 32,
396	. "ave(X)", zsto,zout)
397	c
398	CALL histdef(nid_mth, "cldt", "Total cloudiness", "-",
399	. iim,jjmp1,nhori, 1,1,1, -99, 32,
400	. "ave(X)", zsto,zout)
401	c
402	CALL histdef(nid_mth,"cldq","Cloud liquid water path","kg/m2",
403	. iim,jjmp1,nhori, 1,1,1, -99, 32,
404	. "ave(X)", zsto,zout)
405	cIM: 071003
406	CALL histdef(nid_mth,"lwp","Cloud water path","kg/m2",
407	. iim,jjmp1,nhori, 1,1,1, -99, 32,
408	. "ave(X)", zsto,zout)
409	c
410	CALL histdef(nid_mth,"iwp","Cloud ice water path","kg/m2",
411	. iim,jjmp1,nhori, 1,1,1, -99, 32,
412	. "ave(X)", zsto,zout)
413	cIM: 071003
414	c
415	CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
416	. iim,jjmp1,nhori, 1,1,1, -99, 32,
417	. "ave(X)", zsto,zout)
418	c
419	CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
420	. iim,jjmp1,nhori, 1,1,1, -99, 32,
421	. "ave(X)", zsto,zout)
422	c
423	CALL histdef(nid_mth, "uq", "Zonal humidity transport", "-",
424	. iim,jjmp1,nhori, 1,1,1, -99, 32,
425	. "ave(X)", zsto,zout)
426	c
427	CALL histdef(nid_mth, "vq", "Merid humidity transport", "-",
428	. iim,jjmp1,nhori, 1,1,1, -99, 32,
429	. "ave(X)", zsto,zout)
430	cKE43
431	IF(iflag_con.GE.3) THEN ! sb
432	c
433	CALL histdef(nid_mth, "cape", "Conv avlbl pot ener", "J/kg",
434	. iim,jjmp1,nhori, 1,1,1, -99, 32,
435	. "ave(X)", zsto,zout)
436	c
437	CALL histdef(nid_mth, "pbase", "Cld base pressure", "hPa",
438	. iim,jjmp1,nhori, 1,1,1, -99, 32,
439	. "ave(X)", zsto,zout)
440	c
441	CALL histdef(nid_mth, "ptop", "Cld top pressure", "hPa",
442	. iim,jjmp1,nhori, 1,1,1, -99, 32,
443	. "ave(X)", zsto,zout)
444	c
445	CALL histdef(nid_mth, "fbase", "Cld base mass flux", "kg/m2/s",
446	. iim,jjmp1,nhori, 1,1,1, -99, 32,
447	. "ave(X)", zsto,zout)
448	c
449	CALL histdef(nid_mth, "prw", "Precipitable water", "kg/m2",
450	. iim,jjmp1,nhori, 1,1,1, -99, 32,
451	. "ave(X)", zsto,zout)
452	c
453	ENDIF !iflag_con .GE. 3
454	c34EK
455	c
456	c Champs interpolles sur des niveaux de pression
457	DO k=1, nlevENS
458	CALL histdef(nid_mth, "u"//clev(k),
459	. "Zonal wind"//clev(k)//"mb","m/s",
460	. iim,jjmp1,nhori, 1,1,1, -99, 32,
461	. "ave(X)", zsto,zout)
462	c
463	CALL histdef(nid_mth, "v"//clev(k),
464	. "Meridional wind"//clev(k)//"mb","m/s",
465	. iim,jjmp1,nhori, 1,1,1, -99, 32,
466	. "ave(X)", zsto,zout)
467	c
468	CALL histdef(nid_mth, "w"//clev(k),
469	. "Vertical wind"//clev(k)//"mb","m/s",
470	. iim,jjmp1,nhori, 1,1,1, -99, 32,
471	. "ave(X)", zsto,zout)
472	c
473	CALL histdef(nid_mth, "phi"//clev(k),
474	. "Geopotential"//clev(k)//"mb","m2/s2",
475	. iim,jjmp1,nhori, 1,1,1, -99, 32,
476	. "ave(X)", zsto,zout)
477	c
478	ENDDO
479	ENDIF !lev_histmth.GE.1
480	c
481	IF(lev_histmth.GE.2) THEN
482	c
483	c Champs 3D:
484	c
485	c
486	CALL histdef(nid_mth,"lwcon","Cloud 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,"iwcon","Cloud ice water content","kg/kg",
491	. iim,jjmp1,nhori, klev,1,klev, nvert, 32,
492	. "ave(X)", zsto,zout)
493	c
494	CALL histdef(nid_mth, "temp", "Air temperature", "K",
495	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
496	. "ave(X)", zsto,zout)
497	c
498	CALL histdef(nid_mth, "ovap", "Specific humidity", "kg/kg",
499	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
500	. "ave(X)", zsto,zout)
501	c
502	c CALL histdef(nid_mth,"wvap","Water vapor mixing ratio","kg/kg",
503	c . iim,jjmp1,nhori, klev,1,klev,nvert, 32,
504	c . "ave(X)", zsto,zout)
505	c
506	CALL histdef(nid_mth, "geop", "Geopotential height", "m",
507	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
508	. "ave(X)", zsto,zout)
509	c
510	CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s",
511	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
512	. "ave(X)", zsto,zout)
513	c
514	CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s",
515	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
516	. "ave(X)", zsto,zout)
517	c
518	CALL histdef(nid_mth, "vitw", "Vertical wind", "m/s",
519	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
520	. "ave(X)", zsto,zout)
521	c
522	CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
523	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
524	. "ave(X)", zsto,zout)
525	c
526	CALL histdef(nid_mth, "rneb", "Cloud fraction", "-",
527	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
528	. "ave(X)", zsto,zout)
529	c
530	CALL histdef(nid_mth, "rnebcon", "Convective Cloud Fraction"
531	. , "-",
532	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
533	. "ave(X)", zsto,zout)
534	c
535	CALL histdef(nid_mth, "rhum", "Relative humidity", "-",
536	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
537	. "ave(X)", zsto,zout)
538	c
539	CALL histdef(nid_mth, "ozone", "Ozone concentration", "-",
540	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
541	. "ave(X)", zsto,zout)
542	c
543	CALL histdef(nid_mth, "upwd", "saturated updraft", "kg/m2/s",
544	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
545	. "ave(X)", zsto,zout)
546	c
547	CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
548	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
549	. "ave(X)", zsto,zout)
550	c
551	CALL histdef(nid_mth, "dqphy", "Physics dQ", "kg/kg/s",
552	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
553	. "ave(X)", zsto,zout)
554	c
555	c#define histmthNMC
556	c#ifdef histmthNMC
557	c#include "ini_histmthNMC.h"
558	c#endif
559	c
560	ENDIF !lev_histmth.GE.2
561	c
562	IF(lev_histmth.GE.3) THEN
563	c
564	DO nsrf=1, nbsrf
565	c
566	call histdef(nid_mth, "albe_"//clnsurf(nsrf),
567	$ "Albedo surf. "//clnsurf(nsrf), "-",
568	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
569	$ "ave(X)", zsto,zout)
570	c
571	call histdef(nid_mth, "rugs_"//clnsurf(nsrf),
572	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
573	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
574	$ "ave(X)", zsto,zout)
575	c
576	CALL histdef(nid_mth, "ages_"//clnsurf(nsrf), "Snow age","day",
577	. iim,jjmp1,nhori, 1,1,1, -99, 32,
578	. "ave(X)", zsto,zout)
579	c
580	ENDDO !nsrf=1, nbsrf
581	c
582	CALL histdef(nid_mth, "albs", "Surface albedo", "-",
583	. iim,jjmp1,nhori, 1,1,1, -99, 32,
584	. "ave(X)", zsto,zout)
585	CALL histdef(nid_mth, "albslw", "Surface albedo LW", "-",
586	. iim,jjmp1,nhori, 1,1,1, -99, 32,
587	. "ave(X)", zsto,zout)
588	c
589	ENDIF !lev_histmth.GE.3
590	c
591	IF(lev_histmth.GE.4) THEN
592	c
593	CALL histdef(nid_mth, "clwcon",
594	. "Convective Cloud Liquid water content"
595	. , "kg/kg",
596	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
597	. "ave(X)", zsto,zout)
598	c
599	CALL histdef(nid_mth,"Ma","undilute adiab updraft","kg/m2/s",
600	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
601	. "ave(X)", zsto,zout)
602	c
603	CALL histdef(nid_mth, "dnwd", "saturated 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, "dnwd0", "unsat. downdraft", "kg/m2/s",
608	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
609	. "ave(X)", zsto,zout)
610	c
611	CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
612	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
613	. "ave(X)", zsto,zout)
614	c
615	CALL histdef(nid_mth, "dqdyn", "Dynamics dQ", "kg/kg/s",
616	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
617	. "ave(X)", zsto,zout)
618	c
619	CALL histdef(nid_mth, "dtcon", "Convection dT", "K/s",
620	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
621	. "ave(X)", zsto,zout)
622	c
623	c CALL histdef(nid_mth, "ducon", "Convection du", "m/s2",
624	c . iim,jjmp1,nhori, klev,1,klev,nvert, 32,
625	c . "ave(X)", zsto,zout)
626	c
627	CALL histdef(nid_mth, "dqcon", "Convection dQ", "kg/kg/s",
628	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
629	. "ave(X)", zsto,zout)
630	c
631	CALL histdef(nid_mth, "dtlsc", "Condensation dT", "K/s",
632	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
633	. "ave(X)", zsto,zout)
634	cIM: 071003
635	CALL histdef(nid_mth, "dtlschr",
636	$ "Large-scale condensational heating rate", "K/s",iim,jjmp1
637	$ ,nhori, klev,1,klev,nvert, 32,"ave(X)", zsto,zout)
638	c
639	CALL histdef(nid_mth, "dqlsc", "Condensation dQ", "kg/kg/s",
640	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
641	. "ave(X)", zsto,zout)
642	c
643	CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s",
644	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
645	. "ave(X)", zsto,zout)
646	c
647	CALL histdef(nid_mth, "dqvdf", "Boundary-layer dQ", "kg/kg/s",
648	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
649	. "ave(X)", zsto,zout)
650	c
651	CALL histdef(nid_mth, "dteva", "Reevaporation dT", "K/s",
652	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
653	. "ave(X)", zsto,zout)
654	c
655	CALL histdef(nid_mth, "dqeva", "Reevaporation dQ", "kg/kg/s",
656	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
657	. "ave(X)", zsto,zout)
658
659	CALL histdef(nid_mth, "ptconv", "POINTS CONVECTIFS"," ",
660	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
661	. "ave(X)", zsto,zout)
662
663	CALL histdef(nid_mth, "ratqs", "RATQS"," ",
664	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
665	. "ave(X)", zsto,zout)
666
667	c
668	CALL histdef(nid_mth, "dtajs", "Dry adjust. dT", "K/s",
669	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
670	. "ave(X)", zsto,zout)
671
672	CALL histdef(nid_mth, "dqajs", "Dry adjust. dQ", "kg/kg/s",
673	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
674	. "ave(X)", zsto,zout)
675	c
676	CALL histdef(nid_mth, "dtswr", "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, "dtsw0", "CS SW 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, "dtlwr", "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,"dtlw0","CS LW radiation dT","K/s",
689	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
690	. "ave(X)", zsto,zout)
691	c
692	CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
693	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
694	. "ave(X)", zsto,zout)
695	c
696	CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
697	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
698	. "ave(X)", zsto,zout)
699	c
700	CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
701	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
702	. "ave(X)", zsto,zout)
703	c
704	IF (ok_orodr) THEN
705	CALL histdef(nid_mth, "duoro", "Orography dU", "m/s2",
706	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
707	. "ave(X)", zsto,zout)
708	c
709	CALL histdef(nid_mth, "dvoro", "Orography dV", "m/s2",
710	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
711	. "ave(X)", zsto,zout)
712	c
713	ENDIF
714	C
715	IF (ok_orolf) THEN
716	CALL histdef(nid_mth, "dulif", "Orography dU", "m/s2",
717	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
718	. "ave(X)", zsto,zout)
719	c
720	CALL histdef(nid_mth, "dvlif", "Orography dV", "m/s2",
721	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
722	. "ave(X)", zsto,zout)
723	ENDIF
724	c
725	c Effets des aerosols
726	c
727	c IF (ok_ade.OR.ok_aie) THEN
728	CALL histdef(nid_mth, "re", "CDR", "um",
729	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
730	. "ave(X)", zsto,zout)
731	c
732	CALL histdef(nid_mth, "redenom", "CDR denominator", "-",
733	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
734	. "ave(X)", zsto,zout)
735	c
736	CALL histdef(nid_mth, "tau", "cloud opt thickness", "-",
737	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
738	. "ave(X)", zsto,zout)
739	c
740	CALL histdef(nid_mth, "taupi", "cloud opt thickn. (pi)", "-",
741	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
742	. "ave(X)", zsto,zout)
743	c endif
744	c
745	c CALL histdef(nid_mth, "ozone", "Ozone concentration", "-",
746	c . iim,jjmp1,nhori, klev,1,klev,nvert, 32,
747	c . "ave(X)", zsto,zout)
748	c
749	if (nqmax.GE.3) THEN
750	DO iq=1,nqmax-2
751	IF (iq.LE.99) THEN
752	WRITE(str2,'(i2.2)') iq
753	CALL histdef(nid_mth, "trac"//str2, "Tracer No."//str2, "-",
754	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
755	. "ave(X)", zsto,zout)
756	ELSE
757	PRINT*, "Trop de traceurs"
758	CALL abort
759	ENDIF
760	ENDDO
761	ENDIF
762	c
763	ENDIF !lev_histmth.GE.4
764	c
765	CALL histend(nid_mth)
766	c
767	ndex2d = 0
768	ndex3d = 0
769	c
770	ENDIF ! fin de test sur ok_mensuel

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

Context Navigation

source: LMDZ4/branches/IPSL-CM4_IPCC_branch/libf/phylmd/ini_histmth.h @ 710

Download in other formats: