AMBER Archive (2006)Subject: AMBER: Ignoring every other frame in a .crd trajectory in Ptraj
From: Xiaojian Deng (xjd002_at_gmail.com)
Date: Tue Aug 01 2006 - 12:55:31 CDT
Suppose I split a 2 ns sander simulation into two parts.
Part 1 has dt = 0.002, nstlim = 50000, ntpr = 100
Part 2 has dt = 0.001, nstlim = 100000, ntpr = 100
The coordinate files for each part, respectively, are crd.1 and crd.2. Even
though each output coordinate file represents 1ns of the simulation, the
second one is twice as large because the dt is half the size. Now if I
decide I want to output the RMSD values into a file (that I'll plot), I
could create an input file for PTRAJ like:
trajin crd.1
trajin crd.2
rmsd first out rmsd.txt
But obviously, when I try to plot rmsd.txt, the scaling would be off,
meaning the last 2/3 of the graph would only represent the second half of
the simulation. Is there any way that I can specify in ptraj to extract
only every other frame in crd.2 to avoid this pitfall? Thanks.
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
|