In [1]:

    

%matplotlib inline
import numpy as np
import matplotlib.pyplot as plt
from ttim import *


    

In [2]:

    

ml = ModelMaq(kaq=10, z=[10, 0], Saq=1e-4, tmin=0.01, tmax=10)
w = Well(ml, 0, 20, rw=0.3, tsandQ=[(0, 100)], layers=0)
ml.solve()


    

    




self.neq  1
solution complete

In [3]:

    

t = np.logspace(-2, 1, 100)
h = w.headinside(t)


    

In [4]:

    

plt.plot(t, h[0])


    

    
Out[4]:





[<matplotlib.lines.Line2D at 0x125fd4b90>]






    

In [5]:

    

ml.contour(win=[-30, 50, -10, 30], ngr=40, t=1, labels=True, decimals=2)


    

In [6]:

    

ml2 = ModelMaq(kaq=10, z=[10, 0], Saq=1e-4, tmin=0.01, tmax=10)
w2 = Well(ml2, 0, 20, rw=0.3, tsandQ=[(0, 100)])
wall = LeakyLineDoubletString(ml2, xy=[(-20, 0), (20, 0), (40, 20)], res='imp', order=5, layers=0)
ml2.solve()


    

    




self.neq  13
solution complete

In [7]:

    

h2 = w2.headinside(t)


    

In [8]:

    

plt.plot(t, h[0], label='no wall')
plt.plot(t, h2[0], label='wall')
plt.legend()


    

    
Out[8]:





<matplotlib.legend.Legend at 0x1288b9ed0>






    

In [9]:

    

ml2.contour(win=[-30, 50, -10, 30], ngr=40, t=1, labels=True, decimals=2)