c----------------------------------------------------------------------- subroutine read_orb(orb,iset) c----------------------------------------------------------------------- implicit double precision (a-h,o-z) dimension orb(*)Reads orbitals (square) from an active dump record. iset (input) is the set number of the orbitals. The matrix is automatically expanded to ntb dimentions, if appropriate.