Abstract
Observations of Lake Ontario during the International Field Year for the Great Lakes are used to develop a. three-dimensional numerical model for calculating temperature and current. The model has a variable grid resolution and a horizontal smoothing which filters out small-scale vertical motion caused by truncation error but has little effect on the strong currents of the coastal boundary layer. Resolution of the shore zones and reduced horizontal smoothing improve simulation of both long-term mean flow and current reversals due to low-frequency waves.