Based on the moving least squares method and the finite difference scheme, a new element\|free method is proposed and used to simulate horizontal two dimension river flow with complicated boundary. Reasonable agreement is found between the measured data and the model results.