In [1]:

    

from timml import *
from pylab import *
%matplotlib notebook


    

In [2]:

    

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

ml.contour([-10, 10, -15, 10], 50, [0, 1], 30)


    

    




Number of elements, Number of equations: 13 , 81
.............
solution complete






    

In [3]:

    

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

ml.contour([-10, 10, -5, 10], 50, [0], 30)
ml.contour([-10, 10, -5, 10], 50, [1], 30)


    

    




Number of elements, Number of equations: 13 , 129
.............
solution complete






    













    

In [4]:

    

for ls in ls1.lslist:
    for i in range(ls.ncp):
        print(ml.head(ls.xc[i], ls.yc[i], ls.layers))


    

    




[-2.28658869]
[-2.28658869]
[-2.28658869]
[-2.28658869]
[-2.28658869]
[-2.28658869]
[-2.28658869]
[-2.28658869]

In [5]:

    

ml = ModelMaq(kaq=[1, 2], z=[10, 5, 4, 0], c=20)
xy = [(-5, 0), (0, 0), (5, 0), (5, 8), (-5, 8)]
p1 = PolygonInhomMaq(ml, xy=xy, kaq=[1, 2], z=[11,10,5,4,0], 
                     c=[2, 2], topboundary='semi', hstar=1.0, order=5, ndeg=2)
ls1 = LineSinkDitchString(ml, [(0, -4), (0, 0), (0, 4)], 100, \
                          order=3, layers=[1])
rf = Constant(ml, xr=0, yr=-100, hr=2)
ml.solve()

ml.contour([-10, 10, -5, 10], 51, [0], 30)
ml.contour([-10, 10, -5, 10], 51, [1], 30)


    

    




Number of elements, Number of equations: 13 , 129
.............
solution complete






    













    

In [ ]: