Implement unwind_positions() subroutine
Unwind positions properly in analyze_md
Calculate and subtract the center of mass
Save the COM coordintes
Plot the COM
Simplify loop
Average the graphs