subroutine load_orb_perm(perm) integer, intent(in) :: perm(:) orbital_permutation = perm end subroutine load_orb_perm