def shelf1(r):
"""
Ocean followed by continental slope, continental shelf, and beach.
The ocean is flat, the slope is cubic, and the shelf and beach are linear.
"""
rad1 = 1500e3 # beginning of continental slope
rad2 = rad1 + 60e3 # end of slope, start of shelf
rad3 = rad2 + 80e3 # end of shelf, start of beach
rad4 = rad3 + 5e3 # radius of shoreline where z=0
z1 = -4000. # depth from r=0 to rad1 (ocean)
z2 = -100. # depth at r=rad2 (start of shelf)
z3 = -100. # depth at r=rad3 (start of beach)
z4 = 0. # depth at shoreline
xi = array([0., rad1, rad2, rad3, rad4])
zl = array([z1, z1, z2, z3, z4])
zr = zl # continuous!
slope_of_shelf = (z3 - z2) / (rad3 - rad2)
slope_of_beach = (z4 - z3) / (rad4 - rad3)
print "Slope of shelf = ", slope_of_shelf
print "Slope of beach = ", slope_of_beach
slopel = array([0., 0., slope_of_shelf, slope_of_shelf, slope_of_beach])
sloper = array([0., 0., slope_of_shelf, slope_of_beach, slope_of_beach])
z = pwcubic(xi, zl, zr, slopel, sloper, r)
return z
def island1(r):
"""
Island created using piecewise cubic with radius rad with zero slope at peak and
base. Raises topo by ztop at the peak.
"""
rad = 30e3
ztop = 120. # 20m above sealevel on the 100m deep shelf.
xi = array([0., rad])
zl = array([ztop, 0.])
zr = zl
slopel = array([0., 0.])
sloper = slopel
z = pwcubic(xi, zl, zr, slopel, sloper, r)
return z