axes1[1].plot([0, time[-1] / 3600], [avg, avg], label='Avg: {}'.format(exp_avg_prot))