Changeset 18991

Timestamp:

01/06/15 13:01:44 (10 years ago)

Author:

bdef

Message:

BUG: adding the case where effective pressure is null

File:

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

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

@@ -195 +195 @@
         /*Check to prevent dividing by zero if vmag==0*/
         if(vmag==0.) alpha2=0.;
+        else if (Neff==0) alpha2=0.0;
         else alpha2= Neff * C_param * pow((Chi/(1 + alpha * pow(Chi,q_exp))),1/n) * 1/vmag;
         _assert_(!xIsNan<IssmDouble>(alpha2));