Index: /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp	(revision 18989)
+++ /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp	(revision 18990)
@@ -191,5 +191,9 @@
 	}
 	Chi=vmag/(pow(C_param,n)*pow(Neff,n)*As);
-	alpha2= Neff * C_param * pow((Chi/(1 + alpha * pow(Chi,q_exp))),1/n) * 1/vmag;
+
+
+	/*Check to prevent dividing by zero if vmag==0*/
+	if(vmag==0.) alpha2=0.;
+	else alpha2= Neff * C_param * pow((Chi/(1 + alpha * pow(Chi,q_exp))),1/n) * 1/vmag;
 	_assert_(!xIsNan<IssmDouble>(alpha2));
 	/*Assign output pointers:*/