diff --git a/common/TwoPhase.cpp b/common/TwoPhase.cpp
index 9f60281e..81c9b9d6 100644
--- a/common/TwoPhase.cpp
+++ b/common/TwoPhase.cpp
@@ -333,6 +333,15 @@ void TwoPhase::UpdateMeshValues()
 	//...........................................................................
 	Dm.CommunicateMeshHalo(SDn_z);
 	//...........................................................................
+	Dm.CommunicateMeshHalo(SDs);
+	pmmc_MeshGradient(SDs,SDs_x,SDs_y,SDs_z,Nx,Ny,Nz);
+	//...........................................................................
+	Dm.CommunicateMeshHalo(SDs_x);
+	//...........................................................................
+	Dm.CommunicateMeshHalo(SDs_y);
+	//...........................................................................
+	Dm.CommunicateMeshHalo(SDs_z);
+	//...........................................................................
 	// Compute the mesh curvature of the phase indicator field
 	pmmc_MeshCurvature(SDn, MeanCurvature, GaussCurvature, Nx, Ny, Nz);
 	//...........................................................................