Averaging over a trajectory
dcd_export.py
# Conversion of an MMTK trajectory to DCD format (CHARMM/XPlor).
# A PDB file with the same atom order is created at the same time,
# because otherwise it would be impossible to interpret the data
# in the DCD file.
#
# Note: In order to create the trajectory file "bala1.nc" which is
# read by this example, you must run the example
# MolecularDynamics/protein.py!
#
from MMTK import *
from MMTK.Trajectory import Trajectory
from MMTK.DCD import writeDCDPDB
trajectory = Trajectory(None, "bala1.nc")
universe = trajectory.universe
writeDCDPDB(trajectory.configuration, 'bala1.dcd', 'bala1.pdb')