Spherical averaging of occupied orbitals on atoms

For my localization scheme, I need to get MO coefficients for the full molecule, and for atomic fragments. The problem is that for the atomic fragments like carbon, oxygen, …, the occupied p orbitals are not spherically averaged. In the code DIRAC, one can specify fractional electron number in those orbitals (like 2 electrons in 6 spin-orbitals for carbon) to get this spherical averaging. It is called “Open shell SCF calculation using Average-of-Configuration”. Basically, instead of optimizing the energy for a single open-shell state, it optimizes the energy for a limited set of open-shell states.
Is there anything like this in psi4 ?

Thanks