diff options
| author | hawke <hawke@b7a48df3-cbbf-4440-997f-b4b717c9f7fc> | 2003-07-09 08:46:15 +0000 |
|---|---|---|
| committer | hawke <hawke@b7a48df3-cbbf-4440-997f-b4b717c9f7fc> | 2003-07-09 08:46:15 +0000 |
| commit | b1b623e51da66ab731dda3246a4f66e16600dd9d (patch) | |
| tree | ece3460bc6fd49d57d39440c8b3286f6773a519b /src | |
| parent | 5874f2056f8e9fe3628b319392aac529d101d965 (diff) | |
Change the local_spatial_order within the loop to cope with boundaries when using fourth order differencing.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/ADMConstraints/trunk@90 b7a48df3-cbbf-4440-997f-b4b717c9f7fc
Diffstat (limited to 'src')
| -rw-r--r-- | src/ADMConstraints.F | 8 |
1 files changed, 8 insertions, 0 deletions
diff --git a/src/ADMConstraints.F b/src/ADMConstraints.F index cb8db2f..f896b37 100644 --- a/src/ADMConstraints.F +++ b/src/ADMConstraints.F @@ -111,6 +111,14 @@ c Calculate constraints. do j=2,ny-1 do i=2,nx-1 + if ( (i<3).or.(i>cctk_lsh(1)-2).or. + . (j<3).or.(j>cctk_lsh(2)-2).or. + . (k<3).or.(k>cctk_lsh(3)-2) ) then + local_spatial_order = 2 + else + local_spatial_order = spatial_order + end if + docalc = .TRUE. if (use_mask .eq. 1) then |