& v_flux_bdry_SI_1
_RL v_flux_bdry_SI_1
& v_flux_bdry_SI_1(i,j,bi,bj)
& v_flux_bdry_SI_1,
& + (1-fac) * v_flux_bdry_SI_1(i,j,bi,bj)