Context Navigation

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

Last change on this file since 894 was 892, checked in by slebonnois, 12 years ago
SL: Important commit ! Adaptation of Venus physics to parallel computation / template for arch on the LMD servers using ifort / documentation for 1D column physics and for parallel computations
File size: 9.2 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	c zsto1: pour des valeurs "instantanees" mensuelles
10	zsto1 = dtime * ecrit_mth
11	c zsto2: pour des flux radiatifs calcules avec un autre pas de temps
12	zsto2 = dtime * radpas
13	PRINT*,' zsto,zsto1,zsto2,zout=',zsto, zsto1, zsto2,zout
14
15	idayref = day_ref
16	CALL ymds2ju(annee_ref, 1, idayref, zero, zjulian)
17
18	call histbeg_phy("histmth.nc",itau_phy,
19	. zjulian,zsto,nhori,nid_mth)
20
21	!$OMP MASTER
22	CALL histvert(nid_mth, "presnivs", "Vertical levels", "Pa",
23	. klev, presnivs, nvert)
24
25	c-------------------------------------------------------
26	IF(lev_histmth.GE.1) THEN
27	c
28	ccccccccccccc 2D fields, basics
29	c
30	CALL histdef(nid_mth, "phis", "Surface geop. height", "-",
31	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
32	. "once", zsto,zout)
33	c
34	CALL histdef(nid_mth, "aire", "Grid area", "-",
35	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
36	. "once", zsto,zout)
37	c
38	CALL histdef(nid_mth, "tsol", "Surface Temperature", "K",
39	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
40	. "ave(X)", zsto,zout)
41	c
42	CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa",
43	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
44	. "ave(X)", zsto,zout)
45	c
46	c CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
47	c . iim,jj_nb,nhori, 1,1,1, nvert, 32,
48	c . "ave(X)", zsto,zout)
49	c
50	c CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
51	c . iim,jj_nb,nhori, 1,1,1, nvert, 32,
52	c . "ave(X)", zsto,zout)
53	c
54	CALL histdef(nid_mth, "cdragh", "Drag coef on T", "-",
55	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
56	. "ave(X)", zsto,zout)
57	c
58	CALL histdef(nid_mth, "cdragm", "Drag coef on U", "-",
59	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
60	. "ave(X)", zsto,zout)
61	c
62	ENDIF !lev_histmth.GE.1
63	c
64	c-------------------------------------------------------
65	IF(lev_histmth.GE.2) THEN
66	c
67	ccccccccccccc 3D fields, basics
68	c
69	CALL histdef(nid_mth, "temp", "Air temperature", "K",
70	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
71	. "ave(X)", zsto,zout)
72	c
73	CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
74	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
75	. "ave(X)", zsto,zout)
76	c
77	CALL histdef(nid_mth, "geop", "Geopotential height", "m",
78	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
79	. "ave(X)", zsto,zout)
80	c
81	CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s",
82	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
83	. "ave(X)", zsto,zout)
84	c
85	CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s",
86	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
87	. "ave(X)", zsto,zout)
88	c
89	CALL histdef(nid_mth, "vitw", "Vertical wind", "Pa/s",
90	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
91	. "ave(X)", zsto,zout)
92	c
93	CALL histdef(nid_mth, "dudyn", "Dynamics dU", "m/s2",
94	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
95	. "ave(X)", zsto,zout)
96	c
97	CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
98	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
99	. "ave(X)", zsto,zout)
100	c
101	c CALL histdef(nid_mth, "mang", "Angular momentum", "kg m2/s",
102	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
103	c . "ave(X)", zsto,zout)
104	c
105	CALL histdef(nid_mth, "Kz", "vertical diffusion coef", "m2/s",
106	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
107	. "ave(X)", zsto,zout)
108	c
109	c plusieurs traceurs
110	if (iflag_trac.eq.1) THEN
111	DO iq=1,nqmax
112	IF (iq.LE.99) THEN
113	WRITE(str2,'(i2.2)') iq
114	CALL histdef(nid_mth, tname(iq), ttext(iq), "ppm",
115	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
116	. "ave(X)", zsto,zout)
117	ELSE
118	PRINT*, "Trop de traceurs"
119	CALL abort
120	ENDIF
121	ENDDO
122	endif
123	c
124	CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2",
125	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
126	. "ave(X)", zsto2,zout)
127	c
128	ENDIF !lev_histmth.GE.2
129	c
130	c-------------------------------------------------------
131	IF(lev_histmth.GE.3) THEN
132	c
133	cccccccccccccccccc Radiative transfer
134	c
135	c 2D
136	c
137	CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2",
138	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
139	. "ave(X)", zsto2,zout)
140	c
141	CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2",
142	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
143	. "ave(X)", zsto2,zout)
144	c
145	CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2",
146	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
147	. "ave(X)", zsto2,zout)
148	c
149	c 3D
150	c
151	CALL histdef(nid_mth, "SWnet", "Net SW flux","W/m2",
152	. iim,jj_nb,nhori, klev,1,klev,nvert,
153	. 32, "ave(X)", zsto2,zout)
154	c
155	CALL histdef(nid_mth, "LWnet", "Net LW flux","W/m2",
156	. iim,jj_nb,nhori, klev,1,klev,nvert,
157	. 32, "ave(X)", zsto2,zout)
158	c
159	CALL histdef(nid_mth, "fluxvdf", "PBL net flux","W/m2",
160	. iim,jj_nb,nhori, klev,1,klev,nvert,
161	. 32, "ave(X)", zsto,zout)
162	c
163	CALL histdef(nid_mth, "fluxdyn", "Dyn. net flux","W/m2",
164	. iim,jj_nb,nhori, klev,1,klev,nvert,
165	. 32, "ave(X)", zsto,zout)
166	c
167	CALL histdef(nid_mth, "fluxajs", "Dry adj. net flux","W/m2",
168	. iim,jj_nb,nhori, klev,1,klev,nvert,
169	. 32, "ave(X)", zsto,zout)
170	c
171	c CALL histdef(nid_mth, "fluxec", "Cin. net flux","W/m2",
172	c . iim,jj_nb,nhori, klev,1,klev,nvert,
173	c . 32, "ave(X)", zsto,zout)
174	c
175	ENDIF !lev_histmth.GE.3
176	c
177	c-------------------------------------------------------
178	IF(lev_histmth.GE.4) THEN
179	c
180	CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
181	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
182	. "ave(X)", zsto,zout)
183	c
184	c CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
185	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
186	c . "ave(X)", zsto,zout)
187	c
188	CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s",
189	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
190	. "ave(X)", zsto,zout)
191	c
192	CALL histdef(nid_mth, "dtajs", "Dry adjust. dT", "K/s",
193	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
194	. "ave(X)", zsto,zout)
195	c
196	CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s",
197	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
198	. "ave(X)", zsto,zout)
199	c
200	CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
201	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
202	. "ave(X)", zsto,zout)
203	c
204	c CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
205	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
206	c . "ave(X)", zsto,zout)
207	c
208	CALL histdef(nid_mth, "duajs", "Dry convection dU", "m/s2",
209	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
210	. "ave(X)", zsto,zout)
211	c
212	CALL histdef(nid_mth, "dugwo", "GW oro dU", "m/s2",
213	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
214	. "ave(X)", zsto,zout)
215	c
216	CALL histdef(nid_mth, "dugwno", "GW non-oro dU", "m/s2",
217	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
218	. "ave(X)", zsto,zout)
219	c
220	c CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
221	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
222	c . "ave(X)", zsto,zout)
223	c
224	ENDIF !lev_histmth.GE.4
225	c
226	c-------------------------------------------------------
227	IF(lev_histmth.GE.5) THEN
228	c
229	c call histdef(nid_mth, "taux",
230	c $ "Zonal wind stress", "Pa",
231	c $ iim,jj_nb,nhori, 1,1,1, nvert, 32,
232	c $ "ave(X)", zsto,zout)
233	c
234	c call histdef(nid_mth, "tauy",
235	c $ "Meridional xind stress", "Pa",
236	c $ iim,jj_nb,nhori, 1,1,1, nvert, 32,
237	c $ "ave(X)", zsto,zout)
238	c
239	c CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
240	c . iim,jj_nb,nhori, 1,1,1, nvert, 32,
241	c . "ave(X)", zsto,zout)
242	c
243	c CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
244	c . iim,jj_nb,nhori, 1,1,1, nvert, 32,
245	c . "ave(X)", zsto,zout)
246	c
247	ENDIF !lev_histmth.GE.5
248	c-------------------------------------------------------
249
250	CALL histend(nid_mth)
251	!$OMP END MASTER
252
253	ENDIF ! fin de test sur ok_mensuel

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

Download in other formats: