initialize w_lorentz as well as the other declared quantities - especially as the velocities

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

1 files changed, 18 insertions, 15 deletions

diff --git a/src/Initialisation.c b/src/Initialisation.c
index 1149897..4ed861f 100644
--- a/src/Initialisation.c
+++ b/src/Initialisation.c
@@ -13,11 +13,12 @@ void HydroBase_Zero (CCTK_ARGUMENTS)
   
 #pragma omp parallel for
   for (int i=0; i<np; ++i) {
-    rho[i     ] = 0.0;
-    vel[i     ] = 0.0;
-    vel[i+  np] = 0.0;
-    vel[i+2*np] = 0.0;
-    eps[i     ] = 0.0;
+    rho[i     ]  = 0.0;
+    vel[i     ]  = 0.0;
+    vel[i+  np]  = 0.0;
+    vel[i+2*np]  = 0.0;
+    w_lorentz[i] = 1.0;
+    eps[i     ]  = 0.0;
   }
   
   if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::temperature") >= 1) {
@@ -45,11 +46,12 @@ void HydroBase_Zero (CCTK_ARGUMENTS)
     if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::rho") >= 2) {
 #pragma omp parallel for
       for (int i=0; i<np; ++i) {
-        rho_p[i     ] = 0.0;
-        vel_p[i     ] = 0.0;
-        vel_p[i+  np] = 0.0;
-        vel_p[i+2*np] = 0.0;
-        eps_p[i     ] = 0.0;
+        rho_p[i      ] = 0.0;
+        vel_p[i      ] = 0.0;
+        vel_p[i+  np ] = 0.0;
+        vel_p[i+2*np ] = 0.0;
+        w_lorentz_p[i] = 1.0;
+        eps_p[i      ] = 0.0;
       }
     }
 
@@ -70,11 +72,12 @@ void HydroBase_Zero (CCTK_ARGUMENTS)
     if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::rho") >= 3) {
 #pragma omp parallel for
       for (int i=0; i<np; ++i) {
-        rho_p_p[i     ] = 0.0;
-        vel_p_p[i     ] = 0.0;
-        vel_p_p[i+  np] = 0.0;
-        vel_p_p[i+2*np] = 0.0;
-        eps_p_p[i     ] = 0.0;
+        rho_p_p[i      ] = 0.0;
+        vel_p_p[i      ] = 0.0;
+        vel_p_p[i+  np ] = 0.0;
+        vel_p_p[i+2*np ] = 0.0;
+        w_lorentz_p_p[i] = 1.0;
+        eps_p_p[i      ] = 0.0;
       }
     }