added script for testing particle energies -- need to fix prn