diff options
| -rw-r--r-- | src/RobinBoundary.F | 87 | ||||
| -rw-r--r-- | src/make.code.defn | 1 |
2 files changed, 0 insertions, 88 deletions
diff --git a/src/RobinBoundary.F b/src/RobinBoundary.F deleted file mode 100644 index fda5e5f..0000000 --- a/src/RobinBoundary.F +++ /dev/null @@ -1,87 +0,0 @@ -#include "cctk.h" - - subroutine RobinBCi(sw, doBC, nxyz, var, r, finf, npow) - - implicit none - - integer sw(3) - integer doBC(6) - integer nxyz(3) - - CCTK_REAL var(nxyz(1),nxyz(2),nxyz(3)) - CCTK_REAL r(nxyz(1),nxyz(2),nxyz(3)) - - CCTK_REAL finf - CCTK_INT npow - - integer nx,ny,nz - -c Grid parameters. - - nx = nxyz(1) - ny = nxyz(2) - nz = nxyz(3) - -c Lower x-bound. - - if (doBC(1).eq.1) then - var(1,:,:) = finf + ( - . r(4,:,:)**npow*( var(4,:,:) - finf ) - . -3.0d0*r(3,:,:)**npow*( var(3,:,:) - finf ) - . +3.0d0*r(2,:,:)**npow*( var(2,:,:) - finf ) - . )/r(1,:,:)**npow - end if - -c Upper x-bound. - - if (doBC(2).eq.1) then - var(nx,:,:) = finf + ( - . r(nx-3,:,:)**npow*( var(nx-3,:,:) - finf ) - . -3.0d0*r(nx-2,:,:)**npow*( var(nx-2,:,:) - finf ) - . +3.0d0*r(nx-1,:,:)**npow*( var(nx-1,:,:) - finf ) - . )/r(nx,:,:)**npow - end if - -c Lower y-bound. - - if ((doBC(3).eq.1)) then - var(:,1,:) = finf + ( - . r(:,4,:)**npow*( var(:,4,:) - finf ) - . -3.0d0*r(:,3,:)**npow*( var(:,3,:) - finf ) - . +3.0d0*r(:,2,:)**npow*( var(:,2,:) - finf ) - . )/r(:,1,:)**npow - end if - -c Upper y-bound. - - if ((doBC(4).eq.1)) then - var(:,ny,:) = finf + ( - . r(:,ny-3,:)**npow*( var(:,ny-3,:) - finf ) - . -3.0d0*r(:,ny-2,:)**npow*( var(:,ny-2,:) - finf ) - . +3.0d0*r(:,ny-1,:)**npow*( var(:,ny-1,:) - finf ) - . )/r(:,ny,:)**npow - end if - -c Lower z-bound. - - if (doBC(5).eq.1) then - var(:,:,1) = finf + ( - . r(:,:,4)**npow*( var(:,:,4) - finf ) - . -3.0d0*r(:,:,3)**npow*( var(:,:,3) - finf ) - . +3.0d0*r(:,:,2)**npow*( var(:,:,2) - finf ) - . )/r(:,:,1)**npow - end if - -c Upper z-bound. - - if (doBC(6).eq.1) then - var(:,:,nz) = finf + ( - . r(:,:,nz-3)**npow*( var(:,:,nz-3) - finf ) - . -3.0d0*r(:,:,nz-2)**npow*( var(:,:,nz-2) - finf ) - . +3.0d0*r(:,:,nz-1)**npow*( var(:,:,nz-1) - finf ) - . )/r(:,:,nz)**npow - end if - - return - end - diff --git a/src/make.code.defn b/src/make.code.defn index da181f3..31e07ca 100644 --- a/src/make.code.defn +++ b/src/make.code.defn @@ -8,7 +8,6 @@ SRCS = ConstantBoundary.c\ RadiationBoundary.c\ RadiationBoundaryWrappers.c \ RobinBoundary.c \ - RobinBoundary.F \ RobinBoundaryWrappers.c \ ScalarBoundary.c |