Abstract.We report results gained with a three-dimensional, semi-implicit, semi-spectral model of the shallow water equations on the rotating Earth that allowed one to compute the wind-induced motion in lakes. The barotropic response to unidirectional, uniform winds, Heaviside in time, is determined in a rectangular basin with constant depth, and in Lake Constance, for different values and vertical distributions of the vertical eddy viscosities. It is computationally demonstrated that both the transitory oscillating, as well as the steady state current distribution, depends strongly upon the absolute value and vertical shape of the vertical eddy viscosity. In particular, the excitation and attenuation in time of the inertial waves, the structure of the Ekman spiral, the thickness of the Ekman layer, and the exact distribution and magnitude of the upwelling and downwelling zones are all significantly affected by the eddy viscosities. Observations indicate that the eddy viscosities must be sufficiently small so that the oscillatory behaviour can be adequately modelled.Comparison of the measured current-time series at depth in one position of Lake Constance with those computed on the basis of the measured wind demonstrates fair agreement, including the rotation-induced inertial oscillation.