Wave propagation is changed by the effect of shoaling, current-depth refraction and shelter-ing etc. To solve these problems. numerous models have been developed. In the present study, a coordinate system is proposed based on the wave ray equation with the wave number equation including diffraction effects . The governing equation for the study was derived from the mild slope wave equation in non-steady state, including current effects (Kirby, 1986a) and trans-formed into an orthogonal curvilinear coordinate system on the basis of the wave ray equation. To obtain a numerical solution, an explicit finite difference scheme was used, and solved by the relaxation method. This model was tested for various cases: Firstly a submersed circular shoal and a constant unit depth. Secondly a submerged elliptic shoal on a slope, and finally a breakwater harbour with obliquely incident waves on a slope. The model was found to simulate the experimental results and other theoretical results in wave height and wave angle fairy well, and the applicability of the model around an arbitrary shaped coastal structure was also verified. To demonstrate the general usefullness of the present approach , the model is to be applied to a field situation with a complex bed topography.