Context Navigation

ini_histmth.h @ 5456

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

Context Navigation

source: LMDZ4/branches/LMDZ4_par_0/libf/phylmd/ini_histmth.h @ 5456

Download in other formats: