A finite‐volume, incompressible Navier Stokes model for studies of the ocean on parallel computers

Abstract

Abstract. The numerical implementation of an ocean model based on the incompressible Navier Stokes equations which is designed for studies of the ocean circulation on horizontal scales less than the depth of the ocean right up to global scale is described. A "pressure correction" method is used which is solved as a Poisson equation for the pressure field with Neumann boundary conditions in a geometry as complicated as that of the ocean basins. A major objective of the study is to make this inversion, and hence nonhydrostatic ocean modeling, efficient on parallel computers. The pressure field is separated into surface, hydrostatic, and nonhydrostatic components. First, as in hydrostatic models, a two-dimensional problem is inverted for the surface pressure which is then made use of in the three-dimensional inversion for the nonhydrostatic pressure. Preconditioned conjugate-gradient iteration is used to invert symmetric elliptic operators in both two and three dimensions. Physically motivated preconditioners are designed which are efficient at reducing computation and minimizing communication between processors. Our method exploits the fact that as the horizontal scale of the motion becomes very much larger than the vertical scale, the motion becomes more and more hydrostatic and the threedimensional Poisson operator becomes increasingly anisotropic and dominated by the vertical axis. Accordingly, a preconditioner is used which, in the hydrostatic limit, is an exact integral of the Poisson operator and so leads to a single algorithm that seamlessly moves from nonhydrostatic to hydrostatic limits. Thus in the hydrostatic limit the model is "fast," competitive with the fastest ocean climate models in use today based on the hydrostatic primitive equations. But as the resolution is increased, the model dynamics asymptote smoothly to the Navier Stokes equations and so can be used to address smallscale processes. A "finite-volume" approach is employed to discretize the model in space in which property fluxes are defined normal to faces that delineate the volumes. The method makes possible a novel treatment of the boundary in which cells abutting the bottom or coast may take on irregular shapes and be "shaved" to fit the boundary. The algorithm can conveniently exploit massively parallel computers and suggests a domain decomposition which allocates vertical columns of ocean to each processing unit. The resulting model, which can handle arbitrarily complex geometry, is efficient and scalable and has been mapped on to massively parallel multiprocessors such as the Connection Machine (CM5) using data-parallel FORTRAN and the Massachusetts Institute of Technology data-flow machine MONSOON using the implicitly parallel language Id. IntroductionDetails of the numerical implementation of a model which has been designed for the study of dynamical processes in the ocean from the convective, through the geostrophic eddy, up to global scale are set out. The "kernel" algorithm solves the incompressible Navier Stokes equations on the...

Dynamics of a Snowball Earth ocean

“…Our model couples the state--of--the--art ocean general circulation model, MITgcm 28 , which can simulate thick ice shelves 29 , with a recently developed two--dimensional (2D, latitude--longitude) ice flow model 15 , extending previous 1D flow models 6,7 . The model is run here in …”

“…The model used here couples the Massachusetts Institute of Technology general circulation ocean model (MITgcm) 28 with its ice--shelf package 29 to a recently--developed two dimensional (latitude--longitude) model of thick ice flow over a Snowball ocean 15 which is an extension of the 1D ice flow model of 6,7 and similar ones [31][32][33] . The ice flow model compensates for melting and sublimation at low latitudes and freezing/ precipitation at high latitudes, with an equatorward ice flow calculated from ice thickness gradients based on Glenn's law.…”

The Deep-Ocean Heat Uptake in Transient Climate Change

“…We use the MIT OGCM (Marshall et al 1997) and its adjoint (Giering 1999) with an idealized Pacific and Atlantic connected by an idealized Drake Passage (Huang et al 2002;Kamenkovich et al 2002b). The longitudinal resolution is 1°near the western and eastern boundaries, but 4°in the central ocean.…”

Modelling ocean-colour-derived chlorophyll <i>a</i>

“…The model resolves the cycling of carbon, phosphorus, nitrogen, silica, iron and oxygen through inorganic, living, dissolved and particulate organic phases (including CDOM). The biogeochemical and biological tracers are transported and mixed by the MIT general circulation model (MITgcm; Marshall et al, 1997), constrained to be consistent with altimetric and hydrographic observations (the ECCO-GODAE state estimates; Wunsch and Heimbach, 2007). This three-dimensional configuration has a coarse resolution (1 • ×1 • horizontally) and 23 levels ranging from 10 m at the surface to 500 m at depth.…”