From bb2a1f5af1ff444f76031aa12ee6cdc309e6d33a Mon Sep 17 00:00:00 2001 From: schnetter Date: Thu, 30 Oct 2003 17:32:02 +0000 Subject: Fix typo in comment git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/ADMConstraints/trunk@95 b7a48df3-cbbf-4440-997f-b4b717c9f7fc --- src/ADMConstraints.F | 26 +++++++++++++++++++++++--- 1 file changed, 23 insertions(+), 3 deletions(-) (limited to 'src') diff --git a/src/ADMConstraints.F b/src/ADMConstraints.F index 8bf3116..87a27a5 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 (emask(i,j,k) .ne. 1) then + if (abs(emask(i,j,k) - 1) < 0.001) then docalc = .FALSE. end if end if @@ -261,13 +261,33 @@ c = - (T_{i0} - beta^j T_{ij})/alpha c Apply symmetry boundary conditions. call CartSymGN(ierr,cctkGH,"admconstraints::hamiltonian") + call CartSymGN(ierr,cctkGH,"admconstraints::normalized_hamiltonian") call CartSymGN(ierr,cctkGH,"admconstraints::momentum") c Apply flat boundary conditions at outer boundaries. if (CCTK_Equals(bound,"flat") == 1) then - call BndFlatGN(ierr,cctkGH,sw,"admconstraints::hamiltonian") - call BndFlatGN(ierr,cctkGH,sw,"admconstraints::momentum") + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 0, -1, + $ "admconstraints::hamiltonian", "Flat") + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 0, -1, + $ "admconstraints::normalized_hamiltonian", "Flat") + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 0, -1, + $ "admconstraints::momentum", "Flat") + else + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 0, -1, + $ "admconstraints::hamiltonian", "None") + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 0, -1, + $ "admconstraints::normalized_hamiltonian", "None") + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 0, -1, + $ "admconstraints::momentum", "None") + end if + +c Synchronize. + + if (constraint_communication.eq.1) then + call CCTK_SyncGroup(ierr,cctkGH,"admconstraints::hamiltonian") + call CCTK_SyncGroup(ierr,cctkGH,"admconstraints::normalized_hamiltonian") + call CCTK_SyncGroup(ierr,cctkGH,"admconstraints::momentum") end if c Cartoon. -- cgit v1.2.3