Index: src/Initialisation.c
===================================================================
--- src/Initialisation.c	(revision 48)
+++ src/Initialisation.c	(working copy)
@@ -13,11 +13,12 @@
   
 #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 @@
     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 @@
     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;
       }
     }