Abstract
The elements of an eigenfunction expansion for time-dependent currents as a function of depth are worked out for viscosity that is given as a parabolic function of depth that goes to zero at both the bottom and top of the water. This yields currents with logarithmic behavior characteristic of turbulent boundary layers at both the bottom and top. Also, solutions are obtained for the two viscosity functions that are half a parabola, going to zero at either the bottom or top but not both. In all cases the solutions are Legendre functions. In some cases the eigenfunctions are Legendre polynomials.