https://hgpu.org/?p=6282
An MPI-CUDA implementation of an improved Roe method for two-layer shallow water systems