get_sym Function

    pure function get_sym(this, orb) result(sym)
        ! gives the symmetry label associated with the k-vector of
        ! spatial orbital (orb)
        class(lattice), intent(in) ::  this
        integer, intent(in) :: orb
        integer :: sym
        unused_var(this)

        sym = this%sites(orb)%k_sym

    end function get_sym