def main():
# Based on observed uplift
# landward_start = ufloat(0.97, 9.8) # in Ma
landward_max_start = 1.04
landward_min_start = 0.9
landward_max_end = 0.5
landward_min_end = 0
# From Strasser, et al, 2009
seaward_start = 1.95
seaward_end = 1.24
fig = plt.figure()
gs = GridSpec(1, 4)
ax1 = fig.add_subplot(gs[:, 0])
ax2 = fig.add_subplot(gs[:, 1:], sharey=ax1)
# plt.setp(ax2.get_yticklabels(), visible=False)
ax2.yaxis.tick_right()
plot_range(ax1, 1.0, 2, calc.total_oost_percentage())
ax1.axis([1, 2, 0, 1])
ax1.set(xticks=[], yticks=[0, 0.25, 0.5, 0.75, 1.0], yticklabels=["0%", "25%", "50%", "75%", "100%"])
ax1.grid(True, axis="y")
plot_range(ax2, seaward_start, seaward_end, calc.seaward_percentage())
plot_range(ax2, landward_min_start, landward_max_end, calc.landward_percentage())
plot_range(ax2, landward_max_start, landward_min_start, calc.landward_percentage(), alpha=0.5, error=False)
plot_range(ax2, landward_max_end, landward_min_end, calc.landward_percentage(), alpha=0.5, error=False)
ax2.axis([calc.age()[-1], 0, 0, 1])
ax2.set_xlabel("Age (Ma)")
ax2.grid(True, axis="y")
plt.show()
def line_balancing_plot():
heaves = calc.bed_length_balancing()
balancing_shortening = calc.bed_length_shortening()
plate_shortening, smin, smax = calc.total_convergence()
yticklabels = []
plot_uncertain(4, heaves)
yticklabels.append(template('Bed Length Balancing\n{} km', heaves))
plot_uncertain(3, balancing_shortening)
yticklabels.append(
template('Total Shortening from Balancing\n{} km', balancing_shortening)
)
plot_uncertain(2, plate_shortening, hard_min=smin, hard_max=smax)
yticklabels.append(
template('Convergence over {} myr\n{} km', [calc.age()[0], plate_shortening])
)
oost = calc.oost_shortening()
plot_uncertain(1, oost)
yticklabels.append(
template('Predicted Shortening on OOSTS\n{} km', oost)
)
plt.yticks([4, 3, 2, 1], yticklabels)
plt.draw()
plt.tight_layout()
plt.xlabel('Shortening Parallel to Section (km)')
return oost, plate_shortening
def main():
# Based on observed uplift
#landward_start = ufloat(0.97, 9.8) # in Ma
landward_max_start = 1.04
landward_min_start = 0.9
landward_max_end = 0.5
landward_min_end = 0
# From Strasser, et al, 2009
seaward_start = 1.95
seaward_end = 1.24
fig = plt.figure()
gs = GridSpec(1, 4)
ax1 = fig.add_subplot(gs[:, 0])
ax2 = fig.add_subplot(gs[:, 1:], sharey=ax1)
#plt.setp(ax2.get_yticklabels(), visible=False)
ax2.yaxis.tick_right()
plot_range(ax1, 1.0, 2, calc.total_oost_percentage())
ax1.axis([1, 2, 0, 1])
ax1.set(xticks=[],
yticks=[0, 0.25, 0.5, 0.75, 1.0],
yticklabels=['0%', '25%', '50%', '75%', '100%'])
ax1.grid(True, axis='y')
plot_range(ax2, seaward_start, seaward_end, calc.seaward_percentage())
plot_range(ax2, landward_min_start, landward_max_end,
calc.landward_percentage())
plot_range(ax2,
landward_max_start,
landward_min_start,
calc.landward_percentage(),
alpha=0.5,
error=False)
plot_range(ax2,
landward_max_end,
landward_min_end,
calc.landward_percentage(),
alpha=0.5,
error=False)
ax2.axis([calc.age()[-1], 0, 0, 1])
ax2.set_xlabel('Age (Ma)')
ax2.grid(True, axis='y')
plt.show()
