aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorallen <allen@4451c3c6-1034-4891-99ea-21147727ccdf>2002-05-31 10:04:56 +0000
committerallen <allen@4451c3c6-1034-4891-99ea-21147727ccdf>2002-05-31 10:04:56 +0000
commit4be592803b4c4941d42fa4ea61bc4b5c7d60e758 (patch)
tree058fee3380880a0bac8d09098d5e265d5742917f
parente1d5b02108a3eb785e3dce737f56235fd9e5c168 (diff)
Standardizing all boundary conditions across wavetoy thorns
git-svn-id: http://svn.cactuscode.org/arrangements/CactusWave/WaveToyF77/trunk@86 4451c3c6-1034-4891-99ea-21147727ccdf
-rw-r--r--param.ccl1
-rw-r--r--src/WaveToy.F7712
2 files changed, 9 insertions, 4 deletions
diff --git a/param.ccl b/param.ccl
index 4d915b8..dd0f4d3 100644
--- a/param.ccl
+++ b/param.ccl
@@ -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