Added a divergence correction for the boundary, which physically corresponds to a vertical outflow where the divergence is non-zero