get_highest_pop_node Subroutine

    subroutine get_highest_pop_node(n_states, largest_dets, all_norm)
        ! routine to give the n_states most populated states largest_dets
        ! per node
        integer, intent(in) :: n_states
        integer(n_int), intent(out) :: largest_dets(0:niftot, n_states)
        real(dp), intent(out), optional :: all_norm
#ifdef DEBUG_
        character(*), parameter :: this_routine = "get_highest_pop_node"
#endif
        real(dp) :: norm

        call proc_most_populated_states(n_symmetry_states, GLOBAL_RUN, &
                                        largest_dets, CurrentDets, TotWalkers, norm)

        if (present(all_norm)) then
            call MpiSum(norm, all_norm)
        end if

    end subroutine get_highest_pop_node