From 3ed9afd03e1b1d2628d0eb5521048e58902a2574 Mon Sep 17 00:00:00 2001
From: James E McClure <mcclurej@vt.edu>
Date: Thu, 17 Sep 2015 15:36:48 -0400
Subject: [PATCH] Try to prevent identification fo isolated blobs near the
 solid

---
 analysis/analysis.cpp        | 2 +-
 tests/lbpm_color_simulator.h | 3 +--
 2 files changed, 2 insertions(+), 3 deletions(-)

diff --git a/analysis/analysis.cpp b/analysis/analysis.cpp
index 85925e0d..bc6f40f2 100644
--- a/analysis/analysis.cpp
+++ b/analysis/analysis.cpp
@@ -145,7 +145,7 @@ int ComputeLocalBlobIDs( const DoubleArray& Phase, const DoubleArray& SignDist,
     Array<bool> isPhase(Nx,Ny,Nz);
     memset(isPhase.get(),0,Nx*Ny*Nz*sizeof(bool));
     for (size_t i=0; i<N; i++) {
-        if ( SignDist(i) < 0.0) {
+        if ( SignDist(i) <= vS) {
             // Solid phase 
             LocalBlobID(i) = -2;
         } else {
diff --git a/tests/lbpm_color_simulator.h b/tests/lbpm_color_simulator.h
index 8f2133c2..91a96de7 100644
--- a/tests/lbpm_color_simulator.h
+++ b/tests/lbpm_color_simulator.h
@@ -56,7 +56,7 @@ public:
         MPI_Comm newcomm;
         MPI_Comm_dup(MPI_COMM_WORLD,&newcomm);
         double vF = 0.0;
-        double vS = 0.0;
+        double vS = -1.0; // one voxel buffer region around solid
         IntArray& ids = new_index->second;
         new_index->first = ComputeGlobalBlobIDs(Nx-2,Ny-2,Nz-2,rank_info,*phase,dist,vF,vS,ids,newcomm);
         MPI_Comm_free(&newcomm);
@@ -168,7 +168,6 @@ private:
     double beta;
 };
 
-
 // Function to start the analysis
 void run_analysis( int timestep, int restart_interval, 
     const RankInfoStruct& rank_info, TwoPhase& Averages,