Algorithm for non-hydrostatic dynamics in the Regional Oceanic Modeling System

A non-hydrostatic algorithm for the Regional Oceanic Modeling System (ROMS) is proposed. It is based on a decomposition technique for hydrostatic and non-hydrostatic pressure. The algorithm has a pressure-correction scheme with split-explicit time-stepping for baroclinic and barotropic vertical modes with a free surface. The algorithm implementation requires solving a Poisson equation for a non-hydrostatic pressure that has a non-symmetric matrix in discrete form. The efficiency of a different class of solvers and preconditioners were tested. The algorithm is successfully implemented with several examples where non-hydrostatic effects are important. These include standing external gravity waves; strongly nonlinear internal wave generation and transformation; stratified shear instability and its associated mixing; and nonlinear internal tidal generation over a ridge. The corresponding changes in the pre-processing and post-processing infrastructure in the existing hydrostatic ROMS code were performed to implement parallel elliptic solvers and a new set of dynamical equations.