Set small number to 1e-3 so the model is stable when w=wsat

classmodel · Mar 10, 2015 · dafad15 · dafad15
1 parent 5870f4b
commit dafad15
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/model.cpp b/model.cpp
@@ -1187,7 +1187,7 @@ void model::runlsmodel()
     d1         = 0.1;
 
     C1         = C1sat * pow(wsat / wg, b/2. + 1.);
-    const double wsmall = 1e-6;
+    const double wsmall = 1e-3;
     C2         = C2ref * (w2 / (wsat - w2 + wsmall));
 
     wgeq       = w2 - wsat * a * ( pow(w2 / wsat, p) * (1. - pow(w2 / wsat,8.*p)) );