diff --git a/analysis/Minkowski.cpp b/analysis/Minkowski.cpp
index 3b22bfcd..6d90ec52 100644
--- a/analysis/Minkowski.cpp
+++ b/analysis/Minkowski.cpp
@@ -139,7 +139,7 @@ void Minkowski::MeasureObject(){
 }
 
 
-void Minkowski::MeasureConnectedPathway(){
+int Minkowski::MeasureConnectedPathway(){
 	/*
 	 * compute the connected pathway for object with LABEL in id field
 	 * compute the labels for connected components
@@ -163,13 +163,14 @@ void Minkowski::MeasureConnectedPathway(){
 	
 	// Extract only the connected part of NWP
 	double vF=0.0; 
-	ComputeGlobalBlobIDs(Nx-2,Ny-2,Nz-2,Dm->rank_info,distance,distance,vF,vF,label,Dm->Comm);
+	int n_components = ComputeGlobalBlobIDs(Nx-2,Ny-2,Nz-2,Dm->rank_info,distance,distance,vF,vF,label,Dm->Comm);
+//	int n_connected_components = ComputeGlobalPhaseComponent(Nx-2,Ny-2,Nz-2,Dm->rank_info,const IntArray &PhaseID, int &VALUE, BlobIDArray &GlobalBlobID, Dm->Comm )
 	MPI_Barrier(Dm->Comm);
 	
 	for (int k=0; k<Nz; k++){
 		for (int j=0; j<Ny; j++){
 			for (int i=0; i<Nx; i++){
-				if (id(i,j,k) == LABEL && label(i,j,k) == 0){
+				if ( label(i,j,k) == 0){
 					id(i,j,k) = 0;
 				}
 				else{
@@ -179,6 +180,7 @@ void Minkowski::MeasureConnectedPathway(){
 		}
 	}
 	MeasureObject();
+	return n_connected_components; 
 }
 
 
diff --git a/analysis/Minkowski.h b/analysis/Minkowski.h
index d29feb67..0d36e33d 100644
--- a/analysis/Minkowski.h
+++ b/analysis/Minkowski.h
@@ -41,7 +41,7 @@ public:
 	double Ai,Ji,Xi,Vi;
 	// Global averages (all processes)
 	double Ai_global,Ji_global,Xi_global,Vi_global;
-
+	int n_connected_components;
 	//...........................................................................
 	int Nx,Ny,Nz;
 	double V(){
@@ -62,7 +62,7 @@ public:
 	Minkowski(std::shared_ptr <Domain> Dm);
 	~Minkowski();
 	void MeasureObject();
-	void MeasureConnectedPathway();
+	int MeasureConnectedPathway();
 	void ComputeScalar(const DoubleArray& Field, const double isovalue);
 
 	void PrintAll();