Context Navigation

ini_histins.h @ 3556

Last change on this file since 3556 was 1661, checked in by slebonnois, 8 years ago
SL: Cloud model for Venus. Not validated yet.
File size: 13.5 KB

Line
1	!
2	! $Header: /home/cvsroot/LMDZ4/libf/phylmd/ini_histins.h,v 1.1.1.1 2004/05/19 12:53:08 lmdzadmin Exp $
3	!
4	IF (ok_instan) THEN
5
6	zsto = dtime * ecrit_ins
7	zout = dtime * ecrit_ins
8
9	idayref = day_ref
10	CALL ymds2ju(annee_ref, 1, idayref, zero, zjulian)
11
12	call histbeg_phy("histins.nc",itau_phy,
13	. zjulian,dtime,nhori,nid_ins)
14
15	!$OMP MASTER
16	CALL histvert(nid_ins, "presnivs", "Vertical levels", "Pa",
17	. klev, presnivs, nvert)
18
19	c-------------------------------------------------------
20	IF(lev_histins.GE.1) THEN
21	c
22	ccccccccccccc 2D fields, basics
23	c
24	CALL histdef(nid_ins, "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_ins, "aire", "Grid area", "-",
29	. nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
30	. "once", zsto,zout)
31	c
32	CALL histdef(nid_ins, "tsol", "Surface Temperature", "K",
33	. nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
34	. "ins(X)", zsto,zout)
35	c
36	CALL histdef(nid_ins, "psol", "Surface Pressure", "Pa",
37	. nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
38	. "ins(X)", zsto,zout)
39	c
40	c CALL histdef(nid_ins, "ue", "Zonal energy transport", "-",
41	c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
42	c . "ins(X)", zsto,zout)
43	c
44	c CALL histdef(nid_ins, "ve", "Merid energy transport", "-",
45	c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
46	c . "ins(X)", zsto,zout)
47	c
48	c CALL histdef(nid_ins, "cdragh", "Drag coef on T", "-",
49	c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
50	c . "ins(X)", zsto,zout)
51	c
52	c CALL histdef(nid_ins, "cdragm", "Drag coef on U", "-",
53	c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
54	c . "ins(X)", zsto,zout)
55	c
56	ENDIF !lev_histins.GE.1
57	c
58	c-------------------------------------------------------
59	IF(lev_histins.GE.2) THEN
60	c
61	ccccccccccccc 3D fields, basics
62	c
63	CALL histdef(nid_ins, "temp", "Air temperature", "K",
64	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
65	. "ins(X)", zsto,zout)
66	c
67	CALL histdef(nid_ins, "pres", "Air pressure", "Pa",
68	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
69	. "ins(X)", zsto,zout)
70	c
71	CALL histdef(nid_ins, "geop", "Geopotential height", "m",
72	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
73	. "ins(X)", zsto,zout)
74	c
75	CALL histdef(nid_ins, "vitu", "Zonal wind", "m/s",
76	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
77	. "ins(X)", zsto,zout)
78	c
79	CALL histdef(nid_ins, "vitv", "Meridional wind", "m/s",
80	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
81	. "ins(X)", zsto,zout)
82	c
83	CALL histdef(nid_ins, "vitw", "Vertical wind", "Pa/s",
84	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
85	. "ins(X)", zsto,zout)
86	c
87	CALL histdef(nid_ins, "dudyn", "Dynamics dU", "m/s2",
88	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
89	. "ins(X)", zsto,zout)
90	c
91	CALL histdef(nid_ins, "duvdf", "Boundary-layer dU", "m/s2",
92	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
93	. "ins(X)", zsto,zout)
94	c
95	c CALL histdef(nid_ins, "mang", "Angular momentum", "kg m2/s",
96	c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
97	c . "ins(X)", zsto,zout)
98	c
99	c CALL histdef(nid_ins, "Kz", "vertical diffusion coef", "m2/s",
100	c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
101	c . "ins(X)", zsto,zout)
102	c
103	CALL histdef(nid_ins, "mmean", "Mean molecular mass", "g/mol",
104	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
105	. "ins(X)", zsto,zout)
106
107	CALL histdef(nid_ins, "rho", "Air density [mass/Vol]", "kg/m3",
108	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
109	. "ins(X)", zsto,zout)
110
111
112	c plusieurs traceurs !!!outputs in [vmr]
113	IF (iflag_trac.eq.1) THEN
114	DO iq=1,nqmax
115	IF (iq.LE.99) THEN
116	WRITE(str2,'(i2.2)') iq
117	CALL histdef(nid_ins, tname(iq), ttext(iq), "mol/mol",
118	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
119	. "ins(X)", zsto,zout)
120	ELSE
121	PRINT*, "Trop de traceurs"
122	CALL abort
123	ENDIF
124	ENDDO
125	ENDIF
126
127	IF (callthermos .and. ok_chem) THEN
128	CALL histdef(nid_ins, "d_qmoldif CO2", "Dif molec" , "kg/kg",
129	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
130	. "ins(X)", zsto,zout)
131	CALL histdef(nid_ins, "d_qmoldif O3p", "Dif molec" , "kg/kg",
132	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
133	. "ins(X)", zsto,zout)
134	CALL histdef(nid_ins, "d_qmoldif N2", "Dif molec" , "kg/kg",
135	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
136	. "ins(X)", zsto,zout)
137	ENDIF
138	c
139	CALL histdef(nid_ins, "tops", "Solar rad. at TOA", "W/m2",
140	. nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
141	. "ins(X)", zsto,zout)
142	c
143	if (ok_cloud.and.(cl_scheme.eq.1)) THEN
144
145	if (nb_mode.GE.1) THEN
146
147	c
148	CALL histdef(nid_ins, "NBRTOTm1", "Nbr total droplet",
149	. "#/cm3",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
150	. "ins(X)", zsto,zout)
151	c
152
153	c
154	c CALL histdef(nid_ins, "R_MEDIANm1", "Median radius
155	c . for log normal distribution" ,
156	c . "fraction",
157	c . nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
158	c . "ins(X)", zsto,zout)
159	c
160
161	c
162	c CALL histdef(nid_ins, "STDDEVm1", "Std Deviation
163	c . for log normal distribution",
164	c . "fraction",
165	c . nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
166	c . "ins(X)", zsto,zout)
167	c
168
169	if (nb_mode.GE.2) THEN
170
171	c
172	CALL histdef(nid_ins, "NBRTOTm2", "Nbr total droplet",
173	. "#/cm3",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
174	. "ins(X)", zsto,zout)
175	c
176
177	c
178	c CALL histdef(nid_ins, "R_MEDIANm2", "Median radius
179	c . for log normal distribution" ,
180	c . "fraction",
181	c . nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
182	c . "ins(X)", zsto,zout)
183	c
184
185	c
186	c CALL histdef(nid_ins, "STDDEVm2", "Std Deviation
187	c . for log normal distribution",
188	c . "fraction",
189	c . nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
190	c . "ins(X)", zsto,zout)
191	c
192
193	if (nb_mode.GE.3) THEN
194
195	c
196	CALL histdef(nid_ins, "NBRTOTm3", "Nbr total droplet", "#/cm3",
197	. nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
198	. "ins(X)", zsto,zout)
199	c
200
201	c
202	c CALL histdef(nid_ins, "R_MEDIANm3", "Median radius
203	c . for log normal distribution" ,
204	c . "fraction",
205	c . nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
206	c . "ins(X)", zsto,zout)
207	c
208
209	c
210	c CALL histdef(nid_ins, "STDDEVm3", "Std Deviation
211	c . for log normal distribution",
212	c . "fraction",
213	c . nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
214	c . "ins(X)", zsto,zout)
215	c
216
217	ENDIF
218	ENDIF
219	ENDIF
220
221	c
222	CALL histdef(nid_ins, "WH2SO4", "Weight fraction H2SO4",
223	. "fraction",
224	. nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
225	. "ins(X)", zsto,zout)
226	c
227
228	c
229	CALL histdef(nid_ins, "rho_droplet", "density cloud droplet",
230	. "kg.m-3",
231	. nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
232	. "ins(X)", zsto,zout)
233	c
234
235	ENDIF
236
237	if (ok_sedim.and.(cl_scheme.eq.1)) THEN
238	c
239	CALL histdef(nid_ins, "d_tr_sed_H2SO4", "var mmr from sedim",
240	. "kg/kg",
241	. nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
242	. "ins(X)", zsto,zout)
243	c
244
245	c
246	CALL histdef(nid_ins, "d_tr_sed_H2O", "var mmr from sedim",
247	. "kg/kg",
248	. nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
249	. "ins(X)", zsto,zout)
250	c
251
252	c
253	CALL histdef(nid_ins, "F_sedim", "tendency from sedim",
254	. "kg.m-2.s-1",
255	. nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
256	. "ins(X)", zsto,zout)
257	c
258	ENDIF
259
260	ENDIF !lev_histins.GE.2
261	c
262	c-------------------------------------------------------
263	IF(lev_histins.GE.3) THEN
264	c
265	cccccccccccccccccc Radiative transfer
266	c
267	c 2D
268	c
269	CALL histdef(nid_ins, "topl", "IR rad. at TOA", "W/m2",
270	. nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
271	. "ins(X)", zsto,zout)
272	c
273	CALL histdef(nid_ins, "sols", "Solar rad. at surf.", "W/m2",
274	. nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
275	. "ins(X)", zsto,zout)
276	c
277	CALL histdef(nid_ins, "soll", "IR rad. at surface", "W/m2",
278	. nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
279	. "ins(X)", zsto,zout)
280	c
281	c 3D
282	c
283	CALL histdef(nid_ins, "SWnet", "Net SW flux","W/m2",
284	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
285	. 32, "ins(X)", zsto,zout)
286	c
287	CALL histdef(nid_ins, "LWnet", "Net LW flux","W/m2",
288	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
289	. 32, "ins(X)", zsto,zout)
290	c
291	CALL histdef(nid_ins, "fluxvdf", "PBL net flux","W/m2",
292	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
293	. 32, "ins(X)", zsto,zout)
294	c
295	CALL histdef(nid_ins, "fluxdyn", "Dyn. net flux","W/m2",
296	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
297	. 32, "ins(X)", zsto,zout)
298	c
299	CALL histdef(nid_ins, "fluxajs", "Dry adj. net flux","W/m2",
300	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
301	. 32, "ins(X)", zsto,zout)
302	c
303	c CALL histdef(nid_ins, "fluxec", "Cin. net flux","W/m2",
304	c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
305	c . 32, "ins(X)", zsto,zout)
306	c
307	ENDIF !lev_histins.GE.3
308	c
309	c-------------------------------------------------------
310	IF(lev_histins.GE.4) THEN
311	c
312	CALL histdef(nid_ins, "dtdyn", "Dynamics dT", "K/s",
313	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
314	. "ins(X)", zsto,zout)
315	c
316	c CALL histdef(nid_ins, "dtphy", "Physics dT", "K/s",
317	c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
318	c . "ins(X)", zsto,zout)
319	c
320	CALL histdef(nid_ins, "dtvdf", "Boundary-layer dT", "K/s",
321	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
322	. "ins(X)", zsto,zout)
323	c
324	CALL histdef(nid_ins, "dtajs", "Dry adjust. dT", "K/s",
325	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
326	. "ins(X)", zsto,zout)
327	c
328	CALL histdef(nid_ins, "dtswr", "SW radiation dT", "K/s",
329	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
330	. "ins(X)", zsto,zout)
331	c
332	CALL histdef(nid_ins, "dtlwr", "LW radiation dT", "K/s",
333	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
334	. "ins(X)", zsto,zout)
335	c
336	c CALL histdef(nid_ins, "dtec", "Cinetic dissip dT", "K/s",
337	c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
338	c . "ins(X)", zsto,zout)
339	c
340	CALL histdef(nid_ins, "duajs", "Dry convection dU", "m/s2",
341	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
342	. "ins(X)", zsto,zout)
343	c
344	CALL histdef(nid_ins, "dugwo", "GW oro dU", "m/s2",
345	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
346	. "ins(X)", zsto,zout)
347	c
348	CALL histdef(nid_ins, "dugwno", "GW non-oro dU", "m/s2",
349	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
350	. "ins(X)", zsto,zout)
351	c
352	c CALL histdef(nid_ins, "dvvdf", "Boundary-layer dV", "m/s2",
353	c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
354	c . "ins(X)", zsto,zout)
355	c
356	ENDIF !lev_histins.GE.4
357	c
358	c-------------------------------------------------------
359	IF(lev_histins.GE.5) THEN
360	c
361	c call histdef(nid_ins, "taux",
362	c $ "Zonal wind stress", "Pa",
363	c $ nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
364	c $ "ins(X)", zsto,zout)
365	c
366	c call histdef(nid_ins, "tauy",
367	c $ "Meridional xind stress", "Pa",
368	c $ nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
369	c $ "ins(X)", zsto,zout)
370	c
371	c CALL histdef(nid_ins, "cdrm", "Momentum drag coef.", "-",
372	c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
373	c . "ins(X)", zsto,zout)
374	c
375	c CALL histdef(nid_ins, "cdrh", "Heat drag coef.", "-",
376	c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
377	c . "ins(X)", zsto,zout)
378	c
379	ENDIF !lev_histins.GE.5
380	c-------------------------------------------------------
381	c
382	CALL histend(nid_ins)
383	!$OMP END MASTER
384
385	ENDIF

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

Context Navigation

source: trunk/LMDZ.VENUS/libf/phyvenus/ini_histins.h @ 3556

Download in other formats: