diff options
| author | schnetter <schnetter@9ea3c693-0e4f-0410-9b72-c04bdc00281e> | 2006-08-14 14:24:36 +0000 |
|---|---|---|
| committer | schnetter <schnetter@9ea3c693-0e4f-0410-9b72-c04bdc00281e> | 2006-08-14 14:24:36 +0000 |
| commit | 8f5ca9fb6e3705e0c2710eaefb8d99d2df916df8 (patch) | |
| tree | 19d0eb52c3777ff0250e71fc1e062dcc86814143 | |
| parent | b7361223e356c0a5f6622e190e06998860b6eb9c (diff) | |
Remove unused source file. I assume that Initial.F was renamed to
Initial.F77, but Initial.F was for some reason not removed. At the
moment, Initial.F seems to be an old (and incorrect) version of
Initial.F77.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDConstraintViolate/trunk@10 9ea3c693-0e4f-0410-9b72-c04bdc00281e
| -rw-r--r-- | src/Initial.F | 71 |
1 files changed, 0 insertions, 71 deletions
diff --git a/src/Initial.F b/src/Initial.F deleted file mode 100644 index 19a1eb4..0000000 --- a/src/Initial.F +++ /dev/null @@ -1,71 +0,0 @@ - /*@@ - @file Initial.F - @date May 2002 - @author Gabrielle Allen - @desc - Constraint violating initial data - @enddesc - @version $Header$ - @@*/ - -#include "cctk.h" -#include "cctk_Arguments.h" -#include "cctk_Parameters.h" - -/******************************************************************** - ********************* External Routines ********************** - ********************************************************************/ - - /*@@ - @routine IDConstraintViolate_Initial - @date May 2002 - @author Gabrielle Allen - @desc - Set up initial data which violates the constraints - @enddesc - @calls -@@*/ - - subroutine IDConstraintViolate_Initial(CCTK_ARGUMENTS) - - implicit none - - integer i,j,k - CCTK_REAL xp,yp,zp,rp - - DECLARE_CCTK_ARGUMENTS - DECLARE_CCTK_PARAMETERS - - do k=1, cctk_lsh(3) - do j=1, cctk_lsh(2) - do i=1, cctk_lsh(1) - - xp=x(i,j,k) - yp=y(i,j,k) - zp=z(i,j,k) - - rp = sqrt(xp*xp+yp*yp+zp*zp) - - gxx(i,j,k) = amplitude*exp(-(rp-radius)**2/sigma**2) - gxy(i,j,k) = amplitude*exp(-(rp-radius)**2/sigma**2) - gxz(i,j,k) = amplitude*exp(-(rp-radius)**2/sigma**2) - gyy(i,j,k) = amplitude*exp(-(rp-radius)**2/sigma**2) - gyz(i,j,k) = amplitude*exp(-(rp-radius)**2/sigma**2) - gzz(i,j,k) = amplitude*exp(-(rp-radius)**2/sigma**2) - - kxx(i,j,k) = 0.0d0 - kxy(i,j,k) = 0.0d0 - kxz(i,j,k) = 0.0d0 - kyy(i,j,k) = 0.0d0 - kyz(i,j,k) = 0.0d0 - kzz(i,j,k) = 0.0d0 - - end do - end do - end do - - end - -/******************************************************************** - ********************* Local Routines ************************* - ********************************************************************/ |