Home

Context Navigation

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

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

Rev	Line
[1682]	1	!> @file sum_up_3d_data.f90
[2000]	2	!------------------------------------------------------------------------------!
[2696]	3	! This file is part of the PALM model system.
[1036]	4	!
[2000]	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.
[1036]	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	!
[2718]	17	! Copyright 1997-2018 Leibniz Universitaet Hannover
[2000]	18	!------------------------------------------------------------------------------!
[1036]	19	!
[484]	20	! Current revisions:
[3569]	21	! ------------------
[1360]	22	!
[3589]	23	!
[1321]	24	! Former revisions:
	25	! -----------------
	26	! $Id: sum_up_3d_data.f90 3589 2018-11-30 15:09:51Z kanani $
[3589]	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
[3569]	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
[3553]	35	! variables documented
	36	!
	37	! 3552 2018-11-22 10:28:35Z suehring
[3525]	38	! Changes related to clean-up of biometeorology (dom_dwd_user)
	39	!
	40	! 3467 2018-10-30 19:05:21Z suehring
[3467]	41	! Implementation of a new aerosol module salsa.
	42	!
	43	! 3448 2018-10-29 18:14:31Z kanani
[3448]	44	! Adjustment of biometeorology calls
	45	!
	46	! 3421 2018-10-24 18:39:32Z gronemeier
[3421]	47	! Renamed output variables
	48	!
	49	! 3337 2018-10-12 15:17:09Z kanani
[3337]	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
[3294]	56	! changes concerning modularization of ocean option
	57	!
	58	! 3291 2018-09-28 11:33:03Z scharf
[3291]	59	! corrected previous commit for 3D topography
	60	!
	61	! 3285 2018-09-27 17:16:52Z scharf
[3285]	62	! bugfix for shf_av and qsws_av
	63	!
	64	! 3274 2018-09-24 15:42:55Z knoop
[3274]	65	! Modularization of all bulk cloud physics code components
	66	!
	67	! 3241 2018-09-12 15:02:00Z raasch
[3241]	68	! unused variables removed
	69	!
	70	! 3176 2018-07-26 17:12:48Z suehring
[3176]	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
[3173]	75	! Bugfix for last commit
	76	!
	77	! 3170 2018-07-25 15:19:37Z suehring
[3170]	78	! Revise output of surface quantities in case of overhanging structures
	79	!
	80	! 3151 2018-07-19 08:45:38Z raasch
[3151]	81	! Remaining preprocessor directive __chem removed
	82	!
	83	! 3004 2018-04-27 12:33:25Z Giersch
[3004]	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
[2963]	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
[2894]	92	! Changed comment
	93	!
	94	! 2817 2018-02-19 16:32:21Z suehring
[2817]	95	! Preliminary gust module interface implemented
	96	!
	97	! 2798 2018-02-09 17:16:39Z suehring
[2798]	98	! Consider also default-type surfaces for surface temperature output.
	99	!
	100	! 2797 2018-02-08 13:24:35Z suehring
[2797]	101	! Enable output of ground-heat flux also at urban surfaces.
	102	!
	103	! 2790 2018-02-06 11:57:19Z suehring
[2790]	104	! Bugfix in summation of surface sensible and latent heat flux
	105	!
	106	! 2766 2018-01-22 17:17:47Z kanani
[2766]	107	! Removed preprocessor directive __chem
	108	!
	109	! 2743 2018-01-12 16:03:39Z suehring
[2743]	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
[2742]	113	! Enable output of surface temperature
	114	!
	115	! 2735 2018-01-11 12:01:27Z suehring
[2735]	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
[2716]	119	! Corrected "Former revisions" section
	120	!
	121	! 2696 2017-12-14 17:12:51Z kanani
	122	! - Change in file header (GPL part)
[2696]	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)
[3582]	128	!
[2696]	129	! 2292 2017-06-20 09:51:42Z schwenkel
[2292]	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
[1321]	135	!
[2233]	136	! 2232 2017-05-30 17:47:52Z suehring
	137	! Adjustments to new surface concept
	138	!
[2032]	139	! 2031 2016-10-21 15:11:58Z knoop
	140	! renamed variable rho to rho_ocean and rho_av to rho_ocean_av
	141	!
