add mask for hydro excision to HydroBase

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinBase/HydroBase/trunk@21 57fe0bb3-ccba-405f-9b23-de0201f165b7

1 files changed, 14 insertions, 1 deletions

diff --git a/src/Initialisation.c b/src/Initialisation.c
index fcd71f5..3391fde 100644
--- a/src/Initialisation.c
+++ b/src/Initialisation.c
@@ -1,7 +1,7 @@
 #include <cctk.h>
 #include <cctk_Arguments.h>
 #include <cctk_Parameters.h>
-
+#include "HydroBase.h"
 
 
 void HydroBase_Zero (CCTK_ARGUMENTS)
@@ -166,3 +166,16 @@ void HydroBase_Bvec_zero (CCTK_ARGUMENTS)
                "Unsupported parameter value for InitBase::initial_data_setup_method");
   }
 }
+
+void HydroBase_InitExcisionMask (CCTK_ARGUMENTS)
+{
+  DECLARE_CCTK_ARGUMENTS;
+  DECLARE_CCTK_PARAMETERS;
+  
+  int const np = cctk_lsh[0] * cctk_lsh[1] * cctk_lsh[2];
+  
+  #pragma omp parallel for
+  for (int i=0; i<np; ++i) {
+    hydro_excision_mask[i] = HYDRO_EXCISION_NORMAL;
+  }
+}