Test well near lake

import timml as tml

ml = tml.ModelMaq(kaq=[1, 2], z=[10, 5, 4, 0], c=20)
xy = [(-5, 0), (5, 0), (5, 8), (-5, 8)]
p1 = tml.PolygonInhomMaq(
    ml,
    xy=[(-5, 0), (5, 0), (5, 8), (-5, 8)],
    kaq=[0.2, 8],
    z=[11, 10, 5, 4, 0],
    c=[2, 20],
    topboundary="semi",
    hstar=1.0,
    order=3,
    ndeg=1,
)
w = tml.Well(ml, xw=0, yw=-10, Qw=100, layers=1)
rf = tml.Constant(ml, xr=0, yr=-100, hr=2)
ml.solve()

Number of elements, Number of equations: 11 , 65
.
.
.
.
.
.
.
.
.
.
.

solution complete

ml.plots.contour([-10, 10, -10, 10], 50, [0, 1], 20)

[<matplotlib.contour.QuadContourSet at 0x73700030b850>,
 <matplotlib.contour.QuadContourSet at 0x736ffe1ac590>]