diff --git a/tests/TestColorBubble.cpp b/tests/TestColorBubble.cpp
index dbfed2c8..9323a712 100644
--- a/tests/TestColorBubble.cpp
+++ b/tests/TestColorBubble.cpp
@@ -397,7 +397,7 @@ int main(int argc, char **argv)
 		// Initialize the phase field and variables
 		ScaLBL_D3Q19_Init(fq, Np);
 		if (rank==0)	printf ("Initializing phase field \n");
-		ScaLBL_PhaseField_Init(dvcMap, Phi, Den, Aq, Bq, Np);
+		ScaLBL_PhaseField_Init(dvcMap, Phi, Den, Aq, Bq, 0, ScaLBL_Comm.last_interior, Np);
 
 		//************ MAIN ITERATION LOOP (timing communications)***************************************
 
diff --git a/tests/TestColorMassBounceback.cpp b/tests/TestColorMassBounceback.cpp
index 4cf8c67a..5c2150a2 100644
--- a/tests/TestColorMassBounceback.cpp
+++ b/tests/TestColorMassBounceback.cpp
@@ -280,7 +280,7 @@ int main(int argc, char **argv)
 		// Initialize the phase field and variables
 		ScaLBL_D3Q19_Init(fq, Np);
 		if (rank==0)	printf ("Initializing phase field \n");
-		ScaLBL_PhaseField_Init(dvcMap, Phi, Den, Aq, Bq, Np);
+		ScaLBL_PhaseField_Init(dvcMap, Phi, Den, Aq, Bq, 0, ScaLBL_Comm.last_interior, Np);
 
 		// *************ODD TIMESTEP*************
 		// Compute the Phase indicator field
@@ -419,7 +419,7 @@ int main(int argc, char **argv)
 		// Initialize the phase field and variables
 		ScaLBL_D3Q19_Init(fq, Np);
 		if (rank==0)	printf ("Initializing phase field \n");
-		ScaLBL_PhaseField_Init(dvcMap, Phi, Den, Aq, Bq, Np);
+		ScaLBL_PhaseField_Init(dvcMap, Phi, Den, Aq, Bq, 0, ScaLBL_Comm.last_interior, Np);
 
 		// *************ODD TIMESTEP*************
 		// Compute the Phase indicator field
diff --git a/tests/TestColorSquareTube.cpp b/tests/TestColorSquareTube.cpp
index 8823da7c..4150eba0 100644
--- a/tests/TestColorSquareTube.cpp
+++ b/tests/TestColorSquareTube.cpp
@@ -321,7 +321,8 @@ int main(int argc, char **argv)
 
 		if (rank==0)	printf ("Initializing distributions \n");
 		// Initialize the phase field and variables
-		ScaLBL_D3Q19_Init(fq, Np);
+		ScaLBL_PhaseField_Init(dvcMap, Phi, Den, Aq, Bq, 0, ScaLBL_Comm.last_interior, Np);
+
 		if (rank==0)	printf ("Initializing phase field \n");
 		ScaLBL_PhaseField_Init(dvcMap, Phi, Den, Aq, Bq, Np);
 		if (Dm.kproc==0){