In [1]:

    

from pylab import *
%matplotlib inline


    

In [2]:

    

# ModelMaq figure
figure()
axes(frameon = 0)
grey = [.9, .9, .9]
plot([-1, 1], [0, 0], 'k', lw=2)
axhspan(-5, -10, color=grey)
axhspan(-20, -25, color=grey)
plot([-1, 1], [-35, -35], 'k', lw=2)
text(0, -2.5, '$k$ = 10 m/d', ha='center', va='center')
text(0, -15, '$k$ = 30 m/d', ha='center', va='center')
text(0, -30, '$k$ = 20 m/d', ha='center', va='center')
text(0, -7.5, '$c$ = 2000 d', ha='center', va='center')
text(0, -22.5, '$c$ = 5000 d', ha='center', va='center')
xlim(-1, 1)
yticks([0, -5, -10, -20, -25, -35])
ylabel('elevation (m)')
xticks([])
savefig('../docs/models/modelmaq.png', bbox_inches='tight')
#ModelMaq(kaq=[10, 30, 20], z=[0, -5, -10, -20, -25, -35], c=[2000, 5000])


    

In [3]:

    

# ModelMaq figure
figure()
axes(frameon = 0)
grey = [.9, .9, .9]
plot([-1, 1], [0, 0], 'k', lw=2)
axhspan(-5, -5, color=grey)
axhspan(-10, -10, color=grey)
axhspan(-20, -20, color=grey)
axhspan(-25, -25, color=grey)
plot([-1, 1], [-35, -35], 'k', lw=2)
text(0, -2.5, '$k$ = 10 m/d', ha='center', va='center')
text(0, -7.5, '$k$ = 0.025 m/d', ha='center', va='center')
text(0, -15, '$k$ = 30 m/d', ha='center', va='center')
text(0, -22.5, '$k$ = 0.01 m/d', ha='center', va='center')
text(0, -30, '$k$ = 20 m/d', ha='center', va='center')
xlim(-1, 1)
yticks([0, -5, -10, -20, -25, -35])
ylabel('elevation (m)')
xticks([])
savefig('../docs/models/model3d.png', bbox_inches='tight')
#Model3D(kaq=[10, 0.0025, 30, 0.001, 20], z=[0, -5, -10, -20, -25, -35], kzoverkh=0.1)


    

In [4]:

    

# Model
figure()
axes(frameon = 0)
grey = [.9, .9, .9]
plot([-1, 1], [0, 0], 'k', lw=2)
axhspan(-5, -5, color=grey)
axhspan(-10, -10, color=grey)
axhspan(-20, -25, color=grey)
plot([-1, 1], [-35, -35], 'k', lw=2)
text(0, -2.5, '$k$ = 10 m/d', ha='center', va='center')
text(0, -7.5, '$k$ = 5 m/d', ha='center', va='center')
text(0, -15, '$k$ = 30 m/d', ha='center', va='center')
text(0, -22.5, '$c$ = 2000 d', ha='center', va='center')
text(0, -30, '$k$ = 20 m/d', ha='center', va='center')
xlim(-1, 1)
yticks([0, -5, -10, -20, -25, -35])
ylabel('elevation (m)')
xticks([])
savefig('../docs/models/model.png', bbox_inches='tight')
#Model(kaq=[10, 5, 30, 20], c=[2, 5, 2000], z=[0, -5, -10, -20, -25, -35], 
# npor=[0.3, 0.3, 0.3, 0.3], ltype=['a', 'a', 'a', 'l', 'a'])


    

In [5]:

    

ls


    

    




Untitled.ipynb                 timml_notebook0_sol.ipynb
figs/                          timml_notebook1_sol.ipynb*
lake_horizontal_well.ipynb     timml_notebook2_sol.ipynb*
test_linesink_discharge.ipynb  timml_notebook3_sol.ipynb*
test_well_near_lake.ipynb      timml_notebook4_sol.ipynb
timml_figures.ipynb            timml_notebook5_sol.ipynb

In [ ]: