store_whole_initiator_space Subroutine

    subroutine store_whole_initiator_space()

        use Parallel_neci, only: MPIAllGatherV

        integer :: ierr
        character(len=*), parameter :: t_r = "store_whole_initiator_space"

        allocate(initiator_space(0:NIfTot, initiator_space_size), stat=ierr)
        initiator_space = 0_n_int

        call MPIAllGatherV(SpawnedParts(0:NIfTot, 1:initiator_sizes(iProcIndex)), &
                           initiator_space, initiator_sizes, initiator_displs)

    end subroutine store_whole_initiator_space