Previously, the code used function ph1() to access the valence
threshold energies (ionization potentials) stored in file phfit.dat.

Define method atmdat.getIonPot() to access the NIST IPs (read from file
ionization_potentials.dat), and use it in lieu of the call to ph1(),
where possible.