In [1]:

    

%matplotlib inline
from timml import *
from pylab import *
ml = ModelMaq(kaq=[1, 2], z=[10, 5, 4, 0], c=20)
xy = [(-5,0), (5,0), (5,8), (-5,8)]
p1 = 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 = Well(ml, xw=0, yw=-10, Qw=100, layers=1)
rf = Constant(ml, xr=0, yr=-100, hr=2)
ml.solve()


    

    




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

In [2]:

    

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


    

In [ ]: