Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |