In this study, we carry out the numerical simulation to trace the yearly shoreline change of Mang-Bang beach, which is suffering from erosion problem. We obtain the basic equation (One Line Model for shoreline) for the numerical simulation by assuming that the amount of shoreline retreat or advance is balanced by the net influx of longshore and cross-shore sediment into the unit discretized shoreline segment. In doing so, the energy flux model for the longshore sediment transport rate is also evoked. For the case of cross sediment transport, the modified Bailard's model (1981) by Cho and Kim (2019) is utilized. At each time step of the numerical simulation, we adjust a closure depth according to pertinent wave conditions based on the Hallermeier's analytical model (1978) having its roots on the Shield's parameter. Numerical results show that from 2017.4.26 to 2017.10.15 during which swells are prevailing, a shoreline advances due to the sustained supply of cross-shore sediment. It is also shown that a shoreline temporarily retreats due to the erosion by the yearly highest waves sequentially occurring from mid-October to the end of October, and is followed by gradual recovery of shoreline as high waves subdue and swells prevail. It is worth mentioning that great yearly circulation of shoreline completes when a shoreline retreats due to the erosion by the higher waves occurring from mid-March to the end of March. The great yearly circulation of shoreline mentioned above can also be found in the measured locations of shoreline on 2017.4.5, 2017.9.7, 2017.11.7, 2018.3.14. However, numerically simulated amount of shoreline retreat or advance is more significant than the physically measured one, and it should be noted that these discrepancies become more substantial for the case of RUN II where a closure depth is sustained to be as in the most morphology models like the Genesis (Hanson and Kraus, 1989).