diff options
| author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-01-14 14:23:14 +0000 |
|---|---|---|
| committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-01-14 14:23:14 +0000 |
| commit | 7fba5304ce3e983755fae6fde3e7d9aec49da881 (patch) | |
| tree | c9ec690a5f3af2fc1a555f22cc2df69069d88ae4 /src | |
| parent | 9d1f0811063763206bbcbdbab92abc097c983ffd (diff) | |
GRHydro: reset points inside the excision region to atmosphere in
InitialAtmosphereReset(M) so that the ID thorns do not have to know
about he excised region.
From: Roland Haas <roland.haas@physics.gatech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@447 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src')
| -rw-r--r-- | src/GRHydro_UpdateMask.F90 | 4 | ||||
| -rw-r--r-- | src/GRHydro_UpdateMaskM.F90 | 4 |
2 files changed, 6 insertions, 2 deletions
diff --git a/src/GRHydro_UpdateMask.F90 b/src/GRHydro_UpdateMask.F90 index b5a84a5..67a42cb 100644 --- a/src/GRHydro_UpdateMask.F90 +++ b/src/GRHydro_UpdateMask.F90 @@ -452,7 +452,9 @@ subroutine GRHydro_InitialAtmosphereReset(CCTK_ARGUMENTS) do j = 1, cctk_lsh(2) do i = 1, cctk_lsh(1) - if (rho(i,j,k) .le. rho_min) then + if (rho(i,j,k) .le. rho_min .or. & + GRHydro_enable_internal_excision /= 0 .and. & + hydro_excision_mask(i,j,k) .ne. 0) then rho(i,j,k) = rho_min det = SPATIAL_DETERMINANT(g11(i,j,k), g12(i,j,k), g13(i,j,k), \ g22(i,j,k), g23(i,j,k), g33(i,j,k)) diff --git a/src/GRHydro_UpdateMaskM.F90 b/src/GRHydro_UpdateMaskM.F90 index b5066e2..a965f53 100644 --- a/src/GRHydro_UpdateMaskM.F90 +++ b/src/GRHydro_UpdateMaskM.F90 @@ -323,7 +323,9 @@ subroutine GRHydro_InitialAtmosphereResetM(CCTK_ARGUMENTS) do j = 1, cctk_lsh(2) do i = 1, cctk_lsh(1) - if (rho(i,j,k) .le. rho_min) then + if (rho(i,j,k) .le. rho_min .or. & + (GRHydro_enable_internal_excision .ne. 0 .and. & + hydro_excision_mask(i,j,k) .ne. 0) ) then rho(i,j,k) = rho_min vup(i,j,k,1) = 0.0d0 vup(i,j,k,2) = 0.0d0 |