An approach for calculating turbulent flows in a wave-current boundary layer over a slowly varying bed is presented. Waves are periodic in time with several harmonics. In this paper, we adopt a time invariant eddy viscosity model, in which the eddy viscos