diff options
| author | allen <allen@4451c3c6-1034-4891-99ea-21147727ccdf> | 1999-09-18 17:50:46 +0000 |
|---|---|---|
| committer | allen <allen@4451c3c6-1034-4891-99ea-21147727ccdf> | 1999-09-18 17:50:46 +0000 |
| commit | 5b29a4033d38fd82482d3cbe0f96b95c3d679e05 (patch) | |
| tree | c2f2ee53c4c2e678f0cec695bd21b370ef38caf5 /src | |
| parent | 327bd7574b3b1eb30b9ae90e620c43a91a4c2506 (diff) | |
Fixing radiation boundary conditions
git-svn-id: http://svn.cactuscode.org/arrangements/CactusWave/WaveToyF77/trunk@25 4451c3c6-1034-4891-99ea-21147727ccdf
Diffstat (limited to 'src')
| -rw-r--r-- | src/WaveToy.F77 | 9 |
1 files changed, 6 insertions, 3 deletions
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 |