diff options
| author | schnetter <schnetter@4ec1db94-0e4f-0410-ada3-8bed251432c9> | 2004-07-06 16:24:46 +0000 |
|---|---|---|
| committer | schnetter <schnetter@4ec1db94-0e4f-0410-ada3-8bed251432c9> | 2004-07-06 16:24:46 +0000 |
| commit | 5d921c7d6457f0392fec9ecee61b7e4372937cf3 (patch) | |
| tree | fb3f1d0a9f69f784b29b35ce87893a20b393f6bd /src | |
| parent | 16725b4c9678c22bf1c23b9c40970a0798735a80 (diff) | |
Give people more choice for what kind of Minkowski they want.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/NoExcision/trunk@3 4ec1db94-0e4f-0410-ada3-8bed251432c9
Diffstat (limited to 'src')
| -rw-r--r-- | src/overwrite.F90 | 14 |
1 files changed, 7 insertions, 7 deletions
diff --git a/src/overwrite.F90 b/src/overwrite.F90 index c7949ae..87cd534 100644 --- a/src/overwrite.F90 +++ b/src/overwrite.F90 @@ -21,7 +21,7 @@ subroutine NoExcision_Overwrite (CCTK_ARGUMENTS) cz = centre_z(n) rad = radius(n) - if (overwrite_geometry /= 0) then + if (overwrite_geometry(n) /= 0) then if (conformal_state >= 1) then where ((x - cx)**2 + (y - cy)**2 + (z - cz)**2 <= rad**2) @@ -47,12 +47,12 @@ subroutine NoExcision_Overwrite (CCTK_ARGUMENTS) end if where ((x - cx)**2 + (y - cy)**2 + (z - cz)**2 <= rad**2) - gxx = 1 + gxx = Minkowski_scale(n) gxy = 0 gxz = 0 - gyy = 1 + gyy = Minkowski_scale(n) gyz = 0 - gzz = 1 + gzz = Minkowski_scale(n) kxx = 0 kxy = 0 kxz = 0 @@ -63,15 +63,15 @@ subroutine NoExcision_Overwrite (CCTK_ARGUMENTS) end if - if (overwrite_lapse /= 0) then + if (overwrite_lapse(n) /= 0) then where ((x - cx)**2 + (y - cy)**2 + (z - cz)**2 <= rad**2) - alp = 1 + alp = lapse_scale(n) end where end if - if (overwrite_shift /= 0) then + if (overwrite_shift(n) /= 0) then if (shift_state /= 0) then where ((x - cx)**2 + (y - cy)**2 + (z - cz)**2 <= rad**2) |