diff --git a/common/Domain.h b/common/Domain.h
index 73b3eee1..8305b4cc 100755
--- a/common/Domain.h
+++ b/common/Domain.h
@@ -727,12 +727,12 @@ inline void SignedDistance(double *Distance, int nspheres, double *List_cx, doub
 		r = List_rad[p];
 		// Check if
 		// Range for this sphere in global indexing
-		imin = int ((cx-r)/hx);
-		imax = int ((cx+r)/hx)+2;
-		jmin = int ((cy-r)/hy);
-		jmax = int ((cy+r)/hy)+2;
-		kmin = int ((cz-r)/hz);
-		kmax = int ((cz+r)/hz)+2;
+		imin = int ((cx-2*r)/hx);
+		imax = int ((cx+2*r)/hx)+2;
+		jmin = int ((cy-2*r)/hy);
+		jmax = int ((cy+2*r)/hy)+2;
+		kmin = int ((cz-2*r)/hz);
+		kmax = int ((cz+2*r)/hz)+2;
 		// Obviously we have to do something at the edges
 		if (imin<0)		imin = 0;
 		if (imin>Nx)	imin = Nx;
diff --git a/tests/BlobAnalysis.cpp b/tests/BlobAnalysis.cpp
index 4d181c70..d1ad2cb3 100644
--- a/tests/BlobAnalysis.cpp
+++ b/tests/BlobAnalysis.cpp
@@ -593,7 +593,7 @@ int main(int argc, char **argv)
 	DoubleArray BlobAverages(NUM_AVERAGES,nblobs);
 	
 	// Map the signed distance for the analysis
-	for (i=0; i<Nx*Ny*Nz; i++)	SignDist.data[i] -= (1.0); 
+//	for (i=0; i<Nx*Ny*Nz; i++)	SignDist.data[i] -= (1.0); 
 	
 	// Compute the porosity
 	porosity=0.0;