diff options
Diffstat (limited to 'src/SymmetryCondition.F77')
| -rw-r--r-- | src/SymmetryCondition.F77 | 16 |
1 files changed, 10 insertions, 6 deletions
diff --git a/src/SymmetryCondition.F77 b/src/SymmetryCondition.F77 index 762ac72..20e7aa5 100644 --- a/src/SymmetryCondition.F77 +++ b/src/SymmetryCondition.F77 @@ -24,8 +24,8 @@ INTEGER nxyz(3) CCTK_REAL var(nxyz(1),nxyz(2),nxyz(3)) INTEGER nghostzones - INTEGER sym(3) - INTEGER doSym(3) + INTEGER sym(6) + INTEGER doSym(6) INTEGER sw @@ -34,6 +34,10 @@ c * the grid chunk has a physical boundary (bbox) c * its size in a direction is bigger than one (sh) c * we actually want a symmetry (sx.ne.ESYM_UNSET) + if (doSym(2) .eq. 1 .or. doSym(4) .eq. 1 .or. doSym(6) .eq. 1) then + call CCTKi_NotYetImplemented("Right hand side boundary conditions") + end if + if (doSym(1) .eq. 1) then do sw=1, nghostzones @@ -46,24 +50,24 @@ c * we actually want a symmetry (sx.ne.ESYM_UNSET) end if - if (doSym(2) .eq. 1) then + if (doSym(3) .eq. 1) then do sw=1, nghostzones do i=1,nxyz(1) do k=1,nxyz(3) - var(i,sw,k) = sym(2)*var(i,nghostzones+sw,k) + var(i,sw,k) = sym(3)*var(i,nghostzones+sw,k) end do end do end do end if - if (doSym(3).eq.1) then + if (doSym(5).eq.1) then do sw=1, nghostzones do i=1,nxyz(1) do j=1,nxyz(2) - var(i,j,sw) = sym(3)*var(i,j,nghostzones+sw) + var(i,j,sw) = sym(5)*var(i,j,nghostzones+sw) end do end do end do |