diff options
Diffstat (limited to 'src/WaveToy.F')
| -rw-r--r-- | src/WaveToy.F | 62 |
1 files changed, 33 insertions, 29 deletions
diff --git a/src/WaveToy.F b/src/WaveToy.F index 48189c7..05d847e 100644 --- a/src/WaveToy.F +++ b/src/WaveToy.F @@ -1,5 +1,5 @@ /*@@ - @file WaveToyF90.F + @file WaveToy.F @date @author Tom Goodale @desc @@ -7,7 +7,8 @@ @enddesc @@*/ -#include "cctk.h" +#include "cctk.h" +#include "cctk_Faces.h" #include "cctk_Parameters.h" #include "cctk_Arguments.h" @@ -128,30 +129,43 @@ finf = 1.0d0 npow = 1 +c Apply the symmetry boundary conditions on any coordinate axes +c ------------------------------------------------------------- call CartSymGN(ierr,cctkGH,"wavetoy::scalarevolve") +c Apply the outer boundary conditions +c ----------------------------------- +c Note: In each of the following calls to Boundary_SelectVarForBC, +c default arguments are used, so an invalid table handle of -1 can +c 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,rzero,rone,"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 BndScalarDirVN(ierr,cctkGH,swdir,mione ,sval,"wavetoy::phi") - BOUNDARY_ERROR - call BndScalarDirVN(ierr,cctkGH,swdir,ione ,sval,"wavetoy::phi") - BOUNDARY_ERROR - call BndScalarDirVN(ierr,cctkGH,swdir,mitwo ,sval,"wavetoy::phi") - BOUNDARY_ERROR - call BndScalarDirVN(ierr,cctkGH,swdir,itwo ,sval,"wavetoy::phi") - BOUNDARY_ERROR - call BndScalarDirVN(ierr,cctkGH,swdir,mithree,sval,"wavetoy::phi") - BOUNDARY_ERROR - call BndScalarDirVN(ierr,cctkGH,swdir,ithree ,sval,"wavetoy::phi") - BOUNDARY_ERROR +c Face specific calls are not working yet with the new boundary interface: +c call BndScalarDirVN(ierr,cctkGH,swdir,mione ,sval,"wavetoy::phi") +c BOUNDARY_ERROR +c call BndScalarDirVN(ierr,cctkGH,swdir,ione ,sval,"wavetoy::phi") +c BOUNDARY_ERROR +c call BndScalarDirVN(ierr,cctkGH,swdir,mitwo ,sval,"wavetoy::phi") +c BOUNDARY_ERROR +c call BndScalarDirVN(ierr,cctkGH,swdir,itwo ,sval,"wavetoy::phi") +c BOUNDARY_ERROR +c call BndScalarDirVN(ierr,cctkGH,swdir,mithree,sval,"wavetoy::phi") +c BOUNDARY_ERROR +c call BndScalarDirVN(ierr,cctkGH,swdir,ithree ,sval,"wavetoy::phi") +c BOUNDARY_ERROR + call Boundary_SelectVarForBC(ierr, cctkGH, CCTK_ALL_FACES, -1, + $ "wavetoy::phi", "Scalar"); else if (.NOT. CCTK_EQUALS(bound,"none")) then call CCTK_WARN(0,"Unrecognized boundary condition") end if @@ -161,13 +175,3 @@ end if end subroutine WaveToyF90_Boundaries - - - - - - - - - - |