Context Navigation

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

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

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

Download in other formats: