Turning the code into nmd code.

Right now stuck on how to efficently collect the velocities from lammps into one array.