& v_flux_bdry_SI_0,
_RL v_flux_bdry_SI_0
v_flux_bdry_SI_0(i,j,bi,bj) =
& fac * v_flux_bdry_SI_0(i,j,bi,bj)