Home

Context Navigation

source: palm/trunk/SOURCE/average_3d_data.f90 @ 1322

Last change on this file since 1322 was 1322, checked in by raasch, 10 years ago
REAL functions and a lot of REAL constants provided with KIND-attribute, some small bugfixes
Property svn:keywords set to `Id`
File size: 10.9 KB

Line
1	SUBROUTINE average_3d_data
2
3	!--------------------------------------------------------------------------------!
4	! This file is part of PALM.
5	!
6	! PALM is free software: you can redistribute it and/or modify it under the terms
7	! of the GNU General Public License as published by the Free Software Foundation,
8	! either version 3 of the License, or (at your option) any later version.
9	!
10	! PALM is distributed in the hope that it will be useful, but WITHOUT ANY
11	! WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
12	! A PARTICULAR PURPOSE. See the GNU General Public License for more details.
13	!
14	! You should have received a copy of the GNU General Public License along with
15	! PALM. If not, see <http://www.gnu.org/licenses/>.
16	!
17	! Copyright 1997-2014 Leibniz Universitaet Hannover
18	!--------------------------------------------------------------------------------!
19	!
20	! Current revisions:
21	! -----------------
22	! REAL functions provided with KIND-attribute
23	!
24	! Former revisions:
25	! -----------------
26	! $Id: average_3d_data.f90 1322 2014-03-20 16:38:49Z raasch $
27	!
28	! 1320 2014-03-20 08:40:49Z raasch
29	! ONLY-attribute added to USE-statements,
30	! kind-parameters added to all INTEGER and REAL declaration statements,
31	! kinds are defined in new module kinds,
32	! revision history before 2012 removed,
33	! comment fields (!:) to be used for variable explanations added to
34	! all variable declaration statements
35	!
36	! 1318 2014-03-17 13:35:16Z raasch
37	! barrier argument removed from cpu_log,
38	! module interfaces removed
39	!
40	! 1115 2013-03-26 18:16:16Z hoffmann
41	! +qc
42	!
43	! 1053 2012-11-13 17:11:03Z hoffmann
44	! averaging of nr, qr added
45	!
46	! 1036 2012-10-22 13:43:42Z raasch
47	! code put under GPL (PALM 3.9)
48	!
49	! 978 2012-08-09 08:28:32Z fricke
50	! +z0h_av
51	!
52	! Revision 1.1 2006/02/23 09:48:58 raasch
53	! Initial revision
54	!
55	!
56	! Description:
57	! ------------
58	! Time-averaging of 3d-data-arrays.
59	!------------------------------------------------------------------------------!
60
61	USE averaging
62
63	USE control_parameters, &
64	ONLY: average_count_3d, doav, doav_n
65
66	USE cpulog, &
67	ONLY: cpu_log, log_point
68
69	USE indices, &
70	ONLY: nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, nzt
71
72	USE kinds
73
74
75	IMPLICIT NONE
76
77	INTEGER(iwp) :: i !:
78	INTEGER(iwp) :: ii !:
79	INTEGER(iwp) :: j !:
80	INTEGER(iwp) :: k !:
81
82
83	CALL cpu_log (log_point(35),'average_3d_data','start')
84
85	!
86	!-- Check, if averaging is necessary
87	IF ( average_count_3d <= 1 ) RETURN
88
89	!
90	!-- Loop of all variables to be averaged.
91	DO ii = 1, doav_n
92
93	!
94	!-- Store the array chosen on the temporary array.
95	SELECT CASE ( TRIM( doav(ii) ) )
96
97	CASE ( 'e' )
98	DO i = nxlg, nxrg
99	DO j = nysg, nyng
100	DO k = nzb, nzt+1
101	e_av(k,j,i) = e_av(k,j,i) / REAL( average_count_3d, KIND=wp )
102	ENDDO
103	ENDDO
104	ENDDO
105
106	CASE ( 'qsws*' )
107	DO i = nxlg, nxrg
108	DO j = nysg, nyng
109	qsws_av(j,i) = qsws_av(j,i) / REAL( average_count_3d, KIND=wp )
110	ENDDO
111	ENDDO
112
113	CASE ( 'lpt' )
114	DO i = nxlg, nxrg
115	DO j = nysg, nyng
116	DO k = nzb, nzt+1
117	lpt_av(k,j,i) = lpt_av(k,j,i) / REAL( average_count_3d, KIND=wp )
118	ENDDO
119	ENDDO
120	ENDDO
121
122	CASE ( 'lwp*' )
123	DO i = nxlg, nxrg
124	DO j = nysg, nyng
125	lwp_av(j,i) = lwp_av(j,i) / REAL( average_count_3d, KIND=wp )
126	ENDDO
127	ENDDO
128
129	CASE ( 'nr' )
130	DO i = nxlg, nxrg
131	DO j = nysg, nyng
132	DO k = nzb, nzt+1
133	nr_av(k,j,i) = nr_av(k,j,i) / REAL( average_count_3d, KIND=wp )
134	ENDDO
135	ENDDO
136	ENDDO
137
138	CASE ( 'p' )
139	DO i = nxlg, nxrg
140	DO j = nysg, nyng
141	DO k = nzb, nzt+1
142	p_av(k,j,i) = p_av(k,j,i) / REAL( average_count_3d, KIND=wp )
143	ENDDO
144	ENDDO
145	ENDDO
146
147	CASE ( 'pc' )
148	DO i = nxl, nxr
149	DO j = nys, nyn
150	DO k = nzb, nzt+1
151	pc_av(k,j,i) = pc_av(k,j,i) / REAL( average_count_3d, KIND=wp )
152	ENDDO
153	ENDDO
154	ENDDO
155
156	CASE ( 'pr' )
157	DO i = nxl, nxr
158	DO j = nys, nyn
159	DO k = nzb, nzt+1
160	pr_av(k,j,i) = pr_av(k,j,i) / REAL( average_count_3d, KIND=wp )
161	ENDDO
162	ENDDO
163	ENDDO
164
165	CASE ( 'prr*' )
166	DO i = nxlg, nxrg
167	DO j = nysg, nyng
168	precipitation_rate_av(j,i) = precipitation_rate_av(j,i) / &
169	REAL( average_count_3d, KIND=wp )
170	ENDDO
171	ENDDO
172
173	CASE ( 'pt' )
174	DO i = nxlg, nxrg
175	DO j = nysg, nyng
176	DO k = nzb, nzt+1
177	pt_av(k,j,i) = pt_av(k,j,i) / REAL( average_count_3d, KIND=wp )
178	ENDDO
179	ENDDO
180	ENDDO
181
182	CASE ( 'q' )
183	DO i = nxlg, nxrg
184	DO j = nysg, nyng
185	DO k = nzb, nzt+1
186	q_av(k,j,i) = q_av(k,j,i) / REAL( average_count_3d, KIND=wp )
187	ENDDO
188	ENDDO
189	ENDDO
190
191	CASE ( 'qc' )
192	DO i = nxlg, nxrg
193	DO j = nysg, nyng
194	DO k = nzb, nzt+1
195	qc_av(k,j,i) = qc_av(k,j,i) / REAL( average_count_3d, KIND=wp )
196	ENDDO
197	ENDDO
198	ENDDO
199
200	CASE ( 'ql' )
201	DO i = nxlg, nxrg
202	DO j = nysg, nyng
203	DO k = nzb, nzt+1
204	ql_av(k,j,i) = ql_av(k,j,i) / REAL( average_count_3d, KIND=wp )
205	ENDDO
206	ENDDO
207	ENDDO
208
209	CASE ( 'ql_c' )
210	DO i = nxlg, nxrg
211	DO j = nysg, nyng
212	DO k = nzb, nzt+1
213	ql_c_av(k,j,i) = ql_c_av(k,j,i) / REAL( average_count_3d, KIND=wp )
214	ENDDO
215	ENDDO
216	ENDDO
217
218	CASE ( 'ql_v' )
219	DO i = nxlg, nxrg
220	DO j = nysg, nyng
221	DO k = nzb, nzt+1
222	ql_v_av(k,j,i) = ql_v_av(k,j,i) / REAL( average_count_3d, KIND=wp )
223	ENDDO
224	ENDDO
225	ENDDO
226
227	CASE ( 'ql_vp' )
228	DO i = nxlg, nxrg
229	DO j = nysg, nyng
230	DO k = nzb, nzt+1
231	ql_vp_av(k,j,i) = ql_vp_av(k,j,i) / &
232	REAL( average_count_3d, KIND=wp )
233	ENDDO
234	ENDDO
235	ENDDO
236
237	CASE ( 'qr' )
238	DO i = nxlg, nxrg
239	DO j = nysg, nyng
240	DO k = nzb, nzt+1
241	qr_av(k,j,i) = qr_av(k,j,i) / REAL( average_count_3d, KIND=wp )
242	ENDDO
243	ENDDO
244	ENDDO
245
246	CASE ( 'qv' )
247	DO i = nxlg, nxrg
248	DO j = nysg, nyng
249	DO k = nzb, nzt+1
250	qv_av(k,j,i) = qv_av(k,j,i) / REAL( average_count_3d, KIND=wp )
251	ENDDO
252	ENDDO
253	ENDDO
254
255	CASE ( 'rho' )
256	DO i = nxlg, nxrg
257	DO j = nysg, nyng
258	DO k = nzb, nzt+1
259	rho_av(k,j,i) = rho_av(k,j,i) / REAL( average_count_3d, KIND=wp )
260	ENDDO
261	ENDDO
262	ENDDO
263
264	CASE ( 's' )
265	DO i = nxlg, nxrg
266	DO j = nysg, nyng
267	DO k = nzb, nzt+1
268	s_av(k,j,i) = s_av(k,j,i) / REAL( average_count_3d, KIND=wp )
269	ENDDO
270	ENDDO
271	ENDDO
272
273	CASE ( 'sa' )
274	DO i = nxlg, nxrg
275	DO j = nysg, nyng
276	DO k = nzb, nzt+1
277	sa_av(k,j,i) = sa_av(k,j,i) / REAL( average_count_3d, KIND=wp )
278	ENDDO
279	ENDDO
280	ENDDO
281
282	CASE ( 'shf*' )
283	DO i = nxlg, nxrg
284	DO j = nysg, nyng
285	shf_av(j,i) = shf_av(j,i) / REAL( average_count_3d, KIND=wp )
286	ENDDO
287	ENDDO
288
289	CASE ( 't*' )
290	DO i = nxlg, nxrg
291	DO j = nysg, nyng
292	ts_av(j,i) = ts_av(j,i) / REAL( average_count_3d, KIND=wp )
293	ENDDO
294	ENDDO
295
296	CASE ( 'u' )
297	DO i = nxlg, nxrg
298	DO j = nysg, nyng
299	DO k = nzb, nzt+1
300	u_av(k,j,i) = u_av(k,j,i) / REAL( average_count_3d, KIND=wp )
301	ENDDO
302	ENDDO
303	ENDDO
304
305	CASE ( 'u*' )
306	DO i = nxlg, nxrg
307	DO j = nysg, nyng
308	us_av(j,i) = us_av(j,i) / REAL( average_count_3d, KIND=wp )
309	ENDDO
310	ENDDO
311
312	CASE ( 'v' )
313	DO i = nxlg, nxrg
314	DO j = nysg, nyng
315	DO k = nzb, nzt+1
316	v_av(k,j,i) = v_av(k,j,i) / REAL( average_count_3d, KIND=wp )
317	ENDDO
318	ENDDO
319	ENDDO
320
321	CASE ( 'vpt' )
322	DO i = nxlg, nxrg
323	DO j = nysg, nyng
324	DO k = nzb, nzt+1
325	vpt_av(k,j,i) = vpt_av(k,j,i) / REAL( average_count_3d, KIND=wp )
326	ENDDO
327	ENDDO
328	ENDDO
329
330	CASE ( 'w' )
331	DO i = nxlg, nxrg
332	DO j = nysg, nyng
333	DO k = nzb, nzt+1
334	w_av(k,j,i) = w_av(k,j,i) / REAL( average_count_3d, KIND=wp )
335	ENDDO
336	ENDDO
337	ENDDO
338
339	CASE ( 'z0*' )
340	DO i = nxlg, nxrg
341	DO j = nysg, nyng
342	z0_av(j,i) = z0_av(j,i) / REAL( average_count_3d, KIND=wp )
343	ENDDO
344	ENDDO
345
346	CASE ( 'z0h*' )
347	DO i = nxlg, nxrg
348	DO j = nysg, nyng
349	z0h_av(j,i) = z0h_av(j,i) / REAL( average_count_3d, KIND=wp )
350	ENDDO
351	ENDDO
352
353	CASE DEFAULT
354	!
355	!-- User-defined quantity
356	CALL user_3d_data_averaging( 'average', doav(ii) )
357
358	END SELECT
359
360	ENDDO
361
362	!
363	!-- Reset the counter
364	average_count_3d = 0.0
365
366	CALL cpu_log( log_point(35), 'average_3d_data', 'stop' )
367
368
369	END SUBROUTINE average_3d_data

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |