get_neighbors_site Function

    function get_neighbors_site(this) result(neighbors)
        ! this is a generic routine to get the neighbors of a site
        class(site) :: this
        ! i need assumed array shape and size here or? check how i do that!
        ! can i use the stored number of neighbors in the type?
        ! i can use n_neighbors directly but not the function which gets me
        ! n_neighbors.. strange and unfortunate..
        integer :: neighbors(this%n_neighbors)

        neighbors = this%neighbors

    end function get_neighbors_site