Index: /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp	(revision 26584)
+++ /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp	(revision 26585)
@@ -198,5 +198,5 @@
 
 	/*Check to prevent dividing by zero if vmag==0*/
-	if(vmag==0. && (s-1.)<0.) alpha_complement=0.;
+	if(vmag==0. && (s-1.)<=0.) alpha_complement=0.;
 	else alpha_complement=pow(Neff,r)*pow(vmag,(s-1));
 
@@ -345,5 +345,5 @@
 
 	/*Check to prevent dividing by zero if vmag==0*/
-	if(vmag==0. && (s-1.)<0.){
+	if(vmag==0. && (s-1.)<=0.){
 		alpha2=0.;
 	}
@@ -541,5 +541,5 @@
 
 	/*Check to prevent dividing by zero if vmag==0*/
-	if(vmag==0. && (s-1.)<0.) alpha2=0.;
+	if(vmag==0. && (s-1.)<=0.) alpha2=0.;
 	else alpha2=drag_coefficient*drag_coefficient*pow(Neff,r)*pow(vmag,(s-1.));
 
@@ -587,5 +587,5 @@
 
 	/*Check to prevent dividing by zero if vmag==0*/
-	if(vmag==0. && (s-1.)<0.) alpha2=0.;
+	if(vmag==0. && (s-1.)<=0.) alpha2=0.;
 	else alpha2=drag_coefficient*drag_coefficient*pow(Neff,r)*pow(vmag,(s-1.));