In [1]:

    

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


    

In [2]:

    

ml = ModelMaq(tmin=0.01, tmax=10)
ls1 = LineSink(ml, -1, 0, 1, 0, tsandQ=[(0, 10)])
ml.solve()


    

    




self.neq  0
No unknowns. Solution complete

In [3]:

    

ml2 = ModelMaq(tmin=0.01, tmax=10)
x = np.arange(-0.9, 1, 0.2)
nwells = len(x)
Qtot = 10
for i in range(nwells):
    Well(ml2, x[i], 0, tsandQ=[(0, Qtot / nwells)])
ml2.solve()


    

    




self.neq  10
solution complete

In [4]:

    

t = np.logspace(-2, 1, 100)
h1a = ml.head(0.5, 1, t)
h1b = ml.head(0.5, 10, t)
h2a = ml2.head(0.5, 1, t)
h2b = ml2.head(0.5, 10, t)


    

In [5]:

    

plt.semilogx(t, h1a[0], 'b', label='line-sink')
plt.semilogx(t, h2a[0], 'r--', label='wells')
plt.semilogx(t, h1b[0], 'b')
plt.semilogx(t, h2b[0], 'r--')
plt.xlabel('time')
plt.ylabel('head')
plt.legend();


    

In [6]:

    

plt.figure(figsize=(12, 4))
ml.contour([-2, 2, -2, 2], [40, 40], t=5)
ml2.contour([-2, 2, -2, 2], [40, 40], t=5);


    

    






<Figure size 864x288 with 0 Axes>






    













    

In [ ]: