The development of a simplified 2-D numerical model was described for wind-driven circulation in reservoir using standard k-ε turbulence model to specify eddy viscosity distribution.The goving equations are transformed and solved on variable vertical grids, which allows refinement at the surface and bottom boundaries.The results of the model simulation for flow are compared with analytical solutions for laminar and turbulent flows, experimental data in a wind-flume and wind wave tank.The sensitivity analysis results show that use of large number of depth layers increases the accuracy for the bottom counter-current flow.Prediction of surface drift was not very sensitive to surface grid refinement.The model was also applied to Baisha reservoir for an assumed wind condition and showed to be able to simulate the general features of surface drift and ret flow under variable flow depth.The model can serve as altative means of studying wind-driven flow beside experiments.It also reduced the problem complexity associated with 3-D circulation models while faithfully reproducing the drift and near bottom ret currents.