Numerical methods for the three-dimensional shallow water equations on supercomputers