AMBER Archive (2003)

Subject: Re: AMBER: question about ptraj

From: Thomas E. Cheatham, III (cheatham_at_chpc.utah.edu)
Date: Wed Jul 30 2003 - 17:38:52 CDT


Yes, you outline two excellent problems with ptraj namely (1) that the
mask parser is funny and very particular with respect to deciphering
mask specifications (i.e. you cannot quote or AND or OR)
and (2) it is often difficult to decipher my cryptic manual writing.

> diffusion :SPM :801-26516 0.2 average DiffSPMKz15
> diffusion :K+ :799 0.2 average DiffKz15

Diffusion only allows specification of a SINGLE mask string. If you
really mean to calculate the diffusion, averaged for all residues named
SPM and residues 801-26516, I think you can specify:

diffusion :SPM,801-26516 0.2 average DiffSPMKz15

A nice way to figure out if a mask string works and is picking what you
think it should be picking, use rdparm and the "checkmask" command

        rdparm prmtop

        checkmask :SPM,801-26516

> radial RDF_Kz15_SPM 0.05 8.0 :SPM :799
> radial RDF_NVEk_OO 0.05 8.0 :801-26516_at_O
> radial RDF_MTS_HH 0.05 8.0 :801-26516_at_H
> radial RDF_MTS_OH 0.05 8.0 :801-26516@* go
>
> I do not understand why my job stops with the message
> diffusion:solvent mask not specified

Now, the above is a BUG. I built the radial function using the diffusion
routine as a prototype; I did not change the error message which SHOULD
read "radial: no solvent mask specified" as your error is occuring not in
the diffusion command but in the radial command...

Beyond this, I cannot see where the error is, nor can I reproduce it.
Which version of ptraj are you using? Looking at the code, the only time you
should get this error is if the mask string is missing (which it appears
not to be) so I am a bit confused... Alternatively, perhaps one of your
atom selections is not selecting any atoms? This would also break it.

An example of a complicated mask strings that works:

  :1-100_at_O,:101-200_at_H?,:201-299_at_H1,:A,@CA

key is no spaces, comma's between entries, ...

We will have an updated parser soon that should make life a tad bit
easier (but still support the current syntax).

\ Thomas E. Cheatham, III (Assistant Professor) College of Pharmacy
| Departments of Medicinal Chemistry and of University of Utah
| Pharmaceutics and Pharmaceutical Chemistry 30 South 2000 East, Room 201
| & Center for High Performance Computing Salt Lake City, Utah 84112
|
| e-mail: tec3_at_utah.edu phone: (801) 587-9652 FAX: (801) 585-9119
\ http://www.chpc.utah.edu/~cheatham Offices: BPRP295A / INSCC 418

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu