diff options
| -rw-r--r-- | src/ADMConstraints.F | 14 |
1 files changed, 7 insertions, 7 deletions
diff --git a/src/ADMConstraints.F b/src/ADMConstraints.F index 6dc45a7..35d2f83 100644 --- a/src/ADMConstraints.F +++ b/src/ADMConstraints.F @@ -118,7 +118,7 @@ c Calculate constraints. docalc = .TRUE. if (use_mask .eq. 1) then - if (abs(emask(i,j,k)) < 0.001) then + if (abs(emask(i,j,k)-1) > 0.001) then docalc = .FALSE. end if end if @@ -267,18 +267,18 @@ c Apply symmetry boundary conditions. c Apply flat boundary conditions at outer boundaries. if (CCTK_Equals(bound,"flat") == 1) then - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 0, -1, + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, $ "admconstraints::hamiltonian", "Flat") - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 0, -1, + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, $ "admconstraints::normalized_hamiltonian", "Flat") - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 0, -1, + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, $ "admconstraints::momentum", "Flat") else - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 0, -1, + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, $ "admconstraints::hamiltonian", "None") - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 0, -1, + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, $ "admconstraints::normalized_hamiltonian", "None") - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 0, -1, + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, $ "admconstraints::momentum", "None") end if |