Hi All I was wondering. I assume gromacs writes folded coordinates in the trajectory files. If so does it use rx(i) = rx(i) -boxl *nint(rx(i) / boxl) to fold the coordinates at each step? and where does it define the origin (0,0,0) Cheers Gavin