source: issm/oecreview/Archive/26740-27031/ISSM-26865-26866.diff@ 27032

../trunk-jpl/src/c/classes/Materials/Matice.cpp

@@ -393 +393 @@
                 viscosity=(1.-D)*B/(2.*E);
         }
         else{
-
-                /*if no strain rate, return maximum viscosity*/
-                if(eps_eff==0.){
-                        viscosity = 1.e+14/2.;
-                        //viscosity=2.5*pow(10.,17);
-                }
-
+                /*if strain rate is 0., it is probably our initial guess, use an average
+                 * viscosity instead of a large one*/
+                if(eps_eff==0.) viscosity = 1.e+14/2.;
                 else{
+                        /*if no strain rate, return maximum viscosity*/
+                        //if(eps_eff<1.e-6) eps_eff = 1e-6;
                         viscosity=(1.-D)*B/(2.*pow(E,1./n)*pow(eps_eff,(n-1.)/n));
                 }
         }