[2025]	142	! 2024 2016-10-12 16:42:37Z kanani
	143	! Added missing CASE for ssws*
	144	!
[2012]	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	!
[2008]	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	!
[2001]	155	! 2000 2016-08-20 18:09:15Z knoop
	156	! Forced header and separation lines into 80 columns
	157	!
[1993]	158	! 1992 2016-08-12 15:14:59Z suehring
	159	! Bugfix in summation of passive scalar
	160	!
[1977]	161	! 1976 2016-07-27 13:28:04Z maronga
	162	! Radiation actions are now done directly in the respective module
	163	!
[1973]	164	! 1972 2016-07-26 07:52:02Z maronga
	165	! Land surface actions are now done directly in the respective module
	166	!
[1961]	167	! 1960 2016-07-12 16:34:24Z suehring
	168	! Scalar surface flux added
	169	!
[1950]	170	! 1949 2016-06-17 07:19:16Z maronga
	171	! Bugfix: calculation of lai_av, c_veg_av and c_liq_av.
	172	!
[1851]	173	! 1849 2016-04-08 11:33:18Z hoffmann
	174	! precipitation_rate moved to arrays_3d
[1852]	175	!
[1789]	176	! 1788 2016-03-10 11:01:04Z maronga
	177	! Added z0q and z0q_av
	178	!
[1694]	179	! 1693 2015-10-27 08:35:45Z maronga
	180	! Last revision text corrected
	181	!
[1692]	182	! 1691 2015-10-26 16:17:44Z maronga
	183	! Added output of Obukhov length and radiative heating rates for RRTMG.
[1693]	184	! Corrected output of liquid water path.
[1692]	185	!
[1683]	186	! 1682 2015-10-07 23:56:08Z knoop
	187	! Code annotations made doxygen readable
	188	!
[1586]	189	! 1585 2015-04-30 07:05:52Z maronga
	190	! Adapted for RRTMG
	191	!
[1556]	192	! 1555 2015-03-04 17:44:27Z maronga
	193	! Added output of r_a and r_s
	194	!
[1552]	195	! 1551 2015-03-03 14:18:16Z maronga
	196	! Added support for land surface model and radiation model data.
	197	!
[1360]	198	! 1359 2014-04-11 17:15:14Z hoffmann
	199	! New particle structure integrated.
	200	!
[1354]	201	! 1353 2014-04-08 15:21:23Z heinze
	202	! REAL constants provided with KIND-attribute
	203	!
[1321]	204	! 1320 2014-03-20 08:40:49Z raasch
[1320]	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
[1]	212	!
[1319]	213	! 1318 2014-03-17 13:35:16Z raasch
	214	! barrier argument removed from cpu_log,
	215	! module interfaces removed
	216	!
[1116]	217	! 1115 2013-03-26 18:16:16Z hoffmann
	218	! ql is calculated by calc_liquid_water_content
	219	!
[1054]	220	! 1053 2012-11-13 17:11:03Z hoffmann
	221	! +nr, prr, qr
	222	!
[1037]	223	! 1036 2012-10-22 13:43:42Z raasch
	224	! code put under GPL (PALM 3.9)
	225	!
[1008]	226	! 1007 2012-09-19 14:30:36Z franke
	227	! Bugfix in calculation of ql_vp
	228	!
[979]	229	! 978 2012-08-09 08:28:32Z fricke
	230	! +z0h*
	231	!
[1]	232	! Revision 1.1 2006/02/23 12:55:23 raasch
	233	! Initial revision
	234	!
	235	!
	236	! Description:
	237	! ------------
[1682]	238	!> Sum-up the values of 3d-arrays. The real averaging is later done in routine
	239	!> average_3d_data.
[1]	240	!------------------------------------------------------------------------------!
[1682]	241	SUBROUTINE sum_up_3d_data
	242
[1]	243
[1320]	244	USE arrays_3d, &
[3294]	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
[1]	248
[1320]	249	USE averaging, &
[3294]	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
[3241]	254
[3274]	255	USE basic_constants_and_equations_mod, &
	256	ONLY: c_p, lv_d_cp, l_v
	257
[3448]	258	USE biometeorology_mod, &
[3525]	259	ONLY: bio_3d_data_averaging
[3448]	260
[3274]	261	USE bulk_cloud_model_mod, &
	262	ONLY: bulk_cloud_model, bcm_3d_data_averaging
	263
[2696]	264	USE chemistry_model_mod, &
[3241]	265	ONLY: chem_3d_data_averaging
[1320]	266
	267	USE control_parameters, &
[3448]	268	ONLY: air_chemistry, average_count_3d, biometeorology, doav, doav_n, &
[3582]	269	land_surface, ocean_mode, rho_surface, salsa, urban_surface, &
[3569]	270	varnamelength
[1320]	271
	272	USE cpulog, &
	273	ONLY: cpu_log, log_point
	274
[2817]	275	USE gust_mod, &
	276	ONLY: gust_3d_data_averaging, gust_module_enabled
	277
[1320]	278	USE indices, &
	279	ONLY: nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, nzt
	280
	281	USE kinds
	282
[1551]	283	USE land_surface_model_mod, &
[2232]	284	ONLY: lsm_3d_data_averaging
[1551]	285
[3294]	286	USE ocean_mod, &
	287	ONLY: ocean_3d_data_averaging
	288
[1320]	289	USE particle_attributes, &
[1359]	290	ONLY: grid_particles, number_of_particles, particles, prt_count
[1320]	291
[1551]	292	USE radiation_model_mod, &
[1976]	293	ONLY: radiation, radiation_3d_data_averaging
[3467]	294
	295	USE salsa_mod, &
[3582]	296	ONLY: salsa_3d_data_averaging
[1551]	297
[2232]	298	USE surface_mod, &
[2963]	299	ONLY: ind_pav_green, ind_veg_wall, ind_wat_win, &
	300	surf_def_h, surf_lsm_h, surf_usm_h
[2232]	301
[2696]	302	USE turbulence_closure_mod, &
	303	ONLY: tcm_3d_data_averaging
	304
[2007]	305	USE urban_surface_mod, &
[2011]	306	ONLY: usm_average_3d_data
[1691]	307
[2007]	308
[1]	309	IMPLICIT NONE
	310
[3170]	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
[2232]	315	INTEGER(iwp) :: i !< grid index x direction
[2007]	316	INTEGER(iwp) :: ii !< running index
[2232]	317	INTEGER(iwp) :: j !< grid index y direction
	318	INTEGER(iwp) :: k !< grid index x direction
[3552]	319	INTEGER(iwp) :: m !< running index over surfacle elements
	320	INTEGER(iwp) :: n !< running index over number of particles per grid box
[1]	321
[3552]	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
[1]	325
[2011]	326	CHARACTER (LEN=varnamelength) :: trimvar !< TRIM of output-variable string
[2007]	327
	328
[1]	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
[2894]	335	!-- in rrd_local)
[1]	336	IF ( average_count_3d == 0 ) THEN
	337
	338	DO ii = 1, doav_n
[3337]	339
[2007]	340	trimvar = TRIM( doav(ii) )
[3337]	341
[2007]	342	SELECT CASE ( trimvar )
[1]	343
[2797]	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
[1]	350	CASE ( 'e' )
	351	IF ( .NOT. ALLOCATED( e_av ) ) THEN
