A numerical scheme solving the quasi-3D wave-induced circulation is presented. The governing equations have been solve implicitly using a fractional step method in conjunction with the approximate factroization techniques. The equation of each step was discretized by the finite volume scheme which yields more accurate and conservative approximations than the schemes based on finite differences. Examples of computed nearshore current patterns are presented to demonstrate the validity of the model for typical situations through comparison with laboratory experimental data (Gourlay. 1974).