diff options
Diffstat (limited to 'src/WaveToy.F')
| -rw-r--r-- | src/WaveToy.F | 39 |
1 files changed, 20 insertions, 19 deletions
diff --git a/src/WaveToy.F b/src/WaveToy.F index eb9bc0f..1c9bc8f 100644 --- a/src/WaveToy.F +++ b/src/WaveToy.F @@ -32,20 +32,20 @@ @@*/ - subroutine wavetoyf90_evolution(CCTK_FARGUMENTS) + subroutine WaveToyf90_Evolution(CCTK_FARGUMENTS) implicit none DECLARE_CCTK_FARGUMENTS DECLARE_CCTK_PARAMETERS -c Declare local variables +! Declare local variables INTEGER :: i,j,k INTEGER :: istart, jstart, kstart, iend, jend, kend CCTK_REAL :: dx,dy,dz,dt -c Set up shorthands -c ----------------- +! Set up shorthands +! ----------------- dx = CCTK_DELTA_SPACE(1) dy = CCTK_DELTA_SPACE(2) dz = CCTK_DELTA_SPACE(3) @@ -59,8 +59,8 @@ c ----------------- jend = cctk_lsh(2)-1 kend = cctk_lsh(3)-1 -c Do the evolution -c ---------------- +! Do the evolution +! ---------------- do k = kstart, kend do j = jstart, jend do i = istart, iend @@ -77,20 +77,21 @@ c ---------------- end do end do +! Synchronize +! ----------- + call CCTK_SyncGroup(cctkGH,"wavetoyf90::temps") + +! Apply boundary conditions +! ------------------------- + call WaveToyF90_Boundaries(CCTK_PASS_FTOF) + ! Update timeslices ! ----------------- phi_old = phi phi = phi_tmp -! Apply boundary conditions -! ------------------------- - call WaveToyF90_Boundaries(CCTK_FARGUMENTS) - -! Synchronize -! ----------- - call CCTK_SyncGroup(cctkGH,"wavetoy::scalarevolve") - end subroutine wavetoyf90_evolution + end subroutine WaveToyF90_Evolution /*@@ @@ -121,20 +122,20 @@ c ---------------- CCTK_REAL,parameter :: zero = 0.0 CCTK_REAL,parameter :: one = 1.0 - call ApplySymmetry(cctkGH,"wavetoy::scalarevolve") + call ApplySymmetry(cctkGH,"wavetoyf90::temps") if (CCTK_EQUALS(bound,"flat")) then - call ApplyFlatBC(ierr,cctkGH,sw,"wavetoy::phi") + call ApplyFlatBC(ierr,cctkGH,sw,"wavetoyf90::phi_tmp") else if (CCTK_EQUALS(bound,"radiation")) then call ApplyRadiativeBC(ierr,cctkGH,zero,one,sw, - & "wavetoy::phi","wavetoy::phi_old") + & "wavetoyf90::phi_tmp","wavetoy::phi") end if if (ierr < 0) then - call CCTK_WARN(0,"Boundary conditions not applied - giving up!"); + call CCTK_WARN(0,"Boundary conditions not applied - giving up!") end if - end subroutine wavetoyf90_boundaries + end subroutine WaveToyF90_Boundaries |