From 5b29a4033d38fd82482d3cbe0f96b95c3d679e05 Mon Sep 17 00:00:00 2001 From: allen Date: Sat, 18 Sep 1999 17:50:46 +0000 Subject: Fixing radiation boundary conditions git-svn-id: http://svn.cactuscode.org/arrangements/CactusWave/WaveToyF77/trunk@25 4451c3c6-1034-4891-99ea-21147727ccdf --- src/WaveToy.F77 | 9 ++++++--- 1 file changed, 6 insertions(+), 3 deletions(-) (limited to 'src') diff --git a/src/WaveToy.F77 b/src/WaveToy.F77 index c6036a0..e0a9453 100644 --- a/src/WaveToy.F77 +++ b/src/WaveToy.F77 @@ -39,7 +39,7 @@ c Using Cactus arguments lists c Declare variables in argument list DECLARE_CCTK_FARGUMENTS - INTEGER i,j,k + INTEGER i,j,k,ierr INTEGER istart, jstart, kstart, iend, jend, kend CCTK_REAL dx,dy,dz,dt @@ -50,6 +50,8 @@ c ----------------- dz = CCTK_DELTA_SPACE(3) dt = CCTK_DELTA_TIME + test = cctk_time + call CCTK_OutputVarAsByMethod(ierr,cctkGH,"WaveToyF77::test","IOScalar","Gab") istart = 2 jstart = 2 kstart = 2 @@ -126,7 +128,7 @@ c Declare parameters DECLARE_CCTK_PARAMETERS c Local declarations - CCTK_REAL zero + CCTK_REAL zero,one integer ierr integer sw(3) @@ -134,6 +136,7 @@ c Cactus declarations integer CCTK_Equals zero = 0.0 + one = 1.0 c Set the stencil width c --------------------- @@ -152,7 +155,7 @@ c ----------------------------------- else if (CCTK_EQUALS(bound,"zero")) then call ApplyConstantBC(ierr,cctkGH,zero,sw,"wavetoyf77::phi_next") else if (CCTK_Equals(bound,"radiation").eq.1) then - call ApplyRadiativeBC(ierr,cctkGH,zero,zero,sw,"wavetoyf77::phi_next", + call ApplyRadiativeBC(ierr,cctkGH,zero,one,sw,"wavetoyf77::phi_next", & "wavetoy::phi") end if -- cgit v1.2.3