Would be nice, at least for debugging, to add a new routine the exports the full list of kpoints, but partitioned into orbits.