[Commits] [svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 447)

[rhaas at tapir.caltech.edu]

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 at 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