Context Navigation

source: trunk/LMDZ.TITAN/libf/phytitan/ini_histday.h @ 1243

Last change on this file since 1243 was 1056, checked in by slebonnois, 11 years ago
SL: Titan runs ! see DOC/chantiers/commit_importants.log
File size: 13.4 KB

Line
1	IF (ok_journe) THEN
2
3	zsto = dtime
4	zout = dtime * REAL(ecrit_day)
5	c zsto1: pour des flux radiatifs calcules tous les radpas appels physiq
6	zsto1= dtime * REAL(radpas)
7
8	idayref = day_ref
9	CALL ymds2ju(annee_ref, 1, idayref, zero, zjulian)
10
11	CALL histbeg_phy("histday.nc", itau_phy, zjulian, dtime,
12	. nhori, nid_day)
13
14	!$OMP MASTER
15	CALL histvert(nid_day, "presnivs", "Vertical levels", "Pa",
16	. klev, presnivs, nvert)
17
18	c-------------------------------------------------------
19	IF(lev_histday.GE.1) THEN
20
21	ccccccccccccc 2D fields, invariables
22
23	CALL histdef(nid_day, "phis", "Surface geop. height", "-",
24	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
25	. "once", zsto,zout)
26
27	CALL histdef(nid_day, "aire", "Grid area", "-",
28	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
29	. "once", zsto,zout)
30
31	ccccccc axe Ls
32	CALL histdef(nid_day, "ls", "Solar longitude", "degrees",
33	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
34	. "ave(X)", zsto,zout)
35
36	ccccccccccccc 2D fields, variables
37
38	CALL histdef(nid_day, "tsol", "Surface Temperature", "K",
39	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
40	. "ave(X)", zsto,zout)
41
42	CALL histdef(nid_day, "psol", "Surface Pressure", "Pa",
43	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
44	. "ave(X)", zsto,zout)
45
46	c CALL histdef(nid_day, "ue", "Zonal energy transport", "-",
47	c . iim,jj_nb,nhori, 1,1,1, nvert, 32,
48	c . "ave(X)", zsto,zout)
49
50	c CALL histdef(nid_day, "ve", "Merid energy transport", "-",
51	c . iim,jj_nb,nhori, 1,1,1, nvert, 32,
52	c . "ave(X)", zsto,zout)
53
54	ENDIF !lev_histday.GE.1
55
56	c-------------------------------------------------------
57	IF(lev_histday.GE.2) THEN
58
59	ccccccccccccc 3D fields, basics
60
61	CALL histdef(nid_day, "temp", "Air temperature", "K",
62	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
63	. "ave(X)", zsto,zout)
64
65	CALL histdef(nid_day, "pres", "Air pressure", "Pa",
66	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
67	. "ave(X)", zsto,zout)
68
69	CALL histdef(nid_day, "geop", "Geopotential height", "m",
70	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
71	. "ave(X)", zsto,zout)
72
73	CALL histdef(nid_day, "vitu", "Zonal wind", "m/s",
74	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
75	. "ave(X)", zsto,zout)
76
77	CALL histdef(nid_day, "vitv", "Meridional wind", "m/s",
78	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
79	. "ave(X)", zsto,zout)
80
81	CALL histdef(nid_day, "vitw", "Vertical wind", "Pa/s",
82	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
83	. "ave(X)", zsto,zout)
84
85	CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",
86	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
87	. "ave(X)", zsto1,zout)
88
89	CALL histdef(nid_day, "duvdf", "Boundary-layer dU", "m/s2",
90	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
91	. "ave(X)", zsto,zout)
92
93	CALL histdef(nid_day, "dudyn", "Dynamics dU", "m/s2",
94	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
95	. "ave(X)", zsto,zout)
96
97	cccccccccccccccccc Tracers
98
99	if (iflag_trac.eq.1) THEN
100	if (microfi.ge.1) then
101	c DO iq=1,nmicro
102	c CALL histdef(nid_day, tname(iq), ttext(iq), "n/m2",
103	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
104	c . "ave(X)", zsto,zout)
105	c ENDDO
106	CALL histdef(nid_day, "qaer","nb tot aer" , "n/m2",
107	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
108	. "ave(X)", zsto,zout)
109
110	if (clouds.eq.1) then
111	CALL histdef(nid_day, "qnoy","nb tot noy" , "n/m2",
112	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
113	. "ave(X)", zsto,zout)
114	CALL histdef(nid_day, "qgl1","V tot gl1" , "m3/m2",
115	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
116	. "ave(X)", zsto,zout)
117	CALL histdef(nid_day, "qgl2","V tot gl2" , "m3/m2",
118	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
119	. "ave(X)", zsto,zout)
120	CALL histdef(nid_day, "qgl3","V tot gl3" , "m3/m2",
121	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
122	. "ave(X)", zsto,zout)
123	c--------------
124	c ----- SATURATION ESP NUAGES
125	CALL histdef(nid_day,"ch4sat", "saturation CH4", "--",
126	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
127	. "ave(X)", zsto,zout)
128	CALL histdef(nid_day,"c2h6sat", "saturation C2H6", "--",
129	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
130	. "ave(X)", zsto,zout)
131	CALL histdef(nid_day,"c2h2sat", "saturation C2H2", "--",
132	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
133	. "ave(X)", zsto,zout)
134	c --------------
135	c ----- RESERVOIR DE SURFACE
136	CALL histdef(nid_day, "reserv", "Reservoir surface","m",
137	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
138	. "ave(X)", zsto,zout)
139	c --------------
140	c ----- ECHANGE GAZ SURF/ATM (evaporation)
141	CALL histdef(nid_day, "evapch4", "Evaporation CH4","m",
142	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
143	. "ave(X)", zsto,zout)
144	c --------------
145	c ----- PRECIPITATIONS (precipitations cumulatives)
146	CALL histdef(nid_day,"prech4","Precip CH4","m",
147	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
148	. "ave(X)", zsto,zout)
149	CALL histdef(nid_day,"prec2h6","Precip C2H6",
150	. "m",iim,jj_nb,nhori, 1,1,1, nvert, 32,
151	. "ave(X)", zsto,zout)
152	CALL histdef(nid_day,"prec2h2","Precip C2H2",
153	. "m",iim,jj_nb,nhori, 1,1,1, nvert, 32,
154	. "ave(X)", zsto,zout)
155	CALL histdef(nid_day,"prenoy","Precip NOY",
156	. "um/s",iim,jj_nb,nhori, 1,1,1, nvert, 32,
157	. "ave(X)", zsto,zout)
158	CALL histdef(nid_day,"preaer","Precip AER",
159	. "um/s",iim,jj_nb,nhori, 1,1,1, nvert, 32,
160	. "ave(X)", zsto,zout)
161	c --------------
162	c ----- FLUX GLACE
163	CALL histdef(nid_day,"flxgl1", "flux gl CH4",
164	. "kg/m2/s",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
165	. "ave(X)", zsto,zout)
166	CALL histdef(nid_day,"flxgl2", "flux gl C2H6",
167	. "kg/m2/s",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
168	. "ave(X)", zsto,zout)
169	CALL histdef(nid_day,"flxgl3", "flux gl C2H2",
170	. "kg/m2/s",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
171	. "ave(X)", zsto,zout)
172	c --------------
173	c ----- RAYON DES GOUTTES
174	CALL histdef(nid_day,"rcldbar", "rayon moyen goutte",
175	. "m",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
176	. "ave(X)", zsto,zout)
177	endif
178	endif
179	c --------------
180	c ----- TRACEURS CHIMIQUES
181	if (nmicro.lt.nqmax) then
182	DO iq=nmicro+1,nqmax
183	CALL histdef(nid_day, tname(iq), ttext(iq), "ppm",
184	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
185	. "ave(X)", zsto,zout)
186	ENDDO
187	endif
188	endif
189
190	ENDIF !lev_histday.GE.2
191
192	c-------------------------------------------------------
193	IF(lev_histday.GE.3) THEN
194
195	cccccccccccccccccc Radiative transfer
196
197	c 2D
198
199	CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2",
200	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
201	. "ave(X)", zsto1,zout)
202
203	CALL histdef(nid_day, "sols", "Solar rad. at surf.", "W/m2",
204	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
205	. "ave(X)", zsto1,zout)
206
207	CALL histdef(nid_day, "soll", "IR rad. at surface", "W/m2",
208	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
209	. "ave(X)", zsto1,zout)
210
211	c 3D
212
213	CALL histdef(nid_day, "SWnet", "Net SW flux","W/m2",
214	. iim,jj_nb,nhori, klev,1,klev,nvert,
215	. 32, "ave(X)", zsto1,zout)
216
217	CALL histdef(nid_day, "LWnet", "Net LW flux","W/m2",
218	. iim,jj_nb,nhori, klev,1,klev,nvert,
219	. 32, "ave(X)", zsto1,zout)
220
221	c --------------
222	c ----- OPACITE BRUME
223	DO k=7,NSPECV,10
224	write(str2,'(i2.2)') k
225	CALL histdef(nid_day,"thv"//str2,"Haze Opa Vis",
226	. "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
227	. "ave(X)",zsto1,zout)
228	ENDDO
229
230	DO k=8,NSPECI,10
231	write(str2,'(i2.2)') k
232	CALL histdef(nid_day,"thi"//str2,"Haze Opa IR",
233	. "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
234	. "ave(X)",zsto1,zout)
235	ENDDO
236
237	c --------------
238	c ----- EXTINCTION BRUME
239	DO k=7,NSPECV,10
240	write(str2,'(i2.2)') k
241	CALL histdef(nid_day,"khv"//str2,"Haze ext Vis ",
242	. "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
243	. "ave(X)",zsto1,zout)
244	ENDDO
245
246	DO k=8,NSPECI,10
247	write(str2,'(i2.2)') k
248	CALL histdef(nid_day,"khi"//str2,"Haze ext IR ",
249	. "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
250	. "ave(X)",zsto1,zout)
251	ENDDO
252
253	c --------------
254	c ----- OPACITE GAZ
255	DO k=7,NSPECV,10
256	write(str2,'(i2.2)') k
257	CALL histdef(nid_day,"tgv"//str2,"Gas Opa Vis",
258	. "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
259	. "ave(X)",zsto1,zout)
260	ENDDO
261
262	DO k=8,NSPECI,10
263	write(str2,'(i2.2)') k
264	CALL histdef(nid_day,"tgi"//str2,"Gas Opa IR",
265	. "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
266	. "ave(X)",zsto1,zout)
267	ENDDO
268
269	c --------------
270	c ----- EXTINCTION GAZ
271	DO k=7,NSPECV,10
272	write(str2,'(i2.2)') k
273	CALL histdef(nid_day,"kgv"//str2,"Gas ext Vis ",
274	. "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
275	. "ave(X)",zsto1,zout)
276	ENDDO
277
278	DO k=8,NSPECI,10
279	write(str2,'(i2.2)') k
280	CALL histdef(nid_day,"kgi"//str2,"Gas ext IR ",
281	. "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
282	. "ave(X)",zsto1,zout)
283	ENDDO
284
285	c --------------
286	c ----- OPACITE NUAGES
287	if (clouds.eq.1) then
288	CALL histdef(nid_day,"tcld","Cld Opa proxy",
289	. "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
290	. "ave(X)",zsto,zout)
291
292	c --------------
293	c ----- EXTINCTION NUAGES
294	CALL histdef(nid_day,"kcld","Cld Ext proxy",
295	. "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
296	. "ave(X)",zsto,zout)
297	endif
298
299	ENDIF !lev_histday.GE.3
300
301	c-------------------------------------------------------
302	IF(lev_histday.GE.4) THEN
303
304	CALL histdef(nid_day, "dtdyn", "Dynamics dT", "K/s",
305	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
306	. "ave(X)", zsto,zout)
307
308	CALL histdef(nid_day, "dtphy", "Physics dT", "K/s",
309	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
310	. "ave(X)", zsto,zout)
311
312	CALL histdef(nid_day, "dtvdf", "Boundary-layer dT", "K/s",
313	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
314	. "ave(X)", zsto,zout)
315
316	CALL histdef(nid_day, "dtajs", "Dry adjust. dT", "K/s",
317	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
318	. "ave(X)", zsto,zout)
319
320	CALL histdef(nid_day, "dtswr", "SW radiation dT", "K/s",
321	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
322	. "ave(X)", zsto,zout)
323
324	CALL histdef(nid_day, "dtlwr", "LW radiation dT", "K/s",
325	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
326	. "ave(X)", zsto,zout)
327
328	c CALL histdef(nid_day, "dtec", "Cinetic dissip dT", "K/s",
329	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
330	c . "ave(X)", zsto,zout)
331
332	ENDIF !lev_histday.GE.4
333
334	c-------------------------------------------------------
335	IF(lev_histday.GE.5) THEN
336
337
338	c call histdef(nid_day, "taux",
339	c $ "Zonal wind stress", "Pa",
340	c $ iim,jj_nb,nhori, 1,1,1, nvert, 32,
341	c $ "ave(X)", zsto,zout)
342
343	c call histdef(nid_day, "tauy",
344	c $ "Meridional xind stress", "Pa",
345	c $ iim,jj_nb,nhori, 1,1,1, nvert, 32,
346	c $ "ave(X)", zsto,zout)
347
348	c CALL histdef(nid_day, "cdrm", "Momentum drag coef.", "-",
349	c . iim,jj_nb,nhori, 1,1,1, nvert, 32,
350	c . "ave(X)", zsto,zout)
351
352	c CALL histdef(nid_day, "cdrh", "Heat drag coef.", "-",
353	c . iim,jj_nb,nhori, 1,1,1, nvert, 32,
354	c . "ave(X)", zsto,zout)
355
356	ENDIF !lev_histday.GE.5
357	c-------------------------------------------------------
358
359	CALL histend(nid_day)
360
361	ENDIF ! fin de test sur ok_journe

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

Download in other formats: