Two three-dimensional primitive equation numerical ocean models are applied to the problem of internal Kelvin waves and coastal upwelling in the Great Lakes. One is the Princeton Ocean Model (POM) with a terrain-following (sigma) vertical coordinate, and the other is the Dietrich/Center for Air Sea Technology (DIECAST) model with constant z-level coordinates. The sigma coordinate system is part...