def forearc_plot(): from process_bootstrap_results import shortening_parallel_to_section from process_bootstrap_results import heave_parallel_to_section # Shortening from fault kinematics is in meters... landward_shortening = shortening_parallel_to_section() / 1000 min_bound = heave_parallel_to_section() / 1000 total_oost = oost_shortening() yticklabels = [] plot_uncertain(3, total_oost) yticklabels.append(template('Total OOST Shortening\n{} km', total_oost)) plot_uncertain(2, landward_shortening, hard_min=min_bound) yticklabels.append( template('Shortening on Landward Branch\nFrom Forearc Uplift {} km', landward_shortening) ) seaward_shortening = total_oost - landward_shortening # Hard minimum from Strasser, et al plot_uncertain(1, seaward_shortening, hard_min=1.9) yticklabels.append( template('Predicted Shortening on \n Seaward Branch {} km', seaward_shortening) ) xmin, xmax = plt.xlim() plt.xlim([0, xmax]) plt.yticks([3, 2, 1], yticklabels) plt.tight_layout() plt.xlabel('Shortening Parallel to Section (km)')
def forearc_plot(): from process_bootstrap_results import shortening_parallel_to_section from process_bootstrap_results import heave_parallel_to_section # Shortening from fault kinematics is in meters... landward_shortening = shortening_parallel_to_section() / 1000 min_bound = heave_parallel_to_section() / 1000 total_oost = calc.oost_shortening() yticklabels = [] plot_uncertain(3, total_oost) yticklabels.append(template('Total OOST Shortening\n{} km', total_oost)) plot_uncertain(2, landward_shortening, hard_min=min_bound) yticklabels.append( template('Shortening on Landward Branch\nFrom Forearc Uplift {} km', landward_shortening) ) seaward_shortening = total_oost - landward_shortening # Hard minimum from Strasser, et al plot_uncertain(1, seaward_shortening, hard_min=1.9) yticklabels.append( template('Predicted Shortening on \n Seaward Branch {} km', seaward_shortening) ) xmin, xmax = plt.xlim() plt.xlim([0, xmax]) plt.yticks([3, 2, 1], yticklabels) plt.tight_layout() plt.xlabel('Shortening Parallel to Section (km)')
def seaward_shortening(): """Shortening accomodated on the seaward branch of the OOSTS based on comparing the total (`oost_shortening()`) shortening with the shortening predicted on the landward branch from forearc uplift. Returns: -------- shortening : a ufloat with 2 sigma error in kilometers. """ from process_bootstrap_results import shortening_parallel_to_section landward_shortening = shortening_parallel_to_section() / 1000 return oost_shortening() - landward_shortening
def landward_percentage(): """ Maximum percentage of total plate convergence accomodated by the landward branch of the OOSTS during its period of activity. Returns: -------- percentage : A ufloat with a 2 sigma error representing a unitless ratio (e.g. multiply by 100 to get percentage). """ from process_bootstrap_results import shortening_parallel_to_section landward_shortening = shortening_parallel_to_section() / 1000 duration = 1.04 - 0.5 # Probably too short... rate = plate_motion() total = duration * rate return landward_shortening / total
def landward_percentage(): """ Maximum percentage of total plate convergence accomodated by the landward branch of the OOSTS during its period of activity. Returns: -------- percentage : A ufloat with a 2 sigma error representing a unitless ratio (e.g. multiply by 100 to get percentage). """ from process_bootstrap_results import shortening_parallel_to_section landward_shortening = shortening_parallel_to_section() / 1000 duration = ufloat(0.97, 0.07) - ufloat(0.25, 0.25) rate = plate_motion() total = duration * rate return landward_shortening / total