Index: /issm/trunk/src/c/objects/Loads/Friction.cpp
===================================================================
--- /issm/trunk/src/c/objects/Loads/Friction.cpp	(revision 7921)
+++ /issm/trunk/src/c/objects/Loads/Friction.cpp	(revision 7922)
@@ -110,4 +110,7 @@
 	else _error_("element_type %s not supported yet",element_type);
 
+	/*Checks that s-1>0 if v=0*/
+	if(vmag==0 && (s-1)<0) _error_("velocity is 0 ans (s-1)=%g<0, alpha_complement is Inf",s-1);
+
 	alpha2=pow(drag_coefficient,2)*pow(Neff,r)*pow(vmag,(s-1));
 	_assert_(!isnan(alpha2));
@@ -171,4 +174,7 @@
 	}
 	else _error_("element_type %s not supported yet",element_type);
+
+	/*Checks that s-1>0 if v=0*/
+	if(vmag==0 && (s-1)<0) _error_("velocity is 0 ans (s-1)=%g<0, alpha_complement is Inf",s-1);
 
 	alpha2=pow(drag_coefficient,2)*pow(Neff,r)*pow(vmag,(s-1));
@@ -230,5 +236,5 @@
 	vmag=sqrt(pow(vx,2)+pow(vy,2));
 
-	/*Checks*/
+	/*Checks that s-1>0 if v=0*/
 	if(vmag==0 && (s-1)<0) _error_("velocity is 0 ans (s-1)=%g<0, alpha_complement is Inf",s-1);