Test well near lake

import timflow.steady as tfs

ml = tfs.ModelMaq(kaq=[1, 2], z=[10, 5, 4, 0], c=20)
xy = [(-5, 0), (5, 0), (5, 8), (-5, 8)]
p1 = tfs.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 = tfs.Well(ml, xw=0, yw=-10, Qw=100, layers=1)
rf = tfs.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);