diff --git a/analysis/runAnalysis.cpp b/analysis/runAnalysis.cpp
index 43e46b82..27c7d5b4 100644
--- a/analysis/runAnalysis.cpp
+++ b/analysis/runAnalysis.cpp
@@ -492,10 +492,13 @@ void runAnalysis::run( int timestep, TwoPhase& Averages, const double *Phi,
          matches(type,AnalysisType::IdentifyBlobs) )
     {
         phase = std::shared_ptr<DoubleArray>(new DoubleArray(d_N[0],d_N[1],d_N[2]));
-        ScaLBL_CopyToHost(phase->data(),Phi,N*sizeof(double));
+        //ScaLBL_CopyToHost(phase->data(),Phi,N*sizeof(double));
+		ScaLBL_Comm.RegularLayout(Map,Phi,phase.get());
+
     }
     if ( matches(type,AnalysisType::CopyPhaseIndicator) ) {
         memcpy(Averages.Phase_tplus.data(),phase->data(),N*sizeof(double));
+
         //Averages.ColorToSignedDistance(d_beta,Averages.Phase,Averages.Phase_tplus);
     }
     if ( matches(type,AnalysisType::ComputeAverages) ) {
diff --git a/tests/lbpm_dfh_simulator.cpp b/tests/lbpm_dfh_simulator.cpp
index a62da5fc..d4daf942 100644
--- a/tests/lbpm_dfh_simulator.cpp
+++ b/tests/lbpm_dfh_simulator.cpp
@@ -751,7 +751,7 @@ int main(int argc, char **argv)
 			PROFILE_STOP("Update");
 
 			// Run the analysis
-			//            analysis.run( timestep, *Averages, Phi, Pressure, Velocity, fq, Den );
+			analysis.run( timestep, *Averages, Phi, Pressure, Velocity, fq, Den );
 
 		}
         analysis.finish();