In this study, the numerical model for the determination of transformations of waves while propagating has been presented. This numerical model was developed to solve the extended mild slope equation that is applicable to the rapidly varying topographies. It includes the effects of wave refraction, diffraction, shoaling, reflection, harbor resonance, higher order bottom configurations; dissipative terms due to wave breaking and bottom friction. Nonlinear wave celerity and group velocity were introduced in the solution to obtain results that are more accurate. Mac Cormack Method and Point Gauss Seidel Method were applied together in the proposed new solution approach. The numerical model was tested on the semicircular shoaling area  and shoreparallel breakwater . The comparison of the numerical model in the current study and the physical experiments that are present in the literature shows the reliability of the model for wave transformations and dissipations over uneven bottoms.