diff options
Diffstat (limited to 'src')
| -rw-r--r-- | src/AHFinder_exp.F | 66 |
1 files changed, 66 insertions, 0 deletions
diff --git a/src/AHFinder_exp.F b/src/AHFinder_exp.F index ac251d1..8a74a32 100644 --- a/src/AHFinder_exp.F +++ b/src/AHFinder_exp.F @@ -312,6 +312,72 @@ end do end do +! Zero out the edges before doing the boundaries + ahf_exp( 1, 1, :) = 0 + ahf_exp( 1,ny, :) = 0 + ahf_exp( 1, :, 1) = 0 + ahf_exp( 1, :,nz) = 0 + ahf_exp(nx, 1, :) = 0 + ahf_exp(nx,ny, :) = 0 + ahf_exp(nx, :, 1) = 0 + ahf_exp(nx, :,nz) = 0 + ahf_exp( :, 1, 1) = 0 + ahf_exp( :,ny, 1) = 0 + ahf_exp( :, 1,nz) = 0 + ahf_exp( :,ny,nz) = 0 + + ahfgradx( 1, 1, :) = 0 + ahfgradx( 1,ny, :) = 0 + ahfgradx( 1, :, 1) = 0 + ahfgradx( 1, :,nz) = 0 + ahfgradx(nx, 1, :) = 0 + ahfgradx(nx,ny, :) = 0 + ahfgradx(nx, :, 1) = 0 + ahfgradx(nx, :,nz) = 0 + ahfgradx( :, 1, 1) = 0 + ahfgradx( :,ny, 1) = 0 + ahfgradx( :, 1,nz) = 0 + ahfgradx( :,ny,nz) = 0 + + ahfgrady( 1, 1, :) = 0 + ahfgrady( 1,ny, :) = 0 + ahfgrady( 1, :, 1) = 0 + ahfgrady( 1, :,nz) = 0 + ahfgrady(nx, 1, :) = 0 + ahfgrady(nx,ny, :) = 0 + ahfgrady(nx, :, 1) = 0 + ahfgrady(nx, :,nz) = 0 + ahfgrady( :, 1, 1) = 0 + ahfgrady( :,ny, 1) = 0 + ahfgrady( :, 1,nz) = 0 + ahfgrady( :,ny,nz) = 0 + + ahfgradz( 1, 1, :) = 0 + ahfgradz( 1,ny, :) = 0 + ahfgradz( 1, :, 1) = 0 + ahfgradz( 1, :,nz) = 0 + ahfgradz(nx, 1, :) = 0 + ahfgradz(nx,ny, :) = 0 + ahfgradz(nx, :, 1) = 0 + ahfgradz(nx, :,nz) = 0 + ahfgradz( :, 1, 1) = 0 + ahfgradz( :,ny, 1) = 0 + ahfgradz( :, 1,nz) = 0 + ahfgradz( :,ny,nz) = 0 + + ahfgradn( 1, 1, :) = 0 + ahfgradn( 1,ny, :) = 0 + ahfgradn( 1, :, 1) = 0 + ahfgradn( 1, :,nz) = 0 + ahfgradn(nx, 1, :) = 0 + ahfgradn(nx,ny, :) = 0 + ahfgradn(nx, :, 1) = 0 + ahfgradn(nx, :,nz) = 0 + ahfgradn( :, 1, 1) = 0 + ahfgradn( :,ny, 1) = 0 + ahfgradn( :, 1,nz) = 0 + ahfgradn( :,ny,nz) = 0 + ! Boundaries on x direction. ahf_exp(1,:,:) = 2.0D0*ahf_exp(2,:,:) - ahf_exp(3,:,:) |