The calmness inside a harbor plays an important role in the appropriate disposition of harbor structures. However, it is not easy to acquire accurate computational results because these are affected by many factors concerned with wave transformation. Recently, numerical model tests, which are quicker and more economical than hydraulic model experiments, were carried out for the purpose of analyzing wave height distributions in harbors. This paper presents a numerical model that is able to calculate wave heights inside a harbor. It is based on a time-dependent mild slope involving wave refraction, diffraction, shoaling effect, and reflection. In particular, arbitrary reflectivity is used at the boundary in order to simulate the real harbor reflection condition. The proposed numerical model is applied to Samcheon new-harbor in order to investigate harbor calmness.