Iterative process to calculate scf open-shell orbitals