diff options
author | shawley <> | 2002-02-27 17:54:00 +0000 |
---|---|---|
committer | shawley <> | 2002-02-27 17:54:00 +0000 |
commit | 3d2984c41611897e009938b476a5d6e5dca73ab3 (patch) | |
tree | ef3812d5d2cb1db51155db6837e31742d31edeeb /CarpetExtra/IDScalarWave | |
parent | 45f371b572c2a08c3d84518e9f536e3dfb5afbc7 (diff) |
added cx, cy & cz for 1/r initial data
darcs-hash:20020227175413-e415b-024c510207e52381e1fee261a973bfa86c4963db.gz
Diffstat (limited to 'CarpetExtra/IDScalarWave')
-rw-r--r-- | CarpetExtra/IDScalarWave/src/InitialData.F77 | 12 |
1 files changed, 8 insertions, 4 deletions
diff --git a/CarpetExtra/IDScalarWave/src/InitialData.F77 b/CarpetExtra/IDScalarWave/src/InitialData.F77 index 5f330c924..61576828f 100644 --- a/CarpetExtra/IDScalarWave/src/InitialData.F77 +++ b/CarpetExtra/IDScalarWave/src/InitialData.F77 @@ -1,5 +1,5 @@ c -*-Fortran-*- -c $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/IDScalarWave/src/InitialData.F77,v 1.9 2002/01/08 14:20:57 schnetter Exp $ +c $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/IDScalarWave/src/InitialData.F77,v 1.10 2002/02/27 18:54:13 shawley Exp $ /*@@ @file InitialData.F77 @@ -138,9 +138,13 @@ c Use kx,ky,kz as number of modes in each direction. do j=1,cctk_lsh(2) do i=1,cctk_lsh(1) - phi(i,j,k) = 1 / spher3d_r(i,j,k) - phi_p(i,j,k) = 1 / spher3d_r(i,j,k) - phi_p_p(i,j,k) = 1 / spher3d_r(i,j,k) + x = cart3d_x(i,j,k) + y = cart3d_y(i,j,k) + z = cart3d_z(i,j,k) + + phi(i,j,k) = 1 / sqrt((x-cx)**2 + (y-cy)**2 + (z-cz)**2) + phi_p(i,j,k) = 1 / phi(i,j,k) + phi_p_p(i,j,k) = 1 / phi(i,j,k) end do end do |