subroutine setup_UMatInd() use SystemData, only: nBasis nBI_umat = numBasisIndices(nBasis) end subroutine setup_UMatInd