diff options
-rw-r--r-- | src/ADMConstraints.F | 17 |
1 files changed, 10 insertions, 7 deletions
diff --git a/src/ADMConstraints.F b/src/ADMConstraints.F index 36308c5..803d2d8 100644 --- a/src/ADMConstraints.F +++ b/src/ADMConstraints.F @@ -283,25 +283,28 @@ c = - (T_{i0} - beta^j T_{ij})/alpha DECLARE_CCTK_PARAMETERS DECLARE_CCTK_FUNCTIONS + CCTK_INT, parameter :: izero = 0 + integer, parameter :: ik = kind (izero) + c Return code from Cactus sync routine and boundary conditions. integer ierr c Apply flat boundary conditions at outer boundaries. - if (CCTK_Equals(bound,"flat") == 1) then - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + if (CCTK_EQUALS(bound,"flat")) then + ierr = Boundary_SelectGroupForBC(cctkGH, int(CCTK_ALL_FACES,ik), 1_ik, -1_ik, $ "admconstraints::hamiltonian", "Flat") - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + ierr = Boundary_SelectGroupForBC(cctkGH, int(CCTK_ALL_FACES,ik), 1_ik, -1_ik, $ "admconstraints::normalized_hamiltonian", "Flat") - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + ierr = Boundary_SelectGroupForBC(cctkGH, int(CCTK_ALL_FACES,ik), 1_ik, -1_ik, $ "admconstraints::momentum", "Flat") else - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + ierr = Boundary_SelectGroupForBC(cctkGH, int(CCTK_ALL_FACES,ik), 1_ik, -1_ik, $ "admconstraints::hamiltonian", "None") - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + ierr = Boundary_SelectGroupForBC(cctkGH, int(CCTK_ALL_FACES,ik), 1_ik, -1_ik, $ "admconstraints::normalized_hamiltonian", "None") - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + ierr = Boundary_SelectGroupForBC(cctkGH, int(CCTK_ALL_FACES,ik), 1_ik, -1_ik, $ "admconstraints::momentum", "None") end if |