Added call to TwoPhase::WriteSurfaces to lbpm_color_simulator so that surfaces are written at the end of simulation

tests/lbpm_color_simulator.cpp

@@ -926,6 +926,8 @@ int main(int argc, char **argv)
 	DeviceBarrier();
 	CopyToHost(Averages.Phase.get(),Phi,N*sizeof(double));
 
+	Averages.WriteSurfaces(0);
+
 	sprintf(LocalRankFilename,"%s%s","Phase.",LocalRankString);
 	FILE *PHASE;
 	PHASE = fopen(LocalRankFilename,"wb");