Home

Context Navigation

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

Last change on this file since 1976 was 1976, checked in by maronga, 8 years ago
further modularization of land surface model (2D/3D output and restart data). Bugfix for restart runs without land surface model
Property svn:keywords set to `Id`
File size: 22.9 KB

Line
1	!> @file sum_up_3d_data.f90
2	!--------------------------------------------------------------------------------!
3	! This file is part of PALM.
4	!
5	! PALM is free software: you can redistribute it and/or modify it under the terms
6	! of the GNU General Public License as published by the Free Software Foundation,
7	! either version 3 of the License, or (at your option) any later version.
8	!
9	! PALM is distributed in the hope that it will be useful, but WITHOUT ANY
10	! WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
11	! A PARTICULAR PURPOSE. See the GNU General Public License for more details.
12	!
13	! You should have received a copy of the GNU General Public License along with
14	! PALM. If not, see <http://www.gnu.org/licenses/>.
15	!
16	! Copyright 1997-2016 Leibniz Universitaet Hannover
17	!--------------------------------------------------------------------------------!
18	!
19	! Current revisions:
20	! -----------------
21	! Radiation actions are now done directly in the respective module
22	!
23	! Former revisions:
24	! -----------------
25	! $Id: sum_up_3d_data.f90 1976 2016-07-27 13:28:04Z maronga $
26	!
27	! 1972 2016-07-26 07:52:02Z maronga
28	! Land surface actions are now done directly in the respective module
29	!
30	! 1960 2016-07-12 16:34:24Z suehring
31	! Scalar surface flux added
32	!
33	! 1949 2016-06-17 07:19:16Z maronga
34	! Bugfix: calculation of lai_av, c_veg_av and c_liq_av.
35	!
36	! 1849 2016-04-08 11:33:18Z hoffmann
37	! precipitation_rate moved to arrays_3d
38	!
39	! 1788 2016-03-10 11:01:04Z maronga
40	! Added z0q and z0q_av
41	!
42	! 1693 2015-10-27 08:35:45Z maronga
43	! Last revision text corrected
44	!
45	! 1691 2015-10-26 16:17:44Z maronga
46	! Added output of Obukhov length and radiative heating rates for RRTMG.
47	! Corrected output of liquid water path.
48	!
49	! 1682 2015-10-07 23:56:08Z knoop
50	! Code annotations made doxygen readable
51	!
52	! 1585 2015-04-30 07:05:52Z maronga
53	! Adapted for RRTMG
54	!
55	! 1555 2015-03-04 17:44:27Z maronga
56	! Added output of r_a and r_s
57	!
58	! 1551 2015-03-03 14:18:16Z maronga
59	! Added support for land surface model and radiation model data.
60	!
61	! 1359 2014-04-11 17:15:14Z hoffmann
62	! New particle structure integrated.
63	!
64	! 1353 2014-04-08 15:21:23Z heinze
65	! REAL constants provided with KIND-attribute
66	!
67	! 1320 2014-03-20 08:40:49Z raasch
68	! ONLY-attribute added to USE-statements,
69	! kind-parameters added to all INTEGER and REAL declaration statements,
70	! kinds are defined in new module kinds,
71	! old module precision_kind is removed,
72	! revision history before 2012 removed,
73	! comment fields (!:) to be used for variable explanations added to
74	! all variable declaration statements
75	!
76	! 1318 2014-03-17 13:35:16Z raasch
77	! barrier argument removed from cpu_log,
78	! module interfaces removed
79	!
80	! 1115 2013-03-26 18:16:16Z hoffmann
81	! ql is calculated by calc_liquid_water_content
82	!
83	! 1053 2012-11-13 17:11:03Z hoffmann
84	! +nr, prr, qr
85	!
86	! 1036 2012-10-22 13:43:42Z raasch
87	! code put under GPL (PALM 3.9)
88	!
89	! 1007 2012-09-19 14:30:36Z franke
90	! Bugfix in calculation of ql_vp
91	!
92	! 978 2012-08-09 08:28:32Z fricke
93	! +z0h*
94	!
95	! Revision 1.1 2006/02/23 12:55:23 raasch
96	! Initial revision
97	!
98	!
99	! Description:
100	! ------------
101	!> Sum-up the values of 3d-arrays. The real averaging is later done in routine
102	!> average_3d_data.
103	!------------------------------------------------------------------------------!
104	SUBROUTINE sum_up_3d_data
105
106
107	USE arrays_3d, &
108	ONLY: dzw, e, nr, ol, p, pt, precipitation_rate, q, qc, ql, ql_c, &
109	ql_v, qr, qsws, rho, sa, shf, ssws, ts, u, us, v, vpt, w, z0, &
110	z0h, z0q
111
112	USE averaging, &
113	ONLY: e_av, lpt_av, lwp_av, nr_av, ol_av, p_av, pc_av, pr_av, prr_av, &
114	precipitation_rate_av, pt_av, q_av, qc_av, ql_av, ql_c_av, &
115	ql_v_av, ql_vp_av, qr_av, qsws_av, qv_av, rho_av, s_av, sa_av, &
116	shf_av, ssws_av, ts_av, u_av, us_av, v_av, vpt_av, w_av, z0_av, &
117	z0h_av, z0q_av
118
119	USE cloud_parameters, &
120	ONLY: l_d_cp, pt_d_t
121
122	USE control_parameters, &
123	ONLY: average_count_3d, cloud_physics, doav, doav_n, rho_surface
124
125	USE cpulog, &
126	ONLY: cpu_log, log_point
127
128	USE indices, &
129	ONLY: nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, nzt
130
131	USE kinds
132
133	USE land_surface_model_mod, &
134	ONLY: land_surface, lsm_3d_data_averaging
135
136	USE particle_attributes, &
137	ONLY: grid_particles, number_of_particles, particles, prt_count
138
139	USE radiation_model_mod, &
140	ONLY: radiation, radiation_3d_data_averaging
141
142
143	IMPLICIT NONE
144
145	INTEGER(iwp) :: i !<
146	INTEGER(iwp) :: ii !<
147	INTEGER(iwp) :: j !<
148	INTEGER(iwp) :: k !<
149	INTEGER(iwp) :: n !<
150	INTEGER(iwp) :: psi !<
151
152	REAL(wp) :: mean_r !<
153	REAL(wp) :: s_r2 !<
154	REAL(wp) :: s_r3 !<
155
156	CALL cpu_log (log_point(34),'sum_up_3d_data','start')
157
158	!
159	!-- Allocate and initialize the summation arrays if called for the very first
160	!-- time or the first time after average_3d_data has been called
161	!-- (some or all of the arrays may have been already allocated
162	!-- in read_3d_binary)
163	IF ( average_count_3d == 0 ) THEN
164
165	DO ii = 1, doav_n
166
167	SELECT CASE ( TRIM( doav(ii) ) )
168
169	CASE ( 'e' )
170	IF ( .NOT. ALLOCATED( e_av ) ) THEN
171	ALLOCATE( e_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
172	ENDIF
173	e_av = 0.0_wp
174
175	CASE ( 'lpt' )
176	IF ( .NOT. ALLOCATED( lpt_av ) ) THEN
177	ALLOCATE( lpt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
178	ENDIF
179	lpt_av = 0.0_wp
180
181	CASE ( 'lwp*' )
182	IF ( .NOT. ALLOCATED( lwp_av ) ) THEN
183	ALLOCATE( lwp_av(nysg:nyng,nxlg:nxrg) )
184	ENDIF
185	lwp_av = 0.0_wp
186
187	CASE ( 'nr' )
188	IF ( .NOT. ALLOCATED( nr_av ) ) THEN
189	ALLOCATE( nr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
190	ENDIF
191	nr_av = 0.0_wp
192
193	CASE ( 'ol*' )
194	IF ( .NOT. ALLOCATED( ol_av ) ) THEN
195	ALLOCATE( ol_av(nysg:nyng,nxlg:nxrg) )
196	ENDIF
197	ol_av = 0.0_wp
198
199	CASE ( 'p' )
200	IF ( .NOT. ALLOCATED( p_av ) ) THEN
201	ALLOCATE( p_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
202	ENDIF
203	p_av = 0.0_wp
204
205	CASE ( 'pc' )
206	IF ( .NOT. ALLOCATED( pc_av ) ) THEN
207	ALLOCATE( pc_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
208	ENDIF
209	pc_av = 0.0_wp
210
211	CASE ( 'pr' )
212	IF ( .NOT. ALLOCATED( pr_av ) ) THEN
213	ALLOCATE( pr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
214	ENDIF
215	pr_av = 0.0_wp
216
217	CASE ( 'prr' )
218	IF ( .NOT. ALLOCATED( prr_av ) ) THEN
219	ALLOCATE( prr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
220	ENDIF
221	prr_av = 0.0_wp
222
223	CASE ( 'prr*' )
224	IF ( .NOT. ALLOCATED( precipitation_rate_av ) ) THEN
225	ALLOCATE( precipitation_rate_av(nysg:nyng,nxlg:nxrg) )
226	ENDIF
227	precipitation_rate_av = 0.0_wp
228
229	CASE ( 'pt' )
230	IF ( .NOT. ALLOCATED( pt_av ) ) THEN
231	ALLOCATE( pt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
232	ENDIF
233	pt_av = 0.0_wp
234
235	CASE ( 'q' )
236	IF ( .NOT. ALLOCATED( q_av ) ) THEN
237	ALLOCATE( q_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
238	ENDIF
239	q_av = 0.0_wp
240
241	CASE ( 'qc' )
242	IF ( .NOT. ALLOCATED( qc_av ) ) THEN
243	ALLOCATE( qc_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
244	ENDIF
245	qc_av = 0.0_wp
246
247	CASE ( 'ql' )
248	IF ( .NOT. ALLOCATED( ql_av ) ) THEN
249	ALLOCATE( ql_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
250	ENDIF
251	ql_av = 0.0_wp
252
253	CASE ( 'ql_c' )
254	IF ( .NOT. ALLOCATED( ql_c_av ) ) THEN
255	ALLOCATE( ql_c_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
256	ENDIF
257	ql_c_av = 0.0_wp
258
259	CASE ( 'ql_v' )
260	IF ( .NOT. ALLOCATED( ql_v_av ) ) THEN
261	ALLOCATE( ql_v_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
262	ENDIF
263	ql_v_av = 0.0_wp
264
265	CASE ( 'ql_vp' )
266	IF ( .NOT. ALLOCATED( ql_vp_av ) ) THEN
267	ALLOCATE( ql_vp_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
268	ENDIF
269	ql_vp_av = 0.0_wp
270
271	CASE ( 'qr' )
272	IF ( .NOT. ALLOCATED( qr_av ) ) THEN
273	ALLOCATE( qr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
274	ENDIF
275	qr_av = 0.0_wp
276
277	CASE ( 'qsws*' )
278	IF ( .NOT. ALLOCATED( qsws_av ) ) THEN
279	ALLOCATE( qsws_av(nysg:nyng,nxlg:nxrg) )
280	ENDIF
281	qsws_av = 0.0_wp
282
283	CASE ( 'qv' )
284	IF ( .NOT. ALLOCATED( qv_av ) ) THEN
285	ALLOCATE( qv_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
286	ENDIF
287	qv_av = 0.0_wp
288
289	CASE ( 'rho' )
290	IF ( .NOT. ALLOCATED( rho_av ) ) THEN
291	ALLOCATE( rho_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
292	ENDIF
293	rho_av = 0.0_wp
294
295	CASE ( 's' )
296	IF ( .NOT. ALLOCATED( s_av ) ) THEN
297	ALLOCATE( s_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
298	ENDIF
299	s_av = 0.0_wp
300
301	CASE ( 'sa' )
302	IF ( .NOT. ALLOCATED( sa_av ) ) THEN
303	ALLOCATE( sa_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
304	ENDIF
305	sa_av = 0.0_wp
306
307	CASE ( 'shf*' )
308	IF ( .NOT. ALLOCATED( shf_av ) ) THEN
309	ALLOCATE( shf_av(nysg:nyng,nxlg:nxrg) )
310	ENDIF
311	shf_av = 0.0_wp
312
313	CASE ( 't*' )
314	IF ( .NOT. ALLOCATED( ts_av ) ) THEN
315	ALLOCATE( ts_av(nysg:nyng,nxlg:nxrg) )
316	ENDIF
317	ts_av = 0.0_wp
318
319	CASE ( 'u' )
320	IF ( .NOT. ALLOCATED( u_av ) ) THEN
321	ALLOCATE( u_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
322	ENDIF
323	u_av = 0.0_wp
324
325	CASE ( 'u*' )
326	IF ( .NOT. ALLOCATED( us_av ) ) THEN
327	ALLOCATE( us_av(nysg:nyng,nxlg:nxrg) )
328	ENDIF
329	us_av = 0.0_wp
330
331	CASE ( 'v' )
332	IF ( .NOT. ALLOCATED( v_av ) ) THEN
333	ALLOCATE( v_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
334	ENDIF
335	v_av = 0.0_wp
336
337	CASE ( 'vpt' )
338	IF ( .NOT. ALLOCATED( vpt_av ) ) THEN
339	ALLOCATE( vpt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
340	ENDIF
341	vpt_av = 0.0_wp
342
343	CASE ( 'w' )
344	IF ( .NOT. ALLOCATED( w_av ) ) THEN
345	ALLOCATE( w_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
346	ENDIF
347	w_av = 0.0_wp
348
349	CASE ( 'z0*' )
350	IF ( .NOT. ALLOCATED( z0_av ) ) THEN
351	ALLOCATE( z0_av(nysg:nyng,nxlg:nxrg) )
352	ENDIF
353	z0_av = 0.0_wp
354
355	CASE ( 'z0h*' )
356	IF ( .NOT. ALLOCATED( z0h_av ) ) THEN
357	ALLOCATE( z0h_av(nysg:nyng,nxlg:nxrg) )
358	ENDIF
359	z0h_av = 0.0_wp
360
361	CASE ( 'z0q*' )
362	IF ( .NOT. ALLOCATED( z0q_av ) ) THEN
363	ALLOCATE( z0q_av(nysg:nyng,nxlg:nxrg) )
364	ENDIF
365	z0q_av = 0.0_wp
366
367	CASE DEFAULT
368
369	!
370	!-- Land surface quantity
371	IF ( land_surface ) THEN
372	CALL lsm_3d_data_averaging( 'allocate', doav(ii) )
373	ENDIF
374
375	!
376	!-- Radiation quantity
377	IF ( radiation ) THEN
378	CALL radiation_3d_data_averaging( 'allocate', doav(ii) )
379	ENDIF
380
381	!
382	!-- User-defined quantity
383	CALL user_3d_data_averaging( 'allocate', doav(ii) )
384
385	END SELECT
386
387	ENDDO
388
389	ENDIF
390
391	!
392	!-- Loop of all variables to be averaged.
393	DO ii = 1, doav_n
394
395	!
396	!-- Store the array chosen on the temporary array.
397	SELECT CASE ( TRIM( doav(ii) ) )
398
399	CASE ( 'e' )
400	DO i = nxlg, nxrg
401	DO j = nysg, nyng
402	DO k = nzb, nzt+1
403	e_av(k,j,i) = e_av(k,j,i) + e(k,j,i)
404	ENDDO
405	ENDDO
406	ENDDO
407
408	CASE ( 'lpt' )
409	DO i = nxlg, nxrg
410	DO j = nysg, nyng
411	DO k = nzb, nzt+1
412	lpt_av(k,j,i) = lpt_av(k,j,i) + pt(k,j,i)
413	ENDDO
414	ENDDO
415	ENDDO
416
417	CASE ( 'lwp*' )
418	DO i = nxlg, nxrg
419	DO j = nysg, nyng
420	lwp_av(j,i) = lwp_av(j,i) + SUM( ql(nzb:nzt,j,i) &
421	* dzw(1:nzt+1) ) * rho_surface
422	ENDDO
423	ENDDO
424
425	CASE ( 'nr' )
426	DO i = nxlg, nxrg
427	DO j = nysg, nyng
428	DO k = nzb, nzt+1
429	nr_av(k,j,i) = nr_av(k,j,i) + nr(k,j,i)
430	ENDDO
431	ENDDO
432	ENDDO
433
434	CASE ( 'ol*' )
435	DO i = nxlg, nxrg
436	DO j = nysg, nyng
437	ol_av(j,i) = ol_av(j,i) + ol(j,i)
438	ENDDO
439	ENDDO
440
441	CASE ( 'p' )
442	DO i = nxlg, nxrg
443	DO j = nysg, nyng
444	DO k = nzb, nzt+1
445	p_av(k,j,i) = p_av(k,j,i) + p(k,j,i)
446	ENDDO
447	ENDDO
448	ENDDO
449
450	CASE ( 'pc' )
451	DO i = nxl, nxr
452	DO j = nys, nyn
453	DO k = nzb, nzt+1
454	pc_av(k,j,i) = pc_av(k,j,i) + prt_count(k,j,i)
455	ENDDO
456	ENDDO
457	ENDDO
458
459	CASE ( 'pr' )
460	DO i = nxl, nxr
461	DO j = nys, nyn
462	DO k = nzb, nzt+1
463	number_of_particles = prt_count(k,j,i)
464	IF ( number_of_particles <= 0 ) CYCLE
465	particles => grid_particles(k,j,i)%particles(1:number_of_particles)
466	s_r2 = 0.0_wp
467	s_r3 = 0.0_wp
468
469	DO n = 1, number_of_particles
470	IF ( particles(n)%particle_mask ) THEN
471	s_r2 = s_r2 + particles(n)%radius*2 &
472	particles(n)%weight_factor
473	s_r3 = s_r3 + particles(n)%radius*3 &
474	particles(n)%weight_factor
475	ENDIF
476	ENDDO
477
478	IF ( s_r2 > 0.0_wp ) THEN
479	mean_r = s_r3 / s_r2
480	ELSE
481	mean_r = 0.0_wp
482	ENDIF
483	pr_av(k,j,i) = pr_av(k,j,i) + mean_r
484	ENDDO
485	ENDDO
486	ENDDO
487
488
489	CASE ( 'pr*' )
490	DO i = nxlg, nxrg
491	DO j = nysg, nyng
492	precipitation_rate_av(j,i) = precipitation_rate_av(j,i) + &
493	precipitation_rate(j,i)
494	ENDDO
495	ENDDO
496
497	CASE ( 'pt' )
498	IF ( .NOT. cloud_physics ) THEN
499	DO i = nxlg, nxrg
500	DO j = nysg, nyng
501	DO k = nzb, nzt+1
502	pt_av(k,j,i) = pt_av(k,j,i) + pt(k,j,i)
503	ENDDO
504	ENDDO
505	ENDDO
506	ELSE
507	DO i = nxlg, nxrg
508	DO j = nysg, nyng
509	DO k = nzb, nzt+1
510	pt_av(k,j,i) = pt_av(k,j,i) + pt(k,j,i) + l_d_cp * &
511	pt_d_t(k) * ql(k,j,i)
512	ENDDO
513	ENDDO
514	ENDDO
515	ENDIF
516
517	CASE ( 'q' )
518	DO i = nxlg, nxrg
519	DO j = nysg, nyng
520	DO k = nzb, nzt+1
521	q_av(k,j,i) = q_av(k,j,i) + q(k,j,i)
522	ENDDO
523	ENDDO
524	ENDDO
525
526	CASE ( 'qc' )
527	DO i = nxlg, nxrg
528	DO j = nysg, nyng
529	DO k = nzb, nzt+1
530	qc_av(k,j,i) = qc_av(k,j,i) + qc(k,j,i)
531	ENDDO
532	ENDDO
533	ENDDO
534
535	CASE ( 'ql' )
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
544	CASE ( 'ql_c' )
545	DO i = nxlg, nxrg
546	DO j = nysg, nyng
547	DO k = nzb, nzt+1
548	ql_c_av(k,j,i) = ql_c_av(k,j,i) + ql_c(k,j,i)
549	ENDDO
550	ENDDO
551	ENDDO
552
553	CASE ( 'ql_v' )
554	DO i = nxlg, nxrg
555	DO j = nysg, nyng
556	DO k = nzb, nzt+1
557	ql_v_av(k,j,i) = ql_v_av(k,j,i) + ql_v(k,j,i)
558	ENDDO
559	ENDDO
560	ENDDO
561
562	CASE ( 'ql_vp' )
563	DO i = nxl, nxr
564	DO j = nys, nyn
565	DO k = nzb, nzt+1
566	number_of_particles = prt_count(k,j,i)
567	IF ( number_of_particles <= 0 ) CYCLE
568	particles => grid_particles(k,j,i)%particles(1:number_of_particles)
569	DO n = 1, number_of_particles
570	IF ( particles(n)%particle_mask ) THEN
571	ql_vp_av(k,j,i) = ql_vp_av(k,j,i) + &
572	particles(n)%weight_factor / &
573	number_of_particles
574	ENDIF
575	ENDDO
576	ENDDO
577	ENDDO
578	ENDDO
579
580	CASE ( 'qr' )
581	DO i = nxlg, nxrg
582	DO j = nysg, nyng
583	DO k = nzb, nzt+1
584	qr_av(k,j,i) = qr_av(k,j,i) + qr(k,j,i)
585	ENDDO
586	ENDDO
587	ENDDO
588
589	CASE ( 'qsws*' )
590	DO i = nxlg, nxrg
591	DO j = nysg, nyng
592	qsws_av(j,i) = qsws_av(j,i) + qsws(j,i)
593	ENDDO
594	ENDDO
595
596	CASE ( 'qv' )
597	DO i = nxlg, nxrg
598	DO j = nysg, nyng
599	DO k = nzb, nzt+1
600	qv_av(k,j,i) = qv_av(k,j,i) + q(k,j,i) - ql(k,j,i)
601	ENDDO
602	ENDDO
603	ENDDO
604
605	CASE ( 'rho' )
606	DO i = nxlg, nxrg
607	DO j = nysg, nyng
608	DO k = nzb, nzt+1
609	rho_av(k,j,i) = rho_av(k,j,i) + rho(k,j,i)
610	ENDDO
611	ENDDO
612	ENDDO
613
614	CASE ( 's' )
615	DO i = nxlg, nxrg
616	DO j = nysg, nyng
617	DO k = nzb, nzt+1
618	s_av(k,j,i) = s_av(k,j,i) + q(k,j,i)
619	ENDDO
620	ENDDO
621	ENDDO
622
623	CASE ( 'sa' )
624	DO i = nxlg, nxrg
625	DO j = nysg, nyng
626	DO k = nzb, nzt+1
627	sa_av(k,j,i) = sa_av(k,j,i) + sa(k,j,i)
628	ENDDO
629	ENDDO
630	ENDDO
631
632	CASE ( 'shf*' )
633	DO i = nxlg, nxrg
634	DO j = nysg, nyng
635	shf_av(j,i) = shf_av(j,i) + shf(j,i)
636	ENDDO
637	ENDDO
638
639	CASE ( 'ssws*' )
640	DO i = nxlg, nxrg
641	DO j = nysg, nyng
642	ssws_av(j,i) = ssws_av(j,i) + ssws(j,i)
643	ENDDO
644	ENDDO
645
646	CASE ( 't*' )
647	DO i = nxlg, nxrg
648	DO j = nysg, nyng
649	ts_av(j,i) = ts_av(j,i) + ts(j,i)
650	ENDDO
651	ENDDO
652
653	CASE ( 'u' )
654	DO i = nxlg, nxrg
655	DO j = nysg, nyng
656	DO k = nzb, nzt+1
657	u_av(k,j,i) = u_av(k,j,i) + u(k,j,i)
658	ENDDO
659	ENDDO
660	ENDDO
661
662	CASE ( 'u*' )
663	DO i = nxlg, nxrg
664	DO j = nysg, nyng
665	us_av(j,i) = us_av(j,i) + us(j,i)
666	ENDDO
667	ENDDO
668
669	CASE ( 'v' )
670	DO i = nxlg, nxrg
671	DO j = nysg, nyng
672	DO k = nzb, nzt+1
673	v_av(k,j,i) = v_av(k,j,i) + v(k,j,i)
674	ENDDO
675	ENDDO
676	ENDDO
677
678	CASE ( 'vpt' )
679	DO i = nxlg, nxrg
680	DO j = nysg, nyng
681	DO k = nzb, nzt+1
682	vpt_av(k,j,i) = vpt_av(k,j,i) + vpt(k,j,i)
683	ENDDO
684	ENDDO
685	ENDDO
686
687	CASE ( 'w' )
688	DO i = nxlg, nxrg
689	DO j = nysg, nyng
690	DO k = nzb, nzt+1
691	w_av(k,j,i) = w_av(k,j,i) + w(k,j,i)
692	ENDDO
693	ENDDO
694	ENDDO
695
696	CASE ( 'z0*' )
697	DO i = nxlg, nxrg
698	DO j = nysg, nyng
699	z0_av(j,i) = z0_av(j,i) + z0(j,i)
700	ENDDO
701	ENDDO
702
703	CASE ( 'z0h*' )
704	DO i = nxlg, nxrg
705	DO j = nysg, nyng
706	z0h_av(j,i) = z0h_av(j,i) + z0h(j,i)
707	ENDDO
708	ENDDO
709
710	CASE ( 'z0q*' )
711	DO i = nxlg, nxrg
712	DO j = nysg, nyng
713	z0q_av(j,i) = z0q_av(j,i) + z0q(j,i)
714	ENDDO
715	ENDDO
716
717	CASE DEFAULT
718	!
719	!-- Land surface quantity
720	IF ( land_surface ) THEN
721	CALL lsm_3d_data_averaging( 'sum', doav(ii) )
722	ENDIF
723
724	!
725	!-- Radiation quantity
726	IF ( radiation ) THEN
727	CALL radiation_3d_data_averaging( 'sum', doav(ii) )
728	ENDIF
729
730	!
731	!-- User-defined quantity
732	CALL user_3d_data_averaging( 'sum', doav(ii) )
733
734	END SELECT
735
736	ENDDO
737
738	CALL cpu_log( log_point(34), 'sum_up_3d_data', 'stop' )
739
740
741	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 |