Faster XYZ trajectory output method

[arvk]

The current method simply loops over the snapshot method. This has the overhead of constantly opening and closing the file.