Home

Context Navigation

source: palm/trunk/SOURCE/sum_up_3d_data.f90 @ 3592

Last change on this file since 3592 was 3589, checked in by suehring, 5 years ago
Remove erroneous UTF encoding; last commit documented
Property svn:keywords set to `Id`
File size: 46.9 KB

Line
1	!> @file sum_up_3d_data.f90
2	!------------------------------------------------------------------------------!
3	! This file is part of the PALM model system.
4	!
5	! PALM is free software: you can redistribute it and/or modify it under the
6	! terms of the GNU General Public License as published by the Free Software
7	! Foundation, either version 3 of the License, or (at your option) any later
8	! 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-2018 Leibniz Universitaet Hannover
18	!------------------------------------------------------------------------------!
19	!
20	! Current revisions:
21	! ------------------
22	!
23	!
24	! Former revisions:
25	! -----------------
26	! $Id: sum_up_3d_data.f90 3589 2018-11-30 15:09:51Z suehring $
27	! Move the control parameter "salsa" from salsa_mod to control_parameters
28	! (M. Kurppa)
29	!
30	! 3582 2018-11-29 19:16:36Z suehring
31	! dom_dwd_user, Schrempf:
32	! Remove CALLs to uv exposure model, this is now part of biometeorology_mod
33	!
34	! 3553 2018-11-22 10:30:48Z suehring
35	! variables documented
36	!
37	! 3552 2018-11-22 10:28:35Z suehring
38	! Changes related to clean-up of biometeorology (dom_dwd_user)
39	!
40	! 3467 2018-10-30 19:05:21Z suehring
41	! Implementation of a new aerosol module salsa.
42	!
43	! 3448 2018-10-29 18:14:31Z kanani
44	! Adjustment of biometeorology calls
45	!
46	! 3421 2018-10-24 18:39:32Z gronemeier
47	! Renamed output variables
48	!
49	! 3337 2018-10-12 15:17:09Z kanani
50	! (from branch resler)
51	! Add biometeorology,
52	! fix chemistry output call,
53	! move usm calls
54	!
55	! 3294 2018-10-01 02:37:10Z raasch
56	! changes concerning modularization of ocean option
57	!
58	! 3291 2018-09-28 11:33:03Z scharf
59	! corrected previous commit for 3D topography
60	!
61	! 3285 2018-09-27 17:16:52Z scharf
62	! bugfix for shf_av and qsws_av
63	!
64	! 3274 2018-09-24 15:42:55Z knoop
65	! Modularization of all bulk cloud physics code components
66	!
67	! 3241 2018-09-12 15:02:00Z raasch
68	! unused variables removed
69	!
70	! 3176 2018-07-26 17:12:48Z suehring
71	! Remove output of latent heat flux at urban-surfaces and set fill values
72	! instead
73	!
74	! 3173 2018-07-26 12:55:23Z suehring
75	! Bugfix for last commit
76	!
77	! 3170 2018-07-25 15:19:37Z suehring
78	! Revise output of surface quantities in case of overhanging structures
79	!
80	! 3151 2018-07-19 08:45:38Z raasch
81	! Remaining preprocessor directive __chem removed
82	!
83	! 3004 2018-04-27 12:33:25Z Giersch
84	! prr field added to ONLY-list, prr* case/pr* case/precipitation_rate_av
85	! removed, further allocation checks implemented
86	!
87	! 2963 2018-04-12 14:47:44Z suehring
88	! Introduce index for vegetation/wall, pavement/green-wall and water/window
89	! surfaces, for clearer access of surface fraction, albedo, emissivity, etc. .
90	!
91	! 2894 2018-03-15 09:17:58Z Giersch
92	! Changed comment
93	!
94	! 2817 2018-02-19 16:32:21Z suehring
95	! Preliminary gust module interface implemented
96	!
97	! 2798 2018-02-09 17:16:39Z suehring
98	! Consider also default-type surfaces for surface temperature output.
99	!
100	! 2797 2018-02-08 13:24:35Z suehring
101	! Enable output of ground-heat flux also at urban surfaces.
102	!
103	! 2790 2018-02-06 11:57:19Z suehring
104	! Bugfix in summation of surface sensible and latent heat flux
105	!
106	! 2766 2018-01-22 17:17:47Z kanani
107	! Removed preprocessor directive __chem
108	!
109	! 2743 2018-01-12 16:03:39Z suehring
110	! In case of natural- and urban-type surfaces output surfaces fluxes in W/m2.
111	!
112	! 2742 2018-01-12 14:59:47Z suehring
113	! Enable output of surface temperature
114	!
115	! 2735 2018-01-11 12:01:27Z suehring
116	! output of r_a moved from land-surface to consider also urban-type surfaces
117	!
118	! 2718 2018-01-02 08:49:38Z maronga
119	! Corrected "Former revisions" section
120	!
121	! 2696 2017-12-14 17:12:51Z kanani
122	! - Change in file header (GPL part)
123	! - Implementation of uv exposure model (FK)
124	! - output of diss_av, kh_av, km_av (turbulence_closure_mod) (TG)
125	! - Implementation of chemistry module (FK)
126	! - Workaround for sum-up usm arrays in case of restart runs, to avoid program
127	! crash (MS)
128	!
129	! 2292 2017-06-20 09:51:42Z schwenkel
130	! Implementation of new microphysic scheme: cloud_scheme = 'morrison'
131	! includes two more prognostic equations for cloud drop concentration (nc)
132	! and cloud water content (qc).
133	!
134	! 2233 2017-05-30 18:08:54Z suehring
135	!
136	! 2232 2017-05-30 17:47:52Z suehring
137	! Adjustments to new surface concept
138	!
139	! 2031 2016-10-21 15:11:58Z knoop
140	! renamed variable rho to rho_ocean and rho_av to rho_ocean_av
141	!
142	! 2024 2016-10-12 16:42:37Z kanani
143	! Added missing CASE for ssws*
144	!
145	! 2011 2016-09-19 17:29:57Z kanani
146	! Flag urban_surface is now defined in module control_parameters,
147	! changed prefix for urban surface model output to "usm_",
148	! introduced control parameter varnamelength for LEN of trimvar.
149	!
150	! 2007 2016-08-24 15:47:17Z kanani
151	! Added support for new urban surface model (temporary modifications of
152	! SELECT CASE ( ) necessary, see variable trimvar),
153	! added comments in variable declaration section
154	!
155	! 2000 2016-08-20 18:09:15Z knoop
156	! Forced header and separation lines into 80 columns
157	!
158	! 1992 2016-08-12 15:14:59Z suehring
159	! Bugfix in summation of passive scalar
160	!
161	! 1976 2016-07-27 13:28:04Z maronga
162	! Radiation actions are now done directly in the respective module
163	!
164	! 1972 2016-07-26 07:52:02Z maronga
165	! Land surface actions are now done directly in the respective module
166	!
167	! 1960 2016-07-12 16:34:24Z suehring
168	! Scalar surface flux added
169	!
170	! 1949 2016-06-17 07:19:16Z maronga
171	! Bugfix: calculation of lai_av, c_veg_av and c_liq_av.
172	!
173	! 1849 2016-04-08 11:33:18Z hoffmann
174	! precipitation_rate moved to arrays_3d
175	!
176	! 1788 2016-03-10 11:01:04Z maronga
177	! Added z0q and z0q_av
178	!
179	! 1693 2015-10-27 08:35:45Z maronga
180	! Last revision text corrected
181	!
182	! 1691 2015-10-26 16:17:44Z maronga
183	! Added output of Obukhov length and radiative heating rates for RRTMG.
184	! Corrected output of liquid water path.
185	!
186	! 1682 2015-10-07 23:56:08Z knoop
187	! Code annotations made doxygen readable
188	!
189	! 1585 2015-04-30 07:05:52Z maronga
190	! Adapted for RRTMG
191	!
192	! 1555 2015-03-04 17:44:27Z maronga
193	! Added output of r_a and r_s
194	!
195	! 1551 2015-03-03 14:18:16Z maronga
196	! Added support for land surface model and radiation model data.
197	!
198	! 1359 2014-04-11 17:15:14Z hoffmann
199	! New particle structure integrated.
200	!
201	! 1353 2014-04-08 15:21:23Z heinze
202	! REAL constants provided with KIND-attribute
203	!
204	! 1320 2014-03-20 08:40:49Z raasch
205	! ONLY-attribute added to USE-statements,
206	! kind-parameters added to all INTEGER and REAL declaration statements,
207	! kinds are defined in new module kinds,
208	! old module precision_kind is removed,
209	! revision history before 2012 removed,
210	! comment fields (!:) to be used for variable explanations added to
211	! all variable declaration statements
212	!
213	! 1318 2014-03-17 13:35:16Z raasch
214	! barrier argument removed from cpu_log,
215	! module interfaces removed
216	!
217	! 1115 2013-03-26 18:16:16Z hoffmann
218	! ql is calculated by calc_liquid_water_content
219	!
220	! 1053 2012-11-13 17:11:03Z hoffmann
221	! +nr, prr, qr
222	!
223	! 1036 2012-10-22 13:43:42Z raasch
224	! code put under GPL (PALM 3.9)
225	!
226	! 1007 2012-09-19 14:30:36Z franke
227	! Bugfix in calculation of ql_vp
228	!
229	! 978 2012-08-09 08:28:32Z fricke
230	! +z0h*
231	!
232	! Revision 1.1 2006/02/23 12:55:23 raasch
233	! Initial revision
234	!
235	!
236	! Description:
237	! ------------
238	!> Sum-up the values of 3d-arrays. The real averaging is later done in routine
239	!> average_3d_data.
240	!------------------------------------------------------------------------------!
241	SUBROUTINE sum_up_3d_data
242
243
244	USE arrays_3d, &
245	ONLY: dzw, d_exner, e, heatflux_output_conversion, nc, nr, p, prr, &
246	pt, q, qc, ql, ql_c, ql_v, qr, s, u, v, vpt, w, &
247	waterflux_output_conversion
248
249	USE averaging, &
250	ONLY: e_av, ghf_av, lpt_av, lwp_av, ol_av, p_av, pc_av, pr_av, pt_av, &
251	q_av, ql_av, ql_c_av, ql_v_av, ql_vp_av, qsws_av, qv_av, &
252	r_a_av, s_av, shf_av, ssws_av, ts_av, tsurf_av, u_av, us_av, &
253	v_av, vpt_av, w_av, z0_av, z0h_av, z0q_av
254
255	USE basic_constants_and_equations_mod, &
256	ONLY: c_p, lv_d_cp, l_v
257
258	USE biometeorology_mod, &
259	ONLY: bio_3d_data_averaging
260
261	USE bulk_cloud_model_mod, &
262	ONLY: bulk_cloud_model, bcm_3d_data_averaging
263
264	USE chemistry_model_mod, &
265	ONLY: chem_3d_data_averaging
266
267	USE control_parameters, &
268	ONLY: air_chemistry, average_count_3d, biometeorology, doav, doav_n, &
269	land_surface, ocean_mode, rho_surface, salsa, urban_surface, &
270	varnamelength
271
272	USE cpulog, &
273	ONLY: cpu_log, log_point
274
275	USE gust_mod, &
276	ONLY: gust_3d_data_averaging, gust_module_enabled
277
278	USE indices, &
279	ONLY: nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, nzt
280
281	USE kinds
282
283	USE land_surface_model_mod, &
284	ONLY: lsm_3d_data_averaging
285
286	USE ocean_mod, &
287	ONLY: ocean_3d_data_averaging
288
289	USE particle_attributes, &
290	ONLY: grid_particles, number_of_particles, particles, prt_count
291
292	USE radiation_model_mod, &
293	ONLY: radiation, radiation_3d_data_averaging
294
295	USE salsa_mod, &
296	ONLY: salsa_3d_data_averaging
297
298	USE surface_mod, &
299	ONLY: ind_pav_green, ind_veg_wall, ind_wat_win, &
300	surf_def_h, surf_lsm_h, surf_usm_h
301
302	USE turbulence_closure_mod, &
303	ONLY: tcm_3d_data_averaging
304
305	USE urban_surface_mod, &
306	ONLY: usm_average_3d_data
307
308
309	IMPLICIT NONE
310
311	LOGICAL :: match_def !< flag indicating default-type surface
312	LOGICAL :: match_lsm !< flag indicating natural-type surface
313	LOGICAL :: match_usm !< flag indicating urban-type surface
314
315	INTEGER(iwp) :: i !< grid index x direction
316	INTEGER(iwp) :: ii !< running index
317	INTEGER(iwp) :: j !< grid index y direction
318	INTEGER(iwp) :: k !< grid index x direction
319	INTEGER(iwp) :: m !< running index over surfacle elements
320	INTEGER(iwp) :: n !< running index over number of particles per grid box
321
322	REAL(wp) :: mean_r !< mean-particle radius witin grid box
323	REAL(wp) :: s_r2 !< mean-particle radius witin grid box to the power of two
324	REAL(wp) :: s_r3 !< mean-particle radius witin grid box to the power of three
325
326	CHARACTER (LEN=varnamelength) :: trimvar !< TRIM of output-variable string
327
328
329	CALL cpu_log (log_point(34),'sum_up_3d_data','start')
330
331	!
332	!-- Allocate and initialize the summation arrays if called for the very first
333	!-- time or the first time after average_3d_data has been called
334	!-- (some or all of the arrays may have been already allocated
335	!-- in rrd_local)
336	IF ( average_count_3d == 0 ) THEN
337
338	DO ii = 1, doav_n
339
340	trimvar = TRIM( doav(ii) )
341
342	SELECT CASE ( trimvar )
343
344	CASE ( 'ghf*' )
345	IF ( .NOT. ALLOCATED( ghf_av ) ) THEN
346	ALLOCATE( ghf_av(nysg:nyng,nxlg:nxrg) )
347	ENDIF
348	ghf_av = 0.0_wp
349
350	CASE ( 'e' )
351	IF ( .NOT. ALLOCATED( e_av ) ) THEN
352	ALLOCATE( e_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
353	ENDIF
354	e_av = 0.0_wp
355
356	CASE ( 'thetal' )
357	IF ( .NOT. ALLOCATED( lpt_av ) ) THEN
358	ALLOCATE( lpt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
359	ENDIF
360	lpt_av = 0.0_wp
361
362	CASE ( 'lwp*' )
363	IF ( .NOT. ALLOCATED( lwp_av ) ) THEN
364	ALLOCATE( lwp_av(nysg:nyng,nxlg:nxrg) )
365	ENDIF
366	lwp_av = 0.0_wp
367
368	CASE ( 'ol*' )
369	IF ( .NOT. ALLOCATED( ol_av ) ) THEN
370	ALLOCATE( ol_av(nysg:nyng,nxlg:nxrg) )
371	ENDIF
372	ol_av = 0.0_wp
373
374	CASE ( 'p' )
375	IF ( .NOT. ALLOCATED( p_av ) ) THEN
376	ALLOCATE( p_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
377	ENDIF
378	p_av = 0.0_wp
379
380	CASE ( 'pc' )
381	IF ( .NOT. ALLOCATED( pc_av ) ) THEN
382	ALLOCATE( pc_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
383	ENDIF
384	pc_av = 0.0_wp
385
386	CASE ( 'pr' )
387	IF ( .NOT. ALLOCATED( pr_av ) ) THEN
388	ALLOCATE( pr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
389	ENDIF
390	pr_av = 0.0_wp
391
392	CASE ( 'theta' )
393	IF ( .NOT. ALLOCATED( pt_av ) ) THEN
394	ALLOCATE( pt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
395	ENDIF
396	pt_av = 0.0_wp
397
398	CASE ( 'q' )
399	IF ( .NOT. ALLOCATED( q_av ) ) THEN
400	ALLOCATE( q_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
401	ENDIF
402	q_av = 0.0_wp
403
404	CASE ( 'ql' )
405	IF ( .NOT. ALLOCATED( ql_av ) ) THEN
406	ALLOCATE( ql_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
407	ENDIF
408	ql_av = 0.0_wp
409
410	CASE ( 'ql_c' )
411	IF ( .NOT. ALLOCATED( ql_c_av ) ) THEN
412	ALLOCATE( ql_c_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
413	ENDIF
414	ql_c_av = 0.0_wp
415
416	CASE ( 'ql_v' )
417	IF ( .NOT. ALLOCATED( ql_v_av ) ) THEN
418	ALLOCATE( ql_v_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
419	ENDIF
420	ql_v_av = 0.0_wp
421
422	CASE ( 'ql_vp' )
423	IF ( .NOT. ALLOCATED( ql_vp_av ) ) THEN
424	ALLOCATE( ql_vp_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
425	ENDIF
426	ql_vp_av = 0.0_wp
427
428	CASE ( 'qsws*' )
429	IF ( .NOT. ALLOCATED( qsws_av ) ) THEN
430	ALLOCATE( qsws_av(nysg:nyng,nxlg:nxrg) )
431	ENDIF
432	qsws_av = 0.0_wp
433
434	CASE ( 'qv' )
435	IF ( .NOT. ALLOCATED( qv_av ) ) THEN
436	ALLOCATE( qv_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
437	ENDIF
438	qv_av = 0.0_wp
439
440	CASE ( 'r_a*' )
441	IF ( .NOT. ALLOCATED( r_a_av ) ) THEN
442	ALLOCATE( r_a_av(nysg:nyng,nxlg:nxrg) )
443	ENDIF
444	r_a_av = 0.0_wp
445
446	CASE ( 's' )
447	IF ( .NOT. ALLOCATED( s_av ) ) THEN
448	ALLOCATE( s_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
449	ENDIF
450	s_av = 0.0_wp
451
452	CASE ( 'shf*' )
453	IF ( .NOT. ALLOCATED( shf_av ) ) THEN
454	ALLOCATE( shf_av(nysg:nyng,nxlg:nxrg) )
455	ENDIF
456	shf_av = 0.0_wp
457
458	CASE ( 'ssws*' )
459	IF ( .NOT. ALLOCATED( ssws_av ) ) THEN
460	ALLOCATE( ssws_av(nysg:nyng,nxlg:nxrg) )
461	ENDIF
462	ssws_av = 0.0_wp
463
464	CASE ( 't*' )
465	IF ( .NOT. ALLOCATED( ts_av ) ) THEN
466	ALLOCATE( ts_av(nysg:nyng,nxlg:nxrg) )
467	ENDIF
468	ts_av = 0.0_wp
469
470	CASE ( 'tsurf*' )
471	IF ( .NOT. ALLOCATED( tsurf_av ) ) THEN
472	ALLOCATE( tsurf_av(nysg:nyng,nxlg:nxrg) )
473	ENDIF
474	tsurf_av = 0.0_wp
475
476	CASE ( 'u' )
477	IF ( .NOT. ALLOCATED( u_av ) ) THEN
478	ALLOCATE( u_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
479	ENDIF
480	u_av = 0.0_wp
481
482	CASE ( 'us*' )
483	IF ( .NOT. ALLOCATED( us_av ) ) THEN
484	ALLOCATE( us_av(nysg:nyng,nxlg:nxrg) )
485	ENDIF
486	us_av = 0.0_wp
487
488	CASE ( 'v' )
489	IF ( .NOT. ALLOCATED( v_av ) ) THEN
490	ALLOCATE( v_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
491	ENDIF
492	v_av = 0.0_wp
493
494	CASE ( 'thetav' )
495	IF ( .NOT. ALLOCATED( vpt_av ) ) THEN
496	ALLOCATE( vpt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
497	ENDIF
498	vpt_av = 0.0_wp
499
500	CASE ( 'w' )
501	IF ( .NOT. ALLOCATED( w_av ) ) THEN
502	ALLOCATE( w_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
503	ENDIF
504	w_av = 0.0_wp
505
506	CASE ( 'z0*' )
507	IF ( .NOT. ALLOCATED( z0_av ) ) THEN
508	ALLOCATE( z0_av(nysg:nyng,nxlg:nxrg) )
509	ENDIF
510	z0_av = 0.0_wp
511
512	CASE ( 'z0h*' )
513	IF ( .NOT. ALLOCATED( z0h_av ) ) THEN
514	ALLOCATE( z0h_av(nysg:nyng,nxlg:nxrg) )
515	ENDIF
516	z0h_av = 0.0_wp
517
518	CASE ( 'z0q*' )
519	IF ( .NOT. ALLOCATED( z0q_av ) ) THEN
520	ALLOCATE( z0q_av(nysg:nyng,nxlg:nxrg) )
521	ENDIF
522	z0q_av = 0.0_wp
523
524
525	CASE DEFAULT
526
527	!
528	!-- Allocating and initializing data arrays for other modules
529
530	IF ( air_chemistry .AND. &
531	(trimvar(1:3) == 'kc_' .OR. trimvar(1:3) == 'em_') ) THEN
532	CALL chem_3d_data_averaging( 'allocate', doav(ii) )
533	ENDIF
534
535	IF ( bulk_cloud_model ) THEN
536	CALL bcm_3d_data_averaging( 'allocate', doav(ii) )
537	ENDIF
538
539	IF ( gust_module_enabled ) THEN
540	CALL gust_3d_data_averaging( 'allocate', doav(ii) )
541	ENDIF
542
543	IF ( biometeorology ) THEN
544	CALL bio_3d_data_averaging( 'allocate', doav(ii) )
545	ENDIF
546
547	IF ( land_surface ) THEN
548	CALL lsm_3d_data_averaging( 'allocate', doav(ii) )
549	ENDIF
550
551	IF ( ocean_mode ) THEN
552	CALL ocean_3d_data_averaging( 'allocate', doav(ii) )
553	ENDIF
554
555	IF ( radiation ) THEN
556	CALL radiation_3d_data_averaging( 'allocate', doav(ii) )
557	ENDIF
558
559	IF ( salsa ) THEN
560	CALL salsa_3d_data_averaging( 'allocate', doav(ii) )
561	ENDIF
562
563	CALL tcm_3d_data_averaging( 'allocate', doav(ii) )
564
565	IF ( urban_surface .AND. trimvar(1:4) == 'usm_' ) THEN
566	CALL usm_average_3d_data( 'allocate', doav(ii) )
567	ENDIF
568
569	!
570	!-- User-defined quantities
571	CALL user_3d_data_averaging( 'allocate', doav(ii) )
572
573	END SELECT
574
575	ENDDO
576
577	ENDIF
578
579	!
580	!-- Loop of all variables to be averaged.
581	DO ii = 1, doav_n
582
583	trimvar = TRIM( doav(ii) )
584	!
585	!-- Store the array chosen on the temporary array.
586	SELECT CASE ( trimvar )
587
588	CASE ( 'ghf*' )
589	IF ( ALLOCATED( ghf_av ) ) THEN
590	DO i = nxl, nxr
591	DO j = nys, nyn
592	!
593	!-- Check whether grid point is a natural- or urban-type
594	!-- surface.
595	match_lsm = surf_lsm_h%start_index(j,i) <= &
596	surf_lsm_h%end_index(j,i)
597	match_usm = surf_usm_h%start_index(j,i) <= &
598	surf_usm_h%end_index(j,i)
599	!
600	!-- In order to avoid double-counting of surface properties,
601	!-- always assume that natural-type surfaces are below urban-
602	!-- type surfaces, e.g. in case of bridges.
603	!-- Further, take only the last suface element, i.e. the
604	!-- uppermost surface which would be visible from above
605	IF ( match_lsm .AND. .NOT. match_usm ) THEN
606	m = surf_lsm_h%end_index(j,i)
607	ghf_av(j,i) = ghf_av(j,i) + &
608	surf_lsm_h%ghf(m)
609	ELSEIF ( match_usm ) THEN
610	m = surf_usm_h%end_index(j,i)
611	ghf_av(j,i) = ghf_av(j,i) + &
612	surf_usm_h%frac(ind_veg_wall,m) * &
613	surf_usm_h%wghf_eb(m) + &
614	surf_usm_h%frac(ind_pav_green,m) * &
615	surf_usm_h%wghf_eb_green(m) + &
616	surf_usm_h%frac(ind_wat_win,m) * &
617	surf_usm_h%wghf_eb_window(m)
618	ENDIF
619	ENDDO
620	ENDDO
621	ENDIF
622
623	CASE ( 'e' )
624	IF ( ALLOCATED( e_av ) ) THEN
625	DO i = nxlg, nxrg
626	DO j = nysg, nyng
627	DO k = nzb, nzt+1
628	e_av(k,j,i) = e_av(k,j,i) + e(k,j,i)
629	ENDDO
630	ENDDO
631	ENDDO
632	ENDIF
633
634	CASE ( 'thetal' )
635	IF ( ALLOCATED( lpt_av ) ) THEN
636	DO i = nxlg, nxrg
637	DO j = nysg, nyng
638	DO k = nzb, nzt+1
639	lpt_av(k,j,i) = lpt_av(k,j,i) + pt(k,j,i)
640	ENDDO
641	ENDDO
642	ENDDO
643	ENDIF
644
645	CASE ( 'lwp*' )
646	IF ( ALLOCATED( lwp_av ) ) THEN
647	DO i = nxlg, nxrg
648	DO j = nysg, nyng
649	lwp_av(j,i) = lwp_av(j,i) + SUM( ql(nzb:nzt,j,i) &
650	* dzw(1:nzt+1) ) * rho_surface
651	ENDDO
652	ENDDO
653	ENDIF
654
655	CASE ( 'ol*' )
656	IF ( ALLOCATED( ol_av ) ) THEN
657	DO i = nxl, nxr
658	DO j = nys, nyn
659	match_def = surf_def_h(0)%start_index(j,i) <= &
660	surf_def_h(0)%end_index(j,i)
661	match_lsm = surf_lsm_h%start_index(j,i) <= &
662	surf_lsm_h%end_index(j,i)
663	match_usm = surf_usm_h%start_index(j,i) <= &
664	surf_usm_h%end_index(j,i)
665
666	IF ( match_def ) THEN
667	m = surf_def_h(0)%end_index(j,i)
668	ol_av(j,i) = ol_av(j,i) + &
669	surf_def_h(0)%ol(m)
670	ELSEIF ( match_lsm .AND. .NOT. match_usm ) THEN
671	m = surf_lsm_h%end_index(j,i)
672	ol_av(j,i) = ol_av(j,i) + &
673	surf_lsm_h%ol(m)
674	ELSEIF ( match_usm ) THEN
675	m = surf_usm_h%end_index(j,i)
676	ol_av(j,i) = ol_av(j,i) + &
677	surf_usm_h%ol(m)
678	ENDIF
679	ENDDO
680	ENDDO
681	ENDIF
682
683	CASE ( 'p' )
684	IF ( ALLOCATED( p_av ) ) THEN
685	DO i = nxlg, nxrg
686	DO j = nysg, nyng
687	DO k = nzb, nzt+1
688	p_av(k,j,i) = p_av(k,j,i) + p(k,j,i)
689	ENDDO
690	ENDDO
691	ENDDO
692	ENDIF
693
694	CASE ( 'pc' )
695	IF ( ALLOCATED( pc_av ) ) THEN
696	DO i = nxl, nxr
697	DO j = nys, nyn
698	DO k = nzb, nzt+1
699	pc_av(k,j,i) = pc_av(k,j,i) + prt_count(k,j,i)
700	ENDDO
701	ENDDO
702	ENDDO
703	ENDIF
704
705	CASE ( 'pr' )
706	IF ( ALLOCATED( pr_av ) ) THEN
707	DO i = nxl, nxr
708	DO j = nys, nyn
709	DO k = nzb, nzt+1
710	number_of_particles = prt_count(k,j,i)
711	IF ( number_of_particles <= 0 ) CYCLE
712	particles => &
713	grid_particles(k,j,i)%particles(1:number_of_particles)
714	s_r2 = 0.0_wp
715	s_r3 = 0.0_wp
716
717	DO n = 1, number_of_particles
718	IF ( particles(n)%particle_mask ) THEN
719	s_r2 = s_r2 + particles(n)%radius*2 &
720	particles(n)%weight_factor
721	s_r3 = s_r3 + particles(n)%radius*3 &
722	particles(n)%weight_factor
723	ENDIF
724	ENDDO
725
726	IF ( s_r2 > 0.0_wp ) THEN
727	mean_r = s_r3 / s_r2
728	ELSE
729	mean_r = 0.0_wp
730	ENDIF
731	pr_av(k,j,i) = pr_av(k,j,i) + mean_r
732	ENDDO
733	ENDDO
734	ENDDO
735	ENDIF
736
737	CASE ( 'theta' )
738	IF ( ALLOCATED( pt_av ) ) THEN
739	IF ( .NOT. bulk_cloud_model ) THEN
740	DO i = nxlg, nxrg
741	DO j = nysg, nyng
742	DO k = nzb, nzt+1
743	pt_av(k,j,i) = pt_av(k,j,i) + pt(k,j,i)
744	ENDDO
745	ENDDO
746	ENDDO
747	ELSE
748	DO i = nxlg, nxrg
749	DO j = nysg, nyng
750	DO k = nzb, nzt+1
751	pt_av(k,j,i) = pt_av(k,j,i) + pt(k,j,i) + lv_d_cp * &
752	d_exner(k) * ql(k,j,i)
753	ENDDO
754	ENDDO
755	ENDDO
756	ENDIF
757	ENDIF
758
759	CASE ( 'q' )
760	IF ( ALLOCATED( q_av ) ) THEN
761	DO i = nxlg, nxrg
762	DO j = nysg, nyng
763	DO k = nzb, nzt+1
764	q_av(k,j,i) = q_av(k,j,i) + q(k,j,i)
765	ENDDO
766	ENDDO
767	ENDDO
768	ENDIF
769
770	CASE ( 'ql' )
771	IF ( ALLOCATED( ql_av ) ) THEN
772	DO i = nxlg, nxrg
773	DO j = nysg, nyng
774	DO k = nzb, nzt+1
775	ql_av(k,j,i) = ql_av(k,j,i) + ql(k,j,i)
776	ENDDO
777	ENDDO
778	ENDDO
779	ENDIF
780
781	CASE ( 'ql_c' )
782	IF ( ALLOCATED( ql_c_av ) ) THEN
783	DO i = nxlg, nxrg
784	DO j = nysg, nyng
785	DO k = nzb, nzt+1
786	ql_c_av(k,j,i) = ql_c_av(k,j,i) + ql_c(k,j,i)
787	ENDDO
788	ENDDO
789	ENDDO
790	ENDIF
791
792	CASE ( 'ql_v' )
793	IF ( ALLOCATED( ql_v_av ) ) THEN
794	DO i = nxlg, nxrg
795	DO j = nysg, nyng
796	DO k = nzb, nzt+1
797	ql_v_av(k,j,i) = ql_v_av(k,j,i) + ql_v(k,j,i)
798	ENDDO
799	ENDDO
800	ENDDO
801	ENDIF
802
803	CASE ( 'ql_vp' )
804	IF ( ALLOCATED( ql_vp_av ) ) THEN
805	DO i = nxl, nxr
806	DO j = nys, nyn
807	DO k = nzb, nzt+1
808	number_of_particles = prt_count(k,j,i)
809	IF ( number_of_particles <= 0 ) CYCLE
810	particles => &
811	grid_particles(k,j,i)%particles(1:number_of_particles)
812	DO n = 1, number_of_particles
813	IF ( particles(n)%particle_mask ) THEN
814	ql_vp_av(k,j,i) = ql_vp_av(k,j,i) + &
815	particles(n)%weight_factor / &
816	number_of_particles
817	ENDIF
818	ENDDO
819	ENDDO
820	ENDDO
821	ENDDO
822	ENDIF
823
824	CASE ( 'qsws*' )
825	!
826	!-- In case of default surfaces, clean-up flux by density.
827	!-- In case of land- and urban-surfaces, convert fluxes into
828	!-- dynamic units.
829	IF ( ALLOCATED( qsws_av ) ) THEN
830	DO i = nxl, nxr
831	DO j = nys, nyn
832	match_def = surf_def_h(0)%start_index(j,i) <= &
833	surf_def_h(0)%end_index(j,i)
834	match_lsm = surf_lsm_h%start_index(j,i) <= &
835	surf_lsm_h%end_index(j,i)
836	match_usm = surf_usm_h%start_index(j,i) <= &
837	surf_usm_h%end_index(j,i)
838
839	IF ( match_def ) THEN
840	m = surf_def_h(0)%end_index(j,i)
841	qsws_av(j,i) = qsws_av(j,i) + &
842	surf_def_h(0)%qsws(m) * &
843	waterflux_output_conversion(nzb)
844	ELSEIF ( match_lsm .AND. .NOT. match_usm ) THEN
845	m = surf_lsm_h%end_index(j,i)
846	qsws_av(j,i) = qsws_av(j,i) + &
847	surf_lsm_h%qsws(m) * l_v
848	ENDIF
849	ENDDO
850	ENDDO
851	ENDIF
852
853	CASE ( 'qv' )
854	IF ( ALLOCATED( qv_av ) ) THEN
855	DO i = nxlg, nxrg
856	DO j = nysg, nyng
857	DO k = nzb, nzt+1
858	qv_av(k,j,i) = qv_av(k,j,i) + q(k,j,i) - ql(k,j,i)
859	ENDDO
860	ENDDO
861	ENDDO
862	ENDIF
863
864	CASE ( 'r_a*' )
865	IF ( ALLOCATED( r_a_av ) ) THEN
866	DO i = nxl, nxr
867	DO j = nys, nyn
868	match_lsm = surf_lsm_h%start_index(j,i) <= &
869	surf_lsm_h%end_index(j,i)
870	match_usm = surf_usm_h%start_index(j,i) <= &
871	surf_usm_h%end_index(j,i)
872
873	IF ( match_lsm .AND. .NOT. match_usm ) THEN
874	m = surf_lsm_h%end_index(j,i)
875	r_a_av(j,i) = r_a_av(j,i) + &
876	surf_lsm_h%r_a(m)
877	ELSEIF ( match_usm ) THEN
878	m = surf_usm_h%end_index(j,i)
879	r_a_av(j,i) = r_a_av(j,i) + &
880	surf_usm_h%frac(ind_veg_wall,m) * &
881	surf_usm_h%r_a(m) + &
882	surf_usm_h%frac(ind_pav_green,m) * &
883	surf_usm_h%r_a_green(m) + &
884	surf_usm_h%frac(ind_wat_win,m) * &
885	surf_usm_h%r_a_window(m)
886	ENDIF
887	ENDDO
888	ENDDO
889	ENDIF
890
891	CASE ( 's' )
892	IF ( ALLOCATED( s_av ) ) THEN
893	DO i = nxlg, nxrg
894	DO j = nysg, nyng
895	DO k = nzb, nzt+1
896	s_av(k,j,i) = s_av(k,j,i) + s(k,j,i)
897	ENDDO
898	ENDDO
899	ENDDO
900	ENDIF
901
902	CASE ( 'shf*' )
903	!
904	!-- In case of default surfaces, clean-up flux by density.
905	!-- In case of land- and urban-surfaces, convert fluxes into
906	!-- dynamic units.
907	IF ( ALLOCATED( shf_av ) ) THEN
908	DO i = nxl, nxr
909	DO j = nys, nyn
910	match_def = surf_def_h(0)%start_index(j,i) <= &
911	surf_def_h(0)%end_index(j,i)
912	match_lsm = surf_lsm_h%start_index(j,i) <= &
913	surf_lsm_h%end_index(j,i)
914	match_usm = surf_usm_h%start_index(j,i) <= &
915	surf_usm_h%end_index(j,i)
916
917	IF ( match_def ) THEN
918	m = surf_def_h(0)%end_index(j,i)
919	shf_av(j,i) = shf_av(j,i) + &
920	surf_def_h(0)%shf(m) * &
921	heatflux_output_conversion(nzb)
922	ELSEIF ( match_lsm .AND. .NOT. match_usm ) THEN
923	m = surf_lsm_h%end_index(j,i)
924	shf_av(j,i) = shf_av(j,i) + &
925	surf_lsm_h%shf(m) * c_p
926	ELSEIF ( match_usm ) THEN
927	m = surf_usm_h%end_index(j,i)
928	shf_av(j,i) = shf_av(j,i) + &
929	surf_usm_h%shf(m) * c_p
930	ENDIF
931	ENDDO
932	ENDDO
933	ENDIF
934
935	CASE ( 'ssws*' )
936	IF ( ALLOCATED( ssws_av ) ) THEN
937	DO i = nxl, nxr
938	DO j = nys, nyn
939	match_def = surf_def_h(0)%start_index(j,i) <= &
940	surf_def_h(0)%end_index(j,i)
941	match_lsm = surf_lsm_h%start_index(j,i) <= &
942	surf_lsm_h%end_index(j,i)
943	match_usm = surf_usm_h%start_index(j,i) <= &
944	surf_usm_h%end_index(j,i)
945
946	IF ( match_def ) THEN
947	m = surf_def_h(0)%end_index(j,i)
948	ssws_av(j,i) = ssws_av(j,i) + &
949	surf_def_h(0)%ssws(m)
950	ELSEIF ( match_lsm .AND. .NOT. match_usm ) THEN
951	m = surf_lsm_h%end_index(j,i)
952	ssws_av(j,i) = ssws_av(j,i) + &
953	surf_lsm_h%ssws(m)
954	ELSEIF ( match_usm ) THEN
955	m = surf_usm_h%end_index(j,i)
956	ssws_av(j,i) = ssws_av(j,i) + &
957	surf_usm_h%ssws(m)
958	ENDIF
959	ENDDO
960	ENDDO
961	ENDIF
962
963	CASE ( 't*' )
964	IF ( ALLOCATED( ts_av ) ) THEN
965	DO i = nxl, nxr
966	DO j = nys, nyn
967	match_def = surf_def_h(0)%start_index(j,i) <= &
968	surf_def_h(0)%end_index(j,i)
969	match_lsm = surf_lsm_h%start_index(j,i) <= &
970	surf_lsm_h%end_index(j,i)
971	match_usm = surf_usm_h%start_index(j,i) <= &
972	surf_usm_h%end_index(j,i)
973
974	IF ( match_def ) THEN
975	m = surf_def_h(0)%end_index(j,i)
976	ts_av(j,i) = ts_av(j,i) + &
977	surf_def_h(0)%ts(m)
978	ELSEIF ( match_lsm .AND. .NOT. match_usm ) THEN
979	m = surf_lsm_h%end_index(j,i)
980	ts_av(j,i) = ts_av(j,i) + &
981	surf_lsm_h%ts(m)
982	ELSEIF ( match_usm ) THEN
983	m = surf_usm_h%end_index(j,i)
984	ts_av(j,i) = ts_av(j,i) + &
985	surf_usm_h%ts(m)
986	ENDIF
987	ENDDO
988	ENDDO
989	ENDIF
990
991	CASE ( 'tsurf*' )
992	IF ( ALLOCATED( tsurf_av ) ) THEN
993	DO i = nxl, nxr
994	DO j = nys, nyn
995	match_def = surf_def_h(0)%start_index(j,i) <= &
996	surf_def_h(0)%end_index(j,i)
997	match_lsm = surf_lsm_h%start_index(j,i) <= &
998	surf_lsm_h%end_index(j,i)
999	match_usm = surf_usm_h%start_index(j,i) <= &
1000	surf_usm_h%end_index(j,i)
1001
1002	IF ( match_def ) THEN
1003	m = surf_def_h(0)%end_index(j,i)
1004	tsurf_av(j,i) = tsurf_av(j,i) + &
1005	surf_def_h(0)%pt_surface(m)
1006	ELSEIF ( match_lsm .AND. .NOT. match_usm ) THEN
1007	m = surf_lsm_h%end_index(j,i)
1008	tsurf_av(j,i) = tsurf_av(j,i) + &
1009	surf_lsm_h%pt_surface(m)
1010	ELSEIF ( match_usm ) THEN
1011	m = surf_usm_h%end_index(j,i)
1012	tsurf_av(j,i) = tsurf_av(j,i) + &
1013	surf_usm_h%pt_surface(m)
1014	ENDIF
1015	ENDDO
1016	ENDDO
1017	ENDIF
1018
1019	CASE ( 'u' )
1020	IF ( ALLOCATED( u_av ) ) THEN
1021	DO i = nxlg, nxrg
1022	DO j = nysg, nyng
1023	DO k = nzb, nzt+1
1024	u_av(k,j,i) = u_av(k,j,i) + u(k,j,i)
1025	ENDDO
1026	ENDDO
1027	ENDDO
1028	ENDIF
1029
1030	CASE ( 'us*' )
1031	IF ( ALLOCATED( us_av ) ) THEN
1032	DO i = nxl, nxr
1033	DO j = nys, nyn
1034	match_def = surf_def_h(0)%start_index(j,i) <= &
1035	surf_def_h(0)%end_index(j,i)
1036	match_lsm = surf_lsm_h%start_index(j,i) <= &
1037	surf_lsm_h%end_index(j,i)
1038	match_usm = surf_usm_h%start_index(j,i) <= &
1039	surf_usm_h%end_index(j,i)
1040
1041	IF ( match_def ) THEN
1042	m = surf_def_h(0)%end_index(j,i)
1043	us_av(j,i) = us_av(j,i) + &
1044	surf_def_h(0)%us(m)
1045	ELSEIF ( match_lsm .AND. .NOT. match_usm ) THEN
1046	m = surf_lsm_h%end_index(j,i)
1047	us_av(j,i) = us_av(j,i) + &
1048	surf_lsm_h%us(m)
1049	ELSEIF ( match_usm ) THEN
1050	m = surf_usm_h%end_index(j,i)
1051	us_av(j,i) = us_av(j,i) + &
1052	surf_usm_h%us(m)
1053	ENDIF
1054	ENDDO
1055	ENDDO
1056	ENDIF
1057
1058	CASE ( 'v' )
1059	IF ( ALLOCATED( v_av ) ) THEN
1060	DO i = nxlg, nxrg
1061	DO j = nysg, nyng
1062	DO k = nzb, nzt+1
1063	v_av(k,j,i) = v_av(k,j,i) + v(k,j,i)
1064	ENDDO
1065	ENDDO
1066	ENDDO
1067	ENDIF
1068
1069	CASE ( 'thetav' )
1070	IF ( ALLOCATED( vpt_av ) ) THEN
1071	DO i = nxlg, nxrg
1072	DO j = nysg, nyng
1073	DO k = nzb, nzt+1
1074	vpt_av(k,j,i) = vpt_av(k,j,i) + vpt(k,j,i)
1075	ENDDO
1076	ENDDO
1077	ENDDO
1078	ENDIF
1079
1080	CASE ( 'w' )
1081	IF ( ALLOCATED( w_av ) ) THEN
1082	DO i = nxlg, nxrg
1083	DO j = nysg, nyng
1084	DO k = nzb, nzt+1
1085	w_av(k,j,i) = w_av(k,j,i) + w(k,j,i)
1086	ENDDO
1087	ENDDO
1088	ENDDO
1089	ENDIF
1090
1091	CASE ( 'z0*' )
1092	IF ( ALLOCATED( z0_av ) ) THEN
1093	DO i = nxl, nxr
1094	DO j = nys, nyn
1095	match_def = surf_def_h(0)%start_index(j,i) <= &
1096	surf_def_h(0)%end_index(j,i)
1097	match_lsm = surf_lsm_h%start_index(j,i) <= &
1098	surf_lsm_h%end_index(j,i)
1099	match_usm = surf_usm_h%start_index(j,i) <= &
1100	surf_usm_h%end_index(j,i)
1101
1102	IF ( match_def ) THEN
1103	m = surf_def_h(0)%end_index(j,i)
1104	z0_av(j,i) = z0_av(j,i) + &
1105	surf_def_h(0)%z0(m)
1106	ELSEIF ( match_lsm .AND. .NOT. match_usm ) THEN
1107	m = surf_lsm_h%end_index(j,i)
1108	z0_av(j,i) = z0_av(j,i) + &
1109	surf_lsm_h%z0(m)
1110	ELSEIF ( match_usm ) THEN
1111	m = surf_usm_h%end_index(j,i)
1112	z0_av(j,i) = z0_av(j,i) + &
1113	surf_usm_h%z0(m)
1114	ENDIF
1115	ENDDO
1116	ENDDO
1117	ENDIF
1118
1119	CASE ( 'z0h*' )
1120	IF ( ALLOCATED( z0h_av ) ) THEN
1121	DO i = nxl, nxr
1122	DO j = nys, nyn
1123	match_def = surf_def_h(0)%start_index(j,i) <= &
1124	surf_def_h(0)%end_index(j,i)
1125	match_lsm = surf_lsm_h%start_index(j,i) <= &
1126	surf_lsm_h%end_index(j,i)
1127	match_usm = surf_usm_h%start_index(j,i) <= &
1128	surf_usm_h%end_index(j,i)
1129
1130	IF ( match_def ) THEN
1131	m = surf_def_h(0)%end_index(j,i)
1132	z0h_av(j,i) = z0h_av(j,i) + &
1133	surf_def_h(0)%z0h(m)
1134	ELSEIF ( match_lsm .AND. .NOT. match_usm ) THEN
1135	m = surf_lsm_h%end_index(j,i)
1136	z0h_av(j,i) = z0h_av(j,i) + &
1137	surf_lsm_h%z0h(m)
1138	ELSEIF ( match_usm ) THEN
1139	m = surf_usm_h%end_index(j,i)
1140	z0h_av(j,i) = z0h_av(j,i) + &
1141	surf_usm_h%z0h(m)
1142	ENDIF
1143	ENDDO
1144	ENDDO
1145	ENDIF
1146
1147	CASE ( 'z0q*' )
1148	IF ( ALLOCATED( z0q_av ) ) THEN
1149	DO i = nxl, nxr
1150	DO j = nys, nyn
1151	match_def = surf_def_h(0)%start_index(j,i) <= &
1152	surf_def_h(0)%end_index(j,i)
1153	match_lsm = surf_lsm_h%start_index(j,i) <= &
1154	surf_lsm_h%end_index(j,i)
1155	match_usm = surf_usm_h%start_index(j,i) <= &
1156	surf_usm_h%end_index(j,i)
1157
1158	IF ( match_def ) THEN
1159	m = surf_def_h(0)%end_index(j,i)
1160	z0q_av(j,i) = z0q_av(j,i) + &
1161	surf_def_h(0)%z0q(m)
1162	ELSEIF ( match_lsm .AND. .NOT. match_usm ) THEN
1163	m = surf_lsm_h%end_index(j,i)
1164	z0q_av(j,i) = z0q_av(j,i) + &
1165	surf_lsm_h%z0q(m)
1166	ELSEIF ( match_usm ) THEN
1167	m = surf_usm_h%end_index(j,i)
1168	z0q_av(j,i) = z0q_av(j,i) + &
1169	surf_usm_h%z0q(m)
1170	ENDIF
1171	ENDDO
1172	ENDDO
1173	ENDIF
1174
1175	CASE DEFAULT
1176	!
1177	!-- Summing up data from other modules
1178	IF ( bulk_cloud_model ) THEN
1179	CALL bcm_3d_data_averaging( 'sum', doav(ii) )
1180	ENDIF
1181
1182	IF ( air_chemistry .AND. &
1183	(trimvar(1:3) == 'kc_' .OR. trimvar(1:3) == 'em_') ) THEN
1184	CALL chem_3d_data_averaging( 'sum',doav(ii) )
1185	ENDIF
1186
1187	IF ( gust_module_enabled ) THEN
1188	CALL gust_3d_data_averaging( 'sum', doav(ii) )
1189	ENDIF
1190
1191	IF ( biometeorology ) THEN
1192	CALL bio_3d_data_averaging( 'sum', doav(ii) )
1193	ENDIF
1194
1195	IF ( land_surface ) THEN
1196	CALL lsm_3d_data_averaging( 'sum', doav(ii) )
1197	ENDIF
1198
1199	IF ( ocean_mode ) THEN
1200	CALL ocean_3d_data_averaging( 'sum', doav(ii) )
1201	ENDIF
1202
1203	IF ( radiation ) THEN
1204	CALL radiation_3d_data_averaging( 'sum', doav(ii) )
1205	ENDIF
1206
1207	IF ( salsa ) THEN
1208	CALL salsa_3d_data_averaging( 'sum', doav(ii) )
1209	ENDIF
1210
1211	CALL tcm_3d_data_averaging( 'sum', doav(ii) )
1212
1213	!-- In case of urban surface variables it should be always checked
1214	!-- if respective arrays are allocated, at least in case of a restart
1215	!-- run, as averaged usm arrays are not read from file at the moment.
1216	IF ( urban_surface .AND. trimvar(1:4) == 'usm_' ) THEN
1217	CALL usm_average_3d_data( 'allocate', doav(ii) )
1218	CALL usm_average_3d_data( 'sum', doav(ii) )
1219	ENDIF
1220
1221	!
1222	!-- User-defined quantities
1223	CALL user_3d_data_averaging( 'sum', doav(ii) )
1224
1225	END SELECT
1226
1227	ENDDO
1228
1229	CALL cpu_log( log_point(34), 'sum_up_3d_data', 'stop' )
1230
1231
1232	END SUBROUTINE sum_up_3d_data

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |