Context Navigation

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

Last change on this file since 1314 was 1305, checked in by slebonnois, 11 years ago
SL: VENUS PHOTOCHEMISTRY. Needs Lapack (see arch files...)
File size: 10.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	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
26	. "once", zsto,zout)
27	c
28	CALL histdef(nid_ins, "aire", "Grid area", "-",
29	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
30	. "once", zsto,zout)
31	c
32	CALL histdef(nid_ins, "tsol", "Surface Temperature", "K",
33	. iim,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	. iim,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 . iim,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 . iim,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 . iim,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 . iim,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	. iim,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	. iim,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	. iim,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	. iim,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	. iim,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	. iim,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	. iim,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	. iim,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 . iim,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 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
101	c . "ins(X)", zsto,zout)
102	c
103	c plusieurs traceurs
104	if (iflag_trac.eq.1) THEN
105	DO iq=1,nqmax
106	IF (iq.LE.99) THEN
107	WRITE(str2,'(i2.2)') iq
108	CALL histdef(nid_ins, tname(iq), ttext(iq), "vmr",
109	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
110	. "ins(X)", zsto,zout)
111	ELSE
112	PRINT*, "Trop de traceurs"
113	CALL abort
114	ENDIF
115	ENDDO
116	endif
117	c
118	CALL histdef(nid_ins, "tops", "Solar rad. at TOA", "W/m2",
119	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
120	. "ins(X)", zsto,zout)
121	c
122	if (ok_cloud) THEN
123	CALL histdef(nid_ins, "NBRTOT", "Nbr total droplet", "#/cm3",
124	. iim,jj_nb,nhori, klev,1,klev, nvert, 32,
125	. "ins(X)", zsto,zout)
126	CALL histdef(nid_ins, "WH2SO4", "Weight fraction H2SO4",
127	. "fraction",iim,jj_nb,nhori, klev,1,klev, nvert, 32,
128	. "ins(X)", zsto,zout)
129	CALL histdef(nid_ins, "R_MEDIAN",
130	. "Median radius fo log normal distribution" ,
131	. "fraction",iim,jj_nb,nhori, klev,1,klev, nvert, 32,
132	. "ins(X)", zsto,zout)
133	CALL histdef(nid_ins, "STDDEV",
134	. "Std Deviation for lor normaldistribution",
135	. "fraction",iim,jj_nb,nhori, klev,1,klev, nvert, 32,
136	. "ins(X)", zsto,zout)
137	CALL histdef(nid_ins, "rho_droplet", "density cloud droplet",
138	. "kg.m-3",iim,jj_nb,nhori, klev,1,klev, nvert, 32,
139	. "ins(X)", zsto,zout)
140	endif
141
142	if (ok_sedim) THEN
143	CALL histdef(nid_ins,"d_tr_sed_H2SO4","H2SO4 mmr from sedim",
144	. "kg/kg",iim,jj_nb,nhori, klev,1,klev, nvert, 32,
145	. "ins(X)", zsto,zout)
146	CALL histdef(nid_ins,"d_tr_sed_H2O", "H2O mmr from sedim",
147	. "kg/kg",iim,jj_nb,nhori, klev,1,klev, nvert, 32,
148	. "ins(X)", zsto,zout)
149	CALL histdef(nid_ins, "F_sedim", "tendency from sedim",
150	. "kg.m-2.s-1",iim,jj_nb,nhori, klev,1,klev, nvert, 32,
151	. "ins(X)", zsto,zout)
152	endif
153
154	ENDIF !lev_histins.GE.2
155	c
156	c-------------------------------------------------------
157	IF(lev_histins.GE.3) THEN
158	c
159	cccccccccccccccccc Radiative transfer
160	c
161	c 2D
162	c
163	CALL histdef(nid_ins, "topl", "IR rad. at TOA", "W/m2",
164	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
165	. "ins(X)", zsto,zout)
166	c
167	CALL histdef(nid_ins, "sols", "Solar rad. at surf.", "W/m2",
168	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
169	. "ins(X)", zsto,zout)
170	c
171	CALL histdef(nid_ins, "soll", "IR rad. at surface", "W/m2",
172	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
173	. "ins(X)", zsto,zout)
174	c
175	c 3D
176	c
177	CALL histdef(nid_ins, "SWnet", "Net SW flux","W/m2",
178	. iim,jj_nb,nhori, klev,1,klev,nvert,
179	. 32, "ins(X)", zsto,zout)
180	c
181	CALL histdef(nid_ins, "LWnet", "Net LW flux","W/m2",
182	. iim,jj_nb,nhori, klev,1,klev,nvert,
183	. 32, "ins(X)", zsto,zout)
184	c
185	CALL histdef(nid_ins, "fluxvdf", "PBL net flux","W/m2",
186	. iim,jj_nb,nhori, klev,1,klev,nvert,
187	. 32, "ins(X)", zsto,zout)
188	c
189	CALL histdef(nid_ins, "fluxdyn", "Dyn. net flux","W/m2",
190	. iim,jj_nb,nhori, klev,1,klev,nvert,
191	. 32, "ins(X)", zsto,zout)
192	c
193	CALL histdef(nid_ins, "fluxajs", "Dry adj. net flux","W/m2",
194	. iim,jj_nb,nhori, klev,1,klev,nvert,
195	. 32, "ins(X)", zsto,zout)
196	c
197	c CALL histdef(nid_ins, "fluxec", "Cin. net flux","W/m2",
198	c . iim,jj_nb,nhori, klev,1,klev,nvert,
199	c . 32, "ins(X)", zsto,zout)
200	c
201	ENDIF !lev_histins.GE.3
202	c
203	c-------------------------------------------------------
204	IF(lev_histins.GE.4) THEN
205	c
206	CALL histdef(nid_ins, "dtdyn", "Dynamics dT", "K/s",
207	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
208	. "ins(X)", zsto,zout)
209	c
210	c CALL histdef(nid_ins, "dtphy", "Physics dT", "K/s",
211	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
212	c . "ins(X)", zsto,zout)
213	c
214	CALL histdef(nid_ins, "dtvdf", "Boundary-layer dT", "K/s",
215	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
216	. "ins(X)", zsto,zout)
217	c
218	CALL histdef(nid_ins, "dtajs", "Dry adjust. dT", "K/s",
219	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
220	. "ins(X)", zsto,zout)
221	c
222	CALL histdef(nid_ins, "dtswr", "SW radiation dT", "K/s",
223	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
224	. "ins(X)", zsto,zout)
225	c
226	CALL histdef(nid_ins, "dtlwr", "LW radiation dT", "K/s",
227	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
228	. "ins(X)", zsto,zout)
229	c
230	c CALL histdef(nid_ins, "dtec", "Cinetic dissip dT", "K/s",
231	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
232	c . "ins(X)", zsto,zout)
233	c
234	CALL histdef(nid_ins, "duajs", "Dry convection dU", "m/s2",
235	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
236	. "ins(X)", zsto,zout)
237	c
238	CALL histdef(nid_ins, "dugwo", "GW oro dU", "m/s2",
239	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
240	. "ins(X)", zsto,zout)
241	c
242	CALL histdef(nid_ins, "dugwno", "GW non-oro dU", "m/s2",
243	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
244	. "ins(X)", zsto,zout)
245	c
246	c CALL histdef(nid_ins, "dvvdf", "Boundary-layer dV", "m/s2",
247	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
248	c . "ins(X)", zsto,zout)
249	c
250	ENDIF !lev_histins.GE.4
251	c
252	c-------------------------------------------------------
253	IF(lev_histins.GE.5) THEN
254	c
255	c call histdef(nid_ins, "taux",
256	c $ "Zonal wind stress", "Pa",
257	c $ iim,jj_nb,nhori, 1,1,1, nvert, 32,
258	c $ "ins(X)", zsto,zout)
259	c
260	c call histdef(nid_ins, "tauy",
261	c $ "Meridional xind stress", "Pa",
262	c $ iim,jj_nb,nhori, 1,1,1, nvert, 32,
263	c $ "ins(X)", zsto,zout)
264	c
265	c CALL histdef(nid_ins, "cdrm", "Momentum drag coef.", "-",
266	c . iim,jj_nb,nhori, 1,1,1, nvert, 32,
267	c . "ins(X)", zsto,zout)
268	c
269	c CALL histdef(nid_ins, "cdrh", "Heat drag coef.", "-",
270	c . iim,jj_nb,nhori, 1,1,1, nvert, 32,
271	c . "ins(X)", zsto,zout)
272	c
273	ENDIF !lev_histins.GE.5
274	c-------------------------------------------------------
275	c
276	CALL histend(nid_ins)
277	!$OMP END MASTER
278
279	ENDIF

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

Download in other formats: