GRHydro: copy integer atmosphere mask in Post_Recover_Variables

this is required since the real valued atmosphere mask is not
checkpointed. It is valid since atmosphere_mask_real is non-zero
wherever atmosphere_mask is non-zero due to GRHydroUpdateAtmosphereMask
and atmosphere_mask is non-zero wherever atmosphere_mask_real is
non-zero because of GRHydroPostSyncAtmosphereMask, ie. the two masks are
non-zero/zero at the same points.

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

1 files changed, 40 insertions, 0 deletions

diff --git a/src/GRHydro_UpdateMask.F90 b/src/GRHydro_UpdateMask.F90
index 19cacf3..87247e0 100644
--- a/src/GRHydro_UpdateMask.F90
+++ b/src/GRHydro_UpdateMask.F90
@@ -109,6 +109,45 @@ end subroutine GRHydroPostSyncAtmosphereMask
 
 
  /*@@
+   @routine    GRHydroCopyIntegerMask
+   @date       Wed Jul  4 15:40:16 PDT 2012
+   @author     Roland Haas
+   @desc 
+   Initializes real valued mask with integer valued one.
+   @enddesc 
+   @calls     
+   @calledby   
+   @history 
+ 
+   @endhistory 
+
+@@*/
+
+subroutine GRHydroCopyIntegerMask(CCTK_ARGUMENTS)
+
+  implicit none
+  
+  DECLARE_CCTK_ARGUMENTS
+  DECLARE_CCTK_PARAMETERS
+
+  CCTK_INT :: i,j,k
+
+!! This sets the real atmo mask based on the integer-valued 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)
+	    atmosphere_mask_real(i,j,k) = atmosphere_mask(i,j,k)
+	 end do
+      end do
+   end do
+   !$OMP END PARALLEL DO
+
+end subroutine GRHydroCopyIntegerMask
+
+
+ /*@@
    @routine    GRHydro_SetupMask
    @date       Thu Jun 20 13:27:28 2002
    @author     Ian Hawke
@@ -181,6 +220,7 @@ subroutine GRHydro_InitAtmosMask(CCTK_ARGUMENTS)
       do i = 1, cctk_lsh(1)
         if (rho(i,j,k) .le. GRHydro_rho_min) then
           atmosphere_mask(i,j,k) = 1
+          atmosphere_mask_real(i,j,k) = 1
         end if
       end do
     end do