Context Navigation

ini_histmth.h @ 3461

Last change on this file since 3461 was 3451, checked in by emillour, 7 weeks ago
Venus PCM: Add "Age of Air" scheme from Maureen Cohen. Activated with flag "ok_aoa=y" and requires that there is also an "aoa" tracer in traceur.def MC
File size: 12.4 KB

Line
1	!
2	! $Header: /home/cvsroot/LMDZ4/libf/phylmd/ini_histmth.h,v 1.4 2005/03/09 12:30:16 fairhead Exp $
3	!
4	IF (ok_mensuel) THEN
5
6	zsto = dtime
7	zout = dtime * ecrit_mth
8
9	idayref = day_ref
10	CALL ymds2ju(annee_ref, 1, idayref, zero, zjulian)
11
12	call histbeg_phy("histmth.nc",itau_phy,
13	. zjulian,dtime,nhori,nid_mth)
14
15	!$OMP MASTER
16	CALL histvert(nid_mth, "presnivs", "Vertical levels", "Pa",
17	. klev, presnivs, nvert)
18
19	c-------------------------------------------------------
20	IF(lev_histmth.GE.1) THEN
21	c
22	ccccccccccccc 2D fields, basics
23	c
24	CALL histdef(nid_mth, "phis", "Surface geop. height", "-",
25	. nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
26	. "once", zsto,zout)
27	c
28	CALL histdef(nid_mth, "aire", "Grid area", "-",
29	. nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
30	. "once", zsto,zout)
31	c
32	CALL histdef(nid_mth, "tsol", "Surface Temperature", "K",
33	. nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
34	. "ave(X)", zsto,zout)
35	c
36	CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa",
37	. nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
38	. "ave(X)", zsto,zout)
39	c
40	c CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
41	c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
42	c . "ave(X)", zsto,zout)
43	c
44	c CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
45	c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
46	c . "ave(X)", zsto,zout)
47	c
48	c CALL histdef(nid_mth, "cdragh", "Drag coef on T", "-",
49	c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
50	c . "ave(X)", zsto,zout)
51	c
52	c CALL histdef(nid_mth, "cdragm", "Drag coef on U", "-",
53	c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
54	c . "ave(X)", zsto,zout)
55	c
56	ENDIF !lev_histmth.GE.1
57	c
58	c-------------------------------------------------------
59	IF(lev_histmth.GE.2) THEN
60	c
61	ccccccccccccc 3D fields, basics
62	c
63	CALL histdef(nid_mth, "temp", "Air temperature", "K",
64	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
65	. "ave(X)", zsto,zout)
66	c
67	CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
68	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
69	. "ave(X)", zsto,zout)
70	c
71	CALL histdef(nid_mth, "geop", "Geopotential height", "m",
72	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
73	. "ave(X)", zsto,zout)
74	c
75	CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s",
76	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
77	. "ave(X)", zsto,zout)
78	c
79	CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s",
80	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
81	. "ave(X)", zsto,zout)
82	c
83	CALL histdef(nid_mth, "vitw", "Vertical wind", "Pa/s",
84	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
85	. "ave(X)", zsto,zout)
86	c
87	CALL histdef(nid_mth, "dudyn", "Dynamics dU", "m/s2",
88	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
89	. "ave(X)", zsto,zout)
90	c
91	CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
92	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
93	. "ave(X)", zsto,zout)
94	c
95	c CALL histdef(nid_mth, "mang", "Angular momentum", "kg m2/s",
96	c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
97	c . "ave(X)", zsto,zout)
98	c
99	c CALL histdef(nid_mth, "Kz", "vertical diffusion coef", "m2/s",
100	c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
101	c . "ave(X)", zsto,zout)
102	c
103	CALL histdef(nid_mth, "mmean", "Mean molecular mass", "g/mol",
104	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
105	. "ave(X)", zsto,zout)
106
107	c CALL histdef(nid_mth, "rho", "Air density [mass/Vol]", "kg/m3",
108	c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
109	c . "ave(X)", zsto,zout)
110
111	IF ((iflag_trac.eq.1) .and. ok_aoa) THEN
112	CALL histdef(nid_mth, "age", "Age of air", "s",
113	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
114	. "ave(X)", zsto,zout)
115	ENDIF
116
117	c plusieurs traceurs !!!outputs in [vmr]
118	IF (iflag_trac.eq.1) THEN
119	DO iq=1,nqmax
120	IF (iq.LE.99) THEN
121	WRITE(str2,'(i2.2)') iq
122	CALL histdef(nid_mth, tname(iq), ttext(iq), "mol/mol",
123	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
124	. "ave(X)", zsto,zout)
125	ELSE
126	PRINT*, "Trop de traceurs"
127	CALL abort
128	ENDIF
129	ENDDO
130	ENDIF
131
132	IF (callthermos .and. ok_chem) THEN
133	CALL histdef(nid_mth, "d_qmoldifCO2", "Dif molec" , "kg/kg",
134	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
135	. "ave(X)", zsto,zout)
136	CALL histdef(nid_mth, "d_qmoldifO3p", "Dif molec" , "kg/kg",
137	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
138	. "ave(X)", zsto,zout)
139	CALL histdef(nid_mth, "d_qmoldifN2", "Dif molec" , "kg/kg",
140	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
141	. "ave(X)", zsto,zout)
142	ENDIF
143	c
144	CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2",
145	. nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
146	. "ave(X)", zsto,zout)
147	c
148	ENDIF !lev_histmth.GE.2
149	c
150	c-------------------------------------------------------
151	IF(lev_histmth.GE.3) THEN
152	c
153	cccccccccccccccccc Radiative transfer
154	c
155	c 2D
156	c
157	CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2",
158	. nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
159	. "ave(X)", zsto,zout)
160	c
161	CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2",
162	. nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
163	. "ave(X)", zsto,zout)
164	c
165	CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2",
166	. nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
167	. "ave(X)", zsto,zout)
168	c
169	c 3D
170	c
171	CALL histdef(nid_mth, "SWnet", "Net SW flux","W/m2",
172	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
173	. 32, "ave(X)", zsto,zout)
174	c
175	CALL histdef(nid_mth, "LWnet", "Net LW flux","W/m2",
176	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
177	. 32, "ave(X)", zsto,zout)
178	c
179	c CALL histdef(nid_mth, "fluxvdf", "PBL net flux","W/m2",
180	c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
181	c . 32, "ave(X)", zsto,zout)
182	c
183	c CALL histdef(nid_mth, "fluxdyn", "Dyn. net flux","W/m2",
184	c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
185	c . 32, "ave(X)", zsto,zout)
186	c
187	c CALL histdef(nid_mth, "fluxajs", "Dry adj. net flux","W/m2",
188	c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
189	c . 32, "ave(X)", zsto,zout)
190	c
191	c CALL histdef(nid_mth, "fluxec", "Cin. net flux","W/m2",
192	c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
193	c . 32, "ave(X)", zsto,zout)
194	c
195	ENDIF !lev_histmth.GE.3
196	c
197	c-------------------------------------------------------
198	IF(lev_histmth.GE.4) THEN
199	c
200	CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
201	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
202	. "ave(X)", zsto,zout)
203	c
204	c CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
205	c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
206	c . "ave(X)", zsto,zout)
207	c
208	CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s",
209	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
210	. "ave(X)", zsto,zout)
211	c
212	CALL histdef(nid_mth, "dtajs", "Dry adjust. dT", "K/s",
213	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
214	. "ave(X)", zsto,zout)
215	c
216	CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s",
217	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
218	. "ave(X)", zsto,zout)
219	c
220	CALL histdef(nid_mth,"dtswrNLTE", "SWNLTE radiation dT",
221	. "K/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
222	. "ave(X)", zsto,zout)
223	c
224	CALL histdef(nid_mth,"dtswrDCrisp","SWDCrisp radiation dT",
225	. "K/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
226	. "ave(X)", zsto,zout)
227	c
228	CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
229	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
230	. "ave(X)", zsto,zout)
231	c
232	c CALL histdef(nid_mth,"dtlwrNLTE", "LWNLTE radiation dT",
233	c . "K/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
234	c . "ave(X)", zsto,zout)
235	c
236	c CALL histdef(nid_mth,"dtlwrLTE","LW_LTE radiation dT",
237	c . "K/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
238	c . "ave(X)", zsto,zout)
239	c
240	c CALL histdef(nid_mth, "dteuv", "UV radiation dT", "K/s",
241	c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
242	c . "ave(X)", zsto,zout)
243	c CALL histdef(nid_mth, "dtcond", "Therm conduction", "K/s",
244	c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
245	c . "ave(X)", zsto,zout)
246	c CALL histdef(nid_mth, "dumolvis", "molec viscosity (u)"
247	c . ,"m/s2",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
248	c . "ave(X)", zsto,zout)
249	c CALL histdef(nid_mth, "dvmolvis", "molec viscosity (v)",
250	c . "m/s2",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
251	c . "ave(X)", zsto,zout)
252
253	c
254	c CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
255	c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
256	c . "ave(X)", zsto,zout)
257	c
258	CALL histdef(nid_mth, "duajs", "Dry convection dU", "m/s2",
259	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
260	. "ave(X)", zsto,zout)
261	c
262	CALL histdef(nid_mth, "dugwo", "GW oro dU", "m/s2",
263	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
264	. "ave(X)", zsto,zout)
265	c
266	CALL histdef(nid_mth, "dugwno", "GW non-oro dU", "m/s2",
267	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
268	. "ave(X)", zsto,zout)
269
270	c CALL histdef(nid_mth, "co2_vmr", "volume mixture ratio co2",
271	c . "cm-3",nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
272	c . "ave(X)", zsto,zout)
273
274	c CALL histdef(nid_mth, "o_vmr", "volume mixture ratio o",
275	c . "cm-3",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
276	c . "ave(X)", zsto,zout)
277	c CALL histdef(nid_mth, "co_vmr", "volume mixture ratio co",
278	c . "cm-3",nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
279	c . "ave(X)", zsto,zout)
280
281	c CALL histdef(nid_mth, "n2_vmr", "volume mixture ratio n2",
282	c . "cm-3",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
283	c . "ave(X)", zsto,zout)
284	c
285	c CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
286	c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
287	c . "ave(X)", zsto,zout)
288	c
289	ENDIF !lev_histmth.GE.4
290	c
291	c-------------------------------------------------------
292	IF(lev_histmth.GE.5) THEN
293	c
294	c call histdef(nid_mth, "taux",
295	c $ "Zonal wind stress", "Pa",
296	c $ nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
297	c $ "ave(X)", zsto,zout)
298	c
299	c call histdef(nid_mth, "tauy",
300	c $ "Meridional xind stress", "Pa",
301	c $ nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
302	c $ "ave(X)", zsto,zout)
303	c
304	c CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
305	c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
306	c . "ave(X)", zsto,zout)
307	c
308	c CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
309	c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
310	c . "ave(X)", zsto,zout)
311	c
312	ENDIF !lev_histmth.GE.5
313	c-------------------------------------------------------
314
315	CALL histend(nid_mth)
316	!$OMP END MASTER
317
318	ENDIF ! fin de test sur ok_mensuel
319

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

Context Navigation

source: trunk/LMDZ.VENUS/libf/phyvenus/ini_histmth.h @ 3461

Download in other formats: