Home

Context Navigation

source: palm/trunk/SOURCE/diffusion_w.f90 @ 1682

Last change on this file since 1682 was 1682, checked in by knoop, 9 years ago
Code annotations made doxygen readable
Property svn:keywords set to `Id`
File size: 22.7 KB

Rev	Line
[1682]	1	!> @file diffusion_w.f90
[1036]	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	!
[1310]	16	! Copyright 1997-2014 Leibniz Universitaet Hannover
[1036]	17	!--------------------------------------------------------------------------------!
	18	!
[484]	19	! Current revisions:
[1]	20	! -----------------
[1682]	21	! Code annotations made doxygen readable
[1341]	22	!
[1321]	23	! Former revisions:
	24	! -----------------
	25	! $Id: diffusion_w.f90 1682 2015-10-07 23:56:08Z knoop $
	26	!
[1375]	27	! 1374 2014-04-25 12:55:07Z raasch
	28	! vsws + vswst removed from acc-present-list
	29	!
[1354]	30	! 1353 2014-04-08 15:21:23Z heinze
	31	! REAL constants provided with KIND-attribute
	32	!
[1341]	33	! 1340 2014-03-25 19:45:13Z kanani
	34	! REAL constants defined as wp-kind
	35	!
[1323]	36	! 1322 2014-03-20 16:38:49Z raasch
	37	! REAL constants defined as wp-kind
	38	!
[1321]	39	! 1320 2014-03-20 08:40:49Z raasch
[1320]	40	! ONLY-attribute added to USE-statements,
	41	! kind-parameters added to all INTEGER and REAL declaration statements,
	42	! kinds are defined in new module kinds,
	43	! revision history before 2012 removed,
	44	! comment fields (!:) to be used for variable explanations added to
	45	! all variable declaration statements
[1]	46	!
[1258]	47	! 1257 2013-11-08 15:18:40Z raasch
	48	! openacc loop and loop vector clauses removed, declare create moved after
	49	! the FORTRAN declaration statement
	50	!
[1132]	51	! 1128 2013-04-12 06:19:32Z raasch
	52	! loop index bounds in accelerator version replaced by i_left, i_right, j_south,
	53	! j_north
	54	!
[1037]	55	! 1036 2012-10-22 13:43:42Z raasch
	56	! code put under GPL (PALM 3.9)
	57	!
[1017]	58	! 1015 2012-09-27 09:23:24Z raasch
	59	! accelerator version (*_acc) added
	60	!
[1002]	61	! 1001 2012-09-13 14:08:46Z raasch
	62	! arrays comunicated by module instead of parameter list
	63	!
[979]	64	! 978 2012-08-09 08:28:32Z fricke
	65	! outflow damping layer removed
	66	! kmxm_x/_z and kmxp_x/_z change to kmxm and kmxp
	67	! kmym_y/_z and kmyp_y/_z change to kmym and kmyp
	68	!
[1]	69	! Revision 1.1 1997/09/12 06:24:11 raasch
	70	! Initial revision
	71	!
	72	!
	73	! Description:
	74	! ------------
[1682]	75	!> Diffusion term of the w-component
[1]	76	!------------------------------------------------------------------------------!
[1682]	77	MODULE diffusion_w_mod
	78
[1]	79
[1320]	80	USE wall_fluxes_mod, &
	81	ONLY : wall_fluxes, wall_fluxes_acc
[56]	82
[1]	83	PRIVATE
[1015]	84	PUBLIC diffusion_w, diffusion_w_acc
[1]	85
	86	INTERFACE diffusion_w
	87	MODULE PROCEDURE diffusion_w
	88	MODULE PROCEDURE diffusion_w_ij
	89	END INTERFACE diffusion_w
	90
[1015]	91	INTERFACE diffusion_w_acc
	92	MODULE PROCEDURE diffusion_w_acc
	93	END INTERFACE diffusion_w_acc
	94
[1]	95	CONTAINS
	96
	97
	98	!------------------------------------------------------------------------------!
[1682]	99	! Description:
	100	! ------------
	101	!> Call for all grid points
[1]	102	!------------------------------------------------------------------------------!
[1001]	103	SUBROUTINE diffusion_w
[1]	104
[1320]	105	USE arrays_3d, &
	106	ONLY : ddzu, ddzw, km, tend, u, v, w
	107
	108	USE control_parameters, &
	109	ONLY : topography
	110
	111	USE grid_variables, &
	112	ONLY : ddx, ddy, fwxm, fwxp, fwym, fwyp, wall_w_x, wall_w_y
	113
	114	USE indices, &
	115	ONLY : nxl, nxr, nyn, nys, nzb, nzb_w_inner, nzb_w_outer, nzt
	116
	117	USE kinds
[1]	118
	119	IMPLICIT NONE
	120
[1682]	121	INTEGER(iwp) :: i !<
	122	INTEGER(iwp) :: j !<
	123	INTEGER(iwp) :: k !<
[1320]	124
[1682]	125	REAL(wp) :: kmxm !<
	126	REAL(wp) :: kmxp !<
	127	REAL(wp) :: kmym !<
	128	REAL(wp) :: kmyp !<
[1001]	129
[1682]	130	REAL(wp), DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr) :: wsus !<
	131	REAL(wp), DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr) :: wsvs !<
[1]	132
	133
[56]	134	!
	135	!-- First calculate horizontal momentum flux w'u' and/or w'v' at vertical
	136	!-- walls, if neccessary
	137	IF ( topography /= 'flat' ) THEN
[1320]	138	CALL wall_fluxes( wsus, 0.0_wp, 0.0_wp, 0.0_wp, 1.0_wp, nzb_w_inner, &
[56]	139	nzb_w_outer, wall_w_x )
[1320]	140	CALL wall_fluxes( wsvs, 0.0_wp, 0.0_wp, 1.0_wp, 0.0_wp, nzb_w_inner, &
[56]	141	nzb_w_outer, wall_w_y )
	142	ENDIF
	143
[1]	144	DO i = nxl, nxr
	145	DO j = nys, nyn
	146	DO k = nzb_w_outer(j,i)+1, nzt-1
	147	!
	148	!-- Interpolate eddy diffusivities on staggered gridpoints
[1322]	149	kmxp = 0.25_wp * &
[978]	150	( km(k,j,i)+km(k,j,i+1)+km(k+1,j,i)+km(k+1,j,i+1) )
[1322]	151	kmxm = 0.25_wp * &
[978]	152	( km(k,j,i)+km(k,j,i-1)+km(k+1,j,i)+km(k+1,j,i-1) )
[1322]	153	kmyp = 0.25_wp * &
[978]	154	( km(k,j,i)+km(k+1,j,i)+km(k,j+1,i)+km(k+1,j+1,i) )
[1322]	155	kmym = 0.25_wp * &
[978]	156	( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
[1]	157
	158	tend(k,j,i) = tend(k,j,i) &
[978]	159	& + ( kmxp * ( w(k,j,i+1) - w(k,j,i) ) * ddx &
	160	& + kmxp * ( u(k+1,j,i+1) - u(k,j,i+1) ) * ddzu(k+1) &
	161	& - kmxm * ( w(k,j,i) - w(k,j,i-1) ) * ddx &
	162	& - kmxm * ( u(k+1,j,i) - u(k,j,i) ) * ddzu(k+1) &
[1]	163	& ) * ddx &
[978]	164	& + ( kmyp * ( w(k,j+1,i) - w(k,j,i) ) * ddy &
	165	& + kmyp * ( v(k+1,j+1,i) - v(k,j+1,i) ) * ddzu(k+1) &
	166	& - kmym * ( w(k,j,i) - w(k,j-1,i) ) * ddy &
	167	& - kmym * ( v(k+1,j,i) - v(k,j,i) ) * ddzu(k+1) &
[1]	168	& ) * ddy &
[1322]	169	& + 2.0_wp * ( &
[1]	170	& km(k+1,j,i) * ( w(k+1,j,i) - w(k,j,i) ) * ddzw(k+1) &
	171	& - km(k,j,i) * ( w(k,j,i) - w(k-1,j,i) ) * ddzw(k) &
[1340]	172	& ) * ddzu(k+1)
[1]	173	ENDDO
	174
	175	!
	176	!-- Wall functions at all vertical walls, where necessary
[1340]	177	IF ( wall_w_x(j,i) /= 0.0_wp .OR. wall_w_y(j,i) /= 0.0_wp ) THEN
[51]	178
[1]	179	DO k = nzb_w_inner(j,i)+1, nzb_w_outer(j,i)
	180	!
	181	!-- Interpolate eddy diffusivities on staggered gridpoints
[1322]	182	kmxp = 0.25_wp * &
[978]	183	( km(k,j,i)+km(k,j,i+1)+km(k+1,j,i)+km(k+1,j,i+1) )
[1322]	184	kmxm = 0.25_wp * &
[978]	185	( km(k,j,i)+km(k,j,i-1)+km(k+1,j,i)+km(k+1,j,i-1) )
[1322]	186	kmyp = 0.25_wp * &
[978]	187	( km(k,j,i)+km(k+1,j,i)+km(k,j+1,i)+km(k+1,j+1,i) )
[1322]	188	kmym = 0.25_wp * &
[978]	189	( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
[1]	190
	191	tend(k,j,i) = tend(k,j,i) &
	192	+ ( fwxp(j,i) * ( &
[978]	193	kmxp * ( w(k,j,i+1) - w(k,j,i) ) * ddx &
	194	+ kmxp * ( u(k+1,j,i+1) - u(k,j,i+1) ) * ddzu(k+1) &
[1]	195	) &
	196	- fwxm(j,i) * ( &
[978]	197	kmxm * ( w(k,j,i) - w(k,j,i-1) ) * ddx &
	198	+ kmxm * ( u(k+1,j,i) - u(k,j,i) ) * ddzu(k+1) &
[1]	199	) &
[56]	200	+ wall_w_x(j,i) * wsus(k,j,i) &
[1]	201	) * ddx &
	202	+ ( fwyp(j,i) * ( &
[978]	203	kmyp * ( w(k,j+1,i) - w(k,j,i) ) * ddy &
	204	+ kmyp * ( v(k+1,j+1,i) - v(k,j+1,i) ) * ddzu(k+1) &
[1]	205	) &
	206	- fwym(j,i) * ( &
[978]	207	kmym * ( w(k,j,i) - w(k,j-1,i) ) * ddy &
	208	+ kmym * ( v(k+1,j,i) - v(k,j,i) ) * ddzu(k+1) &
[1]	209	) &
[56]	210	+ wall_w_y(j,i) * wsvs(k,j,i) &
[1]	211	) * ddy &
[1340]	212	+ 2.0_wp * ( &
[1]	213	km(k+1,j,i) * ( w(k+1,j,i) - w(k,j,i) ) * ddzw(k+1) &
	214	- km(k,j,i) * ( w(k,j,i) - w(k-1,j,i) ) * ddzw(k) &
[1340]	215	) * ddzu(k+1)
[1]	216	ENDDO
	217	ENDIF
	218
	219	ENDDO
	220	ENDDO
	221
	222	END SUBROUTINE diffusion_w
	223
	224
	225	!------------------------------------------------------------------------------!
[1682]	226	! Description:
	227	! ------------
	228	!> Call for all grid points - accelerator version
[1015]	229	!------------------------------------------------------------------------------!
	230	SUBROUTINE diffusion_w_acc
	231
[1320]	232	USE arrays_3d, &
	233	ONLY : ddzu, ddzw, km, tend, u, v, w
	234
	235	USE control_parameters, &
	236	ONLY : topography
	237
	238	USE grid_variables, &
	239	ONLY : ddx, ddy, fwxm, fwxp, fwym, fwyp, wall_w_x, wall_w_y
	240
	241	USE indices, &
	242	ONLY : i_left, i_right, j_north, j_south, nxl, nxr, nyn, nys, nzb, &
	243	nzb_w_inner, nzb_w_outer, nzt
	244
	245	USE kinds
[1015]	246
	247	IMPLICIT NONE
	248
[1682]	249	INTEGER(iwp) :: i !<
	250	INTEGER(iwp) :: j !<
	251	INTEGER(iwp) :: k !<
[1320]	252
[1682]	253	REAL(wp) :: kmxm !<
	254	REAL(wp) :: kmxp !<
	255	REAL(wp) :: kmym !<
	256	REAL(wp) :: kmyp !<
[1015]	257
[1682]	258	REAL(wp), DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr) :: wsus !<
	259	REAL(wp), DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr) :: wsvs !<
[1015]	260	!$acc declare create ( wsus, wsvs )
	261
	262	!
	263	!-- First calculate horizontal momentum flux w'u' and/or w'v' at vertical
	264	!-- walls, if neccessary
	265	IF ( topography /= 'flat' ) THEN
[1320]	266	CALL wall_fluxes_acc( wsus, 0.0_wp, 0.0_wp, 0.0_wp, 1.0_wp, &
	267	nzb_w_inner, nzb_w_outer, wall_w_x )
	268	CALL wall_fluxes_acc( wsvs, 0.0_wp, 0.0_wp, 0.0_wp, 1.0_wp, &
	269	nzb_w_inner, nzb_w_outer, wall_w_y )
[1015]	270	ENDIF
	271
[1374]	272	!$acc kernels present ( u, v, w, km, tend ) &
[1015]	273	!$acc present ( ddzu, ddzw, fwxm, fwxp, fwym, fwyp, wall_w_x, wall_w_y ) &
	274	!$acc present ( nzb_w_inner, nzb_w_outer )
[1128]	275	DO i = i_left, i_right
	276	DO j = j_south, j_north
[1015]	277	DO k = 1, nzt
	278	IF ( k > nzb_w_outer(j,i) ) THEN
	279	!
	280	!-- Interpolate eddy diffusivities on staggered gridpoints
[1322]	281	kmxp = 0.25_wp * &
[1015]	282	( km(k,j,i)+km(k,j,i+1)+km(k+1,j,i)+km(k+1,j,i+1) )
[1322]	283	kmxm = 0.25_wp * &
[1015]	284	( km(k,j,i)+km(k,j,i-1)+km(k+1,j,i)+km(k+1,j,i-1) )
[1322]	285	kmyp = 0.25_wp * &
[1015]	286	( km(k,j,i)+km(k+1,j,i)+km(k,j+1,i)+km(k+1,j+1,i) )
[1322]	287	kmym = 0.25_wp * &
[1015]	288	( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
	289
	290	tend(k,j,i) = tend(k,j,i) &
	291	& + ( kmxp * ( w(k,j,i+1) - w(k,j,i) ) * ddx &
	292	& + kmxp * ( u(k+1,j,i+1) - u(k,j,i+1) ) * ddzu(k+1) &
	293	& - kmxm * ( w(k,j,i) - w(k,j,i-1) ) * ddx &
	294	& - kmxm * ( u(k+1,j,i) - u(k,j,i) ) * ddzu(k+1) &
	295	& ) * ddx &
	296	& + ( kmyp * ( w(k,j+1,i) - w(k,j,i) ) * ddy &
	297	& + kmyp * ( v(k+1,j+1,i) - v(k,j+1,i) ) * ddzu(k+1) &
	298	& - kmym * ( w(k,j,i) - w(k,j-1,i) ) * ddy &
	299	& - kmym * ( v(k+1,j,i) - v(k,j,i) ) * ddzu(k+1) &
	300	& ) * ddy &
[1322]	301	& + 2.0_wp * ( &
[1015]	302	& km(k+1,j,i) * ( w(k+1,j,i) - w(k,j,i) ) * ddzw(k+1) &
	303	& - km(k,j,i) * ( w(k,j,i) - w(k-1,j,i) ) * ddzw(k) &
[1340]	304	& ) * ddzu(k+1)
[1015]	305	ENDIF
	306	ENDDO
	307
	308	!
	309	!-- Wall functions at all vertical walls, where necessary
	310	DO k = 1,nzt
	311
	312	IF ( k > nzb_w_inner(j,i) .AND. k <= nzb_w_outer(j,i) .AND. &
[1340]	313	wall_w_x(j,i) /= 0.0_wp .AND. wall_w_y(j,i) /= 0.0_wp ) THEN
[1015]	314	!
	315	!-- Interpolate eddy diffusivities on staggered gridpoints
[1322]	316	kmxp = 0.25_wp * &
[1015]	317	( km(k,j,i)+km(k,j,i+1)+km(k+1,j,i)+km(k+1,j,i+1) )
[1322]	318	kmxm = 0.25_wp * &
[1015]	319	( km(k,j,i)+km(k,j,i-1)+km(k+1,j,i)+km(k+1,j,i-1) )
[1322]	320	kmyp = 0.25_wp * &
[1015]	321	( km(k,j,i)+km(k+1,j,i)+km(k,j+1,i)+km(k+1,j+1,i) )
[1322]	322	kmym = 0.25_wp * &
[1015]	323	( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
	324
	325	tend(k,j,i) = tend(k,j,i) &
	326	+ ( fwxp(j,i) * ( &
	327	kmxp * ( w(k,j,i+1) - w(k,j,i) ) * ddx &
	328	+ kmxp * ( u(k+1,j,i+1) - u(k,j,i+1) ) * ddzu(k+1) &
	329	) &
	330	- fwxm(j,i) * ( &
	331	kmxm * ( w(k,j,i) - w(k,j,i-1) ) * ddx &
	332	+ kmxm * ( u(k+1,j,i) - u(k,j,i) ) * ddzu(k+1) &
	333	) &
	334	+ wall_w_x(j,i) * wsus(k,j,i) &
	335	) * ddx &
	336	+ ( fwyp(j,i) * ( &
	337	kmyp * ( w(k,j+1,i) - w(k,j,i) ) * ddy &
	338	+ kmyp * ( v(k+1,j+1,i) - v(k,j+1,i) ) * ddzu(k+1) &
	339	) &
	340	- fwym(j,i) * ( &
	341	kmym * ( w(k,j,i) - w(k,j-1,i) ) * ddy &
	342	+ kmym * ( v(k+1,j,i) - v(k,j,i) ) * ddzu(k+1) &
	343	) &
	344	+ wall_w_y(j,i) * wsvs(k,j,i) &
	345	) * ddy &
[1340]	346	+ 2.0_wp * ( &
[1015]	347	km(k+1,j,i) * ( w(k+1,j,i) - w(k,j,i) ) * ddzw(k+1) &
	348	- km(k,j,i) * ( w(k,j,i) - w(k-1,j,i) ) * ddzw(k) &
[1340]	349	) * ddzu(k+1)
[1015]	350	ENDIF
	351	ENDDO
	352
	353	ENDDO
	354	ENDDO
	355	!$acc end kernels
	356
	357	END SUBROUTINE diffusion_w_acc
	358
	359
	360	!------------------------------------------------------------------------------!
[1682]	361	! Description:
	362	! ------------
	363	!> Call for grid point i,j
[1]	364	!------------------------------------------------------------------------------!
[1001]	365	SUBROUTINE diffusion_w_ij( i, j )
[1]	366
[1320]	367	USE arrays_3d, &
	368	ONLY : ddzu, ddzw, km, tend, u, v, w
	369
	370	USE control_parameters, &
	371	ONLY : topography
	372
	373	USE grid_variables, &
	374	ONLY : ddx, ddy, fwxm, fwxp, fwym, fwyp, wall_w_x, wall_w_y
	375
	376	USE indices, &
	377	ONLY : nxl, nxr, nyn, nys, nzb, nzb_w_inner, nzb_w_outer, nzt
	378
	379	USE kinds
[1]	380
	381	IMPLICIT NONE
	382
[1682]	383	INTEGER(iwp) :: i !<
	384	INTEGER(iwp) :: j !<
	385	INTEGER(iwp) :: k !<
[1320]	386
[1682]	387	REAL(wp) :: kmxm !<
	388	REAL(wp) :: kmxp !<
	389	REAL(wp) :: kmym !<
	390	REAL(wp) :: kmyp !<
[1]	391
[1320]	392	REAL(wp), DIMENSION(nzb:nzt+1) :: wsus
	393	REAL(wp), DIMENSION(nzb:nzt+1) :: wsvs
[1]	394
[1001]	395
[1]	396	DO k = nzb_w_outer(j,i)+1, nzt-1
	397	!
	398	!-- Interpolate eddy diffusivities on staggered gridpoints
[1322]	399	kmxp = 0.25_wp * ( km(k,j,i)+km(k,j,i+1)+km(k+1,j,i)+km(k+1,j,i+1) )
	400	kmxm = 0.25_wp * ( km(k,j,i)+km(k,j,i-1)+km(k+1,j,i)+km(k+1,j,i-1) )
	401	kmyp = 0.25_wp * ( km(k,j,i)+km(k+1,j,i)+km(k,j+1,i)+km(k+1,j+1,i) )
	402	kmym = 0.25_wp * ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
[1]	403
	404	tend(k,j,i) = tend(k,j,i) &
[978]	405	& + ( kmxp * ( w(k,j,i+1) - w(k,j,i) ) * ddx &
	406	& + kmxp * ( u(k+1,j,i+1) - u(k,j,i+1) ) * ddzu(k+1) &
	407	& - kmxm * ( w(k,j,i) - w(k,j,i-1) ) * ddx &
	408	& - kmxm * ( u(k+1,j,i) - u(k,j,i) ) * ddzu(k+1) &
[1]	409	& ) * ddx &
[978]	410	& + ( kmyp * ( w(k,j+1,i) - w(k,j,i) ) * ddy &
	411	& + kmyp * ( v(k+1,j+1,i) - v(k,j+1,i) ) * ddzu(k+1) &
	412	& - kmym * ( w(k,j,i) - w(k,j-1,i) ) * ddy &
	413	& - kmym * ( v(k+1,j,i) - v(k,j,i) ) * ddzu(k+1) &
[1]	414	& ) * ddy &
[1322]	415	& + 2.0_wp * ( &
[1]	416	& km(k+1,j,i) * ( w(k+1,j,i) - w(k,j,i) ) * ddzw(k+1) &
	417	& - km(k,j,i) * ( w(k,j,i) - w(k-1,j,i) ) * ddzw(k) &
[1340]	418	& ) * ddzu(k+1)
[1]	419	ENDDO
	420
	421	!
	422	!-- Wall functions at all vertical walls, where necessary
[1340]	423	IF ( wall_w_x(j,i) /= 0.0_wp .OR. wall_w_y(j,i) /= 0.0_wp ) THEN
[51]	424
	425	!
	426	!-- Calculate the horizontal momentum fluxes w'u' and/or w'v'
[1340]	427	IF ( wall_w_x(j,i) /= 0.0_wp ) THEN
[1320]	428	CALL wall_fluxes( i, j, nzb_w_inner(j,i)+1, nzb_w_outer(j,i), &
	429	wsus, 0.0_wp, 0.0_wp, 0.0_wp, 1.0_wp )
[51]	430	ELSE
[1340]	431	wsus = 0.0_wp
[51]	432	ENDIF
	433
[1340]	434	IF ( wall_w_y(j,i) /= 0.0_wp ) THEN
[1320]	435	CALL wall_fluxes( i, j, nzb_w_inner(j,i)+1, nzb_w_outer(j,i), &
	436	wsvs, 0.0_wp, 0.0_wp, 1.0_wp, 0.0_wp )
[51]	437	ELSE
[1340]	438	wsvs = 0.0_wp
[51]	439	ENDIF
	440
[1]	441	DO k = nzb_w_inner(j,i)+1, nzb_w_outer(j,i)
	442	!
	443	!-- Interpolate eddy diffusivities on staggered gridpoints
[1322]	444	kmxp = 0.25_wp * ( km(k,j,i)+km(k,j,i+1)+km(k+1,j,i)+km(k+1,j,i+1) )
	445	kmxm = 0.25_wp * ( km(k,j,i)+km(k,j,i-1)+km(k+1,j,i)+km(k+1,j,i-1) )
	446	kmyp = 0.25_wp * ( km(k,j,i)+km(k+1,j,i)+km(k,j+1,i)+km(k+1,j+1,i) )
	447	kmym = 0.25_wp * ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
[1]	448
	449	tend(k,j,i) = tend(k,j,i) &
	450	+ ( fwxp(j,i) * ( &
[978]	451	kmxp * ( w(k,j,i+1) - w(k,j,i) ) * ddx &
	452	+ kmxp * ( u(k+1,j,i+1) - u(k,j,i+1) ) * ddzu(k+1) &
[1]	453	) &
	454	- fwxm(j,i) * ( &
[978]	455	kmxm * ( w(k,j,i) - w(k,j,i-1) ) * ddx &
	456	+ kmxm * ( u(k+1,j,i) - u(k,j,i) ) * ddzu(k+1) &
[1]	457	) &
[51]	458	+ wall_w_x(j,i) * wsus(k) &
[1]	459	) * ddx &
	460	+ ( fwyp(j,i) * ( &
[978]	461	kmyp * ( w(k,j+1,i) - w(k,j,i) ) * ddy &
	462	+ kmyp * ( v(k+1,j+1,i) - v(k,j+1,i) ) * ddzu(k+1) &
[1]	463	) &
	464	- fwym(j,i) * ( &
[978]	465	kmym * ( w(k,j,i) - w(k,j-1,i) ) * ddy &
	466	+ kmym * ( v(k+1,j,i) - v(k,j,i) ) * ddzu(k+1) &
[1]	467	) &
[51]	468	+ wall_w_y(j,i) * wsvs(k) &
[1]	469	) * ddy &
[1340]	470	+ 2.0_wp * ( &
[1]	471	km(k+1,j,i) * ( w(k+1,j,i) - w(k,j,i) ) * ddzw(k+1) &
	472	- km(k,j,i) * ( w(k,j,i) - w(k-1,j,i) ) * ddzw(k) &
[1340]	473	) * ddzu(k+1)
[1]	474	ENDDO
	475	ENDIF
	476
	477	END SUBROUTINE diffusion_w_ij
	478
[1323]	479	END MODULE diffusion_w_mod

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |