gkowal/amun-code
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Remove unnecessary dimension variable fm from bnd_prol().

Browse Source
This commit is contained in:
Grzegorz Kowal 2010-02-12 15:38:05 -02:00
parent 25eee11a94
commit ea00024b65
2 changed files with 5 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
src/boundaries.F90
View File

@ -861,7 +861,7 @@ module boundaries
! local arrays ! local arrays
! !
integer, dimension(3) :: dm, cm, fm, pm integer, dimension(3) :: dm, cm, pm
real , dimension(:,:,:), allocatable :: ux real , dimension(:,:,:), allocatable :: ux
! parameters ! parameters
@ -880,12 +880,10 @@ module boundaries
dm(idir) = ng dm(idir) = ng
cm(:) = dm(:) / 2 + 2 * del cm(:) = dm(:) / 2 + 2 * del
cm(idir) = ng / 2 + 2 * del cm(idir) = ng / 2 + 2 * del
fm(:) = 2 * cm(:)
pm(:) = 2 * dm(:) pm(:) = 2 * dm(:)
#if NDIMS == 2 #if NDIMS == 2
dm(3) = 1 dm(3) = 1
cm(3) = 1 cm(3) = 1
fm(3) = 1
pm(3) = 1 pm(3) = 1
#endif /* NDIMS == 2 */ #endif /* NDIMS == 2 */
@ -1163,20 +1161,20 @@ module boundaries
#ifdef MHD #ifdef MHD
! expand face centered Bx ! expand face centered Bx
! !
call expand(cm, fm, 0, un(ibx,i1:i2,j1:j2,k1:k2), ux(:,:,:), 'c', 'l', 'l') call expand(cm, pm, 0, un(ibx,i1:i2,j1:j2,k1:k2), ux(:,:,:), 'c', 'l', 'l')
pdata%u(ibx,is:it,js:jt,ks:kt) = ux(il:iu,jl:ju,kl:ku) pdata%u(ibx,is:it,js:jt,ks:kt) = ux(il:iu,jl:ju,kl:ku)
! expand face centered By ! expand face centered By
! !
call expand(cm, fm, 0, un(iby,i1:i2,j1:j2,k1:k2), ux(:,:,:), 'l', 'c', 'l') call expand(cm, pm, 0, un(iby,i1:i2,j1:j2,k1:k2), ux(:,:,:), 'l', 'c', 'l')
pdata%u(iby,is:it,js:jt,ks:kt) = ux(il:iu,jl:ju,kl:ku) pdata%u(iby,is:it,js:jt,ks:kt) = ux(il:iu,jl:ju,kl:ku)
#if NDIMS == 3 #if NDIMS == 3
! expand face centered Bz ! expand face centered Bz
! !
call expand(cm, fm, 0, un(ibz,i1:i2,j1:j2,k1:k2), ux(:,:,:), 'l', 'l', 'c') call expand(cm, pm, 0, un(ibz,i1:i2,j1:j2,k1:k2), ux(:,:,:), 'l', 'l', 'c')
pdata%u(ibz,is:it,js:jt,ks:kt) = ux(il:iu,jl:ju,kl:ku) pdata%u(ibz,is:it,js:jt,ks:kt) = ux(il:iu,jl:ju,kl:ku)
#endif /* NDIMS == 3 */ #endif /* NDIMS == 3 */

2
src/hosts/acrux.home
View File

@ -9,7 +9,7 @@
# MPI = Y or N to enable or disable MPI parallelization # MPI = Y or N to enable or disable MPI parallelization
# OMP = Y or N to enable or disable OpenMP parallelization # OMP = Y or N to enable or disable OpenMP parallelization
# #
MPI = Y MPI = N
OMP = N OMP = N
# compilation options: # compilation options: