Context Navigation

source: LMDZ4/trunk/libf/phylmd/ini_histday.h @ 671

Last change on this file since 671 was 644, checked in by Laurent Fairhead, 19 years ago
Synchronisation avec tous les diagnostiques de Ionela IM Inclusion du slab ocean IM LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 22.5 KB

Line
1	c
2	c $Header$
3	c
4	IF (ok_journe) THEN
5	c
6	zstophy = dtime
7	zstoday= dtime * FLOAT(ecrit_day)
8	zout = dtime * FLOAT(ecrit_day)
9	c
10	idayref = day_ref
11	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
12	c
13	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
14	DO i = 1, iim
15	zx_lon(i,1) = rlon(i+1)
16	zx_lon(i,jjmp1) = rlon(i+1)
17	ENDDO
18	DO ll=1,klev
19	znivsig(ll)=float(ll)
20	ENDDO
21	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
22	CALL histbeg("histday", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
23	. 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
24	. nhori, nid_day)
25	CALL histvert(nid_day, "presnivs", "Vertical levels", "mb",
26	. klev, presnivs/100., nvert)
27	c
28	IF(lev_histday.GE.1) THEN
29	c
30	CALL histdef(nid_day, "phis", "Surface geop. height", "m2/s2",
31	. iim,jjmp1,nhori, 1,1,1, -99, 32,
32	. "once", zstophy,zout)
33	c
34	CALL histdef(nid_day, "aire", "Grid area", "-",
35	. iim,jjmp1,nhori, 1,1,1, -99, 32,
36	. "once", zstophy,zout)
37	c
38	CALL histdef(nid_day, "contfracATM","% sfce ter+lic ","-",
39	. iim,jjmp1,nhori, 1,1,1, -99, 32,
40	. "once", zstophy,zout)
41	c
42	CALL histdef(nid_day, "contfracOR","% sfce terre OR", "-",
43	. iim,jjmp1,nhori, 1,1,1, -99, 32,
44	. "once", zstophy,zout)
45	c
46	c Champs 2D:
47	c
48	CALL histdef(nid_day, "tsol", "Surface Temperature", "K",
49	. iim,jjmp1,nhori, 1,1,1, -99, 32,
50	. "ave(X)", zstophy,zout)
51	c
52	CALL histdef(nid_day, "t2m", "Temperature 2m", "K",
53	. iim,jjmp1,nhori, 1,1,1, -99, 32,
54	. "ave(X)", zstophy,zout)
55	c
56	c Champs retires momentannéement en attendant un hypothetique
57	c debugage
58
59	CALL histdef(nid_day, "t2m_min", "Temp. 2m min.",
60	. "K",
61	. iim,jjmp1,nhori, 1,1,1, -99, 32,
62	. t2mincels, zstophy,zout)
63	c
64	CALL histdef(nid_day, "t2m_max", "Temp. 2m max.",
65	. "K",
66	. iim,jjmp1,nhori, 1,1,1, -99, 32,
67	. t2maxcels, zstophy,zout)
68	c
69	CALL histdef(nid_day, "plul", "Large-scale Precip.",
70	. "kg/(s*m2)",
71	. iim,jjmp1,nhori, 1,1,1, -99, 32,
72	. "ave(X)", zstophy,zout)
73	c
74	CALL histdef(nid_day, "pluc", "Convective Precip.",
75	. "kg/(s*m2)",
76	. iim,jjmp1,nhori, 1,1,1, -99, 32,
77	. "ave(X)", zstophy,zout)
78	c
79	CALL histdef(nid_day, "snowl", "Solid Large-scale Precip.",
80	. "kg/(m2*s)",
81	. iim,jjmp1,nhori, 1,1,1, -99, 32,
82	. "ave(X)", zstophy,zout)
83	c
84	IF(1.EQ.0) THEN !snowc=0.
85	CALL histdef(nid_day, "snowc", "Solid Convective Precip.",
86	. "kg/(m2*s)",
87	. iim,jjmp1,nhori, 1,1,1, -99, 32,
88	. "ave(X)", zstophy,zout)
89	ENDIF !snowc=0.
90	c
91	CALL histdef(nid_day, "flat", "Latent heat flux", "W/m2",
92	. iim,jjmp1,nhori, 1,1,1, -99, 32,
93	. "ave(X)", zstophy,zout)
94	c
95	CALL histdef(nid_day, "sicf", "Sea-ice fraction", "-",
96	. iim,jjmp1,nhori, 1,1,1, -99, 32,
97	. "inst(X)", zstoday,zout)
98	c
99	CALL histdef(nid_day, "q2m", "Specific humidity", "kg/kg",
100	. iim,jjmp1,nhori, 1,1,1, -99, 32,
101	. "ave(X)", zstophy,zout)
102	c
103	CALL histdef(nid_day, "u10m", "Vent zonal 10m", "m/s",
104	. iim,jjmp1,nhori, 1,1,1, -99, 32,
105	. "ave(X)", zstophy,zout)
106	c
107	CALL histdef(nid_day, "v10m", "Vent meridien 10m", "m/s",
108	. iim,jjmp1,nhori, 1,1,1, -99, 32,
109	. "ave(X)", zstophy,zout)
110	c
111	CALL histdef(nid_day, "wind10m","10-m wind speed","m/s",
112	. iim,jjmp1,nhori, 1,1,1, -99, 32,
113	. "ave(X)", zstophy,zout)
114	c
115	CALL histdef(nid_day, "wind10max", "10-m wind speed max.",
116	. "m/s",
117	. iim,jjmp1,nhori, 1,1,1, -99, 32,
118	. t2maxcels, zstophy,zout)
119	c
120	CALL histdef(nid_day, "psol", "Surface Pressure", "Pa",
121	. iim,jjmp1,nhori, 1,1,1, -99, 32,
122	. "ave(X)", zstophy,zout)
123	c
124	CALL histdef(nid_day, "precip","Precipitation Totale liq+sol"
125	. , "kg/(s*m2)",
126	. iim,jjmp1,nhori, 1,1,1, -99, 32,
127	. "ave(X)", zstophy,zout)
128	c
129	CALL histdef(nid_day, "snowf", "Snow fall", "kg/(m2*s)",
130	. iim,jjmp1,nhori, 1,1,1, -99, 32,
131	. "ave(X)", zstophy,zout)
132	c
133	c CALL histdef(nid_day, "snow_mass", "Snow Mass", "kg/m2",
134	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
135	c . "ave(X)", zstophy,zout)
136	c
137	CALL histdef(nid_day, "evap", "Evaporation", "kg/(m2*s)",
138	. iim,jjmp1,nhori, 1,1,1, -99, 32,
139	. "ave(X)", zstophy,zout)
140	c
141	CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",
142	. iim,jjmp1,nhori, 1,1,1, -99, 32,
143	. "ave(X)", zstophy,zout)
144	c
145	CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2",
146	. iim,jjmp1,nhori, 1,1,1, -99, 32,
147	. "ave(X)", zstophy,zout)
148	c
149	CALL histdef(nid_day, "sols", "Net Solar rad. at surf.",
150	. "W/m2",
151	. iim,jjmp1,nhori, 1,1,1, -99, 32,
152	. "ave(X)", zstophy,zout)
153	c
154	CALL histdef(nid_day, "soll", "Net IR rad. at surface", "W/m2",
155	. iim,jjmp1,nhori, 1,1,1, -99, 32,
156	. "ave(X)", zstophy,zout)
157	c
158	CALL histdef(nid_day, "cldl", "Low-level cloudiness", "-",
159	. iim,jjmp1,nhori, 1,1,1, -99, 32,
160	. "ave(X)", zstophy,zout)
161	c
162	CALL histdef(nid_day, "cldm", "Mid-level cloudiness", "-",
163	. iim,jjmp1,nhori, 1,1,1, -99, 32,
164	. "ave(X)", zstophy,zout)
165	c
166	CALL histdef(nid_day, "cldh", "High-level cloudiness", "-",
167	. iim,jjmp1,nhori, 1,1,1, -99, 32,
168	. "ave(X)", zstophy,zout)
169	c
170	CALL histdef(nid_day, "cldt", "Total cloudiness", "-",
171	. iim,jjmp1,nhori, 1,1,1, -99, 32,
172	. "ave(X)", zstophy,zout)
173	c
174	CALL histdef(nid_day, "cldq", "Cloud liquid water path",
175	. "kg/m2",
176	. iim,jjmp1,nhori, 1,1,1, -99, 32,
177	. "ave(X)", zstophy,zout)
178	c
179	CALL histdef(nid_day, "prw", "Precipitable water", "kg/m2",
180	. iim,jjmp1,nhori, 1,1,1, -99, 32,
181	. "ave(X)", zstophy,zout)
182	c
183	c Champs dynamiques sur niveaux de pression
184	c
185	c DO k=1, nlevENS
186	DO k=1, nlevSTD
187	c
188	c bb=clevSTD(k)
189	c
190	IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
191	IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
192	c aa=clevSTD(k)
193	c bb2=clevSTD(k)
194	c ENDIF
195	c
196	IF(bb2.EQ."850".OR.bb2.EQ."700"
197	. .OR.bb2.EQ."500".OR.bb2.EQ."200") THEN
198	c
199	CALL histdef(nid_day, "u"//bb2,
200	. "Zonal wind "//bb2//"mb","m/s",
201	. iim,jjmp1,nhori, 1,1,1, -99, 32,
202	. "inst(X)", zout,zout)
203	c
204	CALL histdef(nid_day, "v"//bb2,
205	. "Meridional wind "//bb2//"mb","m/s",
206	. iim,jjmp1,nhori, 1,1,1, -99, 32,
207	. "inst(X)", zout,zout)
208	c
209	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.bb2.EQ."500".OR.bb2.EQ."200") THEN
210	c
211	IF(bb2.EQ."500") THEN
212	c
213	CALL histdef(nid_day, "w"//bb2, "Vertical wind "//bb2//"mb",
214	. "m/s",
215	. iim,jjmp1,nhori, 1,1,1, -99, 32,
216	. "inst(X)", zout,zout)
217	c
218	CALL histdef(nid_day,"phi"//bb2, "Geopotential "//bb2//"mb",
219	. "m",
220	. iim,jjmp1,nhori, 1,1,1, -99, 32,
221	. "inst(X)", zout,zout)
222	c
223	ENDIF !(bb2.EQ."500")
224	c
225	ENDDO !nlevSTD
226	c
227	CALL histdef(nid_day, "slp", "Sea Level Pressure", "Pa",
228	. iim,jjmp1,nhori, 1,1,1, -99, 32,
229	. "ave(X)", zstophy,zout)
230	c
231	CALL histdef(nid_day, "cape_max", "CAPE max.",
232	. "J/kg",
233	. iim,jjmp1,nhori, 1,1,1, -99, 32,
234	. capemaxcels, zstophy,zout)
235	c
236	CALL histdef(nid_day, "solldown", "Down. IR rad. at surface",
237	. "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
238	. "ave(X)", zstophy,zout)
239	c
240	CALL histdef(nid_day, "sens", "Sensible heat flux", "W/m2",
241	. iim,jjmp1,nhori, 1,1,1, -99, 32,
242	. "ave(X)", zstophy,zout)
243	c
244	CALL histdef(nid_day, "SWdnSFC", "SWdn at surface","W/m2",
245	. iim,jjmp1,nhori, 1,1,1,-99,
246	. 32, "ave(X)", zstophy,zout)
247	c
248	cIM diagnostiques flux ocean-atm ou ocean-glace de mer
249	cIM pour utilisation dans un modele de "slab" ocean
250	c
251	IF (OCEAN.EQ.'force ') THEN
252	c
253	CALL histdef(nid_day, "lmt_bils",
254	. "Bilan au sol atmosphere forcee", "W/m2",
255	. iim,jjmp1,nhori, 1,1,1, -99, 32,
256	. "ave(X)", zstophy,zout)
257	c
258	ELSE IF (OCEAN.EQ.'slab ') THEN
259	c
260	CALL histdef(nid_day, "slab_bils",
261	. "Bilan au sol slab", "W/m2",
262	. iim,jjmp1,nhori, 1,1,1, -99, 32,
263	. "ave(X)", zstophy,zout)
264	c
265	CALL histdef(nid_day, "tslab", "Slab SST ", "K",
266	. iim,jjmp1,nhori, 1,1,1, -99, 32,
267	. "ave(X)", zstophy,zout)
268	c
269	CALL histdef(nid_day, "seaice", "Slab seaice", "kg/m2",
270	. iim,jjmp1,nhori, 1,1,1, -99, 32,
271	. "ave(X)", zstophy,zout)
272	c
273	CALL histdef(nid_day, "siceh", "Slab seaice height", "m",
274	. iim,jjmp1,nhori, 1,1,1, -99, 32,
275	. "ave(X)", zstophy,zout)
276	c
277	ENDIF !OCEAN.EQ.force/slab
278	c
279	CALL histdef(nid_day, "fluxo",
280	. "Flux turbulents ocean-atmosphere", "W/m2",
281	. iim,jjmp1,nhori, 1,1,1, -99, 32,
282	. "ave(X)", zstophy,zout)
283	c
284	CALL histdef(nid_day, "fluxg",
285	. "Flux ocean-glace de mer", "W/m2",
286	. iim,jjmp1,nhori, 1,1,1, -99, 32,
287	. "ave(X)", zstophy,zout)
288	c
289	ENDIF !lev_histday.GE.1
290	c
291	IF (lev_histday.GE.2) THEN
292	c
293	CALL histdef(nid_day, "bils", "Surf. total heat flux", "W/m2",
294	. iim,jjmp1,nhori, 1,1,1, -99, 32,
295	. "ave(X)", zstophy,zout)
296	c
297	CALL histdef(nid_day, "fder", "Heat flux derivation", "W/m2/K",
298	. iim,jjmp1,nhori, 1,1,1, -99, 32,
299	. "ave(X)", zstophy,zout)
300	c
301	ENDIF !lev_histday.GE.2
302	c=================================================================
303	IF(lev_histday.GE.3) THEN
304	c=================================================================
305	c INITIALISATION DES CHAMPS 3D
306	c=================================================================
307	c=================================================================
308	c=================================================================
309	c Champs 3D:
310	c
311	CALL histdef(nid_day, "temp", "Air temperature", "K",
312	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
313	. "ave(X)", zstophy,zout)
314	c
315	CALL histdef(nid_day, "ovap", "Specific humidity", "kg/kg",
316	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
317	. "ave(X)", zstophy,zout)
318	c
319	CALL histdef(nid_day, "geop", "Geopotential height", "m2/s2",
320	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
321	. "ave(X)", zstophy,zout)
322	c
323	CALL histdef(nid_day, "vitu", "Zonal wind", "m/s",
324	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
325	. "ave(X)", zstophy,zout)
326	c
327	CALL histdef(nid_day, "vitv", "Meridional wind", "m/s",
328	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
329	. "ave(X)", zstophy,zout)
330	c
331	CALL histdef(nid_day, "vitw", "Vertical wind", "m/s",
332	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
333	. "ave(X)", zstophy,zout)
334	c
335	CALL histdef(nid_day, "pres", "Air pressure", "Pa",
336	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
337	. "ave(X)", zstophy,zout)
338	c
339	c=================================================================
340	c FIN INITIALISATION DES CHAMPS 3D
341	c=================================================================
342	ENDIF ! lev_histday.GE.3
343	c
344	IF (lev_histday.GE.4) THEN
345	c=================================================================
346	c
347	c INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
348	c
349	c=================================================================
350	c
351	CALL histdef(nid_day, "SWupTOA", "SWup at TOA","W/m2",
352	. iim,jjmp1,nhori, 1,1,1,-99,
353	. 32, "ave(X)", zstophy,zout)
354	c
355	CALL histdef(nid_day, "SWupSFC", "SWup at surface","W/m2",
356	. iim,jjmp1,nhori, 1,1,1,-99,
357	. 32, "ave(X)", zstophy,zout)
358	c
359	CALL histdef(nid_day, "SWdnTOA", "SWdn at TOA","W/m2",
360	. iim,jjmp1,nhori, 1,1,1,-99,
361	. 32, "ave(X)", zstophy,zout)
362	c
363	CALL histdef(nid_day, "SWupTOAclr",
364	. "SWup clear sky at TOA","W/m2",
365	. iim,jjmp1,nhori, 1,1,1,-99,
366	. 32, "ave(X)", zstophy,zout)
367
368	CALL histdef(nid_day, "SWupSFCclr",
369	. "SWup clear sky at surface","W/m2",
370	. iim,jjmp1,nhori, 1,1,1,-99,
371	. 32, "ave(X)", zstophy,zout)
372
373	CALL histdef(nid_day, "SWdnTOAclr",
374	. "SWdn clear sky at TOA","W/m2",
375	. iim,jjmp1,nhori, 1,1,1,-99,
376	. 32, "ave(X)", zstophy,zout)
377
378	CALL histdef(nid_day, "SWdnSFCclr",
379	. "SWdn clear sky at surface","W/m2",
380	. iim,jjmp1,nhori, 1,1,1,-99,
381	. 32, "ave(X)", zstophy,zout)
382	c
383	c=================================================================
384	c=================================================================
385	c=================================================================
386	c INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
387	c=================================================================
388	c
389	CALL histdef(nid_day, "tter", "Surface Temperature ter", "K",
390	. iim,jjmp1,nhori, 1,1,1, -99, 32,
391	. "ave(X)", zstophy,zout)
392	c
393	CALL histdef(nid_day, "tlic", "Surface Temperature lic", "K",
394	. iim,jjmp1,nhori, 1,1,1, -99, 32,
395	. "ave(X)", zstophy,zout)
396	c
397	CALL histdef(nid_day, "toce", "Surface Temperature oce", "K",
398	. iim,jjmp1,nhori, 1,1,1, -99, 32,
399	. "ave(X)", zstophy,zout)
400	c
401	CALL histdef(nid_day, "tsic", "Surface Temperature sic", "K",
402	. iim,jjmp1,nhori, 1,1,1, -99, 32,
403	. "ave(X)", zstophy,zout)
404	c
405	CALL histdef(nid_day, "t2mter", "Temp.tere 2m", "K",
406	. iim,jjmp1,nhori, 1,1,1, -99, 32,
407	. "ave(X)", zstophy,zout)
408	c
409	CALL histdef(nid_day, "t2mlic", "Temp.lic 2m", "K",
410	. iim,jjmp1,nhori, 1,1,1, -99, 32,
411	. "ave(X)", zstophy,zout)
412	c
413	CALL histdef(nid_day, "t2moce", "Temp.oce 2m", "K",
414	. iim,jjmp1,nhori, 1,1,1, -99, 32,
415	. "ave(X)", zstophy,zout)
416	c
417	CALL histdef(nid_day, "t2msic", "Temp.sic 2m", "K",
418	. iim,jjmp1,nhori, 1,1,1, -99, 32,
419	. "ave(X)", zstophy,zout)
420	c
421	CALL histdef(nid_day, "t2mter_min", "Temp.terre 2m min.",
422	. "K",
423	. iim,jjmp1,nhori, 1,1,1, -99, 32,
424	. t2mincels, zstophy,zout)
425	c
426	CALL histdef(nid_day, "t2mter_max", "Temp.terre 2m max.",
427	. "K",
428	. iim,jjmp1,nhori, 1,1,1, -99, 32,
429	. t2maxcels, zstophy,zout)
430	c
431	CALL histdef(nid_day, "u10mter", "Vent zonal ter 10m", "m/s",
432	. iim,jjmp1,nhori, 1,1,1, -99, 32,
433	. "ave(X)", zstophy,zout)
434	c
435	CALL histdef(nid_day, "u10mlic", "Vent zonal lic 10m", "m/s",
436	. iim,jjmp1,nhori, 1,1,1, -99, 32,
437	. "ave(X)", zstophy,zout)
438	c
439	CALL histdef(nid_day, "u10moce", "Vent zonal oce 10m", "m/s",
440	. iim,jjmp1,nhori, 1,1,1, -99, 32,
441	. "ave(X)", zstophy,zout)
442	c
443	CALL histdef(nid_day, "u10msic", "Vent zonal sic 10m",
444	. "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
445	. "ave(X)", zstophy,zout)
446	c
447	CALL histdef(nid_day, "v10mter", "Vent meridien ter 10m",
448	. "m/s", iim,jjmp1,nhori, 1,1,1, -99, 32,
449	. "ave(X)", zstophy,zout)
450	c
451	CALL histdef(nid_day, "v10mlic", "Vent meridien lic 10m",
452	. "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
453	. "ave(X)", zstophy,zout)
454	c
455	CALL histdef(nid_day, "v10moce", "Vent meridien oce 10m",
456	. "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
457	. "ave(X)", zstophy,zout)
458	c
459	CALL histdef(nid_day, "v10msic", "Vent meridien sic 10m",
460	. "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
461	. "ave(X)", zstophy,zout)
462	c
463	DO nsrf = 1, nbsrf
464	C
465	call histdef(nid_day, "pourc_"//clnsurf(nsrf),
466	$ "% "//clnsurf(nsrf), "%",
467	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
468	$ "ave(X)", zstophy,zout)
469	C
470	call histdef(nid_day, "fract_"//clnsurf(nsrf),
471	$ "Fraction "//clnsurf(nsrf), "1",
472	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
473	$ "ave(X)", zstophy,zout)
474	C
475	call histdef(nid_day, "tsol_"//clnsurf(nsrf),
476	$ "Temperature "//clnsurf(nsrf), "K",
477	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
478	$ "ave(X)", zstophy,zout)
479	C
480	call histdef(nid_day, "sens_"//clnsurf(nsrf),
481	$ "Sensible heat flux "//clnsurf(nsrf), "W/m2",
482	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
483	$ "ave(X)", zstophy,zout)
484	c
485	call histdef(nid_day, "lat_"//clnsurf(nsrf),
486	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
487	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
488	$ "ave(X)", zstophy,zout)
489	C
490	call histdef(nid_day, "taux_"//clnsurf(nsrf),
491	$ "Zonal wind stress "//clnsurf(nsrf),"Pa",
492	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
493	$ "ave(X)", zstophy,zout)
494
495	call histdef(nid_day, "tauy_"//clnsurf(nsrf),
496	$ "Meridional wind stress "//clnsurf(nsrf), "Pa",
497	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
498	$ "ave(X)", zstophy,zout)
499	C
500	call histdef(nid_day, "albe_"//clnsurf(nsrf),
501	$ "Albedo surf. "//clnsurf(nsrf), "-",
502	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
503	$ "ave(X)", zstophy,zout)
504	C
505	call histdef(nid_day, "rugs_"//clnsurf(nsrf),
506	$ "Rugosite "//clnsurf(nsrf), "m",
507	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
508	$ "ave(X)", zstophy,zout)
509	C
510	END DO
511	C
512	c=================================================================
513	c
514	c FIN INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
515	c
516	c=================================================================
517	ENDIF !lev_histday.GE.4
518	c
519	IF(lev_histday.GE.5) THEN !lev_histday.GE.5
520	c
521	c rajout sorties F. Aires
522	CALL histdef(nid_day, "tops0", "CS Solar rad. at TOA",
523	. "W/m2",
524	. iim,jjmp1,nhori, 1,1,1, -99, 32,
525	. "ave(X)", zstophy,zout)
526	c
527	CALL histdef(nid_day, "topl0", "CS IR rad. at TOA",
528	. "W/m2",
529	. iim,jjmp1,nhori, 1,1,1, -99, 32,
530	. "ave(X)", zstophy,zout)
531	c
532	CALL histdef(nid_day, "sols0", "CS Net Solar rad. at surf.",
533	. "W/m2",
534	. iim,jjmp1,nhori, 1,1,1, -99, 32,
535	. "ave(X)", zstophy,zout)
536	c
537	CALL histdef(nid_day, "soll0", "CS Net IR rad. at surface",
538	. "W/m2",
539	. iim,jjmp1,nhori, 1,1,1, -99, 32,
540	. "ave(X)", zstophy,zout)
541	c
542	CALL histdef(nid_day, "rneb", "Cloud fraction", "-",
543	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
544	. "ave(X)", zstophy,zout)
545	c
546	CALL histdef(nid_day, "rnebcon", "Convective Cloud Fraction"
547	. , "-",
548	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
549	. "ave(X)", zstophy,zout)
550	c
551	CALL histdef(nid_day,"lwcon","Cloud liquid water content",
552	. "kg/kg",
553	. iim,jjmp1,nhori, klev,1,klev, nvert, 32,
554	. "ave(X)", zstophy,zout)
555	c
556	CALL histdef(nid_day,"iwcon","Cloud ice water content","kg/kg",
557	. iim,jjmp1,nhori, klev,1,klev, nvert, 32,
558	. "ave(X)", zstophy,zout)
559	c
560	CALL histdef(nid_day,"lwp","Cloud water path","kg/m2",
561	. iim,jjmp1,nhori, 1,1,1, -99, 32,
562	. "ave(X)", zstophy,zout)
563	c
564	CALL histdef(nid_day,"iwp","Cloud ice water path","kg/m2",
565	. iim,jjmp1,nhori, 1,1,1, -99, 32,
566	. "ave(X)", zstophy,zout)
567	c
568	CALL histdef(nid_day,"meantaucld",
569	. "ISCCP mean cloud optical thickness","1",
570	. iim,jjmp1,nhori, 1,1,1, -99, 32,
571	. "ave(X)", zstophy,zout)
572	c
573	CALL histdef(nid_day,"cldtau","Cloud optical thickness","1",
574	. iim,jjmp1,nhori, klev,1,klev, nvert, 32,
575	. "ave(X)", zstophy,zout)
576	c
577	CALL histdef(nid_day,"cldemi","Cloud optical emissivity","1",
578	. iim,jjmp1,nhori, klev,1,klev, nvert, 32,
579	. "ave(X)", zstophy,zout)
580	c
581	ENDIF !lev_histday.GE.5
582	c=================================================================
583	c
584	CALL histend(nid_day)
585	c
586	ndex2d = 0
587	ndex3d = 0
588	c
589	c=================================================================
590	ENDIF ! fin de test sur ok_journe

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

Download in other formats: