A numerical method is applied to the problem of an incompressible fluid in a slowly rotating rectangular basin for the simulation of wave propagation in shallow water. The present work is a complete study of the wave motion through evaluation of the wave height and the velocity components. The results are found by the application of a fractional step method and illustrated graphically. The technique is applied by splitting the shallow water equations and successive integration in every direction along the characteristics using the Riemann invariants associated with cubic spline interpolation. It has the advantage of reducing the multidimensional matrix inversion problem into an equivalent one-dimensional problem. Numerical results are represented in three dimensions for the velocity components at different times. The distribution of temperature and concentration are also calculated and plotted.