Context Navigation

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

Last change on this file since 937 was 902, checked in by slebonnois, 12 years ago
SL: small modifs on Venus outputs
File size: 9.0 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	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
26	. "once", zsto,zout)
27	c
28	CALL histdef(nid_mth, "aire", "Grid area", "-",
29	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
30	. "once", zsto,zout)
31	c
32	CALL histdef(nid_mth, "tsol", "Surface Temperature", "K",
33	. iim,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	. iim,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 . iim,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 . iim,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 . iim,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 . iim,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	. iim,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	. iim,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	. iim,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	. iim,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	. iim,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	. iim,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	. iim,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	. iim,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 . iim,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 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
101	c . "ave(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_mth, tname(iq), ttext(iq), "ppm",
109	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
110	. "ave(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_mth, "tops", "Solar rad. at TOA", "W/m2",
119	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
120	. "ave(X)", zsto,zout)
121	c
122	ENDIF !lev_histmth.GE.2
123	c
124	c-------------------------------------------------------
125	IF(lev_histmth.GE.3) THEN
126	c
127	cccccccccccccccccc Radiative transfer
128	c
129	c 2D
130	c
131	CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2",
132	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
133	. "ave(X)", zsto,zout)
134	c
135	CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2",
136	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
137	. "ave(X)", zsto,zout)
138	c
139	CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2",
140	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
141	. "ave(X)", zsto,zout)
142	c
143	c 3D
144	c
145	CALL histdef(nid_mth, "SWnet", "Net SW flux","W/m2",
146	. iim,jj_nb,nhori, klev,1,klev,nvert,
147	. 32, "ave(X)", zsto,zout)
148	c
149	CALL histdef(nid_mth, "LWnet", "Net LW flux","W/m2",
150	. iim,jj_nb,nhori, klev,1,klev,nvert,
151	. 32, "ave(X)", zsto,zout)
152	c
153	CALL histdef(nid_mth, "fluxvdf", "PBL net flux","W/m2",
154	. iim,jj_nb,nhori, klev,1,klev,nvert,
155	. 32, "ave(X)", zsto,zout)
156	c
157	CALL histdef(nid_mth, "fluxdyn", "Dyn. net flux","W/m2",
158	. iim,jj_nb,nhori, klev,1,klev,nvert,
159	. 32, "ave(X)", zsto,zout)
160	c
161	CALL histdef(nid_mth, "fluxajs", "Dry adj. net flux","W/m2",
162	. iim,jj_nb,nhori, klev,1,klev,nvert,
163	. 32, "ave(X)", zsto,zout)
164	c
165	c CALL histdef(nid_mth, "fluxec", "Cin. net flux","W/m2",
166	c . iim,jj_nb,nhori, klev,1,klev,nvert,
167	c . 32, "ave(X)", zsto,zout)
168	c
169	ENDIF !lev_histmth.GE.3
170	c
171	c-------------------------------------------------------
172	IF(lev_histmth.GE.4) THEN
173	c
174	CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
175	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
176	. "ave(X)", zsto,zout)
177	c
178	c CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
179	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
180	c . "ave(X)", zsto,zout)
181	c
182	CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s",
183	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
184	. "ave(X)", zsto,zout)
185	c
186	CALL histdef(nid_mth, "dtajs", "Dry adjust. dT", "K/s",
187	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
188	. "ave(X)", zsto,zout)
189	c
190	CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s",
191	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
192	. "ave(X)", zsto,zout)
193	c
194	CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
195	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
196	. "ave(X)", zsto,zout)
197	c
198	c CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
199	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
200	c . "ave(X)", zsto,zout)
201	c
202	CALL histdef(nid_mth, "duajs", "Dry convection dU", "m/s2",
203	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
204	. "ave(X)", zsto,zout)
205	c
206	CALL histdef(nid_mth, "dugwo", "GW oro dU", "m/s2",
207	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
208	. "ave(X)", zsto,zout)
209	c
210	CALL histdef(nid_mth, "dugwno", "GW non-oro dU", "m/s2",
211	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
212	. "ave(X)", zsto,zout)
213	c
214	c CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
215	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
216	c . "ave(X)", zsto,zout)
217	c
218	ENDIF !lev_histmth.GE.4
219	c
220	c-------------------------------------------------------
221	IF(lev_histmth.GE.5) THEN
222	c
223	c call histdef(nid_mth, "taux",
224	c $ "Zonal wind stress", "Pa",
225	c $ iim,jj_nb,nhori, 1,1,1, nvert, 32,
226	c $ "ave(X)", zsto,zout)
227	c
228	c call histdef(nid_mth, "tauy",
229	c $ "Meridional xind stress", "Pa",
230	c $ iim,jj_nb,nhori, 1,1,1, nvert, 32,
231	c $ "ave(X)", zsto,zout)
232	c
233	c CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
234	c . iim,jj_nb,nhori, 1,1,1, nvert, 32,
235	c . "ave(X)", zsto,zout)
236	c
237	c CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
238	c . iim,jj_nb,nhori, 1,1,1, nvert, 32,
239	c . "ave(X)", zsto,zout)
240	c
241	ENDIF !lev_histmth.GE.5
242	c-------------------------------------------------------
243
244	CALL histend(nid_mth)
245	!$OMP END MASTER
246
247	ENDIF ! fin de test sur ok_mensuel

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

Download in other formats: