Re: AMBER: Ignoring every other frame in a .crd trajectory in Ptraj

From: Carlos Simmerling <carlos.csb.sunysb.edu>
Date: Tue, 01 Aug 2006 14:13:07 -0400

check the ptraj section of the manual, you can
specify range and stride through the frames in the trajin command.

Xiaojian Deng wrote:

> 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.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Wed Aug 02 2006 - 06:07:13 PDT
Custom Search