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.5, -2.5, '$k$ = 10 m/d', ha='center', va='center')
text(-0.5, -15, '$k$ = 30 m/d', ha='center', va='center')
text(-0.5, -30, '$k$ = 20 m/d', ha='center', va='center')
text(-0.5, -7.5, '$c$ = 2000 d', ha='center', va='center')
text(-0.5, -22.5, '$c$ = 5000 d', ha='center', va='center')
text(0.5, -2.5, '$S$ = 0.1', ha='center', va='center')
text(0.5, -15, '$S_s$ = 0.0001', ha='center', va='center')
text(0.5, -30, '$S_s$ = 0.0002', ha='center', va='center')
text(0.5, -7.5, '$S_s$ = 0.0001', ha='center', va='center')
text(0.5, -22.5, '$S_s$ = 0.0004', 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.5, -2.5, '$k$ = 10 m/d', ha='center', va='center')
text(-0.5, -7.5, '$k$ = 0.025 m/d', ha='center', va='center')
text(-0.5, -15, '$k$ = 30 m/d', ha='center', va='center')
text(-0.5, -22.5, '$k$ = 0.01 m/d', ha='center', va='center')
text(-0.5, -30, '$k$ = 20 m/d', ha='center', va='center')
text(0.5, -2.5, '$S_s$ = 0.0001', ha='center', va='center')
text(0.5, -7.5, '$S_s$ = 0.0001', ha='center', va='center')
text(0.5, -15, '$S_s$ = 0.0001', ha='center', va='center')
text(0.5, -22.5, '$S_s$ = 0.0001', ha='center', va='center')
text(0.5, -30, '$S_s$ = 0.0001', 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 [ ]: