apply black+isort; add imports and docs

examples/01_call_theis.py

@@ -7,6 +7,7 @@
 
 Reference: `Theis 1935 <https://doi.org/10.1029/TR016i002p00519>`__
 """
+
 import numpy as np
 from matplotlib import pyplot as plt
 

examples/02_call_ext_theis2d.py

@@ -12,6 +12,7 @@
 
 Reference: `Zech et. al. 2016 <http://dx.doi.org/10.1002/2015WR018509>`__
 """
+
 import numpy as np
 from matplotlib import pyplot as plt
 
@@ -51,12 +52,8 @@
     label_TG = "Theis($T_G$)" if i == 0 else None
     label_TH = "Theis($T_H$)" if i == 0 else None
     label_ef = "extended Theis" if i == 0 else None
-    plt.plot(
-        rad, head_TG[i], label=label_TG, color="C" + str(i), linestyle="--"
-    )
-    plt.plot(
-        rad, head_TH[i], label=label_TH, color="C" + str(i), linestyle=":"
-    )
+    plt.plot(rad, head_TG[i], label=label_TG, color="C" + str(i), linestyle="--")
+    plt.plot(rad, head_TH[i], label=label_TH, color="C" + str(i), linestyle=":")
     plt.plot(rad, head_ef[i], label=label_ef, color="C" + str(i))
     time_ticks.append(head_ef[i][-1])
 

examples/03_call_ext_theis3d.py

@@ -14,6 +14,7 @@
 
 Reference: `Müller 2015 <http://dx.doi.org/10.13140/RG.2.2.34074.24002>`__
 """
+
 import numpy as np
 from matplotlib import pyplot as plt
 
@@ -59,12 +60,8 @@
     label_TG = "Theis($K_{efu}$)" if i == 0 else None
     label_TH = "Theis($K_H$)" if i == 0 else None
     label_ef = "extended Theis 3D" if i == 0 else None
-    plt.plot(
-        rad, head_Kefu[i], label=label_TG, color="C" + str(i), linestyle="--"
-    )
-    plt.plot(
-        rad, head_KH[i], label=label_TH, color="C" + str(i), linestyle=":"
-    )
+    plt.plot(rad, head_Kefu[i], label=label_TG, color="C" + str(i), linestyle="--")
+    plt.plot(rad, head_KH[i], label=label_TH, color="C" + str(i), linestyle=":")
     plt.plot(rad, head_ef[i], label=label_ef, color="C" + str(i))
     time_ticks.append(head_ef[i][-1])
 

examples/04_call_ext_theis_tpl.py

@@ -13,6 +13,7 @@
 
 Reference: (not yet published)
 """
+
 import numpy as np
 from matplotlib import pyplot as plt
 
@@ -62,12 +63,8 @@
     label_TG = "Theis($K_G$)" if i == 0 else None
     label_TH = "Theis($K_H$)" if i == 0 else None
     label_ef = "extended Theis TPL 2D" if i == 0 else None
-    plt.plot(
-        rad, head_KG[i], label=label_TG, color="C" + str(i), linestyle="--"
-    )
-    plt.plot(
-        rad, head_KH[i], label=label_TH, color="C" + str(i), linestyle=":"
-    )
+    plt.plot(rad, head_KG[i], label=label_TG, color="C" + str(i), linestyle="--")
+    plt.plot(rad, head_KH[i], label=label_TH, color="C" + str(i), linestyle=":")
     plt.plot(rad, head_ef[i], label=label_ef, color="C" + str(i))
     time_ticks.append(head_ef[i][-1])
 

examples/05_call_neuman2004.py

@@ -14,6 +14,7 @@
 
 Reference: `Neuman 2004 <https://doi.org/10.1029/2003WR002405>`__
 """
+
 import numpy as np
 from matplotlib import pyplot as plt
 
@@ -61,12 +62,8 @@
     label_TG = "Theis($T_G$)" if i == 0 else None
     label_TH = "Theis($T_H$)" if i == 0 else None
     label_ef = "transient Neuman [2004]" if i == 0 else None
-    plt.plot(
-        rad, head_TG[i], label=label_TG, color="C" + str(i), linestyle="--"
-    )
-    plt.plot(
-        rad, head_TH[i], label=label_TH, color="C" + str(i), linestyle=":"
-    )
+    plt.plot(rad, head_TG[i], label=label_TG, color="C" + str(i), linestyle="--")
+    plt.plot(rad, head_TH[i], label=label_TH, color="C" + str(i), linestyle=":")
     plt.plot(rad, head_ef[i], label=label_ef, color="C" + str(i))
     time_ticks.append(head_ef[i][-1])
 

examples/06_compare_extthiem2d_grfsteady.py

@@ -11,6 +11,7 @@
 - `Schneider & Attinger 2008 <https://doi.org/10.1029/2007WR005898>`__
 - `Zech & Attinger 2016 <https://doi.org/10.5194/hess-20-1655-2016>`__
 """
+
 import numpy as np
 from matplotlib import pyplot as plt
 

examples/07_compare_extthiem3d_grfsteady.py

@@ -8,6 +8,7 @@
 
 Reference: `Zech et. al. 2012 <https://doi.org/10.1029/2012WR011852>`__
 """
+
 import numpy as np
 from matplotlib import pyplot as plt
 

examples/08_compare_extthiem2d_neuman.py

@@ -13,6 +13,7 @@
 - `Neuman 2004 <https://doi.org/10.1029/2003WR002405>`__
 - `Zech & Attinger 2016 <https://doi.org/10.5194/hess-20-1655-2016>`__
 """
+
 import numpy as np
 from matplotlib import pyplot as plt
 

examples/09_compare_exttheis2d_neuman.py

@@ -13,6 +13,7 @@
 - `Neuman 2004 <https://doi.org/10.1029/2003WR002405>`__
 - `Zech et. al. 2016 <http://dx.doi.org/10.1002/2015WR018509>`__
 """
+
 import numpy as np
 from matplotlib import pyplot as plt
 
@@ -40,9 +41,7 @@
     plt.plot(rad, head2[i], label=label2, color="C" + str(i), linestyle="--")
     time_ticks.append(head1[i][-1])
 
-plt.title(
-    "$T_G={}$, $\sigma^2={}$, $\ell={}$, $S={}$".format(TG, var, len_scale, S)
-)
+plt.title("$T_G={}$, $\sigma^2={}$, $\ell={}$, $S={}$".format(TG, var, len_scale, S))
 plt.xlabel("r in [m]")
 plt.ylabel("h in [m]")
 plt.legend()

examples/10_convergence_ext_theis_tpl.py

@@ -7,6 +7,7 @@
 
 Reference: (not yet published)
 """
+
 import numpy as np
 from matplotlib import pyplot as plt
 
@@ -25,9 +26,7 @@
 S = 1e-4  # storativity
 rate = -1e-4  # pumping rate
 
-head1 = ext_thiem_tpl(
-    rad, r_ref, KG, len_scale, hurst, var, dim=dim, rate=rate
-)
+head1 = ext_thiem_tpl(rad, r_ref, KG, len_scale, hurst, var, dim=dim, rate=rate)
 head2 = ext_theis_tpl(
     time, rad, S, KG, len_scale, hurst, var, dim=dim, rate=rate, r_bound=r_ref
 )

examples/11_convergence_ext_grf.py

@@ -9,6 +9,7 @@
 
 Reference: `Barker 1988 <https://doi.org/10.1029/WR024i010p01796>`__
 """
+
 import numpy as np
 from matplotlib import pyplot as plt
 
@@ -30,9 +31,7 @@
 
 plt.plot(rad, head1, label="Ext GRF steady")
 plt.plot(rad, head2, label="Ext GRF (t={})".format(time), linestyle="--")
-plt.plot(
-    rad, head3, label="GRF quasi-steady (t={})".format(time), linestyle=":"
-)
+plt.plot(rad, head3, label="GRF quasi-steady (t={})".format(time), linestyle=":")
 
 plt.xlabel("r in [m]")
 plt.ylabel("h in [m]")

examples/12_compare_theis_quasi_steady.py

@@ -5,6 +5,7 @@
 The quasi steady is reached, when the radial shape of the drawdown in not
 changing anymore.
 """
+
 import numpy as np
 from matplotlib import pyplot as plt
 
@@ -21,12 +22,8 @@
     transmissivity=1e-4,
     rate=-1e-4,
 )
-head1 = (
-    theis(time, rad, storage=1e-3, transmissivity=1e-4, rate=-1e-4) - head_ref
-)
-head2 = theis(
-    time, rad, storage=1e-3, transmissivity=1e-4, rate=-1e-4, r_bound=r_ref
-)
+head1 = theis(time, rad, storage=1e-3, transmissivity=1e-4, rate=-1e-4) - head_ref
+head2 = theis(time, rad, storage=1e-3, transmissivity=1e-4, rate=-1e-4, r_bound=r_ref)
 head3 = thiem(rad, r_ref, transmissivity=1e-4, rate=-1e-4)
 
 for i, step in enumerate(time):

examples/13_self_defined_transmissivity.py

@@ -15,6 +15,7 @@
 
 Reference: (not yet published)
 """
+
 import matplotlib.gridspec as gridspec
 import numpy as np
 from matplotlib import pyplot as plt
@@ -78,9 +79,7 @@ def cond(rad, K_far, K_well, len_scale):
 
 rad_lin = np.linspace(rad[0], rad[-1], 1000)
 ax1.plot(rad_lin, step_f(rad_lin, R_part, K_part), label="step Conductivity")
-ax1.plot(
-    rad_lin, cond(rad_lin, K_far, K_well, len_scale), label="Conductivity"
-)
+ax1.plot(rad_lin, cond(rad_lin, K_far, K_well, len_scale), label="Conductivity")
 ax1.set_yticks([K_well, K_far])
 ax1.set_ylabel(r"$K$ in $[\frac{m}{s}]$")
 plt.setp(ax1.get_xticklabels(), visible=False)

examples/14_interval_theis.py

@@ -8,6 +8,7 @@
 Unfortunatly the Stehfest algorithm is not suitable for this kind of solution,
 which is demonstrated in the following script.
 """
+
 import numpy as np
 from matplotlib import pyplot as plt
 

examples/15_accruing_theis.py

@@ -7,6 +7,7 @@
 This could be interpreted as that the water pump needs a certain time to
 reach its constant rate state.
 """
+
 import matplotlib.gridspec as gridspec
 import numpy as np
 from matplotlib import pyplot as plt

pyproject.toml

@@ -79,11 +79,9 @@ fallback_version = "0.0.0.dev0"
 [tool.isort]
 profile = "black"
 multi_line_output = 3
-line_length = 79
 
 [tool.black]
 exclude = "_version.py"
-line-length = 79
 target-version = ["py37"]
 
 [tool.coverage]

src/anaflow/__init__.py

@@ -39,10 +39,14 @@
    ext_thiem_3d
    ext_thiem_tpl
    ext_thiem_tpl_3d
+   ext_thiem_int
+   ext_thiem_int_3d
    ext_theis_2d
    ext_theis_3d
    ext_theis_tpl
    ext_thiem_tpl_3d
+   ext_theis_int
+   ext_thiem_int_3d
    neuman2004
    neuman2004_steady
 
@@ -78,15 +82,20 @@
    specialrange
    specialrange_cut
 """
+
 from anaflow.flow import (
     ext_grf,
     ext_grf_steady,
     ext_theis_2d,
     ext_theis_3d,
+    ext_theis_int,
+    ext_theis_int_3d,
     ext_theis_tpl,
     ext_theis_tpl_3d,
     ext_thiem_2d,
     ext_thiem_3d,
+    ext_thiem_int,
+    ext_thiem_int_3d,
     ext_thiem_tpl,
     ext_thiem_tpl_3d,
     grf,
@@ -95,13 +104,7 @@
     theis,
     thiem,
 )
-from anaflow.tools import (
-    get_lap,
-    get_lap_inv,
-    specialrange,
-    specialrange_cut,
-    step_f,
-)
+from anaflow.tools import get_lap, get_lap_inv, specialrange, specialrange_cut, step_f
 
 try:
     from anaflow._version import __version__
@@ -117,10 +120,14 @@
     "ext_thiem_3d",
     "ext_thiem_tpl",
     "ext_thiem_tpl_3d",
+    "ext_thiem_int",
+    "ext_thiem_int_3d",
     "ext_theis_2d",
     "ext_theis_3d",
     "ext_theis_tpl",
     "ext_theis_tpl_3d",
+    "ext_theis_int",
+    "ext_theis_int_3d",
     "neuman2004",
     "neuman2004_steady",
     "grf",

src/anaflow/flow/__init__.py

@@ -38,10 +38,14 @@
    ext_thiem_3d
    ext_thiem_tpl
    ext_thiem_tpl_3d
+   ext_thiem_int
+   ext_thiem_int_3d
    ext_theis_2d
    ext_theis_3d
    ext_theis_tpl
    ext_theis_tpl_3d
+   ext_theis_int
+   ext_theis_int_3d
    neuman2004
    neuman2004_steady
 
@@ -56,14 +60,19 @@
    ext_grf
    ext_grf_steady
 """
+
 from anaflow.flow.ext_grf_model import ext_grf, ext_grf_steady
 from anaflow.flow.heterogeneous import (
     ext_theis_2d,
     ext_theis_3d,
+    ext_theis_int,
+    ext_theis_int_3d,
     ext_theis_tpl,
     ext_theis_tpl_3d,
     ext_thiem_2d,
     ext_thiem_3d,
+    ext_thiem_int,
+    ext_thiem_int_3d,
     ext_thiem_tpl,
     ext_thiem_tpl_3d,
     neuman2004,
@@ -78,10 +87,14 @@
     "ext_thiem_3d",
     "ext_thiem_tpl",
     "ext_thiem_tpl_3d",
+    "ext_thiem_int",
+    "ext_thiem_int_3d",
     "ext_theis_2d",
     "ext_theis_3d",
     "ext_theis_tpl",
     "ext_theis_tpl_3d",
+    "ext_theis_int",
+    "ext_theis_int_3d",
     "neuman2004",
     "neuman2004_steady",
     "grf",

src/anaflow/flow/ext_grf_model.py

@@ -9,6 +9,7 @@
    ext_grf
    ext_grf_steady
 """
+
 # pylint: disable=C0103
 import numpy as np
 from scipy.integrate import quad as integ
@@ -200,9 +201,7 @@ def integrand(val):
         if np.isclose(dim, 2):
             res = np.log(r_ref / Input.rad) / con
         else:
-            res = (
-                (r_ref ** (2 - dim) - Input.rad ** (2 - dim)) / (2 - dim) / con
-            )
+            res = (r_ref ** (2 - dim) - Input.rad ** (2 - dim)) / (2 - dim) / con
 
     res = Input.reshape(res)
     # rescale by pumping rate