Changeset 18776

Timestamp:

11/12/14 17:25:14 (10 years ago)

Author:

bdef

Message:

CHG:correction in frictionhydro aw computation

File:

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

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

@@ -168 +168 @@
         /*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  **/
@@ -220 +221 @@
         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));