From 4be592803b4c4941d42fa4ea61bc4b5c7d60e758 Mon Sep 17 00:00:00 2001 From: allen Date: Fri, 31 May 2002 10:04:56 +0000 Subject: Standardizing all boundary conditions across wavetoy thorns git-svn-id: http://svn.cactuscode.org/arrangements/CactusWave/WaveToyF77/trunk@86 4451c3c6-1034-4891-99ea-21147727ccdf --- param.ccl | 1 + src/WaveToy.F77 | 12 ++++++++---- 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 -- cgit v1.2.3