Revert "* improve handling of cs2 < 0 in the "standard" (non-MHD, non-hot EOS) eigenvalue routine."

actually this introduce the real valued atmosphere mask.

This reverts commit deccf528b463820b615dfa73e46e137cd59c65be.

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@372 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45

1 files changed, 2 insertions, 35 deletions

diff --git a/src/GRHydro_UpdateMask.F90 b/src/GRHydro_UpdateMask.F90
index c11d218..1808df3 100644
--- a/src/GRHydro_UpdateMask.F90
+++ b/src/GRHydro_UpdateMask.F90
@@ -48,9 +48,7 @@ subroutine GRHydroUpdateAtmosphereMask(CCTK_ARGUMENTS)
                  densrhs(i,j,k) = 0.0d0
                  srhs(i,j,k,:)   = 0.0d0
                  taurhs(i,j,k)  = 0.0d0
-                 ! Set real-valued mask! This will be sync'ed and right after syncing translated to
-                 ! our standard integer based mask (so that atmosphere_mask is still valid!).
-                 atmosphere_mask_real(i,j,k) = 1
+                 atmosphere_mask(i,j,k) = 1
               end if
            end do
         end do
@@ -69,9 +67,7 @@ subroutine GRHydroUpdateAtmosphereMask(CCTK_ARGUMENTS)
                  densrhs(i,j,k) = 0.0d0
                  srhs(i,j,k,:)   = 0.0d0
                  taurhs(i,j,k)  = 0.0d0
-                 ! Set real-valued mask! This will be sync'ed and right after syncing translated to
-                 ! our standard integer based mask (so that atmosphere_mask is still valid!).
-                 atmosphere_mask_real(i,j,k) = 1
+                 atmosphere_mask(i,j,k) = 1
               end if
            end do
         end do
@@ -81,33 +77,6 @@ subroutine GRHydroUpdateAtmosphereMask(CCTK_ARGUMENTS)
 
 end subroutine GRHydroUpdateAtmosphereMask
 
-
-subroutine GRHydroPostSyncAtmosphereMask(CCTK_ARGUMENTS)
-
-  implicit none
-  
-  DECLARE_CCTK_ARGUMENTS
-  DECLARE_CCTK_PARAMETERS
-
-  CCTK_INT :: i,j,k
-
-!! This sets the integer atmo mask based on the real-valued (and sync'ed) atmo mask
-
-   !$OMP PARALLEL DO PRIVATE(k,j,i)
-   do k = 1, cctk_lsh(3)
-      do j = 1, cctk_lsh(2)
-	 do i = 1, cctk_lsh(1)
-	    if ( atmosphere_mask_real(i,j,k) .ne. 0) then
-	       atmosphere_mask(i,j,k) = 1
-	    end if
-	 end do
-      end do
-   end do
-   !$OMP END PARALLEL DO
-
-end subroutine GRHydroPostSyncAtmosphereMask
-
-
  /*@@
    @routine    GRHydro_SetupMask
    @date       Thu Jun 20 13:27:28 2002
@@ -147,7 +116,6 @@ subroutine GRHydro_SetupMask(CCTK_ARGUMENTS)
   endif
 
   atmosphere_mask = 0
-  atmosphere_mask_real = 0
   
   call CCTK_INFO("Setting up the atmosphere mask: all points are not_atmosphere")
   
@@ -295,7 +263,6 @@ subroutine GRHydro_AtmosphereReset(CCTK_ARGUMENTS)
                   velz(i,j,k), eps(i,j,k), press(i,j,k), w_lorentz(i,j,k))
              if (wk_atmosphere .eq. 0) then
                 atmosphere_mask(i, j, k) = 0
-                atmosphere_mask_real(i, j, k) = 0
              end if
           endif