bugfix: typo for r_well>0 in grf-laplace

GeoStat-Framework · Mar 7, 2019 · c41309e · c41309e
1 parent 99dfcbc
commit c41309e
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/anaflow/flow/laplace.py b/anaflow/flow/laplace.py
@@ -132,7 +132,7 @@ def grf_laplace(
     diff_sr0 = np.sqrt(Spart[0] / Kpart[0])
     # set the general pumping-condtion depending on the well-radius
     if rpart[0] > 0.0:
-        Qs = -s ** (-1.5) / diff_sr0 * rpart[0] ** (1 - nu)
+        Qs = -s ** (-1.5) / diff_sr0 * rpart[0] ** (nu - 1)
     else:
         Qs = (2 / diff_sr0) ** nu / gamma(1 - nu) * s ** (-nu / 2 - 1)