The equations for the wave orthogonals and wave heights are presented along the depth using time as the independent variable, The effect of wave refracion is included. In this paper, a numerical procedure is described for the calculation and two analytical solutions are presneted in case of straight and parallel bottom contours and horizontal bottom in order to check the computed results. The computer outputs are in the form of automatically plotted wave orthogonals and wave heights along the wave path. In practical calculation of wave heights, Samchuckk Pohang, and Hanlim are selected as plane sloping bottoms.