Home

Context Navigation

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

Last change on this file since 4327 was 4182, checked in by scharf, 5 years ago
corrected "Former revisions" section minor formatting in "Former revisions" section added "Author" section
Property svn:keywords set to `Id`
File size: 39.7 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-2019 Leibniz Universitaet Hannover
18	!------------------------------------------------------------------------------!
19	!
20	! Current revisions:
21	! ------------------
22	!
23	!
24	! Former revisions:
25	! -----------------
26	! $Id: sum_up_3d_data.f90 4182 2019-08-22 15:20:23Z Giersch $
27	! Corrected "Former revisions" section
28	!
29	! 4048 2019-06-21 21:00:21Z knoop
30	! Moved tcm_3d_data_averaging to module_interface
31	!
32	! 4039 2019-06-18 10:32:41Z suehring
33	! Modularize diagnostic output
34	!
35	! 3994 2019-05-22 18:08:09Z suehring
36	! output of turbulence intensity added
37	!
38	! 3943 2019-05-02 09:50:41Z maronga
39	! Added output of qsws_av for green roofs.
40	!
41	! 3933 2019-04-25 12:33:20Z kanani
42	! Formatting
43	!
44	! 3773 2019-03-01 08:56:57Z maronga
45	! Added output of theta_2m*_xy_av
46	!
47	! 3761 2019-02-25 15:31:42Z raasch
48	! unused variables removed
49	!
50	! 3655 2019-01-07 16:51:22Z knoop
51	! Implementation of the PALM module interface
52	!
53	! Revision 1.1 2006/02/23 12:55:23 raasch
54	! Initial revision
55	!
56	!
57	! Description:
58	! ------------
59	!> Sum-up the values of 3d-arrays. The real averaging is later done in routine
60	!> average_3d_data.
61	!------------------------------------------------------------------------------!
62	SUBROUTINE sum_up_3d_data
63
64
65	USE arrays_3d, &
66	ONLY: dzw, d_exner, e, heatflux_output_conversion, p, &
67	pt, q, ql, ql_c, ql_v, s, u, v, vpt, w, &
68	waterflux_output_conversion
69
70	USE averaging, &
71	ONLY: e_av, ghf_av, lpt_av, lwp_av, ol_av, p_av, pc_av, pr_av, pt_av, &
72	pt_2m_av, q_av, ql_av, ql_c_av, ql_v_av, ql_vp_av, qsws_av, &
73	qv_av, r_a_av, s_av, shf_av, ssws_av, ts_av, tsurf_av, u_av, &
74	us_av, v_av, vpt_av, w_av, z0_av, z0h_av, z0q_av
75
76	USE basic_constants_and_equations_mod, &
77	ONLY: c_p, lv_d_cp, l_v
78
79	USE bulk_cloud_model_mod, &
80	ONLY: bulk_cloud_model
81
82	USE control_parameters, &
83	ONLY: average_count_3d, doav, doav_n, rho_surface, urban_surface, &
84	varnamelength
85
86	USE cpulog, &
87	ONLY: cpu_log, log_point
88
89	USE indices, &
90	ONLY: nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, nzt
91
92	USE kinds
93
94	USE module_interface, &
95	ONLY: module_interface_3d_data_averaging
96
97	USE particle_attributes, &
98	ONLY: grid_particles, number_of_particles, particles, prt_count
99
100	USE surface_mod, &
101	ONLY: ind_pav_green, ind_veg_wall, ind_wat_win, &
102	surf_def_h, surf_lsm_h, surf_usm_h
103
104	USE urban_surface_mod, &
105	ONLY: usm_3d_data_averaging
106
107
108	IMPLICIT NONE
109
110	LOGICAL :: match_def !< flag indicating default-type surface
111	LOGICAL :: match_lsm !< flag indicating natural-type surface
112	LOGICAL :: match_usm !< flag indicating urban-type surface
113
114	INTEGER(iwp) :: i !< grid index x direction
115	INTEGER(iwp) :: ii !< running index
116	INTEGER(iwp) :: j !< grid index y direction
117	INTEGER(iwp) :: k !< grid index x direction
118	INTEGER(iwp) :: m !< running index over surfacle elements
119	INTEGER(iwp) :: n !< running index over number of particles per grid box
120
121	REAL(wp) :: mean_r !< mean-particle radius witin grid box
122	REAL(wp) :: s_r2 !< mean-particle radius witin grid box to the power of two
123	REAL(wp) :: s_r3 !< mean-particle radius witin grid box to the power of three
124
125	CHARACTER (LEN=varnamelength) :: trimvar !< TRIM of output-variable string
126
127
128	CALL cpu_log (log_point(34),'sum_up_3d_data','start')
129
130	!
131	!-- Allocate and initialize the summation arrays if called for the very first
132	!-- time or the first time after average_3d_data has been called
133	!-- (some or all of the arrays may have been already allocated
134	!-- in rrd_local)
135	IF ( average_count_3d == 0 ) THEN
136
137	DO ii = 1, doav_n
138
139	trimvar = TRIM( doav(ii) )
140
141	SELECT CASE ( trimvar )
142
143	CASE ( 'ghf*' )
144	IF ( .NOT. ALLOCATED( ghf_av ) ) THEN
145	ALLOCATE( ghf_av(nysg:nyng,nxlg:nxrg) )
146	ENDIF
147	ghf_av = 0.0_wp
148
149	CASE ( 'e' )
150	IF ( .NOT. ALLOCATED( e_av ) ) THEN
151	ALLOCATE( e_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
152	ENDIF
153	e_av = 0.0_wp
154
155	CASE ( 'thetal' )
156	IF ( .NOT. ALLOCATED( lpt_av ) ) THEN
157	ALLOCATE( lpt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
158	ENDIF
159	lpt_av = 0.0_wp
160
161	CASE ( 'lwp*' )
162	IF ( .NOT. ALLOCATED( lwp_av ) ) THEN
163	ALLOCATE( lwp_av(nysg:nyng,nxlg:nxrg) )
164	ENDIF
165	lwp_av = 0.0_wp
166
167	CASE ( 'ol*' )
168	IF ( .NOT. ALLOCATED( ol_av ) ) THEN
169	ALLOCATE( ol_av(nysg:nyng,nxlg:nxrg) )
170	ENDIF
171	ol_av = 0.0_wp
172
173	CASE ( 'p' )
174	IF ( .NOT. ALLOCATED( p_av ) ) THEN
175	ALLOCATE( p_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
176	ENDIF
177	p_av = 0.0_wp
178
179	CASE ( 'pc' )
180	IF ( .NOT. ALLOCATED( pc_av ) ) THEN
181	ALLOCATE( pc_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
182	ENDIF
183	pc_av = 0.0_wp
184
185	CASE ( 'pr' )
186	IF ( .NOT. ALLOCATED( pr_av ) ) THEN
187	ALLOCATE( pr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
188	ENDIF
189	pr_av = 0.0_wp
190
191	CASE ( 'theta' )
192	IF ( .NOT. ALLOCATED( pt_av ) ) THEN
193	ALLOCATE( pt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
194	ENDIF
195	pt_av = 0.0_wp
196
197	CASE ( 'q' )
198	IF ( .NOT. ALLOCATED( q_av ) ) THEN
199	ALLOCATE( q_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
200	ENDIF
201	q_av = 0.0_wp
202
203	CASE ( 'ql' )
204	IF ( .NOT. ALLOCATED( ql_av ) ) THEN
205	ALLOCATE( ql_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
206	ENDIF
207	ql_av = 0.0_wp
208
209	CASE ( 'ql_c' )
210	IF ( .NOT. ALLOCATED( ql_c_av ) ) THEN
211	ALLOCATE( ql_c_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
212	ENDIF
213	ql_c_av = 0.0_wp
214
215	CASE ( 'ql_v' )
216	IF ( .NOT. ALLOCATED( ql_v_av ) ) THEN
217	ALLOCATE( ql_v_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
218	ENDIF
219	ql_v_av = 0.0_wp
220
221	CASE ( 'ql_vp' )
222	IF ( .NOT. ALLOCATED( ql_vp_av ) ) THEN
223	ALLOCATE( ql_vp_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
224	ENDIF
225	ql_vp_av = 0.0_wp
226
227	CASE ( 'qsws*' )
228	IF ( .NOT. ALLOCATED( qsws_av ) ) THEN
229	ALLOCATE( qsws_av(nysg:nyng,nxlg:nxrg) )
230	ENDIF
231	qsws_av = 0.0_wp
232
233	CASE ( 'qv' )
234	IF ( .NOT. ALLOCATED( qv_av ) ) THEN
235	ALLOCATE( qv_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
236	ENDIF
237	qv_av = 0.0_wp
238
239	CASE ( 'r_a*' )
240	IF ( .NOT. ALLOCATED( r_a_av ) ) THEN
241	ALLOCATE( r_a_av(nysg:nyng,nxlg:nxrg) )
242	ENDIF
243	r_a_av = 0.0_wp
244
245	CASE ( 's' )
246	IF ( .NOT. ALLOCATED( s_av ) ) THEN
247	ALLOCATE( s_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
248	ENDIF
249	s_av = 0.0_wp
250
251	CASE ( 'shf*' )
252	IF ( .NOT. ALLOCATED( shf_av ) ) THEN
253	ALLOCATE( shf_av(nysg:nyng,nxlg:nxrg) )
254	ENDIF
255	shf_av = 0.0_wp
256
257	CASE ( 'ssws*' )
258	IF ( .NOT. ALLOCATED( ssws_av ) ) THEN
259	ALLOCATE( ssws_av(nysg:nyng,nxlg:nxrg) )
260	ENDIF
261	ssws_av = 0.0_wp
262
263	CASE ( 't*' )
264	IF ( .NOT. ALLOCATED( ts_av ) ) THEN
265	ALLOCATE( ts_av(nysg:nyng,nxlg:nxrg) )
266	ENDIF
267	ts_av = 0.0_wp
268
269	CASE ( 'tsurf*' )
270	IF ( .NOT. ALLOCATED( tsurf_av ) ) THEN
271	ALLOCATE( tsurf_av(nysg:nyng,nxlg:nxrg) )
272	ENDIF
273	tsurf_av = 0.0_wp
274
275	CASE ( 'u' )
276	IF ( .NOT. ALLOCATED( u_av ) ) THEN
277	ALLOCATE( u_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
278	ENDIF
279	u_av = 0.0_wp
280
281	CASE ( 'us*' )
282	IF ( .NOT. ALLOCATED( us_av ) ) THEN
283	ALLOCATE( us_av(nysg:nyng,nxlg:nxrg) )
284	ENDIF
285	us_av = 0.0_wp
286
287	CASE ( 'v' )
288	IF ( .NOT. ALLOCATED( v_av ) ) THEN
289	ALLOCATE( v_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
290	ENDIF
291	v_av = 0.0_wp
292
293	CASE ( 'thetav' )
294	IF ( .NOT. ALLOCATED( vpt_av ) ) THEN
295	ALLOCATE( vpt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
296	ENDIF
297	vpt_av = 0.0_wp
298
299	CASE ( 'theta_2m*' )
300	IF ( .NOT. ALLOCATED( pt_2m_av ) ) THEN
301	ALLOCATE( pt_2m_av(nysg:nyng,nxlg:nxrg) )
302	ENDIF
303	pt_2m_av = 0.0_wp
304
305	CASE ( 'w' )
306	IF ( .NOT. ALLOCATED( w_av ) ) THEN
307	ALLOCATE( w_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
308	ENDIF
309	w_av = 0.0_wp
310
311	CASE ( 'z0*' )
312	IF ( .NOT. ALLOCATED( z0_av ) ) THEN
313	ALLOCATE( z0_av(nysg:nyng,nxlg:nxrg) )
314	ENDIF
315	z0_av = 0.0_wp
316
317	CASE ( 'z0h*' )
318	IF ( .NOT. ALLOCATED( z0h_av ) ) THEN
319	ALLOCATE( z0h_av(nysg:nyng,nxlg:nxrg) )
320	ENDIF
321	z0h_av = 0.0_wp
322
323	CASE ( 'z0q*' )
324	IF ( .NOT. ALLOCATED( z0q_av ) ) THEN
325	ALLOCATE( z0q_av(nysg:nyng,nxlg:nxrg) )
326	ENDIF
327	z0q_av = 0.0_wp
328
329
330	CASE DEFAULT
331
332	!
333	!-- Allocating and initializing data arrays for all other modules
334	CALL module_interface_3d_data_averaging( 'allocate', trimvar )
335
336
337	END SELECT
338
339	ENDDO
340
341	ENDIF
342
343	!
344	!-- Loop of all variables to be averaged.
345	DO ii = 1, doav_n
346
347	trimvar = TRIM( doav(ii) )
348	!
349	!-- Store the array chosen on the temporary array.
350	SELECT CASE ( trimvar )
351
352	CASE ( 'ghf*' )
353	IF ( ALLOCATED( ghf_av ) ) THEN
354	DO i = nxl, nxr
355	DO j = nys, nyn
356	!
357	!-- Check whether grid point is a natural- or urban-type
358	!-- surface.
359	match_lsm = surf_lsm_h%start_index(j,i) <= &
360	surf_lsm_h%end_index(j,i)
361	match_usm = surf_usm_h%start_index(j,i) <= &
362	surf_usm_h%end_index(j,i)
363	!
364	!-- In order to avoid double-counting of surface properties,
365	!-- always assume that natural-type surfaces are below urban-
366	!-- type surfaces, e.g. in case of bridges.
367	!-- Further, take only the last suface element, i.e. the
368	!-- uppermost surface which would be visible from above
369	IF ( match_lsm .AND. .NOT. match_usm ) THEN
370	m = surf_lsm_h%end_index(j,i)
371	ghf_av(j,i) = ghf_av(j,i) + &
372	surf_lsm_h%ghf(m)
373	ELSEIF ( match_usm ) THEN
374	m = surf_usm_h%end_index(j,i)
375	ghf_av(j,i) = ghf_av(j,i) + &
376	surf_usm_h%frac(ind_veg_wall,m) * &
377	surf_usm_h%wghf_eb(m) + &
378	surf_usm_h%frac(ind_pav_green,m) * &
379	surf_usm_h%wghf_eb_green(m) + &
380	surf_usm_h%frac(ind_wat_win,m) * &
381	surf_usm_h%wghf_eb_window(m)
382	ENDIF
383	ENDDO
384	ENDDO
385	ENDIF
386
387	CASE ( 'e' )
388	IF ( ALLOCATED( e_av ) ) THEN
389	DO i = nxlg, nxrg
390	DO j = nysg, nyng
391	DO k = nzb, nzt+1
392	e_av(k,j,i) = e_av(k,j,i) + e(k,j,i)
393	ENDDO
394	ENDDO
395	ENDDO
396	ENDIF
397
398	CASE ( 'thetal' )
399	IF ( ALLOCATED( lpt_av ) ) THEN
400	DO i = nxlg, nxrg
401	DO j = nysg, nyng
402	DO k = nzb, nzt+1
403	lpt_av(k,j,i) = lpt_av(k,j,i) + pt(k,j,i)
404	ENDDO
405	ENDDO
406	ENDDO
407	ENDIF
408
409	CASE ( 'lwp*' )
410	IF ( ALLOCATED( lwp_av ) ) THEN
411	DO i = nxlg, nxrg
412	DO j = nysg, nyng
413	lwp_av(j,i) = lwp_av(j,i) + SUM( ql(nzb:nzt,j,i) &
414	* dzw(1:nzt+1) ) * rho_surface
415	ENDDO
416	ENDDO
417	ENDIF
418
419	CASE ( 'ol*' )
420	IF ( ALLOCATED( ol_av ) ) THEN
421	DO i = nxl, nxr
422	DO j = nys, nyn
423	match_def = surf_def_h(0)%start_index(j,i) <= &
424	surf_def_h(0)%end_index(j,i)
425	match_lsm = surf_lsm_h%start_index(j,i) <= &
426	surf_lsm_h%end_index(j,i)
427	match_usm = surf_usm_h%start_index(j,i) <= &
428	surf_usm_h%end_index(j,i)
429
430	IF ( match_def ) THEN
431	m = surf_def_h(0)%end_index(j,i)
432	ol_av(j,i) = ol_av(j,i) + &
433	surf_def_h(0)%ol(m)
434	ELSEIF ( match_lsm .AND. .NOT. match_usm ) THEN
435	m = surf_lsm_h%end_index(j,i)
436	ol_av(j,i) = ol_av(j,i) + &
437	surf_lsm_h%ol(m)
438	ELSEIF ( match_usm ) THEN
439	m = surf_usm_h%end_index(j,i)
440	ol_av(j,i) = ol_av(j,i) + &
441	surf_usm_h%ol(m)
442	ENDIF
443	ENDDO
444	ENDDO
445	ENDIF
446
447	CASE ( 'p' )
448	IF ( ALLOCATED( p_av ) ) THEN
449	DO i = nxlg, nxrg
450	DO j = nysg, nyng
451	DO k = nzb, nzt+1
452	p_av(k,j,i) = p_av(k,j,i) + p(k,j,i)
453	ENDDO
454	ENDDO
455	ENDDO
456	ENDIF
457
458	CASE ( 'pc' )
459	IF ( ALLOCATED( pc_av ) ) THEN
460	DO i = nxl, nxr
461	DO j = nys, nyn
462	DO k = nzb, nzt+1
463	pc_av(k,j,i) = pc_av(k,j,i) + prt_count(k,j,i)
464	ENDDO
465	ENDDO
466	ENDDO
467	ENDIF
468
469	CASE ( 'pr' )
470	IF ( ALLOCATED( pr_av ) ) THEN
471	DO i = nxl, nxr
472	DO j = nys, nyn
473	DO k = nzb, nzt+1
474	number_of_particles = prt_count(k,j,i)
475	IF ( number_of_particles <= 0 ) CYCLE
476	particles => &
477	grid_particles(k,j,i)%particles(1:number_of_particles)
478	s_r2 = 0.0_wp
479	s_r3 = 0.0_wp
480
481	DO n = 1, number_of_particles
482	IF ( particles(n)%particle_mask ) THEN
483	s_r2 = s_r2 + particles(n)%radius*2 &
484	particles(n)%weight_factor
485	s_r3 = s_r3 + particles(n)%radius*3 &
486	particles(n)%weight_factor
487	ENDIF
488	ENDDO
489
490	IF ( s_r2 > 0.0_wp ) THEN
491	mean_r = s_r3 / s_r2
492	ELSE
493	mean_r = 0.0_wp
494	ENDIF
495	pr_av(k,j,i) = pr_av(k,j,i) + mean_r
496	ENDDO
497	ENDDO
498	ENDDO
499	ENDIF
500
501	CASE ( 'theta' )
502	IF ( ALLOCATED( pt_av ) ) THEN
503	IF ( .NOT. bulk_cloud_model ) THEN
504	DO i = nxlg, nxrg
505	DO j = nysg, nyng
506	DO k = nzb, nzt+1
507	pt_av(k,j,i) = pt_av(k,j,i) + pt(k,j,i)
508	ENDDO
509	ENDDO
510	ENDDO
511	ELSE
512	DO i = nxlg, nxrg
513	DO j = nysg, nyng
514	DO k = nzb, nzt+1
515	pt_av(k,j,i) = pt_av(k,j,i) + pt(k,j,i) + lv_d_cp * &
516	d_exner(k) * ql(k,j,i)
517	ENDDO
518	ENDDO
519	ENDDO
520	ENDIF
521	ENDIF
522
523	CASE ( 'q' )
524	IF ( ALLOCATED( q_av ) ) THEN
525	DO i = nxlg, nxrg
526	DO j = nysg, nyng
527	DO k = nzb, nzt+1
528	q_av(k,j,i) = q_av(k,j,i) + q(k,j,i)
529	ENDDO
530	ENDDO
531	ENDDO
532	ENDIF
533
534	CASE ( 'ql' )
535	IF ( ALLOCATED( ql_av ) ) THEN
536	DO i = nxlg, nxrg
537	DO j = nysg, nyng
538	DO k = nzb, nzt+1
539	ql_av(k,j,i) = ql_av(k,j,i) + ql(k,j,i)
540	ENDDO
541	ENDDO
542	ENDDO
543	ENDIF
544
545	CASE ( 'ql_c' )
546	IF ( ALLOCATED( ql_c_av ) ) THEN
547	DO i = nxlg, nxrg
548	DO j = nysg, nyng
549	DO k = nzb, nzt+1
550	ql_c_av(k,j,i) = ql_c_av(k,j,i) + ql_c(k,j,i)
551	ENDDO
552	ENDDO
553	ENDDO
554	ENDIF
555
556	CASE ( 'ql_v' )
557	IF ( ALLOCATED( ql_v_av ) ) THEN
558	DO i = nxlg, nxrg
559	DO j = nysg, nyng
560	DO k = nzb, nzt+1
561	ql_v_av(k,j,i) = ql_v_av(k,j,i) + ql_v(k,j,i)
562	ENDDO
563	ENDDO
564	ENDDO
565	ENDIF
566
567	CASE ( 'ql_vp' )
568	IF ( ALLOCATED( ql_vp_av ) ) THEN
569	DO i = nxl, nxr
570	DO j = nys, nyn
571	DO k = nzb, nzt+1
572	number_of_particles = prt_count(k,j,i)
573	IF ( number_of_particles <= 0 ) CYCLE
574	particles => &
575	grid_particles(k,j,i)%particles(1:number_of_particles)
576	DO n = 1, number_of_particles
577	IF ( particles(n)%particle_mask ) THEN
578	ql_vp_av(k,j,i) = ql_vp_av(k,j,i) + &
579	particles(n)%weight_factor / &
580	number_of_particles
581	ENDIF
582	ENDDO
583	ENDDO
584	ENDDO
585	ENDDO
586	ENDIF
587
588	CASE ( 'qsws*' )
589	!
590	!-- In case of default surfaces, clean-up flux by density.
591	!-- In case of land- and urban-surfaces, convert fluxes into
592	!-- dynamic units.
593	!-- Question (maronga): are the .NOT. statements really required?
594	IF ( ALLOCATED( qsws_av ) ) THEN
595	DO i = nxl, nxr
596	DO j = nys, nyn
597	match_def = surf_def_h(0)%start_index(j,i) <= &
598	surf_def_h(0)%end_index(j,i)
599	match_lsm = surf_lsm_h%start_index(j,i) <= &
600	surf_lsm_h%end_index(j,i)
601	match_usm = surf_usm_h%start_index(j,i) <= &
602	surf_usm_h%end_index(j,i)
603
604	IF ( match_def ) THEN
605	m = surf_def_h(0)%end_index(j,i)
606	qsws_av(j,i) = qsws_av(j,i) + &
607	surf_def_h(0)%qsws(m) * &
608	waterflux_output_conversion(nzb)
609	ELSEIF ( match_lsm .AND. .NOT. match_usm ) THEN
610	m = surf_lsm_h%end_index(j,i)
611	qsws_av(j,i) = qsws_av(j,i) + &
612	surf_lsm_h%qsws(m) * l_v
613	ELSEIF ( match_usm .AND. .NOT. match_lsm ) THEN
614	m = surf_usm_h%end_index(j,i)
615	qsws_av(j,i) = qsws_av(j,i) + &
616	surf_usm_h%qsws(m) * l_v
617	ENDIF
618	ENDDO
619	ENDDO
620	ENDIF
621
622	CASE ( 'qv' )
623	IF ( ALLOCATED( qv_av ) ) THEN
624	DO i = nxlg, nxrg
625	DO j = nysg, nyng
626	DO k = nzb, nzt+1
627	qv_av(k,j,i) = qv_av(k,j,i) + q(k,j,i) - ql(k,j,i)
628	ENDDO
629	ENDDO
630	ENDDO
631	ENDIF
632
633	CASE ( 'r_a*' )
634	IF ( ALLOCATED( r_a_av ) ) THEN
635	DO i = nxl, nxr
636	DO j = nys, nyn
637	match_lsm = surf_lsm_h%start_index(j,i) <= &
638	surf_lsm_h%end_index(j,i)
639	match_usm = surf_usm_h%start_index(j,i) <= &
640	surf_usm_h%end_index(j,i)
641
642	IF ( match_lsm .AND. .NOT. match_usm ) THEN
643	m = surf_lsm_h%end_index(j,i)
644	r_a_av(j,i) = r_a_av(j,i) + &
645	surf_lsm_h%r_a(m)
646	ELSEIF ( match_usm ) THEN
647	m = surf_usm_h%end_index(j,i)
648	r_a_av(j,i) = r_a_av(j,i) + &
649	surf_usm_h%frac(ind_veg_wall,m) * &
650	surf_usm_h%r_a(m) + &
651	surf_usm_h%frac(ind_pav_green,m) * &
652	surf_usm_h%r_a_green(m) + &
653	surf_usm_h%frac(ind_wat_win,m) * &
654	surf_usm_h%r_a_window(m)
655	ENDIF
656	ENDDO
657	ENDDO
658	ENDIF
659
660	CASE ( 's' )
661	IF ( ALLOCATED( s_av ) ) THEN
662	DO i = nxlg, nxrg
663	DO j = nysg, nyng
664	DO k = nzb, nzt+1
665	s_av(k,j,i) = s_av(k,j,i) + s(k,j,i)
666	ENDDO
667	ENDDO
668	ENDDO
669	ENDIF
670
671	CASE ( 'shf*' )
672	!
673	!-- In case of default surfaces, clean-up flux by density.
674	!-- In case of land- and urban-surfaces, convert fluxes into
675	!-- dynamic units.
676	IF ( ALLOCATED( shf_av ) ) THEN
677	DO i = nxl, nxr
678	DO j = nys, nyn
679	match_def = surf_def_h(0)%start_index(j,i) <= &
680	surf_def_h(0)%end_index(j,i)
681	match_lsm = surf_lsm_h%start_index(j,i) <= &
682	surf_lsm_h%end_index(j,i)
683	match_usm = surf_usm_h%start_index(j,i) <= &
684	surf_usm_h%end_index(j,i)
685
686	IF ( match_def ) THEN
687	m = surf_def_h(0)%end_index(j,i)
688	shf_av(j,i) = shf_av(j,i) + &
689	surf_def_h(0)%shf(m) * &
690	heatflux_output_conversion(nzb)
691	ELSEIF ( match_lsm .AND. .NOT. match_usm ) THEN
692	m = surf_lsm_h%end_index(j,i)
693	shf_av(j,i) = shf_av(j,i) + &
694	surf_lsm_h%shf(m) * c_p
695	ELSEIF ( match_usm ) THEN
696	m = surf_usm_h%end_index(j,i)
697	shf_av(j,i) = shf_av(j,i) + &
698	surf_usm_h%shf(m) * c_p
699	ENDIF
700	ENDDO
701	ENDDO
702	ENDIF
703
704	CASE ( 'ssws*' )
705	IF ( ALLOCATED( ssws_av ) ) THEN
706	DO i = nxl, nxr
707	DO j = nys, nyn
708	match_def = surf_def_h(0)%start_index(j,i) <= &
709	surf_def_h(0)%end_index(j,i)
710	match_lsm = surf_lsm_h%start_index(j,i) <= &
711	surf_lsm_h%end_index(j,i)
712	match_usm = surf_usm_h%start_index(j,i) <= &
713	surf_usm_h%end_index(j,i)
714
715	IF ( match_def ) THEN
716	m = surf_def_h(0)%end_index(j,i)
717	ssws_av(j,i) = ssws_av(j,i) + &
718	surf_def_h(0)%ssws(m)
719	ELSEIF ( match_lsm .AND. .NOT. match_usm ) THEN
720	m = surf_lsm_h%end_index(j,i)
721	ssws_av(j,i) = ssws_av(j,i) + &
722	surf_lsm_h%ssws(m)
723	ELSEIF ( match_usm ) THEN
724	m = surf_usm_h%end_index(j,i)
725	ssws_av(j,i) = ssws_av(j,i) + &
726	surf_usm_h%ssws(m)
727	ENDIF
728	ENDDO
729	ENDDO
730	ENDIF
731
732	CASE ( 'theta_2m*' )
733	IF ( ALLOCATED( pt_2m_av ) ) THEN
734	DO i = nxl, nxr
735	DO j = nys, nyn
736	match_def = surf_def_h(0)%start_index(j,i) <= &
737	surf_def_h(0)%end_index(j,i)
738	match_lsm = surf_lsm_h%start_index(j,i) <= &
739	surf_lsm_h%end_index(j,i)
740	match_usm = surf_usm_h%start_index(j,i) <= &
741	surf_usm_h%end_index(j,i)
742
743	IF ( match_def ) THEN
744	m = surf_def_h(0)%end_index(j,i)
745	pt_2m_av(j,i) = pt_2m_av(j,i) + &
746	surf_def_h(0)%pt_2m(m)
747	ELSEIF ( match_lsm .AND. .NOT. match_usm ) THEN
748	m = surf_lsm_h%end_index(j,i)
749	pt_2m_av(j,i) = pt_2m_av(j,i) + &
750	surf_lsm_h%pt_2m(m)
751	ELSEIF ( match_usm ) THEN
752	m = surf_usm_h%end_index(j,i)
753	pt_2m_av(j,i) = pt_2m_av(j,i) + &
754	surf_usm_h%pt_2m(m)
755	ENDIF
756	ENDDO
757	ENDDO
758	ENDIF
759
760
761	CASE ( 't*' )
762	IF ( ALLOCATED( ts_av ) ) THEN
763	DO i = nxl, nxr
764	DO j = nys, nyn
765	match_def = surf_def_h(0)%start_index(j,i) <= &
766	surf_def_h(0)%end_index(j,i)
767	match_lsm = surf_lsm_h%start_index(j,i) <= &
768	surf_lsm_h%end_index(j,i)
769	match_usm = surf_usm_h%start_index(j,i) <= &
770	surf_usm_h%end_index(j,i)
771
772	IF ( match_def ) THEN
773	m = surf_def_h(0)%end_index(j,i)
774	ts_av(j,i) = ts_av(j,i) + &
775	surf_def_h(0)%ts(m)
776	ELSEIF ( match_lsm .AND. .NOT. match_usm ) THEN
777	m = surf_lsm_h%end_index(j,i)
778	ts_av(j,i) = ts_av(j,i) + &
779	surf_lsm_h%ts(m)
780	ELSEIF ( match_usm ) THEN
781	m = surf_usm_h%end_index(j,i)
782	ts_av(j,i) = ts_av(j,i) + &
783	surf_usm_h%ts(m)
784	ENDIF
785	ENDDO
786	ENDDO
787	ENDIF
788
789	CASE ( 'tsurf*' )
790	IF ( ALLOCATED( tsurf_av ) ) THEN
791	DO i = nxl, nxr
792	DO j = nys, nyn
793	match_def = surf_def_h(0)%start_index(j,i) <= &
794	surf_def_h(0)%end_index(j,i)
795	match_lsm = surf_lsm_h%start_index(j,i) <= &
796	surf_lsm_h%end_index(j,i)
797	match_usm = surf_usm_h%start_index(j,i) <= &
798	surf_usm_h%end_index(j,i)
799
800	IF ( match_def ) THEN
801	m = surf_def_h(0)%end_index(j,i)
802	tsurf_av(j,i) = tsurf_av(j,i) + &
803	surf_def_h(0)%pt_surface(m)
804	ELSEIF ( match_lsm .AND. .NOT. match_usm ) THEN
805	m = surf_lsm_h%end_index(j,i)
806	tsurf_av(j,i) = tsurf_av(j,i) + &
807	surf_lsm_h%pt_surface(m)
808	ELSEIF ( match_usm ) THEN
809	m = surf_usm_h%end_index(j,i)
810	tsurf_av(j,i) = tsurf_av(j,i) + &
811	surf_usm_h%pt_surface(m)
812	ENDIF
813	ENDDO
814	ENDDO
815	ENDIF
816
817	CASE ( 'u' )
818	IF ( ALLOCATED( u_av ) ) THEN
819	DO i = nxlg, nxrg
820	DO j = nysg, nyng
821	DO k = nzb, nzt+1
822	u_av(k,j,i) = u_av(k,j,i) + u(k,j,i)
823	ENDDO
824	ENDDO
825	ENDDO
826	ENDIF
827
828	CASE ( 'us*' )
829	IF ( ALLOCATED( us_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	us_av(j,i) = us_av(j,i) + &
842	surf_def_h(0)%us(m)
843	ELSEIF ( match_lsm .AND. .NOT. match_usm ) THEN
844	m = surf_lsm_h%end_index(j,i)
845	us_av(j,i) = us_av(j,i) + &
846	surf_lsm_h%us(m)
847	ELSEIF ( match_usm ) THEN
848	m = surf_usm_h%end_index(j,i)
849	us_av(j,i) = us_av(j,i) + &
850	surf_usm_h%us(m)
851	ENDIF
852	ENDDO
853	ENDDO
854	ENDIF
855
856	CASE ( 'v' )
857	IF ( ALLOCATED( v_av ) ) THEN
858	DO i = nxlg, nxrg
859	DO j = nysg, nyng
860	DO k = nzb, nzt+1
861	v_av(k,j,i) = v_av(k,j,i) + v(k,j,i)
862	ENDDO
863	ENDDO
864	ENDDO
865	ENDIF
866
867	CASE ( 'thetav' )
868	IF ( ALLOCATED( vpt_av ) ) THEN
869	DO i = nxlg, nxrg
870	DO j = nysg, nyng
871	DO k = nzb, nzt+1
872	vpt_av(k,j,i) = vpt_av(k,j,i) + vpt(k,j,i)
873	ENDDO
874	ENDDO
875	ENDDO
876	ENDIF
877
878	CASE ( 'w' )
879	IF ( ALLOCATED( w_av ) ) THEN
880	DO i = nxlg, nxrg
881	DO j = nysg, nyng
882	DO k = nzb, nzt+1
883	w_av(k,j,i) = w_av(k,j,i) + w(k,j,i)
884	ENDDO
885	ENDDO
886	ENDDO
887	ENDIF
888
889	CASE ( 'z0*' )
890	IF ( ALLOCATED( z0_av ) ) THEN
891	DO i = nxl, nxr
892	DO j = nys, nyn
893	match_def = surf_def_h(0)%start_index(j,i) <= &
894	surf_def_h(0)%end_index(j,i)
895	match_lsm = surf_lsm_h%start_index(j,i) <= &
896	surf_lsm_h%end_index(j,i)
897	match_usm = surf_usm_h%start_index(j,i) <= &
898	surf_usm_h%end_index(j,i)
899
900	IF ( match_def ) THEN
901	m = surf_def_h(0)%end_index(j,i)
902	z0_av(j,i) = z0_av(j,i) + &
903	surf_def_h(0)%z0(m)
904	ELSEIF ( match_lsm .AND. .NOT. match_usm ) THEN
905	m = surf_lsm_h%end_index(j,i)
906	z0_av(j,i) = z0_av(j,i) + &
907	surf_lsm_h%z0(m)
908	ELSEIF ( match_usm ) THEN
909	m = surf_usm_h%end_index(j,i)
910	z0_av(j,i) = z0_av(j,i) + &
911	surf_usm_h%z0(m)
912	ENDIF
913	ENDDO
914	ENDDO
915	ENDIF
916
917	CASE ( 'z0h*' )
918	IF ( ALLOCATED( z0h_av ) ) THEN
919	DO i = nxl, nxr
920	DO j = nys, nyn
921	match_def = surf_def_h(0)%start_index(j,i) <= &
922	surf_def_h(0)%end_index(j,i)
923	match_lsm = surf_lsm_h%start_index(j,i) <= &
924	surf_lsm_h%end_index(j,i)
925	match_usm = surf_usm_h%start_index(j,i) <= &
926	surf_usm_h%end_index(j,i)
927
928	IF ( match_def ) THEN
929	m = surf_def_h(0)%end_index(j,i)
930	z0h_av(j,i) = z0h_av(j,i) + &
931	surf_def_h(0)%z0h(m)
932	ELSEIF ( match_lsm .AND. .NOT. match_usm ) THEN
933	m = surf_lsm_h%end_index(j,i)
934	z0h_av(j,i) = z0h_av(j,i) + &
935	surf_lsm_h%z0h(m)
936	ELSEIF ( match_usm ) THEN
937	m = surf_usm_h%end_index(j,i)
938	z0h_av(j,i) = z0h_av(j,i) + &
939	surf_usm_h%z0h(m)
940	ENDIF
941	ENDDO
942	ENDDO
943	ENDIF
944
945	CASE ( 'z0q*' )
946	IF ( ALLOCATED( z0q_av ) ) THEN
947	DO i = nxl, nxr
948	DO j = nys, nyn
949	match_def = surf_def_h(0)%start_index(j,i) <= &
950	surf_def_h(0)%end_index(j,i)
951	match_lsm = surf_lsm_h%start_index(j,i) <= &
952	surf_lsm_h%end_index(j,i)
953	match_usm = surf_usm_h%start_index(j,i) <= &
954	surf_usm_h%end_index(j,i)
955
956	IF ( match_def ) THEN
957	m = surf_def_h(0)%end_index(j,i)
958	z0q_av(j,i) = z0q_av(j,i) + &
959	surf_def_h(0)%z0q(m)
960	ELSEIF ( match_lsm .AND. .NOT. match_usm ) THEN
961	m = surf_lsm_h%end_index(j,i)
962	z0q_av(j,i) = z0q_av(j,i) + &
963	surf_lsm_h%z0q(m)
964	ELSEIF ( match_usm ) THEN
965	m = surf_usm_h%end_index(j,i)
966	z0q_av(j,i) = z0q_av(j,i) + &
967	surf_usm_h%z0q(m)
968	ENDIF
969	ENDDO
970	ENDDO
971	ENDIF
972
973	CASE DEFAULT
974
975	!-- In case of urban surface variables it should be always checked
976	!-- if respective arrays are allocated, at least in case of a restart
977	!-- run, as averaged usm arrays are not read from file at the moment.
978	IF ( urban_surface ) THEN
979	CALL usm_3d_data_averaging( 'allocate', trimvar )
980	ENDIF
981
982	!
983	!-- Summing up data from all other modules
984	CALL module_interface_3d_data_averaging( 'sum', trimvar )
985
986
987	END SELECT
988
989	ENDDO
990
991	CALL cpu_log( log_point(34), 'sum_up_3d_data', 'stop' )
992
993
994	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 |