NH-GOMO, an efficient and easily portable non-hydrostatic ocean model

Numerical simulation of nonlinear gravity internal waves with non-hydrostatic ocean models, especially these which using the terrain-following sigma-coordinate, is challenging. The expensive computation cost, which is caused by the dynamic pressure Poisson solver in cases using fine grid resolution in both directions (horizontal and vertical), is the main reason. A non-hydrostatic ocean model named NH-GOMO is constructed based on a partially implicit finite difference scheme for the dynamic pressure and adopts an idea of “decimation and interpolation”. A significant optimization for the pressure Poisson solver, which brings no obvious accuracy loss, is obtained with these technologies. The automatic parallel operator library named OpenArray is used as the bottom layer of this model and make it easy to transport between different computing platforms. Accuracy and efficiency have been validated by several ideal test cases.