[667]	352	ALLOCATE( e_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	353	ENDIF
[1353]	354	e_av = 0.0_wp
[1]	355
[3421]	356	CASE ( 'thetal' )
[771]	357	IF ( .NOT. ALLOCATED( lpt_av ) ) THEN
	358	ALLOCATE( lpt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
	359	ENDIF
[1353]	360	lpt_av = 0.0_wp
[771]	361
[1]	362	CASE ( 'lwp*' )
	363	IF ( .NOT. ALLOCATED( lwp_av ) ) THEN
[667]	364	ALLOCATE( lwp_av(nysg:nyng,nxlg:nxrg) )
[1]	365	ENDIF
[1353]	366	lwp_av = 0.0_wp
[1]	367
[1691]	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
[1]	374	CASE ( 'p' )
	375	IF ( .NOT. ALLOCATED( p_av ) ) THEN
[667]	376	ALLOCATE( p_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	377	ENDIF
[1353]	378	p_av = 0.0_wp
[1]	379
	380	CASE ( 'pc' )
	381	IF ( .NOT. ALLOCATED( pc_av ) ) THEN
[667]	382	ALLOCATE( pc_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	383	ENDIF
[1353]	384	pc_av = 0.0_wp
[1]	385
	386	CASE ( 'pr' )
	387	IF ( .NOT. ALLOCATED( pr_av ) ) THEN
[667]	388	ALLOCATE( pr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	389	ENDIF
[1353]	390	pr_av = 0.0_wp
[1]	391
[3421]	392	CASE ( 'theta' )
[1]	393	IF ( .NOT. ALLOCATED( pt_av ) ) THEN
[667]	394	ALLOCATE( pt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	395	ENDIF
[1353]	396	pt_av = 0.0_wp
[1]	397
	398	CASE ( 'q' )
	399	IF ( .NOT. ALLOCATED( q_av ) ) THEN
[667]	400	ALLOCATE( q_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	401	ENDIF
[1353]	402	q_av = 0.0_wp
[1]	403
	404	CASE ( 'ql' )
	405	IF ( .NOT. ALLOCATED( ql_av ) ) THEN
[667]	406	ALLOCATE( ql_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	407	ENDIF
[1353]	408	ql_av = 0.0_wp
[1]	409
	410	CASE ( 'ql_c' )
	411	IF ( .NOT. ALLOCATED( ql_c_av ) ) THEN
[667]	412	ALLOCATE( ql_c_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	413	ENDIF
[1353]	414	ql_c_av = 0.0_wp
[1]	415
	416	CASE ( 'ql_v' )
	417	IF ( .NOT. ALLOCATED( ql_v_av ) ) THEN
[667]	418	ALLOCATE( ql_v_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	419	ENDIF
[1353]	420	ql_v_av = 0.0_wp
[1]	421
	422	CASE ( 'ql_vp' )
	423	IF ( .NOT. ALLOCATED( ql_vp_av ) ) THEN
[667]	424	ALLOCATE( ql_vp_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	425	ENDIF
[1353]	426	ql_vp_av = 0.0_wp
[1]	427
[354]	428	CASE ( 'qsws*' )
	429	IF ( .NOT. ALLOCATED( qsws_av ) ) THEN
[667]	430	ALLOCATE( qsws_av(nysg:nyng,nxlg:nxrg) )
[354]	431	ENDIF
[1353]	432	qsws_av = 0.0_wp
[354]	433
[1]	434	CASE ( 'qv' )
	435	IF ( .NOT. ALLOCATED( qv_av ) ) THEN
[667]	436	ALLOCATE( qv_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	437	ENDIF
[1353]	438	qv_av = 0.0_wp
[1]	439
[2735]	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
[1]	446	CASE ( 's' )
	447	IF ( .NOT. ALLOCATED( s_av ) ) THEN
[667]	448	ALLOCATE( s_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	449	ENDIF
[1353]	450	s_av = 0.0_wp
[1]	451
[354]	452	CASE ( 'shf*' )
	453	IF ( .NOT. ALLOCATED( shf_av ) ) THEN
[667]	454	ALLOCATE( shf_av(nysg:nyng,nxlg:nxrg) )
[354]	455	ENDIF
[1353]	456	shf_av = 0.0_wp
[2024]	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
[354]	463
[1]	464	CASE ( 't*' )
	465	IF ( .NOT. ALLOCATED( ts_av ) ) THEN
[667]	466	ALLOCATE( ts_av(nysg:nyng,nxlg:nxrg) )
[1]	467	ENDIF
[1353]	468	ts_av = 0.0_wp
[1]	469
[2742]	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
[1]	476	CASE ( 'u' )
	477	IF ( .NOT. ALLOCATED( u_av ) ) THEN
[667]	478	ALLOCATE( u_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	479	ENDIF
[1353]	480	u_av = 0.0_wp
[1]	481
[3421]	482	CASE ( 'us*' )
[1]	483	IF ( .NOT. ALLOCATED( us_av ) ) THEN
[667]	484	ALLOCATE( us_av(nysg:nyng,nxlg:nxrg) )
[1]	485	ENDIF
[1353]	486	us_av = 0.0_wp
[1]	487
	488	CASE ( 'v' )
	489	IF ( .NOT. ALLOCATED( v_av ) ) THEN
[667]	490	ALLOCATE( v_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	491	ENDIF
[1353]	492	v_av = 0.0_wp
[1]	493
[3421]	494	CASE ( 'thetav' )
[1]	495	IF ( .NOT. ALLOCATED( vpt_av ) ) THEN
[667]	496	ALLOCATE( vpt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	497	ENDIF
[1353]	498	vpt_av = 0.0_wp
[1]	499
	500	CASE ( 'w' )
	501	IF ( .NOT. ALLOCATED( w_av ) ) THEN
[667]	502	ALLOCATE( w_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	503	ENDIF
[1353]	504	w_av = 0.0_wp
[1]	505
[72]	506	CASE ( 'z0*' )
	507	IF ( .NOT. ALLOCATED( z0_av ) ) THEN
[667]	508	ALLOCATE( z0_av(nysg:nyng,nxlg:nxrg) )
[72]	509	ENDIF
[1353]	510	z0_av = 0.0_wp
[72]	511
[978]	512	CASE ( 'z0h*' )
	513	IF ( .NOT. ALLOCATED( z0h_av ) ) THEN
	514	ALLOCATE( z0h_av(nysg:nyng,nxlg:nxrg) )
	515	ENDIF
[1353]	516	z0h_av = 0.0_wp
[978]	517
[1788]	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
[3294]	523
[2007]	524
[1]	525	CASE DEFAULT
[1972]	526
[1]	527	!
[3294]	528	!-- Allocating and initializing data arrays for other modules
[3337]	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
[3274]	535	IF ( bulk_cloud_model ) THEN
	536	CALL bcm_3d_data_averaging( 'allocate', doav(ii) )
	537	ENDIF
	538
[3294]	539	IF ( gust_module_enabled ) THEN
	540	CALL gust_3d_data_averaging( 'allocate', doav(ii) )
	541	ENDIF
	542
[3569]	543	IF ( biometeorology ) THEN
[3525]	544	CALL bio_3d_data_averaging( 'allocate', doav(ii) )
[3448]	545	ENDIF
	546
[1972]	547	IF ( land_surface ) THEN
	548	CALL lsm_3d_data_averaging( 'allocate', doav(ii) )
	549	ENDIF
	550
[3294]	551	IF ( ocean_mode ) THEN
	552	CALL ocean_3d_data_averaging( 'allocate', doav(ii) )
	553	ENDIF
	554
[1976]	555	IF ( radiation ) THEN
	556	CALL radiation_3d_data_averaging( 'allocate', doav(ii) )
	557	ENDIF
[3582]	558
[3467]	559	IF ( salsa ) THEN
	560	CALL salsa_3d_data_averaging( 'allocate', doav(ii) )
	561	ENDIF
[1976]	562
[3337]	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
[2696]	569	!
[3294]	570	!-- User-defined quantities
[1]	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
[3337]	582
	583	trimvar = TRIM( doav(ii) )
[1]	584	!
	585	!-- Store the array chosen on the temporary array.
[2007]	586	SELECT CASE ( trimvar )
[1]	587
[2797]	588	CASE ( 'ghf*' )
[3004]	589	IF ( ALLOCATED( ghf_av ) ) THEN
[3170]	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)
[3173]	607	ghf_av(j,i) = ghf_av(j,i) + &
[3170]	608	surf_lsm_h%ghf(m)
	609	ELSEIF ( match_usm ) THEN
	610	m = surf_usm_h%end_index(j,i)
[3173]	611	ghf_av(j,i) = ghf_av(j,i) + &
[3170]	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
[3004]	620	ENDDO
	621	ENDIF
[2797]	622
[1]	623	CASE ( 'e' )
[3004]	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
[1]	630	ENDDO
	631	ENDDO
[3004]	632	ENDIF
[1]	633
[3421]	634	CASE ( 'thetal' )
[3004]	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
[771]	641	ENDDO
	642	ENDDO
[3004]	643	ENDIF
[771]	644
[1]	645	CASE ( 'lwp*' )
[3004]	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
[1]	652	ENDDO
[3004]	653	ENDIF
[1]	654
[1691]	655	CASE ( 'ol*' )
[3004]	656	IF ( ALLOCATED( ol_av ) ) THEN
[3170]	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)
[3173]	668	ol_av(j,i) = ol_av(j,i) + &
[3170]	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)
[3173]	672	ol_av(j,i) = ol_av(j,i) + &
[3170]	673	surf_lsm_h%ol(m)
	674	ELSEIF ( match_usm ) THEN
	675	m = surf_usm_h%end_index(j,i)
[3173]	676	ol_av(j,i) = ol_av(j,i) + &
[3170]	677	surf_usm_h%ol(m)
	678	ENDIF
	679	ENDDO
[3004]	680	ENDDO
	681	ENDIF
[1691]	682
[1]	683	CASE ( 'p' )
[3004]	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
[1]	690	ENDDO
	691	ENDDO
[3004]	692	ENDIF
[1]	693
	694	CASE ( 'pc' )
[3004]	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
[1]	701	ENDDO
	702	ENDDO
[3004]	703	ENDIF
[1]	704
	705	CASE ( 'pr' )
[3004]	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
[1359]	716
[3004]	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
[1359]	730	ENDIF
[3004]	731	pr_av(k,j,i) = pr_av(k,j,i) + mean_r
[1]	732	ENDDO
	733	ENDDO
	734	ENDDO
[3004]	735	ENDIF
[1]	736
[3421]	737	CASE ( 'theta' )
[3004]	738	IF ( ALLOCATED( pt_av ) ) THEN
[3274]	739	IF ( .NOT. bulk_cloud_model ) THEN
[3004]	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
[1]	745	ENDDO
	746	ENDDO
[3004]	747	ELSE
	748	DO i = nxlg, nxrg
	749	DO j = nysg, nyng
	750	DO k = nzb, nzt+1
[3274]	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)
[3004]	753	ENDDO
[1]	754	ENDDO
	755	ENDDO
[3004]	756	ENDIF
[1]	757	ENDIF
	758
	759	CASE ( 'q' )
[3004]	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
[1]	766	ENDDO
	767	ENDDO
[3004]	768	ENDIF
[402]	769
[1]	770	CASE ( 'ql' )
[3004]	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
[1]	777	ENDDO
	778	ENDDO
[3004]	779	ENDIF
[1]	780
	781	CASE ( 'ql_c' )
[3004]	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
[1]	788	ENDDO
	789	ENDDO
[3004]	790	ENDIF
[1]	791
	792	CASE ( 'ql_v' )
[3004]	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
[1]	799	ENDDO
	800	ENDDO
[3004]	801	ENDIF
[1]	802
	803	CASE ( 'ql_vp' )
[3004]	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
[1007]	819	ENDDO
[1]	820	ENDDO
	821	ENDDO
[3004]	822	ENDIF
[1]	823
[402]	824	CASE ( 'qsws*' )
[2743]	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.
[3004]	829	IF ( ALLOCATED( qsws_av ) ) THEN
[3170]	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)
[3173]	841	qsws_av(j,i) = qsws_av(j,i) + &
[3170]	842	surf_def_h(0)%qsws(m) * &
[3285]	843	waterflux_output_conversion(nzb)
[3170]	844	ELSEIF ( match_lsm .AND. .NOT. match_usm ) THEN
	845	m = surf_lsm_h%end_index(j,i)
[3173]	846	qsws_av(j,i) = qsws_av(j,i) + &
[3170]	847	surf_lsm_h%qsws(m) * l_v
	848	ENDIF
	849	ENDDO
[3004]	850	ENDDO
	851	ENDIF
[402]	852
[1]	853	CASE ( 'qv' )
[3004]	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
[1]	860	ENDDO
	861	ENDDO
[3004]	862	ENDIF
[1]	863
[2735]	864	CASE ( 'r_a*' )
[3004]	865	IF ( ALLOCATED( r_a_av ) ) THEN
[3170]	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)
[3173]	875	r_a_av(j,i) = r_a_av(j,i) + &
[3170]	876	surf_lsm_h%r_a(m)
	877	ELSEIF ( match_usm ) THEN
	878	m = surf_usm_h%end_index(j,i)
[3173]	879	r_a_av(j,i) = r_a_av(j,i) + &
[3170]	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
[3004]	888	ENDDO
	889	ENDIF
[2735]	890
[1]	891	CASE ( 's' )
[3004]	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
[1]	898	ENDDO
	899	ENDDO
[3004]	900	ENDIF
[402]	901
	902	CASE ( 'shf*' )
[2743]	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.
[3004]	907	IF ( ALLOCATED( shf_av ) ) THEN
[3170]	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)
[3173]	919	shf_av(j,i) = shf_av(j,i) + &
[3170]	920	surf_def_h(0)%shf(m) * &
[3285]	921	heatflux_output_conversion(nzb)
[3170]	922	ELSEIF ( match_lsm .AND. .NOT. match_usm ) THEN
	923	m = surf_lsm_h%end_index(j,i)
[3173]	924	shf_av(j,i) = shf_av(j,i) + &
[3274]	925	surf_lsm_h%shf(m) * c_p
[3170]	926	ELSEIF ( match_usm ) THEN
	927	m = surf_usm_h%end_index(j,i)
[3173]	928	shf_av(j,i) = shf_av(j,i) + &
[3274]	929	surf_usm_h%shf(m) * c_p
[3170]	930	ENDIF
	931	ENDDO
[3004]	932	ENDDO
	933	ENDIF
[402]	934
[1960]	935	CASE ( 'ssws*' )
[3004]	936	IF ( ALLOCATED( ssws_av ) ) THEN
[3170]	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)
[3173]	948	ssws_av(j,i) = ssws_av(j,i) + &
[3170]	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)
[3173]	952	ssws_av(j,i) = ssws_av(j,i) + &
[3170]	953	surf_lsm_h%ssws(m)
	954	ELSEIF ( match_usm ) THEN
	955	m = surf_usm_h%end_index(j,i)
[3173]	956	ssws_av(j,i) = ssws_av(j,i) + &
[3170]	957	surf_usm_h%ssws(m)
	958	ENDIF
	959	ENDDO
[3004]	960	ENDDO
	961	ENDIF
[1960]	962
[1]	963	CASE ( 't*' )
[3004]	964	IF ( ALLOCATED( ts_av ) ) THEN
[3170]	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)
[3173]	976	ts_av(j,i) = ts_av(j,i) + &
[3170]	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)
[3173]	980	ts_av(j,i) = ts_av(j,i) + &
[3170]	981	surf_lsm_h%ts(m)
	982	ELSEIF ( match_usm ) THEN
	983	m = surf_usm_h%end_index(j,i)
[3173]	984	ts_av(j,i) = ts_av(j,i) + &
[3170]	985	surf_usm_h%ts(m)
	986	ENDIF
	987	ENDDO
[3004]	988	ENDDO
	989	ENDIF
[1]	990
[2742]	991	CASE ( 'tsurf*' )
[3170]	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)
[2798]	1001
[3170]	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
[3004]	1016	ENDDO
	1017	ENDIF
[2742]	1018
[1]	1019	CASE ( 'u' )
[3004]	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
[1]	1026	ENDDO
	1027	ENDDO
[3004]	1028	ENDIF
[1]	1029
[3421]	1030	CASE ( 'us*' )
[3004]	1031	IF ( ALLOCATED( us_av ) ) THEN
[3170]	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)
[3173]	1043	us_av(j,i) = us_av(j,i) + &
[3170]	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)
[3173]	1047	us_av(j,i) = us_av(j,i) + &
[3170]	1048	surf_lsm_h%us(m)
	1049	ELSEIF ( match_usm ) THEN
	1050	m = surf_usm_h%end_index(j,i)
[3173]	1051	us_av(j,i) = us_av(j,i) + &
[3170]	1052	surf_usm_h%us(m)
	1053	ENDIF
	1054	ENDDO
[3004]	1055	ENDDO
	1056	ENDIF
[1]	1057
	1058	CASE ( 'v' )
[3004]	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
[1]	1065	ENDDO
	1066	ENDDO
[3004]	1067	ENDIF
[1]	1068
[3421]	1069	CASE ( 'thetav' )
[3004]	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
[1]	1076	ENDDO
	1077	ENDDO
[3004]	1078	ENDIF
[1]	1079
	1080	CASE ( 'w' )
[3004]	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
[1]	1087	ENDDO
	1088	ENDDO
[3004]	1089	ENDIF
[1]	1090
[72]	1091	CASE ( 'z0*' )
[3004]	1092	IF ( ALLOCATED( z0_av ) ) THEN
[3170]	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)
[3173]	1104	z0_av(j,i) = z0_av(j,i) + &
[3170]	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)
[3173]	1108	z0_av(j,i) = z0_av(j,i) + &
[3170]	1109	surf_lsm_h%z0(m)
	1110	ELSEIF ( match_usm ) THEN
	1111	m = surf_usm_h%end_index(j,i)
[3173]	1112	z0_av(j,i) = z0_av(j,i) + &
[3170]	1113	surf_usm_h%z0(m)
	1114	ENDIF
	1115	ENDDO
	1116	ENDDO
[3004]	1117	ENDIF
[72]	1118
[978]	1119	CASE ( 'z0h*' )
[3004]	1120	IF ( ALLOCATED( z0h_av ) ) THEN
[3170]	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)
[3173]	1132	z0h_av(j,i) = z0h_av(j,i) + &
[3170]	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)
[3173]	1136	z0h_av(j,i) = z0h_av(j,i) + &
[3170]	1137	surf_lsm_h%z0h(m)
	1138	ELSEIF ( match_usm ) THEN
	1139	m = surf_usm_h%end_index(j,i)
