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