In [1]:

from __future__ import division, print_function
from pylab import *
%matplotlib notebook
from semi_interface import *

In [2]:

schead = SemiCoastHead(k=42.5, H=160, c=100000, x=-1000, h=10,
                       rhof=1000, rhos=1025, Ls=120000,
                       ztop= -150, sealevel=0)
schead.plot()
print('toe of interface at:', schead.toe())
print('tip of interface at:', schead.tip())

toe of interface at: [ 12228.00685973]
tip of interface at: [ 76102.88455042]

In [3]:

schead = SemiCoastHead(k=3, H=16, c=200000, x=-1900, h=8.5,
                       rhof=1000, rhos=1025, Ls=50000,
                       ztop= -60, sealevel=0)
schead.plot()
print('toe of interface at:', schead.toe())
print('tip of interface at:', schead.tip())

toe of interface at: [ 7684.10708229]
tip of interface at: [ 15273.5734667]

In [4]:

schead = SemiCoastHead(k=10, H=400, c=400000, x=-5000, h=20.5,
                       rhof=1000, rhos=1025, Ls=30000,
                       ztop= -400, sealevel=0)
schead.plot()
print('toe of interface at:', schead.toe())
print('tip of interface at:', schead.tip())

toe of interface at: [-2665.60068303]
tip of interface at: 30000.0

In [5]:

schead = SemiCoastHead(k=45, H=150, c=31250, x=-5000, h=3.5,
                       rhof=1000, rhos=1025, Ls=50000,
                       ztop= -50, sealevel=0)
schead.plot()
print('toe of interface at:', schead.toe())
print('tip of interface at:', schead.tip())

toe of interface at: [-24396.13913519]
tip of interface at: [ 23641.31305418]

In [6]:

schead = SemiCoastHead(k=50, H=40, c=750, x=-8000, h=11,
                       rhof=1000, rhos=1025, Ls=13000,
                       ztop= -160, sealevel=0)
schead.plot()
print('toe of interface at:', schead.toe())
print('tip of interface at:', schead.tip())

toe of interface at: [ 106.94813568]
tip of interface at: [ 3106.94813568]