diff options
| -rw-r--r-- | interface.ccl | 4 | ||||
| -rw-r--r-- | src/WaveToy.F90 | 12 |
2 files changed, 8 insertions, 8 deletions
diff --git a/interface.ccl b/interface.ccl index 5b2180a..587729d 100644 --- a/interface.ccl +++ b/interface.ccl @@ -15,6 +15,6 @@ cctk_real scalarevolve type=GF Timelevels=3 CCTK_INT FUNCTION Boundary_SelectVarForBC(CCTK_POINTER IN GH, \ - CCTK_INT IN faces, CCTK_INT IN table_handle, CCTK_STRING IN var_name, \ - CCTK_STRING IN bc_name) + CCTK_INT IN faces, CCTK_INT IN boundary_width, CCTK_INT IN table_handle, \ + CCTK_STRING IN var_name, CCTK_STRING IN bc_name) USES FUNCTION Boundary_SelectVarForBC diff --git a/src/WaveToy.F90 b/src/WaveToy.F90 index 16b1677..4e53edb 100644 --- a/src/WaveToy.F90 +++ b/src/WaveToy.F90 @@ -135,22 +135,22 @@ subroutine WaveToyFreeF90_Boundaries(CCTK_ARGUMENTS) ! be passed if (CCTK_EQUALS(bound,"flat")) then - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, -1, & + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, & "wavetoy::phi", "Flat"); else if (CCTK_EQUALS(bound,"static")) then - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, -1, & + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, & "wavetoy::phi", "Static"); else if (CCTK_EQUALS(bound,"radiation")) then - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, -1, & + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, & "wavetoy::phi", "Radiation"); else if (CCTK_EQUALS(bound,"robin")) then - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, -1, & + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, & "wavetoy::phi", "Robin"); else if (CCTK_EQUALS(bound,"zero")) then - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, -1, & + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, & "wavetoy::phi", "Scalar"); else if (CCTK_EQUALS(bound,"none")) then - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, -1, & + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, & "wavetoy::phi", "None"); else call CCTK_WARN(0,"Unrecognized boundary condition") |