A numerical model of low frequency waves is presented. The model is based on that of Roelvink (1993), but the numerical techniques used in the solution are based on the so-called Weighted-Average Flux (WAF) method withTime-Operator-Splitting (TOS) used for the treatment of the source terms. This method allows a small number ofcomputational points to be used, and is particularly efficient in modeling wave setup. The short wave (or primary wave)energy equation is solved with a traditional Lax-Wendroff technique. A nonlinear wave theory is introduced. The modeldescribed in this paper is found to be satisfactory in modeling low frequency waves associated with incident bichromaticwaves.