Salmon describes in this paper a new method to model the ocean based partly on Lattice-Boltzmann model and operator splitting. The resulting code is highly parallel and composed of many short routines that make the computation efficient. Although the code does preserve the main conservation laws,
“our model equations do not conserve potential vorticity”.
The code can handle large-scale processes as well as sound waves which makes it particularly useful to study nonhydrostatic processes.
Salmon R. (2009): An Ocean Circulation Model Based on Operator-Splitting, Hamiltonian Brackets, and the Inclusion of Sound Waves, J. Phys. Ocean., 39, 1615–1633.