Changeset 18990

Timestamp:

01/06/15 10:06:09 (10 years ago)

Author:

bdef

Message:

BUG:adding zero velocity check in friction hydro

File:

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

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

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