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