Re: [AMBER] ptraj closest water script problem

From: Daniel Roe <daniel.r.roe.gmail.com>
Date: Fri, 4 May 2012 06:39:28 -0400

Hi,

On Fri, May 4, 2012 at 5:39 AM, Thomas Cheatham <tec3.utah.edu> wrote:
> The key is how you view with VMD.  You cannot load the original prmtop
> file since this will have all the waters; there is no information in the
> trajectory file to specify what atom is what, only a list of x, y, z for
> atoms 1->N. Therefore you need some way to map and this requires either
> loading a corresponding PDB file (with no ions and only 10 waters) or a
> corresponding prmtop file (with no ions and only 10 waters).

Note that AT12 cpptraj can now generate a matching prmtop for you automatically:

trajin md_reimaged.mdcrd 201 400
trajout closest_waters.mdcrd
strip :Na+
closest 10 :13 oxygen outprefix strip

This will create a prmtop that matches the output from 'closest'
called strip.<parmname> where <parmname> is your original parm
filename. If you use cpptraj though make sure the latest bugfixes are
applied since AT12 bugfix 5 applies to cpptraj 'closest'.

-Dan

_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri May 04 2012 - 04:00:05 PDT
Custom Search