Implementation of a direct nonhydrostatic pressure gradient discretisation into a layered ocean model

► Detailed straight-forward nonhydrostatic extension of a full 3D explicit mode-splitting model. ► Demonstration of consistency and convergence. ► Discussion of stability issues and implications for efficiency and practical applicability.