A model is developed to predict the long-term beach evolution near the long groin considering the combined effects of variation of sea level, wave refraction and diffraction. A numerical solution for this problem is solved by considering the equation as a system subject to the boundary condition for longshore transport rate. One possible method is the centered Crank-Nicolson type implicit scheme. The results which ard obtained by applying this numerical model at Songdo beach, Pohang are as follows. Owing to the approximation used in the calculation of the refraction and diffraction coefficients, the discrepancy between the predicted and actual shoreline occurs to the interior of long groin. However, the shape of shoreline at the exterier of long groins agrees well.