1 files changed, 55 insertions, 0 deletions
diff --git a/src/Initialisation.c b/src/Initialisation.c
index a739091..fcd71f5 100644
--- a/src/Initialisation.c
+++ b/src/Initialisation.c
@@ -111,3 +111,58 @@ void HydroBase_Y_e_one (CCTK_ARGUMENTS)
                "Unsupported parameter value for InitBase::initial_data_setup_method");
   }
 }
+
+
+
+void HydroBase_Bvec_zero (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) {
+    Bvec[i]      = 0.0;
+    Bvec[i+  np] = 0.0;
+    Bvec[i+2*np] = 0.0;
+  }
+  
+  if (CCTK_EQUALS (initial_data_setup_method, "init_some_levels") ||
+      CCTK_EQUALS (initial_data_setup_method, "init_single_levels"))
+  {
+    /* do nothing */
+  }
+  else if (CCTK_EQUALS (initial_data_setup_method, "init_all_levels"))
+  {
+    
+    if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::Bvec") >= 2) {
+#pragma omp parallel for
+      for (int i=0; i<np; ++i) {
+        Bvec_p[i]      = 0.0;
+        Bvec_p[i+  np] = 0.0;
+        Bvec_p[i+2*np] = 0.0;
+      }
+    }
+    
+    if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::Bvec") >= 3) {
+#pragma omp parallel for
+      for (int i=0; i<np; ++i) {
+        Bvec_p_p[i]      = 0.0;
+        Bvec_p_p[i+  np] = 0.0;
+        Bvec_p_p[i+2*np] = 0.0;
+      }
+    }
+    
+    if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::Bvec") >= 4) {
+      CCTK_WARN (CCTK_WARN_ABORT,
+                 "Too many active time levels for HydroBase::Bvec");
+    }
+    
+  }
+  else
+  {
+    CCTK_WARN (CCTK_WARN_ABORT,
+               "Unsupported parameter value for InitBase::initial_data_setup_method");
+  }
+}