|
MITgcm is discretized by carving up the fluid in to
"finite volumes" with six sides. In the interior of
the fluid this creates a grid of box shaped elements. Where the
topography (solid bottom) intersect the grid we shape the volume
to fit the topography. This can be done in varying degrees of
approximation, such as piecewise constant or discontinuous
partial steps ("lopping"), piecewise linear or even
higher order.
In contrast, the finite difference method used in ocean
models until recently would fit the topography to the grid;
there is clearly a large representation error.

Our use of topography using finite volumes also
opens up as yet unexplored possibilities, as suggested by the
figures below.

Representing porous barriers


Representing thin barriers


More
details can be found on the web here.
The paper describing this work is:
Adcroft, A.J., Hill, C.N. and J. Marshall, (1997)
Representation of topography by shaved cells in a height
coordinate ocean model
Mon Wea Rev, vol 125, 2293-2315
|