diff options
| -rw-r--r-- | interface.ccl | 2 | ||||
| -rw-r--r-- | src/ADMConstraints.F | 8 |
2 files changed, 9 insertions, 1 deletions
diff --git a/interface.ccl b/interface.ccl index 4421e12..6429fd8 100644 --- a/interface.ccl +++ b/interface.ccl @@ -2,7 +2,7 @@ # $Header$ implements: admconstraints -inherits: ADMBase, StaticConformal, SpaceMask, grid +inherits: ADMBase, StaticConformal, SpaceMask, grid, ADMMacros friend: ADMCoupling USES INCLUDE: CalcTmunu.inc 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 |