[3173]	1140	z0h_av(j,i) = z0h_av(j,i) + &
[3170]	1141	surf_usm_h%z0h(m)
	1142	ENDIF
	1143	ENDDO
[3004]	1144	ENDDO
	1145	ENDIF
	1146
[1788]	1147	CASE ( 'z0q*' )
[3004]	1148	IF ( ALLOCATED( z0q_av ) ) THEN
[3170]	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)
[3173]	1160	z0q_av(j,i) = z0q_av(j,i) + &
[3170]	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)
[3173]	1164	z0q_av(j,i) = z0q_av(j,i) + &
[3170]	1165	surf_lsm_h%z0q(m)
	1166	ELSEIF ( match_usm ) THEN
	1167	m = surf_usm_h%end_index(j,i)
[3173]	1168	z0q_av(j,i) = z0q_av(j,i) + &
[3170]	1169	surf_usm_h%z0q(m)
	1170	ENDIF
	1171	ENDDO
[3004]	1172	ENDDO
	1173	ENDIF
[3294]	1174
[1]	1175	CASE DEFAULT
	1176	!
[3294]	1177	!-- Summing up data from other modules
[3274]	1178	IF ( bulk_cloud_model ) THEN
	1179	CALL bcm_3d_data_averaging( 'sum', doav(ii) )
	1180	ENDIF
	1181
[3337]	1182	IF ( air_chemistry .AND. &
	1183	(trimvar(1:3) == 'kc_' .OR. trimvar(1:3) == 'em_') ) THEN
[3294]	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
[3569]	1191	IF ( biometeorology ) THEN
[3525]	1192	CALL bio_3d_data_averaging( 'sum', doav(ii) )
[3448]	1193	ENDIF
	1194
[1972]	1195	IF ( land_surface ) THEN
	1196	CALL lsm_3d_data_averaging( 'sum', doav(ii) )
	1197	ENDIF
	1198
[3294]	1199	IF ( ocean_mode ) THEN
	1200	CALL ocean_3d_data_averaging( 'sum', doav(ii) )
	1201	ENDIF
	1202
[1976]	1203	IF ( radiation ) THEN
	1204	CALL radiation_3d_data_averaging( 'sum', doav(ii) )
	1205	ENDIF
[3582]	1206
[3467]	1207	IF ( salsa ) THEN
	1208	CALL salsa_3d_data_averaging( 'sum', doav(ii) )
	1209	ENDIF
[1976]	1210
[3294]	1211	CALL tcm_3d_data_averaging( 'sum', doav(ii) )
[2817]	1212
[3337]	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
[2696]	1221	!
[3294]	1222	!-- User-defined quantities
[1]	1223	CALL user_3d_data_averaging( 'sum', doav(ii) )
	1224
	1225	END SELECT
	1226
	1227	ENDDO
	1228
[1318]	1229	CALL cpu_log( log_point(34), 'sum_up_3d_data', 'stop' )
[1]	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 |