Abstract

In this work, we take a modern high-resolution finite-volume scheme for solving the rotational shallow-water equations and extend it with features required to run real-world ocean simulations. Our contributions include a spatially varying north vector and Coriolis term required for large scale domains, moving wet-dry fronts, a static land mask, bottom shear stress, wind forcing, boundary conditions for nesting in a global model, and an efficient model reformulation that makes it well-suited for massively parallel implementations. Our model order is verified using a grid convergence test, and we show numerical experiments using three different sections along the coast of Norway based on data originating from operational forecasts run at the Norwegian Meteorological Institute. Our simulation framework shows perfect weak scaling on a modern P100 GPU, and is capable of providing tidal wave forecasts that are very close to the operational model at a fraction of the cost. All source code and data used in this work are publicly available under open licenses.

摘要

在这项工作中，我们采用了一种现代的高分辨率有限体积方案来求解浅水旋转方程，并将其扩展为运行真实世界海洋模拟所需的功能。我们的贡献包括空间变化的北向矢量和大范围域所需的科里奥利项，移动的干湿锋面，静态土地遮罩，底部切应力，风强迫，在全局模型中嵌套的边界条件以及有效的模型重构，使它非常适合大规模并行实现。我们的模型阶数已通过网格收敛测试进行了验证，并且我们基于挪威气象研究所运行预报得出的数据，使用了挪威沿海三个不同断面的数值实验。我们的仿真框架显示了在现代P100 GPU上的完美弱缩放，并且能够以非常低的成本提供非常接近于运营模型的海潮预报。本工作中使用的所有源代码和数据都可以在开放许可下公开获得。