A three-dimensional hydrostatic model for coastal and ocean modelling using a generalised topography following co-ordinate system

A three-dimensional hydrostatic model is presented that combines a generalised vertical co-ordinate system with an efficient implicit solution technique for the free surface. The model is capable of maintaining high resolution in the surface and/or bottom boundary layers as well as dealing with steep topography. Horizontal diffusion is calculated using the Smagorinsky formulation and a k-epsilon turbulence model is used in the vertical. In addition the model uses higher-order advection routines. An important aspect in three-dimensional models is the choice of vertical discretisation. If one is mostly interested in problems which are governed by boundary layer flows, a terrain following or sigma co-ordinate system seems attractive. This paper focuses on the development of a generalised sigma-type grid in a three-dimensional hydrostatic model. The generalised grid offers a wide range of possibilities including grid refinement toward the bed or surface, a mixed layer transformation, and a constant layer transformation where the lowermost or uppermost grid cells can be specified to have a constant height above the bed or below the surface. A number of tests are presented which show that the model is capable of simulating both shallow nearshore, estuarine flows as well as large-scale geophysical flows. These include an extreme flooding event in the shallow North Sea and the Odden ice tongue formation in the Greenland Sea. (C) 2002 Elsevier Science Ltd. All rights reserved.