Move calculation of a particles Ua to a method

This moves the calculation of the universal anisotropy
index to a method in the particle class. In doing this
we fix a bug in the calculation of Ua for plotting
(we were always using particle 1 for this).
Allow multiple files to be read

Also, make optons work using argpars