Safer check for cut off value

GeoStat-Framework · Jul 30, 2019 · 3f28a20 · 3f28a20
1 parent 46592ec
commit 3f28a20
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/anaflow/flow/heterogeneous.py b/anaflow/flow/heterogeneous.py
@@ -396,7 +396,7 @@ def ext_theis2D(
     # genearte rlast from a given relativ-error to farfield-transmissivity
     rlast = T_CG_error(T_err, TG, sig2, corr, prop, Twell)
     # generate the partition points
-    if rlast > 0.0:
+    if rlast > rwell:
         rpart = specialrange_cut(rwell, rinf, parts + 1, rlast)
     else:
         rpart = np.array([rwell, rinf])
@@ -580,7 +580,7 @@ def ext_theis3D(
     # genearte rlast from a given relativ-error to farfield-conductivity
     rlast = K_CG_error(K_err, KG, sig2, corr, e, prop, Kwell=Kwell)
     # generate the partition points
-    if rlast > 0.0:
+    if rlast > rwell:
         rpart = specialrange_cut(rwell, rinf, parts + 1, rlast)
     else:
         rpart = np.array([rwell, rinf])