From 5f02c8e3db7ee5bb7a8c6a6f8c97cdd044bcd21b Mon Sep 17 00:00:00 2001
From: James E McClure <jemcclur@gmail.com>
Date: Mon, 2 Feb 2015 17:09:21 -0500
Subject: [PATCH] Fixed the bug from tests/lb2_Color_wia_mpi in
 tests/lbpm_color_simulator

---
 tests/lbpm_color_simulator.cpp | 6 +++++-
 1 file changed, 5 insertions(+), 1 deletion(-)

diff --git a/tests/lbpm_color_simulator.cpp b/tests/lbpm_color_simulator.cpp
index cffc4adf..f014f9e7 100644
--- a/tests/lbpm_color_simulator.cpp
+++ b/tests/lbpm_color_simulator.cpp
@@ -275,10 +275,10 @@ int main(int argc, char **argv)
 		printf("********************************************************\n");
 	}
 
+	// Initialized domain and averaging framework for Two-Phase Flow
 	// Initialized domain and averaging framework for Two-Phase Flow
 	int BC=pBC;
 	Domain Dm(Nx,Ny,Nz,rank,nprocx,nprocy,nprocz,Lx,Ly,Lz,BC);
-	Dm.CommInit(MPI_COMM_WORLD);
 	TwoPhase Averages(Dm);
 
 	InitializeRanks( rank, nprocx, nprocy, nprocz, iproc, jproc, kproc,
@@ -423,6 +423,9 @@ int main(int argc, char **argv)
 		if (rank == 0)	printf("Resetting phi_s = %f, das = %f, dbs = %f \n", phi_s, das, dbs);
 		FlipID(id,Nx*Ny*Nz);
 	}
+	// Initialize communication structures in averaging domain
+	for (i=0; i<Dm.Nx*Dm.Ny*Dm.Nz; i++) Dm.id[i] = 1;
+	Dm.CommInit(MPI_COMM_WORLD);
 
 	// Set up MPI communication structurese
 	if (rank==0)	printf ("Setting up communication control structures \n");
@@ -1583,6 +1586,7 @@ int main(int argc, char **argv)
 			Averages.PrintAll(timestep);
 			//....................................................................
 		}
+
 		if (timestep%RESTART_INTERVAL == 0){
 			if (pBC){
 				//err = fabs(sat_w - sat_w_previous);