User: rhaas Date: 2013/01/14 08:23 AM

Modified: /trunk/src/ GRHydro_UpdateMask.F90, GRHydro_UpdateMaskM.F90

Log: 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

File Changes:

Directory: /trunk/src/ ======================

File [modified]: GRHydro_UpdateMask.F90 Delta lines: +3 -1 =================================================================== --- trunk/src/GRHydro_UpdateMask.F90 2013-01-14 14:23:12 UTC (rev 446) +++ trunk/src/GRHydro_UpdateMask.F90 2013-01-14 14:23:14 UTC (rev 447) @@ -452,7 +452,9 @@ 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))

File [modified]: GRHydro_UpdateMaskM.F90 Delta lines: +3 -1 =================================================================== --- trunk/src/GRHydro_UpdateMaskM.F90 2013-01-14 14:23:12 UTC (rev 446) +++ trunk/src/GRHydro_UpdateMaskM.F90 2013-01-14 14:23:14 UTC (rev 447) @@ -323,7 +323,9 @@ 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