Index: ../trunk-jpl/src/c/classes/Loads/Friction.cpp
===================================================================
--- ../trunk-jpl/src/c/classes/Loads/Friction.cpp	(revision 18775)
+++ ../trunk-jpl/src/c/classes/Loads/Friction.cpp	(revision 18776)
@@ -167,7 +167,8 @@
 
 	/*This routine calculates the basal friction coefficient 
 		Based on Gagliardini 2007, needs a good effective pressure computation
-	  alpha2= C*Neff*[|vel|/{(C^n*Neff^n*As)*(1+alpha*Chi^q)}]^(1/n)  with
+		Not tested so far so use at your own risks
+	  alpha2= [vel^(1-n)/{(As)*(1+alpha*Chi^q)}]^(1/n)  with
 		-Chi=|vel|/(C^n*Neff^n*As)
 		-alpha=(q-1)^(q-1)/q^q  **/
 
@@ -219,7 +220,7 @@
 	alpha=(pow(q_exp-1,q_exp-1))/pow(q_exp,q_exp);
 	Chi=vmag/(pow(C_param,n)*pow(Neff,n)*As);
 
-	alpha2=C_param*Neff*pow((Chi/(1+alpha*pow(Chi,q_exp))),1./n);
+	alpha2=pow((pow(vmag,1.-n)/(As*(1+alpha*pow(Chi,q_exp)))),1./n);
 	_assert_(!xIsNan<IssmDouble>(alpha2));
 	
 	/*Assign output pointers:*/
