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

Last change on this file since 1320 was 1320, checked in by raasch, 7 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.