From 6dc7e7259e1283f9ccd52cd81953439a71398005 Mon Sep 17 00:00:00 2001 From: allen Date: Wed, 17 Jan 2001 15:30:18 +0000 Subject: Making the call to do zero boundary conditions use new directional boundary conditions, and new testsuite to test it. git-svn-id: http://svn.cactuscode.org/arrangements/CactusWave/WaveToyF90/trunk@55 f80f6fb6-8356-4fd4-90bc-d84ad503c100 --- src/WaveToy.F | 21 ++++++++++++++++----- 1 file changed, 16 insertions(+), 5 deletions(-) (limited to 'src/WaveToy.F') diff --git a/src/WaveToy.F b/src/WaveToy.F index bd434d8..0ae3478 100644 --- a/src/WaveToy.F +++ b/src/WaveToy.F @@ -111,17 +111,28 @@ integer i integer :: ierr=-1 integer,dimension(3):: sw=1 - CCTK_REAL,parameter :: zero = 0.0 - CCTK_REAL,parameter :: one = 1.0 - + CCTK_REAL,parameter :: sval = 0.0 + CCTK_REAL,parameter :: rzero = 0.0 + CCTK_REAL,parameter :: rone = 1.0 + integer,parameter :: swdir = 1 + integer,parameter :: ione = 1,mione=-1 + integer,parameter :: itwo = 2,mitwo=-2 + integer,parameter :: ithree = 3,mithree=-3 + integer,parameter :: ifour = 4,izero = 0 call CartSymGN(ierr,cctkGH,"wavetoy::scalarevolve") 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") +# call BndScalarVN(ierr,cctkGH,sw,sval,"wavetoy::phi") + call BndScalarDirVN(ierr,cctkGH,swdir,mione ,sval,"wavetoy::phi") + call BndScalarDirVN(ierr,cctkGH,swdir,ione ,sval,"wavetoy::phi") + call BndScalarDirVN(ierr,cctkGH,swdir,mitwo ,sval,"wavetoy::phi") + call BndScalarDirVN(ierr,cctkGH,swdir,itwo ,sval,"wavetoy::phi") + call BndScalarDirVN(ierr,cctkGH,swdir,mithree,sval,"wavetoy::phi") + call BndScalarDirVN(ierr,cctkGH,swdir,ithree ,sval,"wavetoy::phi") else if (CCTK_EQUALS(bound,"radiation")) then - call BndRadiativeVN(ierr,cctkGH,sw,zero,one, + call BndRadiativeVN(ierr,cctkGH,sw,rzero,rone, & "wavetoy::phi","wavetoy::phi") end if -- cgit v1.2.3