Changeset 27679

Timestamp:

04/05/23 07:51:21 (2 years ago)

Author:

Mathieu Morlighem

Message:

CHG: minor

File:

• issm/trunk-jpl/src/c/classes/Loads/Friction.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/classes/Loads/Friction.cpp

@@ -901 +901 @@
         /*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 +941 @@
         }
         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);
         }