This study describes the application of a two dimensional depth integrated numerical model. The explict scheme of finite difference method had been applied to the model of circulation. The nonlinear terms showed a slight difference for the variations of water elevation when calculated grid was small. They were also found to be minor when calculated grid size was increased. For verification of the numerical model, numerical results were compared with predicted values and field data. In the model, the effect of nonlinear advective terms proved not to be significant.