diff options
| -rw-r--r-- | param.ccl | 1 | ||||
| -rw-r--r-- | src/WaveToy.F77 | 12 |
2 files changed, 9 insertions, 4 deletions
@@ -7,6 +7,7 @@ KEYWORD bound "Type of boundary condition to use" { "none" :: "No boundary condition" "flat" :: "Flat boundary condition" + "static" :: "Static boundary condition" "radiation" :: "Radiation boundary condition" "robin" :: "Robin boundary condition" "zero" :: "Zero boundary condition" diff --git a/src/WaveToy.F77 b/src/WaveToy.F77 index ea38c8f..4b287bc 100644 --- a/src/WaveToy.F77 +++ b/src/WaveToy.F77 @@ -137,13 +137,17 @@ c Apply the outer boundary conditions c ----------------------------------- if (CCTK_EQUALS(bound,"flat")) then call BndFlatVN(ierr,cctkGH,sw,"wavetoy::phi") - else if (CCTK_EQUALS(bound,"zero")) then - call BndScalarVN(ierr,cctkGH,zero,sw,"wavetoy::phi") - else if (CCTK_Equals(bound,"radiation").eq.1) then + else if (CCTK_EQUALS(bound,"static")) then + call BndStaticVN(ierr,cctkGH,sw,"wavetoy::phi") + else if (CCTK_EQUALS(bound,"radiation")) then call BndRadiativeVN(ierr,cctkGH,sw,zero,one,"wavetoy::phi", & "wavetoy::phi") - else if (CCTK_Equals(bound,"robin").eq.1) then + else if (CCTK_EQUALS(bound,"robin")) then call BndRobinVN(ierr,cctkGH, sw, finf, npow,"wavetoy::phi") + 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") end if if (ierr < 0) then |