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