er_w_time_array_A[:,i] = pt.make_slip_rate_w_time(initA, accel, erA1, erA2, time_vector = time_vector) er_w_time_array_B[:,i] = pt.make_slip_rate_w_time(initB, accel, erB1, erB2, time_vector = time_vector) cum_ext_w_time_array_A[:,i] = pt.get_cum_vector(er_w_time_array_A[:,i], time_step) cum_ext_w_time_array_B[:,i] = pt.get_cum_vector(er_w_time_array_B[:,i], time_step) er_w_time_array = er_w_time_array_A + er_w_time_array_B cum_ext_w_time_array = cum_ext_w_time_array_A + cum_ext_w_time_array_B # make plots fig1 = plt.figure(1) pt.make_fault_histograms(fig1, filtered_runs, initA_col, srA1_col, accel_col, srA2_col) fig2 = plt.figure(2) pt.make_fault_histograms(fig2, filtered_runs, initB_col, srB1_col, accel_col, srB2_col) fig3 = plt.figure(3) pt.make_ext_histories(fig3, 'nlrT4', time_vector, er_w_time_array, cum_ext_w_time_array)