Home

Context Navigation

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

Last change on this file since 1359 was 1359, checked in by hoffmann, 10 years ago
new Lagrangian particle structure integrated
Property svn:keywords set to `Id`
File size: 19.9 KB

Rev	Line
[1]	1	SUBROUTINE sum_up_3d_data
	2
[1036]	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	!
[1310]	17	! Copyright 1997-2014 Leibniz Universitaet Hannover
[1036]	18	!--------------------------------------------------------------------------------!
	19	!
[484]	20	! Current revisions:
[1]	21	! -----------------
[1359]	22	! New particle structure integrated.
	23	!
[1321]	24	! Former revisions:
	25	! -----------------
	26	! $Id: sum_up_3d_data.f90 1359 2014-04-11 17:15:14Z hoffmann $
	27	!
[1354]	28	! 1353 2014-04-08 15:21:23Z heinze
	29	! REAL constants provided with KIND-attribute
	30	!
[1321]	31	! 1320 2014-03-20 08:40:49Z raasch
[1320]	32	! ONLY-attribute added to USE-statements,
	33	! kind-parameters added to all INTEGER and REAL declaration statements,
	34	! kinds are defined in new module kinds,
	35	! old module precision_kind is removed,
	36	! revision history before 2012 removed,
	37	! comment fields (!:) to be used for variable explanations added to
	38	! all variable declaration statements
[1]	39	!
[1319]	40	! 1318 2014-03-17 13:35:16Z raasch
	41	! barrier argument removed from cpu_log,
	42	! module interfaces removed
	43	!
[1116]	44	! 1115 2013-03-26 18:16:16Z hoffmann
	45	! ql is calculated by calc_liquid_water_content
	46	!
[1054]	47	! 1053 2012-11-13 17:11:03Z hoffmann
	48	! +nr, prr, qr
	49	!
[1037]	50	! 1036 2012-10-22 13:43:42Z raasch
	51	! code put under GPL (PALM 3.9)
	52	!
[1008]	53	! 1007 2012-09-19 14:30:36Z franke
	54	! Bugfix in calculation of ql_vp
	55	!
[979]	56	! 978 2012-08-09 08:28:32Z fricke
	57	! +z0h*
	58	!
[1]	59	! Revision 1.1 2006/02/23 12:55:23 raasch
	60	! Initial revision
	61	!
	62	!
	63	! Description:
	64	! ------------
	65	! Sum-up the values of 3d-arrays. The real averaging is later done in routine
	66	! average_3d_data.
	67	!------------------------------------------------------------------------------!
	68
[1320]	69	USE arrays_3d, &
	70	ONLY: dzw, e, nr, p, pt, q, qc, ql, ql_c, ql_v, qr, qsws, rho, sa, &
	71	shf, ts, u, us, v, vpt, w, z0, z0h
[1]	72
[1320]	73	USE averaging, &
	74	ONLY: e_av, lpt_av, lwp_av, nr_av, p_av, pc_av, pr_av, prr_av, &
	75	precipitation_rate_av, pt_av, q_av, qc_av, ql_av, ql_c_av, &
	76	ql_v_av, ql_vp_av, qr_av, qsws_av, qv_av, rho_av, s_av, sa_av, &
	77	shf_av, ts_av, u_av, us_av, v_av, vpt_av, w_av, z0_av, z0h_av
	78
	79	USE cloud_parameters, &
	80	ONLY: l_d_cp, precipitation_rate, pt_d_t
	81
	82	USE control_parameters, &
	83	ONLY: average_count_3d, cloud_physics, doav, doav_n
	84
	85	USE cpulog, &
	86	ONLY: cpu_log, log_point
	87
	88	USE indices, &
	89	ONLY: nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, nzt
	90
	91	USE kinds
	92
	93	USE particle_attributes, &
[1359]	94	ONLY: grid_particles, number_of_particles, particles, prt_count
[1320]	95
[1]	96	IMPLICIT NONE
	97
[1320]	98	INTEGER(iwp) :: i !:
	99	INTEGER(iwp) :: ii !:
	100	INTEGER(iwp) :: j !:
	101	INTEGER(iwp) :: k !:
	102	INTEGER(iwp) :: n !:
	103	INTEGER(iwp) :: psi !:
[1]	104
[1320]	105	REAL(wp) :: mean_r !:
[1359]	106	REAL(wp) :: s_r2 !:
[1320]	107	REAL(wp) :: s_r3 !:
[1]	108
	109	CALL cpu_log (log_point(34),'sum_up_3d_data','start')
	110
	111	!
	112	!-- Allocate and initialize the summation arrays if called for the very first
	113	!-- time or the first time after average_3d_data has been called
	114	!-- (some or all of the arrays may have been already allocated
	115	!-- in read_3d_binary)
	116	IF ( average_count_3d == 0 ) THEN
	117
	118	DO ii = 1, doav_n
	119
	120	SELECT CASE ( TRIM( doav(ii) ) )
	121
	122	CASE ( 'e' )
	123	IF ( .NOT. ALLOCATED( e_av ) ) THEN
[667]	124	ALLOCATE( e_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	125	ENDIF
[1353]	126	e_av = 0.0_wp
[1]	127
[771]	128	CASE ( 'lpt' )
	129	IF ( .NOT. ALLOCATED( lpt_av ) ) THEN
	130	ALLOCATE( lpt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
	131	ENDIF
[1353]	132	lpt_av = 0.0_wp
[771]	133
[1]	134	CASE ( 'lwp*' )
	135	IF ( .NOT. ALLOCATED( lwp_av ) ) THEN
[667]	136	ALLOCATE( lwp_av(nysg:nyng,nxlg:nxrg) )
[1]	137	ENDIF
[1353]	138	lwp_av = 0.0_wp
[1]	139
[1053]	140	CASE ( 'nr' )
	141	IF ( .NOT. ALLOCATED( nr_av ) ) THEN
	142	ALLOCATE( nr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
	143	ENDIF
[1353]	144	nr_av = 0.0_wp
[1053]	145
[1]	146	CASE ( 'p' )
	147	IF ( .NOT. ALLOCATED( p_av ) ) THEN
[667]	148	ALLOCATE( p_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	149	ENDIF
[1353]	150	p_av = 0.0_wp
[1]	151
	152	CASE ( 'pc' )
	153	IF ( .NOT. ALLOCATED( pc_av ) ) THEN
[667]	154	ALLOCATE( pc_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	155	ENDIF
[1353]	156	pc_av = 0.0_wp
[1]	157
	158	CASE ( 'pr' )
	159	IF ( .NOT. ALLOCATED( pr_av ) ) THEN
[667]	160	ALLOCATE( pr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	161	ENDIF
[1353]	162	pr_av = 0.0_wp
[1]	163
[1053]	164	CASE ( 'prr' )
	165	IF ( .NOT. ALLOCATED( prr_av ) ) THEN
	166	ALLOCATE( prr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
	167	ENDIF
[1353]	168	prr_av = 0.0_wp
[1053]	169
[72]	170	CASE ( 'prr*' )
	171	IF ( .NOT. ALLOCATED( precipitation_rate_av ) ) THEN
[667]	172	ALLOCATE( precipitation_rate_av(nysg:nyng,nxlg:nxrg) )
[72]	173	ENDIF
[1353]	174	precipitation_rate_av = 0.0_wp
[72]	175
[1]	176	CASE ( 'pt' )
	177	IF ( .NOT. ALLOCATED( pt_av ) ) THEN
[667]	178	ALLOCATE( pt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	179	ENDIF
[1353]	180	pt_av = 0.0_wp
[1]	181
	182	CASE ( 'q' )
	183	IF ( .NOT. ALLOCATED( q_av ) ) THEN
[667]	184	ALLOCATE( q_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	185	ENDIF
[1353]	186	q_av = 0.0_wp
[1]	187
[1115]	188	CASE ( 'qc' )
	189	IF ( .NOT. ALLOCATED( qc_av ) ) THEN
	190	ALLOCATE( qc_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
	191	ENDIF
[1353]	192	qc_av = 0.0_wp
[1115]	193
[1]	194	CASE ( 'ql' )
	195	IF ( .NOT. ALLOCATED( ql_av ) ) THEN
[667]	196	ALLOCATE( ql_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	197	ENDIF
[1353]	198	ql_av = 0.0_wp
[1]	199
	200	CASE ( 'ql_c' )
	201	IF ( .NOT. ALLOCATED( ql_c_av ) ) THEN
[667]	202	ALLOCATE( ql_c_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	203	ENDIF
[1353]	204	ql_c_av = 0.0_wp
[1]	205
	206	CASE ( 'ql_v' )
	207	IF ( .NOT. ALLOCATED( ql_v_av ) ) THEN
[667]	208	ALLOCATE( ql_v_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	209	ENDIF
[1353]	210	ql_v_av = 0.0_wp
[1]	211
	212	CASE ( 'ql_vp' )
	213	IF ( .NOT. ALLOCATED( ql_vp_av ) ) THEN
[667]	214	ALLOCATE( ql_vp_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	215	ENDIF
[1353]	216	ql_vp_av = 0.0_wp
[1]	217
[1053]	218	CASE ( 'qr' )
	219	IF ( .NOT. ALLOCATED( qr_av ) ) THEN
	220	ALLOCATE( qr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
	221	ENDIF
[1353]	222	qr_av = 0.0_wp
[1053]	223
[354]	224	CASE ( 'qsws*' )
	225	IF ( .NOT. ALLOCATED( qsws_av ) ) THEN
[667]	226	ALLOCATE( qsws_av(nysg:nyng,nxlg:nxrg) )
[354]	227	ENDIF
[1353]	228	qsws_av = 0.0_wp
[354]	229
[1]	230	CASE ( 'qv' )
	231	IF ( .NOT. ALLOCATED( qv_av ) ) THEN
[667]	232	ALLOCATE( qv_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	233	ENDIF
[1353]	234	qv_av = 0.0_wp
[1]	235
[96]	236	CASE ( 'rho' )
	237	IF ( .NOT. ALLOCATED( rho_av ) ) THEN
[667]	238	ALLOCATE( rho_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[96]	239	ENDIF
[1353]	240	rho_av = 0.0_wp
[96]	241
[1]	242	CASE ( 's' )
	243	IF ( .NOT. ALLOCATED( s_av ) ) THEN
[667]	244	ALLOCATE( s_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	245	ENDIF
[1353]	246	s_av = 0.0_wp
[1]	247
[96]	248	CASE ( 'sa' )
	249	IF ( .NOT. ALLOCATED( sa_av ) ) THEN
[667]	250	ALLOCATE( sa_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[96]	251	ENDIF
[1353]	252	sa_av = 0.0_wp
[96]	253
[354]	254	CASE ( 'shf*' )
	255	IF ( .NOT. ALLOCATED( shf_av ) ) THEN
[667]	256	ALLOCATE( shf_av(nysg:nyng,nxlg:nxrg) )
[354]	257	ENDIF
[1353]	258	shf_av = 0.0_wp
[354]	259
[1]	260	CASE ( 't*' )
	261	IF ( .NOT. ALLOCATED( ts_av ) ) THEN
[667]	262	ALLOCATE( ts_av(nysg:nyng,nxlg:nxrg) )
[1]	263	ENDIF
[1353]	264	ts_av = 0.0_wp
[1]	265
	266	CASE ( 'u' )
	267	IF ( .NOT. ALLOCATED( u_av ) ) THEN
[667]	268	ALLOCATE( u_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	269	ENDIF
[1353]	270	u_av = 0.0_wp
[1]	271
	272	CASE ( 'u*' )
	273	IF ( .NOT. ALLOCATED( us_av ) ) THEN
[667]	274	ALLOCATE( us_av(nysg:nyng,nxlg:nxrg) )
[1]	275	ENDIF
[1353]	276	us_av = 0.0_wp
[1]	277
	278	CASE ( 'v' )
	279	IF ( .NOT. ALLOCATED( v_av ) ) THEN
[667]	280	ALLOCATE( v_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	281	ENDIF
[1353]	282	v_av = 0.0_wp
[1]	283
	284	CASE ( 'vpt' )
	285	IF ( .NOT. ALLOCATED( vpt_av ) ) THEN
[667]	286	ALLOCATE( vpt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	287	ENDIF
[1353]	288	vpt_av = 0.0_wp
[1]	289
	290	CASE ( 'w' )
	291	IF ( .NOT. ALLOCATED( w_av ) ) THEN
[667]	292	ALLOCATE( w_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
[1]	293	ENDIF
[1353]	294	w_av = 0.0_wp
[1]	295
[72]	296	CASE ( 'z0*' )
	297	IF ( .NOT. ALLOCATED( z0_av ) ) THEN
[667]	298	ALLOCATE( z0_av(nysg:nyng,nxlg:nxrg) )
[72]	299	ENDIF
[1353]	300	z0_av = 0.0_wp
[72]	301
[978]	302	CASE ( 'z0h*' )
	303	IF ( .NOT. ALLOCATED( z0h_av ) ) THEN
	304	ALLOCATE( z0h_av(nysg:nyng,nxlg:nxrg) )
	305	ENDIF
[1353]	306	z0h_av = 0.0_wp
[978]	307
[1]	308	CASE DEFAULT
	309	!
	310	!-- User-defined quantity
	311	CALL user_3d_data_averaging( 'allocate', doav(ii) )
	312
	313	END SELECT
	314
	315	ENDDO
	316
	317	ENDIF
	318
	319	!
	320	!-- Loop of all variables to be averaged.
	321	DO ii = 1, doav_n
	322
	323	!
	324	!-- Store the array chosen on the temporary array.
	325	SELECT CASE ( TRIM( doav(ii) ) )
	326
	327	CASE ( 'e' )
[667]	328	DO i = nxlg, nxrg
	329	DO j = nysg, nyng
[1]	330	DO k = nzb, nzt+1
	331	e_av(k,j,i) = e_av(k,j,i) + e(k,j,i)
	332	ENDDO
	333	ENDDO
	334	ENDDO
	335
[771]	336	CASE ( 'lpt' )
	337	DO i = nxlg, nxrg
	338	DO j = nysg, nyng
	339	DO k = nzb, nzt+1
	340	lpt_av(k,j,i) = lpt_av(k,j,i) + pt(k,j,i)
	341	ENDDO
	342	ENDDO
	343	ENDDO
	344
[1]	345	CASE ( 'lwp*' )
[667]	346	DO i = nxlg, nxrg
	347	DO j = nysg, nyng
[1]	348	lwp_av(j,i) = lwp_av(j,i) + SUM( ql(nzb:nzt,j,i) * &
	349	dzw(1:nzt+1) )
	350	ENDDO
	351	ENDDO
	352
[1053]	353	CASE ( 'nr' )
	354	DO i = nxlg, nxrg
	355	DO j = nysg, nyng
	356	DO k = nzb, nzt+1
	357	nr_av(k,j,i) = nr_av(k,j,i) + nr(k,j,i)
	358	ENDDO
	359	ENDDO
	360	ENDDO
	361
[1]	362	CASE ( 'p' )
[667]	363	DO i = nxlg, nxrg
	364	DO j = nysg, nyng
[1]	365	DO k = nzb, nzt+1
	366	p_av(k,j,i) = p_av(k,j,i) + p(k,j,i)
	367	ENDDO
	368	ENDDO
	369	ENDDO
	370
	371	CASE ( 'pc' )
	372	DO i = nxl, nxr
	373	DO j = nys, nyn
	374	DO k = nzb, nzt+1
	375	pc_av(k,j,i) = pc_av(k,j,i) + prt_count(k,j,i)
	376	ENDDO
	377	ENDDO
	378	ENDDO
	379
	380	CASE ( 'pr' )
	381	DO i = nxl, nxr
	382	DO j = nys, nyn
	383	DO k = nzb, nzt+1
[1359]	384	number_of_particles = prt_count(k,j,i)
	385	IF ( number_of_particles <= 0 ) CYCLE
	386	particles => grid_particles(k,j,i)%particles(1:number_of_particles)
	387	s_r2 = 0.0_wp
[1353]	388	s_r3 = 0.0_wp
[1359]	389
	390	DO n = 1, number_of_particles
	391	IF ( particles(n)%particle_mask ) THEN
	392	s_r2 = s_r2 + particles(n)%radius*2 &
	393	particles(n)%weight_factor
	394	s_r3 = s_r3 + particles(n)%radius*3 &
	395	particles(n)%weight_factor
	396	ENDIF
[1]	397	ENDDO
[1359]	398
	399	IF ( s_r2 > 0.0_wp ) THEN
	400	mean_r = s_r3 / s_r2
[1]	401	ELSE
[1353]	402	mean_r = 0.0_wp
[1]	403	ENDIF
	404	pr_av(k,j,i) = pr_av(k,j,i) + mean_r
	405	ENDDO
	406	ENDDO
	407	ENDDO
	408
[1359]	409
[72]	410	CASE ( 'pr*' )
[667]	411	DO i = nxlg, nxrg
	412	DO j = nysg, nyng
[72]	413	precipitation_rate_av(j,i) = precipitation_rate_av(j,i) + &
	414	precipitation_rate(j,i)
	415	ENDDO
	416	ENDDO
	417
[1]	418	CASE ( 'pt' )
	419	IF ( .NOT. cloud_physics ) THEN
[667]	420	DO i = nxlg, nxrg
	421	DO j = nysg, nyng
	422	DO k = nzb, nzt+1
[1]	423	pt_av(k,j,i) = pt_av(k,j,i) + pt(k,j,i)
	424	ENDDO
	425	ENDDO
	426	ENDDO
	427	ELSE
[667]	428	DO i = nxlg, nxrg
	429	DO j = nysg, nyng
	430	DO k = nzb, nzt+1
[1]	431	pt_av(k,j,i) = pt_av(k,j,i) + pt(k,j,i) + l_d_cp * &
	432	pt_d_t(k) * ql(k,j,i)
	433	ENDDO
	434	ENDDO
	435	ENDDO
	436	ENDIF
	437
	438	CASE ( 'q' )
[667]	439	DO i = nxlg, nxrg
	440	DO j = nysg, nyng
[1]	441	DO k = nzb, nzt+1
	442	q_av(k,j,i) = q_av(k,j,i) + q(k,j,i)
	443	ENDDO
	444	ENDDO
	445	ENDDO
[402]	446
[1115]	447	CASE ( 'qc' )
	448	DO i = nxlg, nxrg
	449	DO j = nysg, nyng
	450	DO k = nzb, nzt+1
	451	qc_av(k,j,i) = qc_av(k,j,i) + qc(k,j,i)
	452	ENDDO
	453	ENDDO
	454	ENDDO
	455
[1]	456	CASE ( 'ql' )
[667]	457	DO i = nxlg, nxrg
	458	DO j = nysg, nyng
[1]	459	DO k = nzb, nzt+1
	460	ql_av(k,j,i) = ql_av(k,j,i) + ql(k,j,i)
	461	ENDDO
	462	ENDDO
	463	ENDDO
	464
	465	CASE ( 'ql_c' )
[667]	466	DO i = nxlg, nxrg
	467	DO j = nysg, nyng
[1]	468	DO k = nzb, nzt+1
	469	ql_c_av(k,j,i) = ql_c_av(k,j,i) + ql_c(k,j,i)
	470	ENDDO
	471	ENDDO
	472	ENDDO
	473
	474	CASE ( 'ql_v' )
[667]	475	DO i = nxlg, nxrg
	476	DO j = nysg, nyng
[1]	477	DO k = nzb, nzt+1
	478	ql_v_av(k,j,i) = ql_v_av(k,j,i) + ql_v(k,j,i)
	479	ENDDO
	480	ENDDO
	481	ENDDO
	482
	483	CASE ( 'ql_vp' )
[1007]	484	DO i = nxl, nxr
	485	DO j = nys, nyn
[1]	486	DO k = nzb, nzt+1
[1359]	487	number_of_particles = prt_count(k,j,i)
	488	IF ( number_of_particles <= 0 ) CYCLE
	489	particles => grid_particles(k,j,i)%particles(1:number_of_particles)
	490	DO n = 1, number_of_particles
	491	IF ( particles(n)%particle_mask ) THEN
	492	ql_vp_av(k,j,i) = ql_vp_av(k,j,i) + &
	493	particles(n)%weight_factor / &
	494	number_of_particles
	495	ENDIF
[1007]	496	ENDDO
[1]	497	ENDDO
	498	ENDDO
	499	ENDDO
	500
[1053]	501	CASE ( 'qr' )
	502	DO i = nxlg, nxrg
	503	DO j = nysg, nyng
	504	DO k = nzb, nzt+1
	505	qr_av(k,j,i) = qr_av(k,j,i) + qr(k,j,i)
	506	ENDDO
	507	ENDDO
	508	ENDDO
	509
[402]	510	CASE ( 'qsws*' )
[667]	511	DO i = nxlg, nxrg
	512	DO j = nysg, nyng
[402]	513	qsws_av(j,i) = qsws_av(j,i) + qsws(j,i)
	514	ENDDO
	515	ENDDO
	516
[1]	517	CASE ( 'qv' )
[667]	518	DO i = nxlg, nxrg
	519	DO j = nysg, nyng
[1]	520	DO k = nzb, nzt+1
	521	qv_av(k,j,i) = qv_av(k,j,i) + q(k,j,i) - ql(k,j,i)
	522	ENDDO
	523	ENDDO
	524	ENDDO
	525
[96]	526	CASE ( 'rho' )
[667]	527	DO i = nxlg, nxrg
	528	DO j = nysg, nyng
[96]	529	DO k = nzb, nzt+1
	530	rho_av(k,j,i) = rho_av(k,j,i) + rho(k,j,i)
	531	ENDDO
	532	ENDDO
	533	ENDDO
[402]	534
[1]	535	CASE ( 's' )
[667]	536	DO i = nxlg, nxrg
	537	DO j = nysg, nyng
[1]	538	DO k = nzb, nzt+1
	539	s_av(k,j,i) = s_av(k,j,i) + q(k,j,i)
	540	ENDDO
	541	ENDDO
	542	ENDDO
[402]	543
[96]	544	CASE ( 'sa' )
[667]	545	DO i = nxlg, nxrg
	546	DO j = nysg, nyng
[96]	547	DO k = nzb, nzt+1
	548	sa_av(k,j,i) = sa_av(k,j,i) + sa(k,j,i)
	549	ENDDO
	550	ENDDO
	551	ENDDO
[402]	552
	553	CASE ( 'shf*' )
[667]	554	DO i = nxlg, nxrg
	555	DO j = nysg, nyng
[402]	556	shf_av(j,i) = shf_av(j,i) + shf(j,i)
	557	ENDDO
	558	ENDDO
	559
[1]	560	CASE ( 't*' )
[667]	561	DO i = nxlg, nxrg
	562	DO j = nysg, nyng
[1]	563	ts_av(j,i) = ts_av(j,i) + ts(j,i)
	564	ENDDO
	565	ENDDO
	566
	567	CASE ( 'u' )
[667]	568	DO i = nxlg, nxrg
	569	DO j = nysg, nyng
[1]	570	DO k = nzb, nzt+1
	571	u_av(k,j,i) = u_av(k,j,i) + u(k,j,i)
	572	ENDDO
	573	ENDDO
	574	ENDDO
	575
	576	CASE ( 'u*' )
[667]	577	DO i = nxlg, nxrg
	578	DO j = nysg, nyng
[1]	579	us_av(j,i) = us_av(j,i) + us(j,i)
	580	ENDDO
	581	ENDDO
	582
	583	CASE ( 'v' )
[667]	584	DO i = nxlg, nxrg
	585	DO j = nysg, nyng
[1]	586	DO k = nzb, nzt+1
	587	v_av(k,j,i) = v_av(k,j,i) + v(k,j,i)
	588	ENDDO
	589	ENDDO
	590	ENDDO
	591
	592	CASE ( 'vpt' )
[667]	593	DO i = nxlg, nxrg
	594	DO j = nysg, nyng
[1]	595	DO k = nzb, nzt+1
	596	vpt_av(k,j,i) = vpt_av(k,j,i) + vpt(k,j,i)
	597	ENDDO
	598	ENDDO
	599	ENDDO
	600
	601	CASE ( 'w' )
[667]	602	DO i = nxlg, nxrg
	603	DO j = nysg, nyng
[1]	604	DO k = nzb, nzt+1
	605	w_av(k,j,i) = w_av(k,j,i) + w(k,j,i)
	606	ENDDO
	607	ENDDO
	608	ENDDO
	609
[72]	610	CASE ( 'z0*' )
[667]	611	DO i = nxlg, nxrg
	612	DO j = nysg, nyng
[72]	613	z0_av(j,i) = z0_av(j,i) + z0(j,i)
	614	ENDDO
	615	ENDDO
	616
[978]	617	CASE ( 'z0h*' )
	618	DO i = nxlg, nxrg
	619	DO j = nysg, nyng
	620	z0h_av(j,i) = z0h_av(j,i) + z0h(j,i)
	621	ENDDO
	622	ENDDO
	623
[1]	624	CASE DEFAULT
	625	!
	626	!-- User-defined quantity
	627	CALL user_3d_data_averaging( 'sum', doav(ii) )
	628
	629	END SELECT
	630
	631	ENDDO
	632
[1318]	633	CALL cpu_log( log_point(34), 'sum_up_3d_data', 'stop' )
[1]	634
	635
	636	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 |