Context Navigation

source: LMDZ4/branches/pre_V3/libf/phylmd/ini_histday.h @ 4249

Last change on this file since 4249 was 684, checked in by lmdzadmin, 19 years ago
Ajout nouveaux flags (ip_ebil_phy, ok_slab_sicOBS) ds physiq.def + Lecture freq. ecriture en nombre de jours IM
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 24.5 KB

Line
1	c
2	c $Header$
3	c
4	IF (ok_journe) THEN
5	c
6	zstophy = dtime
7	zstoday= ecrit_day
8	zout = 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, "radsol", "Rayonnement au sol", "W/m2",
159	. iim,jjmp1,nhori, 1,1,1, -99, 32,
160	. "ave(X)", zstophy,zout)
161	c
162	CALL histdef(nid_day, "cldl", "Low-level cloudiness", "-",
163	. iim,jjmp1,nhori, 1,1,1, -99, 32,
164	. "ave(X)", zstophy,zout)
165	c
166	CALL histdef(nid_day, "cldm", "Mid-level cloudiness", "-",
167	. iim,jjmp1,nhori, 1,1,1, -99, 32,
168	. "ave(X)", zstophy,zout)
169	c
170	CALL histdef(nid_day, "cldh", "High-level cloudiness", "-",
171	. iim,jjmp1,nhori, 1,1,1, -99, 32,
172	. "ave(X)", zstophy,zout)
173	c
174	CALL histdef(nid_day, "cldt", "Total cloudiness", "-",
175	. iim,jjmp1,nhori, 1,1,1, -99, 32,
176	. "ave(X)", zstophy,zout)
177	c
178	CALL histdef(nid_day, "cldq", "Cloud liquid water path",
179	. "kg/m2",
180	. iim,jjmp1,nhori, 1,1,1, -99, 32,
181	. "ave(X)", zstophy,zout)
182	c
183	CALL histdef(nid_day, "prw", "Precipitable water", "kg/m2",
184	. iim,jjmp1,nhori, 1,1,1, -99, 32,
185	. "ave(X)", zstophy,zout)
186	c
187	c Champs dynamiques sur niveaux de pression
188	c
189	c DO k=1, nlevENS
190	DO k=1, nlevSTD
191	c
192	c bb=clevSTD(k)
193	c
194	IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
195	IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
196	c aa=clevSTD(k)
197	c bb2=clevSTD(k)
198	c ENDIF
199	c
200	IF(bb2.EQ."850".OR.bb2.EQ."700"
201	. .OR.bb2.EQ."500".OR.bb2.EQ."200") THEN
202	c
203	CALL histdef(nid_day, "u"//bb2,
204	. "Zonal wind "//bb2//"mb","m/s",
205	. iim,jjmp1,nhori, 1,1,1, -99, 32,
206	. "inst(X)", zout,zout)
207	c
208	CALL histdef(nid_day, "v"//bb2,
209	. "Meridional wind "//bb2//"mb","m/s",
210	. iim,jjmp1,nhori, 1,1,1, -99, 32,
211	. "inst(X)", zout,zout)
212	c
213	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.bb2.EQ."500".OR.bb2.EQ."200") THEN
214	c
215	IF(bb2.EQ."500") THEN
216	c
217	CALL histdef(nid_day, "w"//bb2, "Vertical wind "//bb2//"mb",
218	. "m/s",
219	. iim,jjmp1,nhori, 1,1,1, -99, 32,
220	. "inst(X)", zout,zout)
221	c
222	CALL histdef(nid_day,"phi"//bb2, "Geopotential "//bb2//"mb",
223	. "m",
224	. iim,jjmp1,nhori, 1,1,1, -99, 32,
225	. "inst(X)", zout,zout)
226	c
227	ENDIF !(bb2.EQ."500")
228	c
229	ENDDO !nlevSTD
230	c
231	CALL histdef(nid_day, "slp", "Sea Level Pressure", "Pa",
232	. iim,jjmp1,nhori, 1,1,1, -99, 32,
233	. "ave(X)", zstophy,zout)
234	c
235	CALL histdef(nid_day, "cape_max", "CAPE max.",
236	. "J/kg",
237	. iim,jjmp1,nhori, 1,1,1, -99, 32,
238	. capemaxcels, zstophy,zout)
239	c
240	CALL histdef(nid_day, "solldown", "Down. IR rad. at surface",
241	. "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
242	. "ave(X)", zstophy,zout)
243	c
244	CALL histdef(nid_day, "sens", "Sensible heat flux", "W/m2",
245	. iim,jjmp1,nhori, 1,1,1, -99, 32,
246	. "ave(X)", zstophy,zout)
247	c
248	CALL histdef(nid_day, "SWdnSFC", "SWdn at surface","W/m2",
249	. iim,jjmp1,nhori, 1,1,1,-99,
250	. 32, "ave(X)", zstophy,zout)
251	c
252	cIM diagnostiques flux ocean-atm ou ocean-glace de mer
253	cIM pour utilisation dans un modele de "slab" ocean
254	c
255	IF (OCEAN.EQ.'force ') THEN
256	c
257	CALL histdef(nid_day, "lmt_bils",
258	. "Bilan au sol atmosphere forcee", "W/m2",
259	. iim,jjmp1,nhori, 1,1,1, -99, 32,
260	. "ave(X)", zstophy,zout)
261	c
262	ELSE IF (OCEAN.EQ.'slab ') THEN
263	c
264	CALL histdef(nid_day, "slab_bils",
265	. "Bilan au sol slab", "W/m2",
266	. iim,jjmp1,nhori, 1,1,1, -99, 32,
267	. "ave(X)", zstophy,zout)
268	c
269	CALL histdef(nid_day, "tslab", "Slab SST ", "K",
270	. iim,jjmp1,nhori, 1,1,1, -99, 32,
271	. "ave(X)", zstophy,zout)
272	c
273	CALL histdef(nid_day, "seaice", "Slab seaice", "kg/m2",
274	. iim,jjmp1,nhori, 1,1,1, -99, 32,
275	. "ave(X)", zstophy,zout)
276	c
277	CALL histdef(nid_day, "siceh", "Slab seaice height", "m",
278	. iim,jjmp1,nhori, 1,1,1, -99, 32,
279	. "ave(X)", zstophy,zout)
280	c
281	ENDIF !OCEAN.EQ.force/slab
282	c
283	CALL histdef(nid_day, "fluxo",
284	. "Flux turbulents ocean-atmosphere", "W/m2",
285	. iim,jjmp1,nhori, 1,1,1, -99, 32,
286	. "ave(X)", zstophy,zout)
287	c
288	CALL histdef(nid_day, "fluxg",
289	. "Flux ocean-glace de mer", "W/m2",
290	. iim,jjmp1,nhori, 1,1,1, -99, 32,
291	. "ave(X)", zstophy,zout)
292	c
293	ENDIF !lev_histday.GE.1
294	c
295	IF (lev_histday.GE.2) THEN
296	c
297	CALL histdef(nid_day, "bils", "Surf. total heat flux", "W/m2",
298	. iim,jjmp1,nhori, 1,1,1, -99, 32,
299	. "ave(X)", zstophy,zout)
300	c
301	CALL histdef(nid_day, "fder", "Heat flux derivation", "W/m2/K",
302	. iim,jjmp1,nhori, 1,1,1, -99, 32,
303	. "ave(X)", zstophy,zout)
304	c
305	ENDIF !lev_histday.GE.2
306	c=================================================================
307	IF(lev_histday.GE.3) THEN
308	c=================================================================
309	c INITIALISATION DES CHAMPS 3D
310	c=================================================================
311	c=================================================================
312	c=================================================================
313	c Champs 3D:
314	c
315	CALL histdef(nid_day, "temp", "Air temperature", "K",
316	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
317	. "ave(X)", zstophy,zout)
318	c
319	CALL histdef(nid_day, "ovap", "Specific humidity", "kg/kg",
320	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
321	. "ave(X)", zstophy,zout)
322	c
323	CALL histdef(nid_day, "geop", "Geopotential height", "m2/s2",
324	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
325	. "ave(X)", zstophy,zout)
326	c
327	CALL histdef(nid_day, "vitu", "Zonal wind", "m/s",
328	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
329	. "ave(X)", zstophy,zout)
330	c
331	CALL histdef(nid_day, "vitv", "Meridional wind", "m/s",
332	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
333	. "ave(X)", zstophy,zout)
334	c
335	CALL histdef(nid_day, "vitw", "Vertical wind", "m/s",
336	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
337	. "ave(X)", zstophy,zout)
338	c
339	CALL histdef(nid_day, "pres", "Air pressure", "Pa",
340	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
341	. "ave(X)", zstophy,zout)
342	c
343	c=================================================================
344	c FIN INITIALISATION DES CHAMPS 3D
345	c=================================================================
346	ENDIF ! lev_histday.GE.3
347	c
348	IF (lev_histday.GE.4) THEN
349	c=================================================================
350	c
351	c INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
352	c
353	c=================================================================
354	c
355	CALL histdef(nid_day, "SWupTOA", "SWup at TOA","W/m2",
356	. iim,jjmp1,nhori, 1,1,1,-99,
357	. 32, "ave(X)", zstophy,zout)
358	c
359	CALL histdef(nid_day, "SWupSFC", "SWup at surface","W/m2",
360	. iim,jjmp1,nhori, 1,1,1,-99,
361	. 32, "ave(X)", zstophy,zout)
362	c
363	CALL histdef(nid_day, "SWdnTOA", "SWdn at TOA","W/m2",
364	. iim,jjmp1,nhori, 1,1,1,-99,
365	. 32, "ave(X)", zstophy,zout)
366	c
367	CALL histdef(nid_day, "SWupTOAclr",
368	. "SWup clear sky at TOA","W/m2",
369	. iim,jjmp1,nhori, 1,1,1,-99,
370	. 32, "ave(X)", zstophy,zout)
371
372	CALL histdef(nid_day, "SWupSFCclr",
373	. "SWup clear sky at surface","W/m2",
374	. iim,jjmp1,nhori, 1,1,1,-99,
375	. 32, "ave(X)", zstophy,zout)
376
377	CALL histdef(nid_day, "SWdnTOAclr",
378	. "SWdn clear sky at TOA","W/m2",
379	. iim,jjmp1,nhori, 1,1,1,-99,
380	. 32, "ave(X)", zstophy,zout)
381
382	CALL histdef(nid_day, "SWdnSFCclr",
383	. "SWdn clear sky at surface","W/m2",
384	. iim,jjmp1,nhori, 1,1,1,-99,
385	. 32, "ave(X)", zstophy,zout)
386	c
387	cIM AMMA-MIP
388	CALL histdef(nid_day, "LWdnSFC",
389	. "LW down at surface","W/m2",
390	. iim,jjmp1,nhori, 1,1,1,-99,
391	. 32, "ave(X)", zstophy,zout)
392	c
393	CALL histdef(nid_day, "LWupSFC",
394	. "LW down at surface","W/m2",
395	. iim,jjmp1,nhori, 1,1,1,-99,
396	. 32, "ave(X)", zstophy,zout)
397	c
398	cIM IF (iflag_con.GE.3) THEN
399	CALL histdef(nid_day, "ptop",
400	c . "Cld top pressure","mb",
401	. "Cld top pressure (cf ISCCP simulator)","Pa",
402	. iim,jjmp1,nhori, 1,1,1,-99,
403	. 32, "ave(X)", zstophy,zout)
404	cIM ENDIF !(iflag_con.GE.3) THEN
405	c
406	c=================================================================
407	c=================================================================
408	c=================================================================
409	c INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
410	c=================================================================
411	c
412	CALL histdef(nid_day, "tter", "Surface Temperature ter", "K",
413	. iim,jjmp1,nhori, 1,1,1, -99, 32,
414	. "ave(X)", zstophy,zout)
415	c
416	CALL histdef(nid_day, "tlic", "Surface Temperature lic", "K",
417	. iim,jjmp1,nhori, 1,1,1, -99, 32,
418	. "ave(X)", zstophy,zout)
419	c
420	CALL histdef(nid_day, "toce", "Surface Temperature oce", "K",
421	. iim,jjmp1,nhori, 1,1,1, -99, 32,
422	. "ave(X)", zstophy,zout)
423	c
424	CALL histdef(nid_day, "tsic", "Surface Temperature sic", "K",
425	. iim,jjmp1,nhori, 1,1,1, -99, 32,
426	. "ave(X)", zstophy,zout)
427	c
428	CALL histdef(nid_day, "t2mter", "Temp.tere 2m", "K",
429	. iim,jjmp1,nhori, 1,1,1, -99, 32,
430	. "ave(X)", zstophy,zout)
431	c
432	CALL histdef(nid_day, "t2mlic", "Temp.lic 2m", "K",
433	. iim,jjmp1,nhori, 1,1,1, -99, 32,
434	. "ave(X)", zstophy,zout)
435	c
436	CALL histdef(nid_day, "t2moce", "Temp.oce 2m", "K",
437	. iim,jjmp1,nhori, 1,1,1, -99, 32,
438	. "ave(X)", zstophy,zout)
439	c
440	CALL histdef(nid_day, "t2msic", "Temp.sic 2m", "K",
441	. iim,jjmp1,nhori, 1,1,1, -99, 32,
442	. "ave(X)", zstophy,zout)
443	c
444	CALL histdef(nid_day, "t2mter_min", "Temp.terre 2m min.",
445	. "K",
446	. iim,jjmp1,nhori, 1,1,1, -99, 32,
447	. t2mincels, zstophy,zout)
448	c
449	CALL histdef(nid_day, "t2mter_max", "Temp.terre 2m max.",
450	. "K",
451	. iim,jjmp1,nhori, 1,1,1, -99, 32,
452	. t2maxcels, zstophy,zout)
453	c
454	CALL histdef(nid_day, "u10mter", "Vent zonal ter 10m", "m/s",
455	. iim,jjmp1,nhori, 1,1,1, -99, 32,
456	. "ave(X)", zstophy,zout)
457	c
458	CALL histdef(nid_day, "u10mlic", "Vent zonal lic 10m", "m/s",
459	. iim,jjmp1,nhori, 1,1,1, -99, 32,
460	. "ave(X)", zstophy,zout)
461	c
462	CALL histdef(nid_day, "u10moce", "Vent zonal oce 10m", "m/s",
463	. iim,jjmp1,nhori, 1,1,1, -99, 32,
464	. "ave(X)", zstophy,zout)
465	c
466	CALL histdef(nid_day, "u10msic", "Vent zonal sic 10m",
467	. "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
468	. "ave(X)", zstophy,zout)
469	c
470	CALL histdef(nid_day, "v10mter", "Vent meridien ter 10m",
471	. "m/s", iim,jjmp1,nhori, 1,1,1, -99, 32,
472	. "ave(X)", zstophy,zout)
473	c
474	CALL histdef(nid_day, "v10mlic", "Vent meridien lic 10m",
475	. "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
476	. "ave(X)", zstophy,zout)
477	c
478	CALL histdef(nid_day, "v10moce", "Vent meridien oce 10m",
479	. "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
480	. "ave(X)", zstophy,zout)
481	c
482	CALL histdef(nid_day, "v10msic", "Vent meridien sic 10m",
483	. "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
484	. "ave(X)", zstophy,zout)
485	c
486	DO nsrf = 1, nbsrf
487	C
488	call histdef(nid_day, "pourc_"//clnsurf(nsrf),
489	$ "% "//clnsurf(nsrf), "%",
490	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
491	$ "ave(X)", zstophy,zout)
492	C
493	call histdef(nid_day, "fract_"//clnsurf(nsrf),
494	$ "Fraction "//clnsurf(nsrf), "1",
495	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
496	$ "ave(X)", zstophy,zout)
497	C
498	call histdef(nid_day, "tsol_"//clnsurf(nsrf),
499	$ "Temperature "//clnsurf(nsrf), "K",
500	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
501	$ "ave(X)", zstophy,zout)
502	C
503	call histdef(nid_day, "sens_"//clnsurf(nsrf),
504	$ "Sensible heat flux "//clnsurf(nsrf), "W/m2",
505	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
506	$ "ave(X)", zstophy,zout)
507	c
508	call histdef(nid_day, "lat_"//clnsurf(nsrf),
509	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
510	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
511	$ "ave(X)", zstophy,zout)
512	C
513	call histdef(nid_day, "taux_"//clnsurf(nsrf),
514	$ "Zonal wind stress "//clnsurf(nsrf),"Pa",
515	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
516	$ "ave(X)", zstophy,zout)
517
518	call histdef(nid_day, "tauy_"//clnsurf(nsrf),
519	$ "Meridional wind stress "//clnsurf(nsrf), "Pa",
520	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
521	$ "ave(X)", zstophy,zout)
522	C
523	call histdef(nid_day, "albe_"//clnsurf(nsrf),
524	$ "Albedo surf. "//clnsurf(nsrf), "-",
525	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
526	$ "ave(X)", zstophy,zout)
527	C
528	call histdef(nid_day, "rugs_"//clnsurf(nsrf),
529	$ "Rugosite "//clnsurf(nsrf), "m",
530	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
531	$ "ave(X)", zstophy,zout)
532	C
533	END DO
534	C
535	c=================================================================
536	c
537	c FIN INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
538	c
539	c=================================================================
540	ENDIF !lev_histday.GE.4
541	c
542	IF(lev_histday.GE.5) THEN !lev_histday.GE.5
543	c
544	c rajout sorties F. Aires
545	CALL histdef(nid_day, "tops0", "CS Solar rad. at TOA",
546	. "W/m2",
547	. iim,jjmp1,nhori, 1,1,1, -99, 32,
548	. "ave(X)", zstophy,zout)
549	c
550	CALL histdef(nid_day, "topl0", "CS IR rad. at TOA",
551	. "W/m2",
552	. iim,jjmp1,nhori, 1,1,1, -99, 32,
553	. "ave(X)", zstophy,zout)
554	c
555	CALL histdef(nid_day, "sols0", "CS Net Solar rad. at surf.",
556	. "W/m2",
557	. iim,jjmp1,nhori, 1,1,1, -99, 32,
558	. "ave(X)", zstophy,zout)
559	c
560	CALL histdef(nid_day, "soll0", "CS Net IR rad. at surface",
561	. "W/m2",
562	. iim,jjmp1,nhori, 1,1,1, -99, 32,
563	. "ave(X)", zstophy,zout)
564	c
565	CALL histdef(nid_day, "rneb", "Cloud fraction", "-",
566	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
567	. "ave(X)", zstophy,zout)
568	c
569	CALL histdef(nid_day, "rnebcon", "Convective Cloud Fraction"
570	. , "-",
571	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
572	. "ave(X)", zstophy,zout)
573	c
574	CALL histdef(nid_day,"lwcon","Cloud liquid water content",
575	. "kg/kg",
576	. iim,jjmp1,nhori, klev,1,klev, nvert, 32,
577	. "ave(X)", zstophy,zout)
578	c
579	CALL histdef(nid_day,"iwcon","Cloud ice water content","kg/kg",
580	. iim,jjmp1,nhori, klev,1,klev, nvert, 32,
581	. "ave(X)", zstophy,zout)
582	c
583	CALL histdef(nid_day,"lwp","Cloud water path","kg/m2",
584	. iim,jjmp1,nhori, 1,1,1, -99, 32,
585	. "ave(X)", zstophy,zout)
586	c
587	CALL histdef(nid_day,"iwp","Cloud ice water path","kg/m2",
588	. iim,jjmp1,nhori, 1,1,1, -99, 32,
589	. "ave(X)", zstophy,zout)
590	c
591	CALL histdef(nid_day,"meantaucld",
592	. "ISCCP mean cloud optical thickness","1",
593	. iim,jjmp1,nhori, 1,1,1, -99, 32,
594	. "ave(X)", zstophy,zout)
595	c
596	CALL histdef(nid_day,"cldtau","Cloud optical thickness","1",
597	. iim,jjmp1,nhori, klev,1,klev, nvert, 32,
598	. "ave(X)", zstophy,zout)
599	c
600	CALL histdef(nid_day,"cldemi","Cloud optical emissivity","1",
601	. iim,jjmp1,nhori, klev,1,klev, nvert, 32,
602	. "ave(X)", zstophy,zout)
603	c
604	cIM rajout AMMA-MIP
605	c
606	CALL histdef(nid_day, "dtcon", "Convection dT", "K/s",
607	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
608	. "ave(X)", zstophy,zout)
609	c
610	CALL histdef(nid_day, "dqcon", "Convection dQ", "(kg/kg)/s",
611	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
612	. "ave(X)", zstophy,zout)
613	c
614	CALL histdef(nid_day, "rh2m", "Relative humidity at 2m", "%",
615	. iim,jjmp1,nhori, 1,1,1, -99, 32,
616	. "ave(X)", zstophy,zout)
617	c
618	CALL histdef(nid_day, "qsat2m", "Saturant humidity at 2m", "%",
619	. iim,jjmp1,nhori, 1,1,1, -99, 32,
620	. "ave(X)", zstophy,zout)
621	c
622	CALL histdef(nid_day, "tpot",
623	. "Surface air potential temperature", "K",
624	. iim,jjmp1,nhori, 1,1,1, -99, 32,
625	. "ave(X)", zstophy,zout)
626	c
627	CALL histdef(nid_day, "tpote",
628	. "Surface air equivalent potential temperature", "K",
629	. iim,jjmp1,nhori, 1,1,1, -99, 32,
630	. "ave(X)", zstophy,zout)
631	c
632	ENDIF !lev_histday.GE.5
633	c=================================================================
634	c
635	CALL histend(nid_day)
636	c
637	ndex2d = 0
638	ndex3d = 0
639	c
640	c=================================================================
641	ENDIF ! fin de test sur ok_journe

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

Download in other formats: