From 4b95038531a96c4f7dc48b25cf44287e90f505cf Mon Sep 17 00:00:00 2001 From: yye00 Date: Fri, 29 Apr 2005 21:25:30 +0000 Subject: more \r to \n conversion git-svn-id: http://svn.cactuscode.org/arrangements/CactusWave/WaveToyFreeF90/trunk@76 daab82bb-f315-4ad1-b6d0-9353ff8b6e27 --- src/InitSymBound.F90 | 39 -------------- src/Startup.F90 | 23 -------- src/WaveToy.F90 | 148 --------------------------------------------------- 3 files changed, 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 -- cgit v1.2.3