diff --git a/common/TwoPhase.h b/common/TwoPhase.h
index 9b356892..4012bdaf 100644
--- a/common/TwoPhase.h
+++ b/common/TwoPhase.h
@@ -549,7 +549,7 @@ void TwoPhase::ComputeLocal(){
 void TwoPhase::ComputeLocalBlob(){
     int i,j,k,n,label;
 	double vF,vS;
-	vF = 0.0; vS= -1.0;
+	vF = 0.0; vS=0.0;
 //    const RankInfoStruct rank_info(Dm.rank,Dm.nprocx,Dm.nprocy,Dm.nprocz);
 
 	int cube[8][3] = {{0,0,0},{1,0,0},{0,1,0},{1,1,0},{0,0,1},{1,0,1},{0,1,1},{1,1,1}};
diff --git a/tests/TestBlobAnalyze.cpp b/tests/TestBlobAnalyze.cpp
index 3f676858..fa280702 100644
--- a/tests/TestBlobAnalyze.cpp
+++ b/tests/TestBlobAnalyze.cpp
@@ -207,6 +207,13 @@ int main(int argc, char **argv)
 				else{
 					Dm.id[n] = 1;
 				}
+				Averages.SDn(i,j,k) = Averages.Phase(i,j,k);
+				Averages.Phase_tplus(i,j,k) = Averages.Phase(i,j,k);
+				Averages.Phase_tminus(i,j,k) = Averages.Phase(i,j,k);
+				Averages.Press(i,j,k) = 0.0;
+				Averages.Vel_x(i,j,k) = 0.0;
+				Averages.Vel_y(i,j,k) = 0.0;
+				Averages.Vel_z(i,j,k) = 0.0;
 			}
 		}
 	}