Add average atomic mass: HydroBase::Abar. Ticket #1649.svn

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

1 files changed, 22 insertions, 1 deletions

diff --git a/src/Initialisation.c b/src/Initialisation.c
index 9aa31ec..762b402 100644
--- a/src/Initialisation.c
+++ b/src/Initialisation.c
@@ -20,7 +20,14 @@ void HydroBase_Zero (CCTK_ARGUMENTS)
     w_lorentz[i] = 1.0;
     eps[i     ]  = 0.0;
   }
-  
+
+  if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::Abar") >= 1) {
+#pragma omp parallel for
+    for (int i=0; i<np; ++i) {
+      Abar[i     ] = 0.0;
+    }
+  }
+
   if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::temperature") >= 1) {
 #pragma omp parallel for
     for (int i=0; i<np; ++i) {
@@ -55,6 +62,13 @@ void HydroBase_Zero (CCTK_ARGUMENTS)
       }
     }
 
+    if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::Abar") >= 2) {
+#pragma omp parallel for
+      for (int i=0; i<np; ++i) {
+        Abar_p[i     ] = 0.0;
+      }
+    }
+
     if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::temperature") >= 2) {
 #pragma omp parallel for
       for (int i=0; i<np; ++i) {
@@ -86,6 +100,13 @@ void HydroBase_Zero (CCTK_ARGUMENTS)
                  "Too many active time levels for HydroBase variables");
     }
 
+    if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::Abar") >= 3) {
+#pragma omp parallel for
+      for (int i=0; i<np; ++i) {
+        Abar_p_p[i     ] = 0.0;
+      }
+    }
+
     if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::temperature") >= 3) {
 #pragma omp parallel for
       for (int i=0; i<np; ++i) {