diff --git a/models/ColorModel.cpp b/models/ColorModel.cpp
index b834ba3a..510e54ed 100644
--- a/models/ColorModel.cpp
+++ b/models/ColorModel.cpp
@@ -947,6 +947,7 @@ double ScaLBL_ColorModel::ImageInit(std::string Filename){
 	ScaLBL_CopyToDevice(Phi, PhaseLabel, Nx*Ny*Nz*sizeof(double));
 	MPI_Barrier(comm);
 	
+	ScaLBL_D3Q19_Init(fq, Np);
 	ScaLBL_PhaseField_Init(dvcMap, Phi, Den, Aq, Bq, 0, ScaLBL_Comm->LastExterior(), Np);
 	ScaLBL_PhaseField_Init(dvcMap, Phi, Den, Aq, Bq, ScaLBL_Comm->FirstInterior(), ScaLBL_Comm->LastInterior(), Np);
 	MPI_Barrier(comm);