Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |