A modified sigma equations' approach to the numerical modeling of Great Lakes hydrodynamics