diff --git a/cpu/Color.cpp b/cpu/Color.cpp
index a117bedb..036453df 100644
--- a/cpu/Color.cpp
+++ b/cpu/Color.cpp
@@ -1299,9 +1299,9 @@ extern "C" void ScaLBL_D3Q19_ColorCollide_gen( char *ID, double *disteven, doubl
 			distodd[7*N+n] = f15;
 			distodd[8*N+n] = f17;
 			//...Store the Velocity..........................
-			Velocity[n] = jx;
-			Velocity[N+n] = jy;
-			Velocity[2*N+n] = jz;
+			Velocity[n] = jx/rho_0;
+			Velocity[N+n] = jy/rho_0;
+			Velocity[2*N+n] = jz/rho_0;
 			//***************************************************************
 		}	// check if n is in the solid
 	} // loop over n
diff --git a/gpu/Color.cu b/gpu/Color.cu
index e352a995..4bdf4f5f 100644
--- a/gpu/Color.cu
+++ b/gpu/Color.cu
@@ -1300,9 +1300,9 @@ __global__  void dvc_ScaLBL_D3Q19_ColorCollide_gen( char *ID, double *disteven,
 			distodd[7*N+n] = f15;
 			distodd[8*N+n] = f17;
 			//...Store the Velocity..........................
-			Velocity[n] = jx;
-			Velocity[N+n] = jy;
-			Velocity[2*N+n] = jz;
+			Velocity[n] = jx/rho0;
+			Velocity[N+n] = jy/rho0;
+			Velocity[2*N+n] = jz/rho0;
 			//***************************************************************
 
 			}// check if n is in the solid