diff options
Diffstat (limited to 'src/WaveToy.F')
| -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 |