For an isolated atom calculation, the occupation for various orbitals needs to be specified in the OCCUPATIONS card, in the input file (or so I believe. It wasn't converging otherwise). But what order should I follow while giving them? Is it using the l quantum number, or the Aufbau principle? Also, should I include all the orbitals used in generating the pseudopotential I'm using?