An orthogonal unstructured computational mesh is required for the implementation to be accurate. Mixed meshes e.g. triangles and squares as in the Roskilde Fjord setup are supported. The velocities in the model are defined on the interfaces between the cells and the scalars in the cell centers, hence it is a finite difference-finite volume formulation. The model uses horizontal layers (z-level).
A theta formulation is used to define the degree of implicitness. When theta equals zero the model is explicit and when theta equals one it is fully implicit. It can be show that the highest accuracy and efficiency is achieved when theta equals 0.5 (Casulli and Cattani 1992).
A predictor-corrector algorithm is used, where first, a provisional solution is calculated where the nonhydrostatic pressure component is neglected, thereby obtaining the hydrostatic (barotropic) solution. Subsequently the provisional solution can optionally be corrected by solving the implicit equation system for the nonhydrostatic pressure component, thereby obtaining the nonhydrostatic solution. Both the nonhydrostatic and hydrostatic solutions are mass conservative.
|