corrected velocity output to work correctly for dim=1,2,3

git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@7594 2fb0f335-1f38-0410-981e-8018bf24f1b0

dumux/decoupled/2p/diffusion/fv/fvvelocity2p.hh

@@ -191,8 +191,11 @@ public:
             }
 
             Dune::FieldVector < Scalar, dim > refVelocity(0);
-            refVelocity[0] = 0.5 * (fluxW[1] - fluxW[0]);
-            refVelocity[1] = 0.5 * (fluxW[3] - fluxW[2]);
+                refVelocity[0] = 0.5 * (fluxW[1] - fluxW[0]);
+            if (dim > 1)
+                refVelocity[1] = 0.5 * (fluxW[3] - fluxW[2]);
+            if (dim == 3)
+                refVelocity[2] = 0.5 * (fluxW[5] - fluxW[4]);
 
             const Dune::FieldVector<Scalar, dim>& localPos =
                     ReferenceElementContainer::general(eIt->geometry().type()).position(0, 0);
@@ -211,7 +214,10 @@ public:
 
             refVelocity = 0;
             refVelocity[0] = 0.5 * (fluxNW[1] - fluxNW[0]);
+            if (dim > 1)
             refVelocity[1] = 0.5 * (fluxNW[3] - fluxNW[2]);
+            if (dim == 3)
+                refVelocity[2] = 0.5 * (fluxNW[5] - fluxNW[4]);
 
             // calculate the element velocity by the Piola transformation
             elementVelocity = 0;