A natural and an artificial harbor can exhibit frequency (or period) dependent water surface oscillations when excited by incident waves. Such oscillations in harbors can cause significant damages to moored ships and adjacent structures. This can also induce undesirable current in harbor. Many previous investigators have studied various aspects of harbor resonance problem. In the present paper, both a localized finite element method(LFEM) which is based on the functional constructed by Chen & Mei(1974) and Bai & Yeung(1974) and an integral equation method which was used by Lee(1969) are applied to harbor resonance problem. The LFEM shows computationally more efficient than the integral equation method. Our test results show a good agreement compared with other results. In the present computations, specifically two harbor geometris are treated here. The present method by LFEM can be extended to a fully three dimensional harbor problem.