A two-dimensional depth-integrated hydrodynamic and a depth-averaged passive scalar transport models were developed by using a Compact Finite Volume Method (CFVM) which can assure a higher order accuracy. A typical wave current interaction experimental data set was compared with the computed results by the proposed CFVM model, and resonable agreements were observed from the comparisons. One and two dimensional scalar advection tests were conducted, and very close agreements were observed with very little numerical diffusion. Finally, a turbulent mixing simulation was done in an open channel flow, and a reasonable similarity with LES data was observed.