Index: /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp	(revision 27678)
+++ /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp	(revision 27679)
@@ -901,20 +901,18 @@
 	/*This routine calculates the basal friction coefficient
 	 *
-	 *               C |u_b|^(m-1)
+	 *               C^2 |u_b|^(m-1)
 	 * alpha2= __________________________
-	 *          (1+(C/(Cmax Neff))^1/m |u_b| )^m
+	 *          (1+(C^2/(Cmax Neff))^1/m |u_b| )^m
 	 *
 	 * */
 
 	/*diverse: */
-	IssmDouble  C,coeff,Cmax,m,alpha2;
+	IssmDouble  C,Cmax,m,alpha2;
 
 	/*Recover parameters: */
 	element->GetInputValue(&Cmax,gauss,FrictionCmaxEnum);
-	element->GetInputValue(&coeff,gauss,FrictionCEnum);
+	element->GetInputValue(&C,gauss,FrictionCEnum);
 	element->GetInputValue(&m,gauss,FrictionMEnum);
 
-	/* scale C for a better inversion */
-	C = coeff*coeff;
 
 	/*Get effective pressure*/
@@ -943,5 +941,5 @@
 	}
 	else{
-		alpha2= (C*pow(ub,m-1.)) / pow(1.+  pow(C/(Cmax*Neff),1./m)*ub,m);
+		alpha2= (C*C*pow(ub,m-1.)) / pow(1.+  pow(C*C/(Cmax*Neff),1./m)*ub,m);
 	}