Changeset 18781

Timestamp:

11/13/14 14:55:38 (10 years ago)

Author:

bdef

Message:

CHG:dealing with q exponent 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

@@ -218 +218 @@
 
         //compute alpha and Chi coefficients: */
-        alpha=(pow(q_exp-1,q_exp-1))/pow(q_exp,q_exp);
+        if (q_exp==1){
+                alpha=1;
+        }
+        else{
+                alpha=(pow(q_exp-1,q_exp-1))/pow(q_exp,q_exp);
+        }
+        
         Chi=vmag/(pow(C_param,n)*pow(Neff,n)*As);