Parallel finite element method utilizing the mode splitting and sigma coordinate for shallow water flows