In [1]:

    

%matplotlib inline
import numpy as np
from ttim import *


    

In [2]:

    

ml = Model3D(kaq=[2, 1, 5, 10, 4], z=[10, 8, 6, 4, 2, 0], \
             Saq=[.1, .0001, .0002, .0002, .0001], \
             phreatictop=True, kzoverkh=0.1, topboundary='conf', \
             tmin=1, tmax=10)
w = Well(ml, -25, 0, rw=.3, tsandQ=[(0, 100)], layers=[2, 3])

dxdy = 30 * np.pi / 100 * np.cos(np.pi)
x1 = np.arange(-150, -100, 10)
y1 = (x1 + 100) * dxdy
x2 = np.linspace(-100, 100, 20)
y2 = 30 * np.sin(x2 * np.pi / 100)
yls = np.hstack((y1, y2))
xls = np.hstack((x1, x2))
xy = np.array([xls, yls]).T

ls1 = HeadLineSinkString(ml, xy=xy, tsandh='fixed', layers=[0, 1])


    

In [3]:

    

ml.solve()


    

    




self.neq  50
solution complete

In [4]:

    

ml.contour(win=[-150, 150, -50, 50], ngr=40, t=1, layers=1)


    

In [5]:

    

ml.contour(win=[-150, 150, -50, 50], ngr=40, t=1, layers=2)