diff options
| author | allen <allen@f80f6fb6-8356-4fd4-90bc-d84ad503c100> | 2001-01-17 15:30:18 +0000 |
|---|---|---|
| committer | allen <allen@f80f6fb6-8356-4fd4-90bc-d84ad503c100> | 2001-01-17 15:30:18 +0000 |
| commit | 6dc7e7259e1283f9ccd52cd81953439a71398005 (patch) | |
| tree | 597c10ad9a8c69b0fff45adc874d9ddbdde3073f /src | |
| parent | df15ceb404a0222c1ed8935efcaf972c2e394ed8 (diff) | |
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
Diffstat (limited to 'src')
| -rw-r--r-- | src/WaveToy.F | 21 |
1 files changed, 16 insertions, 5 deletions
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 |