diff options
| author | schnetter <schnetter@b7a48df3-cbbf-4440-997f-b4b717c9f7fc> | 2003-10-30 17:32:02 +0000 |
|---|---|---|
| committer | schnetter <schnetter@b7a48df3-cbbf-4440-997f-b4b717c9f7fc> | 2003-10-30 17:32:02 +0000 |
| commit | bb2a1f5af1ff444f76031aa12ee6cdc309e6d33a (patch) | |
| tree | c50b44021dd92ff2c43c25316d405370f3903e40 /src | |
| parent | ed954040b8313656e91d9ae38d907a049c9dd2b0 (diff) | |
Fix typo in comment
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/ADMConstraints/trunk@95 b7a48df3-cbbf-4440-997f-b4b717c9f7fc
Diffstat (limited to 'src')
| -rw-r--r-- | src/ADMConstraints.F | 26 |
1 files changed, 23 insertions, 3 deletions
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. |