Index: /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp	(revision 23530)
+++ /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp	(revision 23531)
@@ -316,7 +316,7 @@
 		Based on Gagliardini 2007, needs a good effective pressure computation
 		Not tested so far so use at your own risks
-	  alpha2= NeffC[Chi/(1+alpha*Chi^q)]^(1/n)*1/vel  with
-		-Chi=|vel|/(C^n*Neff^n*As)
-		-alpha=(q-1)^(q-1)/q^q */
+	  alpha2= NeffC[Chi/(1+alpha*Chi^q)]^(1/n)*|vel|^(1/n-1)  with
+		 Chi=|vel|/(C^n*Neff^n*As)
+		 alpha=(q-1)^(q-1)/q^q */
 
 	/*diverse: */
@@ -375,5 +375,5 @@
 	if(vmag==0.) alpha2=0.;
 	else	if (Neff==0) alpha2=0.0;
-	else	alpha2=Neff * C_param * pow(Gamma,1./n) * 1/vmag;
+	else	alpha2=Neff * C_param * pow(Gamma,1./n) * pow(vmag,1/n-1);
 
 	_assert_(!xIsNan<IssmDouble>(alpha2));