aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authoryye00 <yye00@daab82bb-f315-4ad1-b6d0-9353ff8b6e27>2005-04-29 21:25:30 +0000
committeryye00 <yye00@daab82bb-f315-4ad1-b6d0-9353ff8b6e27>2005-04-29 21:25:30 +0000
commit4b95038531a96c4f7dc48b25cf44287e90f505cf (patch)
tree6fc1fdb9bb33b00990c222ac1dcaab0d77718e39
parent63a9307937610ee7cbe380199f906fc67ff95765 (diff)
more \r to \n conversion
git-svn-id: http://svn.cactuscode.org/arrangements/CactusWave/WaveToyFreeF90/trunk@76 daab82bb-f315-4ad1-b6d0-9353ff8b6e27
-rw-r--r--src/InitSymBound.F9039
-rw-r--r--src/Startup.F9023
-rw-r--r--src/WaveToy.F90148
3 files changed, 0 insertions, 210 deletions
diff --git a/src/InitSymBound.F90 b/src/InitSymBound.F90
index 8822a73..e69de29 100644
--- a/src/InitSymBound.F90
+++ b/src/InitSymBound.F90
@@ -1,39 +0,0 @@
- /*@@
- @file InitSymBound.F90
- @date
- @author Gabrielle Allen
- @desc
- Sets the symmetries across the coordinate axes
- @enddesc
- @@*/
-
-#include "cctk.h"
-#include "cctk_Arguments.h"
-
- /*@@
- @routine WaveToyFreeF90_InitSymBound
- @date
- @author Gabrielle Allen
- @desc
- Sets the symmetries for Wave Toy
- @enddesc
- @calls
- @calledby
- @history
-
- @endhistory
-
-@@*/
-
-subroutine WaveToyFreeF90_InitSymBound(CCTK_ARGUMENTS)
-
- implicit none
-
- DECLARE_CCTK_ARGUMENTS
-
- INTEGER :: ierr
- INTEGER, DIMENSION(3) :: sym = 1
-
- call SetCartSymVN(ierr, cctkGH, sym, 'wavetoy::phi')
-
-end subroutine WaveToyFreeF90_InitSymbound
diff --git a/src/Startup.F90 b/src/Startup.F90
index d00270f..e69de29 100644
--- a/src/Startup.F90
+++ b/src/Startup.F90
@@ -1,23 +0,0 @@
- /*@@
- @file Startup.F90
- @date
- @author Gabrielle Allen
- @desc
- Register banner
- @enddesc
- @version $Header$
- @@*/
-
-#include "cctk.h"
-
-integer function WaveToyFreeF90_Startup ()
-
- implicit none
-
- integer ierr
-
- call CCTK_RegisterBanner(ierr, "WaveToyFreeF90: Evolutions of a Scalar Field")
-
- WaveToyFreeF90_Startup = 0
-
-end function WaveToyFreeF90_Startup
diff --git a/src/WaveToy.F90 b/src/WaveToy.F90
index 41a650c..e69de29 100644
--- a/src/WaveToy.F90
+++ b/src/WaveToy.F90
@@ -1,148 +0,0 @@
- /*@@
- @file WaveToy.F90
- @date
- @author Tom Goodale
- @desc
- Evolution routines for the wave equation solver
- @enddesc
- @version $Header$
- @@*/
-
-#include "cctk.h"
-#include "cctk_Parameters.h"
-#include "cctk_Arguments.h"
-#include "cctk_Functions.h"
-
-
- /*@@
- @routine WaveToyFreeF90_Evolution
- @date
- @author Tom Goodale
- @desc
- Evolution for the wave equation
- @enddesc
- @calls CCTK_SyncGroup, wavetoy_boundaries
- @calledby
- @history
-
- @endhistory
-
-@@*/
-
-subroutine WaveToyFreeF90_Evolution(CCTK_ARGUMENTS)
-
- implicit none
-
- DECLARE_CCTK_ARGUMENTS
- DECLARE_CCTK_PARAMETERS
-
- ! Declare local variables
- INTEGER :: i,j,k
- INTEGER :: istart, jstart, kstart, iend, jend, kend
-
- CCTK_REAL :: dx,dy,dz,dt
- CCTK_REAL :: dx2,dy2,dz2,dt2
- CCTK_REAL :: dx2i,dy2i,dz2i
-
- CCTK_REAL :: factor
-
- ! Set up shorthands
- ! -----------------
- dx = CCTK_DELTA_SPACE(1)
- dy = CCTK_DELTA_SPACE(2)
- dz = CCTK_DELTA_SPACE(3)
- dt = CCTK_DELTA_TIME
-
- dx2 = dx*dx
- dy2 = dy*dy
- dz2 = dz*dz
- dt2 = dt*dt
-
- dx2i = 1.0/dx2
- dy2i = 1.0/dy2
- dz2i = 1.0/dz2
-
- istart = 2
- jstart = 2
- kstart = 2
-
- iend = cctk_lsh(1)-1
- jend = cctk_lsh(2)-1
- kend = cctk_lsh(3)-1
-
- factor = 2*(1 - (dt2)*(dx2i + dy2i + dz2i))
-
- ! Do the evolution
- ! ----------------
- do k = kstart, kend
- do j = jstart, jend
- do i = istart, iend
-
- phi(i,j,k) = factor*phi_p(i,j,k) - &
- phi_p_p(i,j,k) + (dt2) * &
- ((phi_p(i+1,j,k)+phi_p(i-1,j,k))*dx2i &
- +(phi_p(i,j+1,k)+phi_p(i,j-1,k))*dy2i &
- +(phi_p(i,j,k+1)+phi_p(i,j,k-1))*dz2i)
-
- end do
- end do
- end do
-
-end subroutine WaveToyFreeF90_Evolution
-
-
- /*@@
- @routine WaveToyFortran_Boundaries
- @date
- @author Tom Goodale
- @desc
- Boundary conditions for the wave equation
- @enddesc
- @calls CartSymGH,FlatBC,RadiativeBC
- @calledby
- @history
-
- @endhistory
-
-@@*/
-
-subroutine WaveToyFreeF90_Boundaries(CCTK_ARGUMENTS)
-
- implicit none
-
- DECLARE_CCTK_ARGUMENTS
- DECLARE_CCTK_PARAMETERS
- DECLARE_CCTK_FUNCTIONS
-
-! Local declarations
- CCTK_INT :: ierr
-
- CHARACTER (len=100) :: boundary
- INTEGER :: length
-
- ierr = 0
-
-! The "bound" parameter needs to be converted into a Fortran string.
- call CCTK_FortranString(length,bound,boundary)
-
-! Apply the outer boundary conditions
-! -----------------------------------
-! Note: In each of the following calls to Boundary_SelectVarForBC,
-! default arguments are used, so an invalid table handle of -1 can
-! be passed
-
- if (CCTK_EQUALS(bound,"flat") .or. CCTK_EQUALS(bound,"static") .or. &
- CCTK_EQUALS(bound,"radiation") .or. CCTK_EQUALS(bound,"robin") .or. &
- CCTK_EQUALS(bound,"none") ) then
- ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, &
- "wavetoy::phi", boundary)
- else if (CCTK_EQUALS(bound,"zero")) then
- ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, &
- "wavetoy::phi", "Scalar")
- end if
-
- if (ierr < 0) then
- call CCTK_WARN(0,"WaveToyFreeF90_Boundaries: Error selecting boundary condition")
- end if
-
-end subroutine WaveToyFreeF90_Boundaries