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

Last change on this file since 1053 was 1053, checked in by hoffmann, 11 years ago

two-moment cloud physics implemented

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