Wave distribution in Jangjeon Harbour is numerically simulated for an optimum design of the harbour facilities. A deep-water design wave is estimated based on stochastic extreme wave analysis of wind data in the vicinity of the harbour, and it is applied to the boundary condition at open sea. Boussinesq wave theory that includes effects of frequency dispersion and nonlinearity is employed for the wave simulation. The porosity and sponge layer are adapted at beach to depict partial reflection and complete absorption of waves, respectively. The design wave for breakwater is computed in global domain with coarse grids and the wave distribution inside of wharf is simulated in local domain with fine grids.