diff options
| author | allen <allen@b7a48df3-cbbf-4440-997f-b4b717c9f7fc> | 2002-06-02 21:06:08 +0000 |
|---|---|---|
| committer | allen <allen@b7a48df3-cbbf-4440-997f-b4b717c9f7fc> | 2002-06-02 21:06:08 +0000 |
| commit | d9d0dfb31bfd95f930b0437505c22e17f1d9802d (patch) | |
| tree | 8dc96d8b50ebf45d173e143f494eafc2b47fef86 /src | |
| parent | 8977fd2cb4471de3f3c3dbe733820417265bf5c9 (diff) | |
Fixed my bad logic for deciding on whether a point is excised or not.
This was causing seg faults on T3E and Alphas.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/ADMConstraints/trunk@75 b7a48df3-cbbf-4440-997f-b4b717c9f7fc
Diffstat (limited to 'src')
| -rw-r--r-- | src/ADMConstraints.F | 13 |
1 files changed, 11 insertions, 2 deletions
diff --git a/src/ADMConstraints.F b/src/ADMConstraints.F index 7715ce9..41e9224 100644 --- a/src/ADMConstraints.F +++ b/src/ADMConstraints.F @@ -41,6 +41,7 @@ c Stencil width used for calculating constraints c (for outer boundary condition) integer sw(3) + logical docalc c Return code from Cactus sync routine and boundary conditions. @@ -109,8 +110,16 @@ c Calculate constraints. do k=2,nz-1 do j=2,ny-1 do i=2,nx-1 - - if ((use_mask .eq. 0) .or. (emask(i,j,k) .eq. 1)) then + + docalc = 1 + + if (use_mask .eq. 1) then + if (emask(i,j,k) .ne. 1) then + docalc = 0 + end if + end if + + if (docalc) then ialp = 1.0D0/alp(i,j,k) ialp2 = ialp**2 |