diff options
Diffstat (limited to 'src/ApplySymmetry.c')
-rw-r--r-- | src/ApplySymmetry.c | 58 |
1 files changed, 0 insertions, 58 deletions
diff --git a/src/ApplySymmetry.c b/src/ApplySymmetry.c deleted file mode 100644 index a0b5ba3..0000000 --- a/src/ApplySymmetry.c +++ /dev/null @@ -1,58 +0,0 @@ - -/*@@ - @routine FortranSym - @date Mon Mar 15 15:51:57 1999 - @author Gerd Lanfermann - @desc - Routine performs the symmetry boundary operations. - @enddesc - @calls - @calledby - @history - - @endhistory -@@*/ - - subroutine FortranSym(nxyz,var,nghostzones,sym,doSym) - - implicit none - - INTEGER nxyz(3) - CCTK_REAL var(nxyz(1),nxyz(2),nxyz(3)) - INTEGER nghostzones - INTEGER sym(3) - INTEGER doSym(3) - - INTEGER sw - -c Apply symmetry if -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(1) == 1) then - - do sw=1, nghostzones - var(sw,:,:) = sym(1)*var(nghostzones+sw,:,:) - end do - - end if - - if (doSym(2) == 1) then - - do sw=1, nghostzones - var(:,sw,:) = sym(2)*var(:,nghostzones+sw,:) - end do - - end if - - if (doSym(3).eq.1) then - - do sw=1, nghostzones - var(:,:,sw) = sym(3)*var(:,:,nghostzones+sw) - end do - - end if - - end subroutine - |