diff --git a/analysis/analysis.cpp b/analysis/analysis.cpp
index e24fb85d..e4c13548 100644
--- a/analysis/analysis.cpp
+++ b/analysis/analysis.cpp
@@ -311,7 +311,7 @@ int ComputeLocalBlobIDs( const DoubleArray& Phase, const DoubleArray& SignDist,
     return nblobs;
 }
 
-int ComputeLocalPhaseComponent(IntArray &PhaseID, char VALUE, IntArray &ComponentLabel, bool periodic )
+int ComputeLocalPhaseComponent(IntArray &PhaseID, int VALUE, IntArray &ComponentLabel, bool periodic )
 {
     PROFILE_START("ComputeLocalPhaseComponent");
     size_t Nx = ComponentLabel.size(0);
diff --git a/analysis/analysis.h b/analysis/analysis.h
index 059f7260..2dd7a512 100644
--- a/analysis/analysis.h
+++ b/analysis/analysis.h
@@ -41,6 +41,18 @@ int ComputeBlob( IntArray &blobs, int &nblobs, int &ncubes, IntArray &indicator,
 int ComputeLocalBlobIDs( const DoubleArray& Phase, const DoubleArray& SignDist, 
     double vF, double vS, IntArray& LocalBlobID, bool periodic=true );
 
+/*!
+ *  @brief Compute blob of an arbitrary phase
+ *  @details Compute the connected components for Phase(i,j,k)=VALUE
+ *  @return the number of connected components of the phase
+ *  @param[in] PhaseID
+ *  @param[in] VALUE
+ *  @param[out] ComponentLabel
+ *  @param[in] periodic
+ */
+
+int ComputeLocalPhaseComponent(IntArray &PhaseID, int VALUE, IntArray &ComponentLabel,
+		bool periodic )
 
 /*!
  * @brief  Compute the blob