Context Navigation

source: trunk/LMDZ.VENUS/libf/phyvenus/ini_histday.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: 7.7 KB

Line
1	!
2	! $Header: /home/cvsroot/LMDZ4/libf/phylmd/ini_histday.h,v 1.2 2005/01/27 10:06:12 fairhead Exp $
3	!
4	IF (ok_journe) THEN
5
6	zsto = dtime
7	zout = dtime * REAL(ecrit_day)
8	zsto1= dtime * REAL(ecrit_day)
9
10	idayref = day_ref
11	CALL ymds2ju(annee_ref, 1, idayref, zero, zjulian)
12
13	call histbeg_phy("histday.nc",itau_phy,
14	. zjulian,zsto,nhori,nid_day)
15
16	!$OMP MASTER
17	CALL histvert(nid_day, "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_day, "phis", "Surface geop. height", "-",
26	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
27	. "once", zsto,zout)
28	c
29	CALL histdef(nid_day, "aire", "Grid area", "-",
30	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
31	. "once", zsto,zout)
32	c
33	CALL histdef(nid_day, "tsol", "Surface Temperature", "K",
34	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
35	. "ave(X)", zsto,zout)
36	c
37	CALL histdef(nid_day, "psol", "Surface Pressure", "Pa",
38	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
39	. "ave(X)", zsto,zout)
40	c
41	CALL histdef(nid_day, "ue", "Zonal energy transport", "-",
42	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
43	. "ave(X)", zsto,zout)
44	c
45	CALL histdef(nid_day, "ve", "Merid energy transport", "-",
46	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
47	. "ave(X)", zsto,zout)
48	c
49	ENDIF !lev_histday.GE.1
50	c
51	c-------------------------------------------------------
52	IF(lev_histday.GE.2) THEN
53	c
54	ccccccccccccc 3D fields, basics
55	c
56	CALL histdef(nid_day, "temp", "Air temperature", "K",
57	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
58	. "ave(X)", zsto,zout)
59	c
60	CALL histdef(nid_day, "pres", "Air pressure", "Pa",
61	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
62	. "ave(X)", zsto,zout)
63	c
64	CALL histdef(nid_day, "geop", "Geopotential height", "m",
65	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
66	. "ave(X)", zsto,zout)
67	c
68	CALL histdef(nid_day, "vitu", "Zonal wind", "m/s",
69	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
70	. "ave(X)", zsto,zout)
71	c
72	CALL histdef(nid_day, "vitv", "Meridional wind", "m/s",
73	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
74	. "ave(X)", zsto,zout)
75	c
76	CALL histdef(nid_day, "vitw", "Vertical wind", "Pa/s",
77	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
78	. "ave(X)", zsto,zout)
79	c
80	ENDIF !lev_histday.GE.2
81	c
82	c-------------------------------------------------------
83	IF(lev_histday.GE.3) THEN
84	c
85	cccccccccccccccccc Radiative transfer
86	c
87	c 2D
88	c
89	CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",
90	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
91	. "ave(X)", zsto1,zout)
92	c
93	CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2",
94	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
95	. "ave(X)", zsto1,zout)
96	c
97	CALL histdef(nid_day, "sols", "Solar rad. at surf.", "W/m2",
98	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
99	. "ave(X)", zsto1,zout)
100	c
101	CALL histdef(nid_day, "soll", "IR rad. at surface", "W/m2",
102	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
103	. "ave(X)", zsto1,zout)
104	c
105	c 3D
106	c
107	CALL histdef(nid_day, "SWnet", "Net SW flux","W/m2",
108	. iim,jj_nb,nhori, klev,1,klev,nvert,
109	. 32, "ave(X)", zsto1,zout)
110	c
111	CALL histdef(nid_day, "LWnet", "Net LW flux","W/m2",
112	. iim,jj_nb,nhori, klev,1,klev,nvert,
113	. 32, "ave(X)", zsto1,zout)
114	c
115	ENDIF !lev_histday.GE.3
116	c
117	c-------------------------------------------------------
118	IF(lev_histday.GE.4) THEN
119	c
120	CALL histdef(nid_day, "dtdyn", "Dynamics dT", "K/s",
121	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
122	. "ave(X)", zsto,zout)
123	c
124	c CALL histdef(nid_day, "dtphy", "Physics dT", "K/s",
125	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
126	c . "ave(X)", zsto,zout)
127	c
128	CALL histdef(nid_day, "dtvdf", "Boundary-layer dT", "K/s",
129	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
130	. "ave(X)", zsto,zout)
131	c
132	CALL histdef(nid_day, "dtajs", "Dry adjust. dT", "K/s",
133	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
134	. "ave(X)", zsto,zout)
135	c
136	CALL histdef(nid_day, "dtswr", "SW radiation dT", "K/s",
137	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
138	. "ave(X)", zsto,zout)
139	c
140	CALL histdef(nid_day, "dtlwr", "LW radiation dT", "K/s",
141	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
142	. "ave(X)", zsto,zout)
143	c
144	c CALL histdef(nid_day, "dtec", "Cinetic dissip dT", "K/s",
145	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
146	c . "ave(X)", zsto,zout)
147	c
148	CALL histdef(nid_mth, "dudyn", "Dynamics dU", "m/s2",
149	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
150	. "ave(X)", zsto,zout)
151	c
152	CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
153	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
154	. "ave(X)", zsto,zout)
155	c
156	CALL histdef(nid_mth, "duajs", "Dry convection dU", "m/s2",
157	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
158	. "ave(X)", zsto,zout)
159	c
160	CALL histdef(nid_mth, "dugwo", "GW oro dU", "m/s2",
161	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
162	. "ave(X)", zsto,zout)
163	c
164	CALL histdef(nid_mth, "dugwno", "GW non-oro dU", "m/s2",
165	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
166	. "ave(X)", zsto,zout)
167	c
168	c CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
169	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
170	c . "ave(X)", zsto,zout)
171	c
172	ENDIF !lev_histday.GE.4
173	c
174	c-------------------------------------------------------
175	IF(lev_histday.GE.5) THEN
176	c
177	if (nqmax.GE.3) THEN
178	DO iq=1,nqmax-2
179	IF (iq.LE.99) THEN
180	WRITE(str2,'(i2.2)') iq
181	CALL histdef(nid_day, "trac"//str2, "Tracer No."//str2, "-",
182	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
183	. "ave(X)", zsto,zout)
184	ELSE
185	PRINT*, "Trop de traceurs"
186	CALL abort
187	ENDIF
188	ENDDO
189	endif
190	c
191	c call histdef(nid_day, "taux",
192	c $ "Zonal wind stress", "Pa",
193	c $ iim,jj_nb,nhori, 1,1,1, nvert, 32,
194	c $ "ave(X)", zsto,zout)
195	c
196	c call histdef(nid_day, "tauy",
197	c $ "Meridional xind stress", "Pa",
198	c $ iim,jj_nb,nhori, 1,1,1, nvert, 32,
199	c $ "ave(X)", zsto,zout)
200	c
201	c CALL histdef(nid_day, "cdrm", "Momentum drag coef.", "-",
202	c . iim,jj_nb,nhori, 1,1,1, nvert, 32,
203	c . "ave(X)", zsto,zout)
204	c
205	c CALL histdef(nid_day, "cdrh", "Heat drag coef.", "-",
206	c . iim,jj_nb,nhori, 1,1,1, nvert, 32,
207	c . "ave(X)", zsto,zout)
208	c
209	ENDIF !lev_histday.GE.5
210	c-------------------------------------------------------
211
212	CALL histend(nid_day)
213	!$OMP END MASTER
214
215	ENDIF ! fin de test sur ok_journe

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

Download in other formats: