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()