diff options
Diffstat (limited to 'src/WaveToy.F90')
-rw-r--r-- | src/WaveToy.F90 | 38 |
1 files changed, 27 insertions, 11 deletions
diff --git a/src/WaveToy.F90 b/src/WaveToy.F90 index f82aed7..9b0d9b2 100644 --- a/src/WaveToy.F90 +++ b/src/WaveToy.F90 @@ -1,5 +1,5 @@ /*@@ - @file WaveToyFortran.F90 + @file WaveToy.F90 @date @author Tom Goodale @desc @@ -7,13 +7,14 @@ @enddesc @@*/ -#include "cctk.h" +#include "cctk.h" +#include "cctk_Faces.h" #include "cctk_Parameters.h" #include "cctk_Arguments.h" /*@@ - @routine WaveToyFortran_Evolution + @routine WaveToyFreeF90_Evolution @date @author Tom Goodale @desc @@ -122,21 +123,36 @@ subroutine WaveToyFreeF90_Boundaries(CCTK_ARGUMENTS) finf = 1.0d0 npow = 1 +! Apply the symmetry boundary conditions on any coordinate axes +! ------------------------------------------------------------- call CartSymGN(ierr,cctkGH,"wavetoy::scalarevolve") +! Apply the outer boundary conditions +! ----------------------------------- +! Note: In each of the following calls to Boundary_SelectVarForBC, +! default arguments are used, so an invalid table handle of -1 can +! be passed + if (CCTK_EQUALS(bound,"flat")) then - call BndFlatVN(ierr,cctkGH,sw,"wavetoy::phi") + call Boundary_SelectVarForBC(ierr, cctkGH, CCTK_ALL_FACES, -1, & + "wavetoy::phi", "Flat"); else if (CCTK_EQUALS(bound,"static")) then - call BndStaticVN(ierr,cctkGH,sw,"wavetoy::phi") + call Boundary_SelectVarForBC(ierr, cctkGH, CCTK_ALL_FACES, -1, & + "wavetoy::phi", "Static"); else if (CCTK_EQUALS(bound,"radiation")) then - call BndRadiativeVN(ierr,cctkGH,sw,zero,one,"wavetoy::phi", & - "wavetoy::phi") + call Boundary_SelectVarForBC(ierr, cctkGH, CCTK_ALL_FACES, -1, & + "wavetoy::phi", "Radiative"); else if (CCTK_EQUALS(bound,"robin")) then - call BndRobinVN(ierr,cctkGH, sw, finf, npow,"wavetoy::phi") + call Boundary_SelectVarForBC(ierr, cctkGH, CCTK_ALL_FACES, -1, & + "wavetoy::phi", "Robin"); else if (CCTK_EQUALS(bound,"zero")) then - call BndScalarVN(ierr,cctkGH,zero,sw,"wavetoy::phi") - else if (.NOT. CCTK_EQUALS(bound,"none")) then - call CCTK_WARN(0,"Unrecognized boundary condition") + call Boundary_SelectVarForBC(ierr, cctkGH, CCTK_ALL_FACES, -1, & + "wavetoy::phi", "Scalar"); + else if (CCTK_EQUALS(bound,"none")) then + call Boundary_SelectVarForBC(ierr, cctkGH, CCTK_ALL_FACES, -1, & + "wavetoy::phi", "None"); + else + call CCTK_WARN(0,"Unrecognized boundary condition") end if if (ierr < 0) then |