Computational simulations of ocean circulation rely on the numerical solution of partial differential equations of fluid dynamics, as applied to a relatively thin layer of stratified fluid on a rotating globe. This paper describes some of the physical and mathematical properties of the solutions being sought, some of the issues that are encountered when the governing equations are solved numerically, and some of the numerical methods that are being used in this area.