Abstract
The distribution of atmospheric motion with height following a change in terrain roughness is studied by numerically integrating the equations of motion. Instead of employing the boundary layer approximation method discussed by Elliott and by Panofsky and Townsend, we specify a functional form for the coefficient of eddy viscosity. For a shallow atmospheric layer, where we may neglect the Coriolis force and assume the eddy viscosity to vary linearly with height, there is a remarkable accordance between our numerical calculations and the results of other theoretical and experimental studies.
