diff options
author | allen <allen@5c0f84ea-6048-4d6e-bfa4-55cd5f2e0dd7> | 1999-06-25 17:09:04 +0000 |
---|---|---|
committer | allen <allen@5c0f84ea-6048-4d6e-bfa4-55cd5f2e0dd7> | 1999-06-25 17:09:04 +0000 |
commit | 40c56ceb52068dbe831fa9723fc32e6957ec1154 (patch) | |
tree | 4f6592d767f1e6f41873b76b9cf81ae697f85319 /src/planewaves.F | |
parent | 02288ab35f484eebd60602e6956a49b63c302213 (diff) |
Implementing new naming convention for variables on the cGH. See cGH.h for
details.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDLinearWaves/trunk@16 5c0f84ea-6048-4d6e-bfa4-55cd5f2e0dd7
Diffstat (limited to 'src/planewaves.F')
-rw-r--r-- | src/planewaves.F | 16 |
1 files changed, 8 insertions, 8 deletions
diff --git a/src/planewaves.F b/src/planewaves.F index ca9ed34..83221e4 100644 --- a/src/planewaves.F +++ b/src/planewaves.F @@ -38,7 +38,7 @@ #include "declare_arguments.h" #include "declare_parameters.h" c Using macro definitions from Einstein -#include "../../packages/CactusEinstein/Einstein/src/Einstein.h" +#include "CactusEinstein/Einstein/src/Einstein.h" subroutine planewaves(CCTK_FARGUMENTS) @@ -112,14 +112,14 @@ c precalc wp = (kxp*kxp+kyp*kyp+kzp*kzp)**0.5 c *************** plane waves ******************** - do k=1,sh(3) - do j=1,sh(2) - do i=1,sh(1) - plus = (kx*x(i,j,k)+ky*y(i,j,k)+kz*z(i,j,k)+w*time) - minus = (kx*x(i,j,k)+ky*y(i,j,k)+kz*z(i,j,k)-w*time) + do k=1,cctk_lsh(3) + do j=1,cctk_lsh(2) + do i=1,cctk_lsh(1) + plus = (kx*x(i,j,k)+ky*y(i,j,k)+kz*z(i,j,k)+w*cctk_time) + minus = (kx*x(i,j,k)+ky*y(i,j,k)+kz*z(i,j,k)-w*cctk_time) - plusp =(kxp*x(i,j,k)+kyp*y(i,j,k)+kzp*z(i,j,k)+wp*(time-ra)) - minusp =(kxp*x(i,j,k)+kyp*y(i,j,k)+kzp*z(i,j,k)-wp*(time-ra)) + plusp =(kxp*x(i,j,k)+kyp*y(i,j,k)+kzp*z(i,j,k)+wp*(cctk_time-ra)) + minusp =(kxp*x(i,j,k)+kyp*y(i,j,k)+kzp*z(i,j,k)-wp*(cctk_time-ra)) ain = iin*amplitude*cos(plus)*exp(-(plusp)**2) aout = iout*amplitude*cos(minus)*exp(-(minusp)**2) |