ax[jj].set_ylim(ymin - 0.11 * delta_y, ymax + 0.11 * delta_y)