Efficient multigrid solvers for mixed finite element discretisations in NWP models