diff options
| author | tradke <tradke@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2002-06-03 12:18:05 +0000 |
|---|---|---|
| committer | tradke <tradke@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2002-06-03 12:18:05 +0000 |
| commit | 9640fac572718259e4b4c850fd6d9ec163f35f42 (patch) | |
| tree | 60ad009271927c3b75ff38029625e7cdff3e24ee | |
| parent | 32f4172947cb9d96bcbc2854bf0c28971bb2209f (diff) | |
Zero out the edges of some grid functions before calculating their boundaries.
This fixes core dumps on Alphas and finally lets the testsuites pass.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinder/trunk@307 89daf98e-ef62-4674-b946-b8ff9de2216c
| -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,:,:) |