Home

Context Navigation

source: palm/trunk/SOURCE/average_3d_data.f90 @ 1320

Last change on this file since 1320 was 1320, checked in by raasch, 10 years ago
ONLY-attribute added to USE-statements, kind-parameters added to all INTEGER and REAL declaration statements, kinds are defined in new module kinds, old module precision_kind is removed, revision history before 2012 removed, comment fields (!:) to be used for variable explanations added to all variable declaration statements
Property svn:keywords set to `Id`
File size: 10.5 KB

Rev	Line
[1]	1	SUBROUTINE average_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	! -----------------
[1320]	22	! ONLY-attribute added to USE-statements,
	23	! kind-parameters added to all INTEGER and REAL declaration statements,
	24	! kinds are defined in new module kinds,
	25	! revision history before 2012 removed,
	26	! comment fields (!:) to be used for variable explanations added to
	27	! all variable declaration statements
[772]	28	!
[1]	29	! Former revisions:
	30	! -----------------
[3]	31	! $Id: average_3d_data.f90 1320 2014-03-20 08:40:49Z raasch $
[1054]	32	!
[1319]	33	! 1318 2014-03-17 13:35:16Z raasch
	34	! barrier argument removed from cpu_log,
	35	! module interfaces removed
	36	!
[1116]	37	! 1115 2013-03-26 18:16:16Z hoffmann
	38	! +qc
	39	!
[1054]	40	! 1053 2012-11-13 17:11:03Z hoffmann
	41	! averaging of nr, qr added
	42	!
[1037]	43	! 1036 2012-10-22 13:43:42Z raasch
	44	! code put under GPL (PALM 3.9)
	45	!
[979]	46	! 978 2012-08-09 08:28:32Z fricke
	47	! +z0h_av
[77]	48	!
[1]	49	! Revision 1.1 2006/02/23 09:48:58 raasch
	50	! Initial revision
	51	!
	52	!
	53	! Description:
	54	! ------------
	55	! Time-averaging of 3d-data-arrays.
	56	!------------------------------------------------------------------------------!
	57
	58	USE averaging
	59
[1320]	60	USE control_parameters, &
	61	ONLY: average_count_3d, doav, doav_n
	62
	63	USE cpulog, &
	64	ONLY: cpu_log, log_point
	65
	66	USE indices, &
	67	ONLY: nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, nzt
	68
	69	USE kinds
	70
	71
[1]	72	IMPLICIT NONE
	73
[1320]	74	INTEGER(iwp) :: i !:
	75	INTEGER(iwp) :: ii !:
	76	INTEGER(iwp) :: j !:
	77	INTEGER(iwp) :: k !:
[1]	78
	79
	80	CALL cpu_log (log_point(35),'average_3d_data','start')
	81
	82	!
	83	!-- Check, if averaging is necessary
	84	IF ( average_count_3d <= 1 ) RETURN
	85
	86	!
	87	!-- Loop of all variables to be averaged.
	88	DO ii = 1, doav_n
	89
	90	!
	91	!-- Store the array chosen on the temporary array.
	92	SELECT CASE ( TRIM( doav(ii) ) )
	93
	94	CASE ( 'e' )
[667]	95	DO i = nxlg, nxrg
	96	DO j = nysg, nyng
[1]	97	DO k = nzb, nzt+1
	98	e_av(k,j,i) = e_av(k,j,i) / REAL( average_count_3d )
	99	ENDDO
	100	ENDDO
	101	ENDDO
	102
[354]	103	CASE ( 'qsws*' )
[667]	104	DO i = nxlg, nxrg
	105	DO j = nysg, nyng
[354]	106	qsws_av(j,i) = qsws_av(j,i) / REAL( average_count_3d )
	107	ENDDO
	108	ENDDO
	109
[771]	110	CASE ( 'lpt' )
	111	DO i = nxlg, nxrg
	112	DO j = nysg, nyng
	113	DO k = nzb, nzt+1
	114	lpt_av(k,j,i) = lpt_av(k,j,i) / REAL( average_count_3d )
	115	ENDDO
	116	ENDDO
	117	ENDDO
	118
[1]	119	CASE ( 'lwp*' )
[667]	120	DO i = nxlg, nxrg
	121	DO j = nysg, nyng
[1]	122	lwp_av(j,i) = lwp_av(j,i) / REAL( average_count_3d )
	123	ENDDO
	124	ENDDO
	125
[1053]	126	CASE ( 'nr' )
	127	DO i = nxlg, nxrg
	128	DO j = nysg, nyng
	129	DO k = nzb, nzt+1
	130	nr_av(k,j,i) = nr_av(k,j,i) / REAL( average_count_3d )
	131	ENDDO
	132	ENDDO
	133	ENDDO
	134
[1]	135	CASE ( 'p' )
[667]	136	DO i = nxlg, nxrg
	137	DO j = nysg, nyng
[1]	138	DO k = nzb, nzt+1
	139	p_av(k,j,i) = p_av(k,j,i) / REAL( average_count_3d )
	140	ENDDO
	141	ENDDO
	142	ENDDO
	143
	144	CASE ( 'pc' )
	145	DO i = nxl, nxr
	146	DO j = nys, nyn
	147	DO k = nzb, nzt+1
	148	pc_av(k,j,i) = pc_av(k,j,i) / REAL( average_count_3d )
	149	ENDDO
	150	ENDDO
	151	ENDDO
	152
	153	CASE ( 'pr' )
	154	DO i = nxl, nxr
	155	DO j = nys, nyn
	156	DO k = nzb, nzt+1
	157	pr_av(k,j,i) = pr_av(k,j,i) / REAL( average_count_3d )
	158	ENDDO
	159	ENDDO
	160	ENDDO
	161
[72]	162	CASE ( 'prr*' )
[667]	163	DO i = nxlg, nxrg
	164	DO j = nysg, nyng
[1320]	165	precipitation_rate_av(j,i) = precipitation_rate_av(j,i) / &
[72]	166	REAL( average_count_3d )
	167	ENDDO
	168	ENDDO
	169
[1]	170	CASE ( 'pt' )
[667]	171	DO i = nxlg, nxrg
	172	DO j = nysg, nyng
[1]	173	DO k = nzb, nzt+1
	174	pt_av(k,j,i) = pt_av(k,j,i) / REAL( average_count_3d )
	175	ENDDO
	176	ENDDO
	177	ENDDO
	178
	179	CASE ( 'q' )
[667]	180	DO i = nxlg, nxrg
	181	DO j = nysg, nyng
[1]	182	DO k = nzb, nzt+1
	183	q_av(k,j,i) = q_av(k,j,i) / REAL( average_count_3d )
	184	ENDDO
	185	ENDDO
	186	ENDDO
[367]	187
[1115]	188	CASE ( 'qc' )
	189	DO i = nxlg, nxrg
	190	DO j = nysg, nyng
	191	DO k = nzb, nzt+1
	192	qc_av(k,j,i) = qc_av(k,j,i) / REAL( average_count_3d )
	193	ENDDO
	194	ENDDO
	195	ENDDO
	196
[1]	197	CASE ( 'ql' )
[667]	198	DO i = nxlg, nxrg
	199	DO j = nysg, nyng
[1]	200	DO k = nzb, nzt+1
	201	ql_av(k,j,i) = ql_av(k,j,i) / REAL( average_count_3d )
	202	ENDDO
	203	ENDDO
	204	ENDDO
	205
	206	CASE ( 'ql_c' )
[667]	207	DO i = nxlg, nxrg
	208	DO j = nysg, nyng
[1]	209	DO k = nzb, nzt+1
	210	ql_c_av(k,j,i) = ql_c_av(k,j,i) / REAL( average_count_3d )
	211	ENDDO
	212	ENDDO
	213	ENDDO
	214
	215	CASE ( 'ql_v' )
[667]	216	DO i = nxlg, nxrg
	217	DO j = nysg, nyng
[1]	218	DO k = nzb, nzt+1
	219	ql_v_av(k,j,i) = ql_v_av(k,j,i) / REAL( average_count_3d )
	220	ENDDO
	221	ENDDO
	222	ENDDO
	223
	224	CASE ( 'ql_vp' )
[667]	225	DO i = nxlg, nxrg
	226	DO j = nysg, nyng
[1]	227	DO k = nzb, nzt+1
[1320]	228	ql_vp_av(k,j,i) = ql_vp_av(k,j,i) / &
[1]	229	REAL( average_count_3d )
	230	ENDDO
	231	ENDDO
	232	ENDDO
	233
[1053]	234	CASE ( 'qr' )
	235	DO i = nxlg, nxrg
	236	DO j = nysg, nyng
	237	DO k = nzb, nzt+1
	238	qr_av(k,j,i) = qr_av(k,j,i) / REAL( average_count_3d )
	239	ENDDO
	240	ENDDO
	241	ENDDO
	242
[1]	243	CASE ( 'qv' )
[667]	244	DO i = nxlg, nxrg
	245	DO j = nysg, nyng
[1]	246	DO k = nzb, nzt+1
	247	qv_av(k,j,i) = qv_av(k,j,i) / REAL( average_count_3d )
	248	ENDDO
	249	ENDDO
	250	ENDDO
[367]	251
[96]	252	CASE ( 'rho' )
[667]	253	DO i = nxlg, nxrg
	254	DO j = nysg, nyng
[96]	255	DO k = nzb, nzt+1
	256	rho_av(k,j,i) = rho_av(k,j,i) / REAL( average_count_3d )
	257	ENDDO
	258	ENDDO
	259	ENDDO
[367]	260
[1]	261	CASE ( 's' )
[667]	262	DO i = nxlg, nxrg
	263	DO j = nysg, nyng
[1]	264	DO k = nzb, nzt+1
	265	s_av(k,j,i) = s_av(k,j,i) / REAL( average_count_3d )
	266	ENDDO
	267	ENDDO
	268	ENDDO
[367]	269
[96]	270	CASE ( 'sa' )
[667]	271	DO i = nxlg, nxrg
	272	DO j = nysg, nyng
[96]	273	DO k = nzb, nzt+1
	274	sa_av(k,j,i) = sa_av(k,j,i) / REAL( average_count_3d )
	275	ENDDO
	276	ENDDO
	277	ENDDO
[367]	278
[354]	279	CASE ( 'shf*' )
[667]	280	DO i = nxlg, nxrg
	281	DO j = nysg, nyng
[354]	282	shf_av(j,i) = shf_av(j,i) / REAL( average_count_3d )
	283	ENDDO
	284	ENDDO
[367]	285
[1]	286	CASE ( 't*' )
[667]	287	DO i = nxlg, nxrg
	288	DO j = nysg, nyng
[1]	289	ts_av(j,i) = ts_av(j,i) / REAL( average_count_3d )
	290	ENDDO
	291	ENDDO
	292
	293	CASE ( 'u' )
[667]	294	DO i = nxlg, nxrg
	295	DO j = nysg, nyng
[1]	296	DO k = nzb, nzt+1
	297	u_av(k,j,i) = u_av(k,j,i) / REAL( average_count_3d )
	298	ENDDO
	299	ENDDO
	300	ENDDO
	301
	302	CASE ( 'u*' )
[667]	303	DO i = nxlg, nxrg
	304	DO j = nysg, nyng
[1]	305	us_av(j,i) = us_av(j,i) / REAL( average_count_3d )
	306	ENDDO
	307	ENDDO
	308
	309	CASE ( 'v' )
[667]	310	DO i = nxlg, nxrg
	311	DO j = nysg, nyng
[1]	312	DO k = nzb, nzt+1
	313	v_av(k,j,i) = v_av(k,j,i) / REAL( average_count_3d )
	314	ENDDO
	315	ENDDO
	316	ENDDO
	317
	318	CASE ( 'vpt' )
[667]	319	DO i = nxlg, nxrg
	320	DO j = nysg, nyng
[1]	321	DO k = nzb, nzt+1
	322	vpt_av(k,j,i) = vpt_av(k,j,i) / REAL( average_count_3d )
	323	ENDDO
	324	ENDDO
	325	ENDDO
	326
	327	CASE ( 'w' )
[667]	328	DO i = nxlg, nxrg
	329	DO j = nysg, nyng
[1]	330	DO k = nzb, nzt+1
	331	w_av(k,j,i) = w_av(k,j,i) / REAL( average_count_3d )
	332	ENDDO
	333	ENDDO
	334	ENDDO
	335
[72]	336	CASE ( 'z0*' )
[667]	337	DO i = nxlg, nxrg
	338	DO j = nysg, nyng
[72]	339	z0_av(j,i) = z0_av(j,i) / REAL( average_count_3d )
	340	ENDDO
	341	ENDDO
	342
[978]	343	CASE ( 'z0h*' )
	344	DO i = nxlg, nxrg
	345	DO j = nysg, nyng
	346	z0h_av(j,i) = z0h_av(j,i) / REAL( average_count_3d )
	347	ENDDO
	348	ENDDO
	349
[1]	350	CASE DEFAULT
	351	!
	352	!-- User-defined quantity
	353	CALL user_3d_data_averaging( 'average', doav(ii) )
	354
	355	END SELECT
	356
	357	ENDDO
	358
	359	!
	360	!-- Reset the counter
	361	average_count_3d = 0.0
	362
[1318]	363	CALL cpu_log( log_point(35), 'average_3d_data', 'stop' )
[1]	364
	365
	366	END SUBROUTINE average_3d_data

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |