Home

Context Navigation

source: palm/trunk/SOURCE/diffusion_s.f90 @ 1873

Last change on this file since 1873 was 1873, checked in by maronga, 8 years ago
revised renaming of modules
Property svn:keywords set to `Id`
File size: 20.1 KB

Rev	Line
[1873]	1	!> @file diffusion_s.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	!
[1818]	16	! Copyright 1997-2016 Leibniz Universitaet Hannover
[1036]	17	!--------------------------------------------------------------------------------!
	18	!
[484]	19	! Current revisions:
[1001]	20	! ------------------
[1873]	21	! Module renamed (removed _mod)
[1341]	22	!
[1692]	23	!
[1321]	24	! Former revisions:
	25	! -----------------
	26	! $Id: diffusion_s.f90 1873 2016-04-18 14:50:06Z maronga $
	27	!
[1851]	28	! 1850 2016-04-08 13:29:27Z maronga
	29	! Module renamed
	30	!
	31	!
[1692]	32	! 1691 2015-10-26 16:17:44Z maronga
	33	! Formatting corrections.
	34	!
[1683]	35	! 1682 2015-10-07 23:56:08Z knoop
	36	! Code annotations made doxygen readable
	37	!
[1375]	38	! 1374 2014-04-25 12:55:07Z raasch
	39	! missing variables added to ONLY list
	40	!
[1341]	41	! 1340 2014-03-25 19:45:13Z kanani
	42	! REAL constants defined as wp-kind
	43	!
[1321]	44	! 1320 2014-03-20 08:40:49Z raasch
[1320]	45	! ONLY-attribute added to USE-statements,
	46	! kind-parameters added to all INTEGER and REAL declaration statements,
	47	! kinds are defined in new module kinds,
	48	! revision history before 2012 removed,
	49	! comment fields (!:) to be used for variable explanations added to
	50	! all variable declaration statements
[1321]	51	!
[1258]	52	! 1257 2013-11-08 15:18:40Z raasch
	53	! openacc loop and loop vector clauses removed
	54	!
[1132]	55	! 1128 2013-04-12 06:19:32Z raasch
	56	! loop index bounds in accelerator version replaced by i_left, i_right, j_south,
	57	! j_north
	58	!
[1093]	59	! 1092 2013-02-02 11:24:22Z raasch
	60	! unused variables removed
	61	!
[1037]	62	! 1036 2012-10-22 13:43:42Z raasch
	63	! code put under GPL (PALM 3.9)
	64	!
[1017]	65	! 1015 2012-09-27 09:23:24Z raasch
	66	! accelerator version (*_acc) added
	67	!
[1011]	68	! 1010 2012-09-20 07:59:54Z raasch
	69	! cpp switch __nopointer added for pointer free version
	70	!
[1002]	71	! 1001 2012-09-13 14:08:46Z raasch
	72	! some arrays comunicated by module instead of parameter list
	73	!
[1]	74	! Revision 1.1 2000/04/13 14:54:02 schroeter
	75	! Initial revision
	76	!
	77	!
	78	! Description:
	79	! ------------
[1682]	80	!> Diffusion term of scalar quantities (temperature and water content)
[1]	81	!------------------------------------------------------------------------------!
[1682]	82	MODULE diffusion_s_mod
	83
[1]	84
	85	PRIVATE
[1015]	86	PUBLIC diffusion_s, diffusion_s_acc
[1]	87
	88	INTERFACE diffusion_s
	89	MODULE PROCEDURE diffusion_s
	90	MODULE PROCEDURE diffusion_s_ij
	91	END INTERFACE diffusion_s
	92
[1015]	93	INTERFACE diffusion_s_acc
	94	MODULE PROCEDURE diffusion_s_acc
	95	END INTERFACE diffusion_s_acc
	96
[1]	97	CONTAINS
	98
	99
	100	!------------------------------------------------------------------------------!
[1682]	101	! Description:
	102	! ------------
	103	!> Call for all grid points
[1]	104	!------------------------------------------------------------------------------!
[1001]	105	SUBROUTINE diffusion_s( s, s_flux_b, s_flux_t, wall_s_flux )
[1]	106
[1320]	107	USE arrays_3d, &
	108	ONLY: ddzu, ddzw, kh, tend
	109
	110	USE control_parameters, &
	111	ONLY: use_surface_fluxes, use_top_fluxes
	112
	113	USE grid_variables, &
	114	ONLY: ddx2, ddy2, fwxm, fwxp, fwym, fwyp, wall_w_x, wall_w_y
	115
	116	USE indices, &
[1374]	117	ONLY: nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, &
[1320]	118	nzb_diff_s_inner, nzb_s_inner, nzb_s_outer, nzt, nzt_diff
	119
	120	USE kinds
[1]	121
	122	IMPLICIT NONE
	123
[1682]	124	INTEGER(iwp) :: i !<
	125	INTEGER(iwp) :: j !<
	126	INTEGER(iwp) :: k !<
	127	REAL(wp) :: wall_s_flux(0:4) !<
	128	REAL(wp), DIMENSION(nysg:nyng,nxlg:nxrg) :: s_flux_b, s_flux_t !<
[1010]	129	#if defined( __nopointer )
[1682]	130	REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) :: s !<
[1010]	131	#else
[1682]	132	REAL(wp), DIMENSION(:,:,:), POINTER :: s !<
[1010]	133	#endif
[1]	134
	135	DO i = nxl, nxr
	136	DO j = nys,nyn
	137	!
	138	!-- Compute horizontal diffusion
[19]	139	DO k = nzb_s_outer(j,i)+1, nzt
[1]	140
[1320]	141	tend(k,j,i) = tend(k,j,i) &
[1340]	142	+ 0.5_wp * ( &
[1320]	143	( kh(k,j,i) + kh(k,j,i+1) ) * ( s(k,j,i+1)-s(k,j,i) ) &
	144	- ( kh(k,j,i) + kh(k,j,i-1) ) * ( s(k,j,i)-s(k,j,i-1) ) &
[1340]	145	) * ddx2 &
	146	+ 0.5_wp * ( &
[1320]	147	( kh(k,j,i) + kh(k,j+1,i) ) * ( s(k,j+1,i)-s(k,j,i) ) &
	148	- ( kh(k,j,i) + kh(k,j-1,i) ) * ( s(k,j,i)-s(k,j-1,i) ) &
[1340]	149	) * ddy2
[1]	150	ENDDO
	151
	152	!
	153	!-- Apply prescribed horizontal wall heatflux where necessary
[1691]	154	IF ( ( wall_w_x(j,i) /= 0.0_wp ) .OR. ( wall_w_y(j,i) /= 0.0_wp ) &
	155	) THEN
[1]	156	DO k = nzb_s_inner(j,i)+1, nzb_s_outer(j,i)
	157
[1320]	158	tend(k,j,i) = tend(k,j,i) &
[1340]	159	+ ( fwxp(j,i) * 0.5_wp * &
[1320]	160	( kh(k,j,i) + kh(k,j,i+1) ) * ( s(k,j,i+1)-s(k,j,i) ) &
[1340]	161	+ ( 1.0_wp - fwxp(j,i) ) * wall_s_flux(1) &
	162	-fwxm(j,i) * 0.5_wp * &
[1320]	163	( kh(k,j,i) + kh(k,j,i-1) ) * ( s(k,j,i)-s(k,j,i-1) ) &
[1340]	164	+ ( 1.0_wp - fwxm(j,i) ) * wall_s_flux(2) &
[1320]	165	) * ddx2 &
[1340]	166	+ ( fwyp(j,i) * 0.5_wp * &
[1320]	167	( kh(k,j,i) + kh(k,j+1,i) ) * ( s(k,j+1,i)-s(k,j,i) ) &
[1340]	168	+ ( 1.0_wp - fwyp(j,i) ) * wall_s_flux(3) &
	169	-fwym(j,i) * 0.5_wp * &
[1320]	170	( kh(k,j,i) + kh(k,j-1,i) ) * ( s(k,j,i)-s(k,j-1,i) ) &
[1340]	171	+ ( 1.0_wp - fwym(j,i) ) * wall_s_flux(4) &
[1]	172	) * ddy2
	173	ENDDO
	174	ENDIF
	175
	176	!
	177	!-- Compute vertical diffusion. In case that surface fluxes have been
[19]	178	!-- prescribed or computed at bottom and/or top, index k starts/ends at
	179	!-- nzb+2 or nzt-1, respectively.
	180	DO k = nzb_diff_s_inner(j,i), nzt_diff
[1]	181
[1320]	182	tend(k,j,i) = tend(k,j,i) &
[1340]	183	+ 0.5_wp * ( &
[1320]	184	( kh(k,j,i) + kh(k+1,j,i) ) * ( s(k+1,j,i)-s(k,j,i) ) * ddzu(k+1) &
	185	- ( kh(k,j,i) + kh(k-1,j,i) ) * ( s(k,j,i)-s(k-1,j,i) ) * ddzu(k) &
[1340]	186	) * ddzw(k)
[1]	187	ENDDO
	188
	189	!
[19]	190	!-- Vertical diffusion at the first computational gridpoint along
[1]	191	!-- z-direction
	192	IF ( use_surface_fluxes ) THEN
	193
	194	k = nzb_s_inner(j,i)+1
	195
[1320]	196	tend(k,j,i) = tend(k,j,i) &
[1340]	197	+ ( 0.5_wp * ( kh(k,j,i)+kh(k+1,j,i) ) &
	198	* ( s(k+1,j,i)-s(k,j,i) ) &
	199	* ddzu(k+1) &
[1320]	200	+ s_flux_b(j,i) &
[1]	201	) * ddzw(k)
	202
	203	ENDIF
	204
[19]	205	!
	206	!-- Vertical diffusion at the last computational gridpoint along
	207	!-- z-direction
	208	IF ( use_top_fluxes ) THEN
	209
	210	k = nzt
	211
[1320]	212	tend(k,j,i) = tend(k,j,i) &
	213	+ ( - s_flux_t(j,i) &
[1340]	214	- 0.5_wp * ( kh(k-1,j,i)+kh(k,j,i) )&
	215	* ( s(k,j,i)-s(k-1,j,i) ) &
	216	* ddzu(k) &
[19]	217	) * ddzw(k)
	218
	219	ENDIF
	220
[1]	221	ENDDO
	222	ENDDO
	223
	224	END SUBROUTINE diffusion_s
	225
	226
	227	!------------------------------------------------------------------------------!
[1682]	228	! Description:
	229	! ------------
	230	!> Call for all grid points - accelerator version
[1015]	231	!------------------------------------------------------------------------------!
	232	SUBROUTINE diffusion_s_acc( s, s_flux_b, s_flux_t, wall_s_flux )
	233
[1320]	234	USE arrays_3d, &
	235	ONLY: ddzu, ddzw, kh, tend
	236
	237	USE control_parameters, &
	238	ONLY: use_surface_fluxes, use_top_fluxes
	239
	240	USE grid_variables, &
	241	ONLY: ddx2, ddy2, fwxm, fwxp, fwym, fwyp, wall_w_x, wall_w_y
	242
	243	USE indices, &
	244	ONLY: i_left, i_right, j_north, j_south, nxlg, nxrg, nyng, nysg, &
[1374]	245	nzb, nzb_diff_s_inner, nzb_s_inner, nzb_s_outer, nzt, nzt_diff
[1320]	246
	247	USE kinds
[1015]	248
	249	IMPLICIT NONE
	250
[1682]	251	INTEGER(iwp) :: i !<
	252	INTEGER(iwp) :: j !<
	253	INTEGER(iwp) :: k !<
	254	REAL(wp) :: wall_s_flux(0:4) !<
	255	REAL(wp), DIMENSION(nysg:nyng,nxlg:nxrg) :: s_flux_b !<
	256	REAL(wp), DIMENSION(nysg:nyng,nxlg:nxrg) :: s_flux_t !<
[1015]	257	#if defined( __nopointer )
[1682]	258	REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) :: s !<
[1015]	259	#else
[1682]	260	REAL(wp), DIMENSION(:,:,:), POINTER :: s !<
[1015]	261	#endif
	262
	263	!$acc kernels present( ddzu, ddzw, fwxm, fwxp, fwym, fwyp, kh ) &
	264	!$acc present( nzb_diff_s_inner, nzb_s_inner, nzb_s_outer, s ) &
	265	!$acc present( s_flux_b, s_flux_t, tend, wall_s_flux ) &
	266	!$acc present( wall_w_x, wall_w_y )
[1128]	267	DO i = i_left, i_right
	268	DO j = j_south, j_north
[1015]	269	!
	270	!-- Compute horizontal diffusion
	271	DO k = 1, nzt
	272	IF ( k > nzb_s_outer(j,i) ) THEN
	273
[1320]	274	tend(k,j,i) = tend(k,j,i) &
[1340]	275	+ 0.5_wp * ( &
[1320]	276	( kh(k,j,i) + kh(k,j,i+1) ) * ( s(k,j,i+1)-s(k,j,i) ) &
	277	- ( kh(k,j,i) + kh(k,j,i-1) ) * ( s(k,j,i)-s(k,j,i-1) ) &
[1340]	278	) * ddx2 &
	279	+ 0.5_wp * ( &
[1320]	280	( kh(k,j,i) + kh(k,j+1,i) ) * ( s(k,j+1,i)-s(k,j,i) ) &
	281	- ( kh(k,j,i) + kh(k,j-1,i) ) * ( s(k,j,i)-s(k,j-1,i) ) &
[1340]	282	) * ddy2
[1015]	283	ENDIF
	284	ENDDO
	285
	286	!
	287	!-- Apply prescribed horizontal wall heatflux where necessary
	288	DO k = 1, nzt
	289	IF ( k > nzb_s_inner(j,i) .AND. k <= nzb_s_outer(j,i) .AND. &
[1340]	290	( wall_w_x(j,i) /= 0.0_wp .OR. wall_w_y(j,i) /= 0.0_wp ) ) &
[1015]	291	THEN
[1320]	292	tend(k,j,i) = tend(k,j,i) &
[1340]	293	+ ( fwxp(j,i) * 0.5_wp * &
[1320]	294	( kh(k,j,i) + kh(k,j,i+1) ) * ( s(k,j,i+1)-s(k,j,i) ) &
[1340]	295	+ ( 1.0_wp - fwxp(j,i) ) * wall_s_flux(1) &
	296	-fwxm(j,i) * 0.5_wp * &
[1320]	297	( kh(k,j,i) + kh(k,j,i-1) ) * ( s(k,j,i)-s(k,j,i-1) ) &
[1340]	298	+ ( 1.0_wp - fwxm(j,i) ) * wall_s_flux(2) &
[1320]	299	) * ddx2 &
[1340]	300	+ ( fwyp(j,i) * 0.5_wp * &
[1320]	301	( kh(k,j,i) + kh(k,j+1,i) ) * ( s(k,j+1,i)-s(k,j,i) ) &
[1340]	302	+ ( 1.0_wp - fwyp(j,i) ) * wall_s_flux(3) &
	303	-fwym(j,i) * 0.5_wp * &
[1320]	304	( kh(k,j,i) + kh(k,j-1,i) ) * ( s(k,j,i)-s(k,j-1,i) ) &
[1340]	305	+ ( 1.0_wp - fwym(j,i) ) * wall_s_flux(4) &
[1015]	306	) * ddy2
	307	ENDIF
	308	ENDDO
	309
	310	!
	311	!-- Compute vertical diffusion. In case that surface fluxes have been
	312	!-- prescribed or computed at bottom and/or top, index k starts/ends at
	313	!-- nzb+2 or nzt-1, respectively.
	314	DO k = 1, nzt_diff
	315	IF ( k >= nzb_diff_s_inner(j,i) ) THEN
[1320]	316	tend(k,j,i) = tend(k,j,i) &
[1340]	317	+ 0.5_wp * ( &
[1320]	318	( kh(k,j,i) + kh(k+1,j,i) ) * ( s(k+1,j,i)-s(k,j,i) ) * ddzu(k+1) &
	319	- ( kh(k,j,i) + kh(k-1,j,i) ) * ( s(k,j,i)-s(k-1,j,i) ) * ddzu(k) &
[1340]	320	) * ddzw(k)
[1015]	321	ENDIF
	322	ENDDO
	323
	324	!
	325	!-- Vertical diffusion at the first computational gridpoint along
	326	!-- z-direction
	327	DO k = 1, nzt
	328	IF ( use_surface_fluxes .AND. k == nzb_s_inner(j,i)+1 ) THEN
[1320]	329	tend(k,j,i) = tend(k,j,i) &
[1340]	330	+ ( 0.5_wp * ( kh(k,j,i)+kh(k+1,j,i) )&
	331	* ( s(k+1,j,i)-s(k,j,i) ) &
	332	* ddzu(k+1) &
[1320]	333	+ s_flux_b(j,i) &
[1015]	334	) * ddzw(k)
	335	ENDIF
	336
	337	!
	338	!-- Vertical diffusion at the last computational gridpoint along
	339	!-- z-direction
	340	IF ( use_top_fluxes .AND. k == nzt ) THEN
	341	tend(k,j,i) = tend(k,j,i) &
	342	+ ( - s_flux_t(j,i) &
[1340]	343	- 0.5_wp * ( kh(k-1,j,i)+kh(k,j,i) )&
	344	* ( s(k,j,i)-s(k-1,j,i) ) &
	345	* ddzu(k) &
[1015]	346	) * ddzw(k)
	347	ENDIF
	348	ENDDO
	349
	350	ENDDO
	351	ENDDO
	352	!$acc end kernels
	353
	354	END SUBROUTINE diffusion_s_acc
	355
	356
	357	!------------------------------------------------------------------------------!
[1682]	358	! Description:
	359	! ------------
	360	!> Call for grid point i,j
[1]	361	!------------------------------------------------------------------------------!
[1001]	362	SUBROUTINE diffusion_s_ij( i, j, s, s_flux_b, s_flux_t, wall_s_flux )
[1]	363
[1320]	364	USE arrays_3d, &
	365	ONLY: ddzu, ddzw, kh, tend
	366
	367	USE control_parameters, &
	368	ONLY: use_surface_fluxes, use_top_fluxes
	369
	370	USE grid_variables, &
	371	ONLY: ddx2, ddy2, fwxm, fwxp, fwym, fwyp, wall_w_x, wall_w_y
	372
	373	USE indices, &
[1374]	374	ONLY: nxlg, nxrg, nyng, nysg, nzb, nzb_diff_s_inner, nzb_s_inner, &
[1320]	375	nzb_s_outer, nzt, nzt_diff
	376
	377	USE kinds
[1]	378
	379	IMPLICIT NONE
	380
[1682]	381	INTEGER(iwp) :: i !<
	382	INTEGER(iwp) :: j !<
	383	INTEGER(iwp) :: k !<
	384	REAL(wp) :: wall_s_flux(0:4) !<
	385	REAL(wp), DIMENSION(nysg:nyng,nxlg:nxrg) :: s_flux_b !<
	386	REAL(wp), DIMENSION(nysg:nyng,nxlg:nxrg) :: s_flux_t !<
[1010]	387	#if defined( __nopointer )
[1682]	388	REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) :: s !<
[1010]	389	#else
[1682]	390	REAL(wp), DIMENSION(:,:,:), POINTER :: s !<
[1010]	391	#endif
[1]	392
	393	!
	394	!-- Compute horizontal diffusion
[19]	395	DO k = nzb_s_outer(j,i)+1, nzt
[1]	396
[1320]	397	tend(k,j,i) = tend(k,j,i) &
[1340]	398	+ 0.5_wp * ( &
[1320]	399	( kh(k,j,i) + kh(k,j,i+1) ) * ( s(k,j,i+1)-s(k,j,i) ) &
	400	- ( kh(k,j,i) + kh(k,j,i-1) ) * ( s(k,j,i)-s(k,j,i-1) ) &
[1340]	401	) * ddx2 &
	402	+ 0.5_wp * ( &
[1320]	403	( kh(k,j,i) + kh(k,j+1,i) ) * ( s(k,j+1,i)-s(k,j,i) ) &
	404	- ( kh(k,j,i) + kh(k,j-1,i) ) * ( s(k,j,i)-s(k,j-1,i) ) &
[1340]	405	) * ddy2
[1]	406	ENDDO
	407
	408	!
	409	!-- Apply prescribed horizontal wall heatflux where necessary
[1691]	410	IF ( ( wall_w_x(j,i) /= 0.0_wp ) .OR. ( wall_w_y(j,i) /= 0.0_wp ) ) &
[1]	411	THEN
	412	DO k = nzb_s_inner(j,i)+1, nzb_s_outer(j,i)
	413
[1320]	414	tend(k,j,i) = tend(k,j,i) &
[1340]	415	+ ( fwxp(j,i) * 0.5_wp * &
[1320]	416	( kh(k,j,i) + kh(k,j,i+1) ) * ( s(k,j,i+1)-s(k,j,i) ) &
[1340]	417	+ ( 1.0_wp - fwxp(j,i) ) * wall_s_flux(1) &
	418	-fwxm(j,i) * 0.5_wp * &
[1320]	419	( kh(k,j,i) + kh(k,j,i-1) ) * ( s(k,j,i)-s(k,j,i-1) ) &
[1340]	420	+ ( 1.0_wp - fwxm(j,i) ) * wall_s_flux(2) &
[1320]	421	) * ddx2 &
[1340]	422	+ ( fwyp(j,i) * 0.5_wp * &
[1320]	423	( kh(k,j,i) + kh(k,j+1,i) ) * ( s(k,j+1,i)-s(k,j,i) ) &
[1340]	424	+ ( 1.0_wp - fwyp(j,i) ) * wall_s_flux(3) &
	425	-fwym(j,i) * 0.5_wp * &
[1320]	426	( kh(k,j,i) + kh(k,j-1,i) ) * ( s(k,j,i)-s(k,j-1,i) ) &
[1340]	427	+ ( 1.0_wp - fwym(j,i) ) * wall_s_flux(4) &
[1]	428	) * ddy2
	429	ENDDO
	430	ENDIF
	431
	432	!
	433	!-- Compute vertical diffusion. In case that surface fluxes have been
[19]	434	!-- prescribed or computed at bottom and/or top, index k starts/ends at
	435	!-- nzb+2 or nzt-1, respectively.
	436	DO k = nzb_diff_s_inner(j,i), nzt_diff
[1]	437
[1320]	438	tend(k,j,i) = tend(k,j,i) &
[1340]	439	+ 0.5_wp * ( &
[1320]	440	( kh(k,j,i) + kh(k+1,j,i) ) * ( s(k+1,j,i)-s(k,j,i) ) * ddzu(k+1) &
	441	- ( kh(k,j,i) + kh(k-1,j,i) ) * ( s(k,j,i)-s(k-1,j,i) ) * ddzu(k) &
[1340]	442	) * ddzw(k)
[1]	443	ENDDO
	444
	445	!
[19]	446	!-- Vertical diffusion at the first computational gridpoint along z-direction
[1]	447	IF ( use_surface_fluxes ) THEN
	448
	449	k = nzb_s_inner(j,i)+1
	450
[1340]	451	tend(k,j,i) = tend(k,j,i) + ( 0.5_wp * ( kh(k,j,i)+kh(k+1,j,i) ) &
	452	* ( s(k+1,j,i)-s(k,j,i) ) &
	453	* ddzu(k+1) &
[1320]	454	+ s_flux_b(j,i) &
[19]	455	) * ddzw(k)
[1]	456
	457	ENDIF
	458
[19]	459	!
	460	!-- Vertical diffusion at the last computational gridpoint along z-direction
	461	IF ( use_top_fluxes ) THEN
	462
	463	k = nzt
	464
[1320]	465	tend(k,j,i) = tend(k,j,i) + ( - s_flux_t(j,i) &
[1340]	466	- 0.5_wp * ( kh(k-1,j,i)+kh(k,j,i) ) &
	467	* ( s(k,j,i)-s(k-1,j,i) ) &
	468	* ddzu(k) &
[19]	469	) * ddzw(k)
	470
	471	ENDIF
	472
[1]	473	END SUBROUTINE diffusion_s_ij
	474
	475	END MODULE diffusion_s_mod

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |