Skip to content
Switch branches/tags
Go to file
Cannot retrieve contributors at this time
executable file 31 lines (22 sloc) 834 Bytes
#!/usr/bin/env python
# coding=utf-8
MDAnalysis example: Convert DCD trajectory into XTC
This example shows how one can use MDAnalysis to convert between
different trajectory formats.
from MDAnalysis.tests.datafiles import PDB_small, DCD
from MDAnalysis import Universe, Writer
import os.path
root, ext = os.path.splitext(os.path.basename(DCD))
xtcname = root + '.xtc' # output format determined by extension
u = Universe(PDB_small, DCD)
# create a writer instance for the output trajectory
w = Writer(xtcname, u.trajectory.numatoms)
# loop through the trajectory and write a frame for every step
for ts in u.trajectory:
print "Converted frame {0:d}".format(ts.frame)
print "Converted {0!r} --> {1!r}".format(DCD, xtcname)