Index: /issm/trunk-jpl/src/c/objects/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/objects/Elements/Tria.cpp	(revision 11557)
+++ /issm/trunk-jpl/src/c/objects/Elements/Tria.cpp	(revision 11558)
@@ -524,5 +524,5 @@
 		if(stabilization==2){
 			/*Streamline upwinding*/
-			vel=sqrt(pow(vx,2.)+pow(vy,2.));
+			vel=sqrt(pow(vx,2.)+pow(vy,2.))+1.e-8;
 			K[0][0]=h/(2*vel)*vx*vx;
 			K[1][0]=h/(2*vel)*vy*vx;
@@ -3667,6 +3667,6 @@
 	//	vy_input->GetInputValue(&vy,gauss);
 	//	grade_g[iv]=-2*drag*alpha_complement*((lambda*vx+mu*vy));
-	//	VecSetValues(gradient,NUMVERTICES,doflist1,(const double*)grade_g,INSERT_VALUES);
 	//}
+	//VecSetValues(gradient,NUMVERTICES,doflist1,(const double*)grade_g,INSERT_VALUES);
 	/*End Analytical gradient*/