AMBER Archive (2009)

Subject: Re: AMBER: RMS plot

From: Carlos Simmerling (carlos.simmerling_at_gmail.com)
Date: Fri Jan 02 2009 - 08:08:44 CST


probably the easiest way is to make a script that creates a ptraj
input for this, like the following. then you'll need to average the
outputs, again a script will do that. maybe there is an easier way or
someone else might have a better program- I don't think ptraj does it
directly. hopefully you can use the example below as a starting place
and modify it for your needs. make sure to watch out for the "nofit"
parts and change the masks to what you want.

#!/bin/csh

set i=1
set last = 100
echo "ptraj prmtop << EOF"
echo "trajin ./mdcrd.gz"
echo "reference ./ref.crd"

echo "rmsd reference mass out rmsd.res.ALL.dat :$i-$last"
while ($i < 10)
        echo "rmsd reference mass out rmsd.res.all.00$i.dat :$i nofit"
        @ i++
end
while ($i < 100)
        echo "rmsd reference mass out rmsd.res.all.0$i.dat :$i nofit"
        @ i++
end
while ($i <= $last)
        echo "rmsd reference mass out rmsd.res.all.$i.dat :$i nofit"
        @ i++
end

echo "EOF"

On Wed, Dec 31, 2008 at 2:11 PM, Beale, John <jbeale_at_stlcop.edu> wrote:
> I need to make a plot of C-alpha RMS deviations (in Angstroms) vs. residue number for a protein that I am studying with MD. Is it possible to do this in ptraj? If so, how do I set it up?
>
> John Beale
> -----------------------------------------------------------------------
> The AMBER Mail Reflector
> To post, send mail to amber_at_scripps.edu
> To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
> to majordomo_at_scripps.edu
>
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
      to majordomo_at_scripps.edu