This thesis presents a new methodology for representing the effect of spiral flow on the bed shear stress magnitude and direction, using a low-vertical-resolution three-dimensional model. Furthermore, it discusses the performance of the new methodology in comparison with conventional high-vertical-resolution three-dimensional and depth-averaged models.