An sph numerical scheme for 2D shallow water equations with bathymetry and open boundaries