diff options
-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 ************************* - ********************************************************************/ |