There might be bugs in ptraj, especially in the mask parser. :P
But using the following command should get the heavy atoms you want:
":1-466.C*,N*,O*,P*"
Hope that helps.
Jianyin
On 1/23/08, Vlad Cojocaru <Vlad.Cojocaru.eml-r.villa-bosch.de> wrote:
>
> Dear amber users,
>
> I was trying tu use the following mask string in ptraj (the version
> distributed with AMBER 9)
>
> ":1-466 & !.H="
>
> to select the heavy atoms of a protein. According to the AMBER 9 manual
> this should work if enclosed in double quotes (which I did). However I
> always get an error saying "the = symbol wasnt recognized"
>
> The "=" symbol is described in the manual as an accepted mask wildcard
> in the new mask parser.
>
> Could somebody tell me why is this happening? Is there a bug that I am
> not aware of ?
>
> Thanks
>
> vlad
>
> --
>
> ----------------------------------------------------------------------------
> Dr. Vlad Cojocaru
>
> EML Research gGmbH
> Schloss-Wolfsbrunnenweg 33
> 69118 Heidelberg
>
> Tel: ++49-6221-533266
> Fax: ++49-6221-533298
>
> e-mail:Vlad.Cojocaru[at]eml-r.villa-bosch.de
>
> http://projects.villa-bosch.de/mcm/people/cojocaru/
>
>
> ----------------------------------------------------------------------------
> EML Research gGmbH
> Amtgericht Mannheim / HRB 337446
> Managing Partner: Dr. h.c. Klaus Tschira
> Scientific and Managing Director: Prof. Dr.-Ing. Andreas Reuter
> http://www.eml-r.org
>
> ----------------------------------------------------------------------------
>
>
> -----------------------------------------------------------------------
> The AMBER Mail Reflector
> To post, send mail to amber.scripps.edu
> To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
>
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Sun Jan 27 2008 - 06:07:20 PST
