Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |