Create a separate MPI operation for the top-level multipole reduction. Closes #485 See merge request !661