Dear Amber,
Please note that I observe some strange behavior for the ptraj parser.
Attached is the ptraj script I use. (yes I'm sure there are better ways
to construct the script)
I have found that if the file name begins with the text "nofit"
then the parser mistakenly picks up this file name as the "nofit" key
word and scrambles the rest of the line.
The easy solution is to name the files "sys.nofit.rms" or "notfit" or
anything that does not begin with "nofit"!
Just a friendly reminder to all ptraj users
"Note that the current parser is not very sophisticated." (amber
manuals, 6,7,8)
Now can anyone suggest how I select only the heavy atoms for a limited
number of residues?
example:
:1-50.~H* applies the "NOT" to the residue selection AND atom selection.
:1-50.C*,N*,O*,S*,P* should pick up all non-hydrogens for resids 1-50
right? or does this miss some atoms?
Other suggestions for getting the rms that I want?
Thanks,
Tom
************************** ptraj script ****************
#!/bin/csh
set BAC = ".CA,C,N,O"
set SYS = ":1-105"
set PRO = ":1-76"
set DNA = ":77-105"
set NA = ":106-116"
set CL = ":117-128"
set PARM = "nowat.parm"
set TRAJ1 = "nowat.eq1-2.dcd"
set BASEDIR = "/data/bishop/sims/"
setenv AMBERHOME /usr/local/amber/amber8
$AMBERHOME/exe/ptraj $PARM << EOF
trajin nowat.eq1-2.dcd 1 5 1
rms first mass out nofit.sys.rms $SYS nofit
rms first mass out fit.sys.rms $SYS
rms first mass out fit.pro.rms $PRO
rms first mass out fit.probac.rms $PRO$BAC
go
EOF
--
http://dna.cbr.tulane.edu
Center for Bioenvironmental Research
Tulane and Xavier Universities
1430 Tulane Avenue SL-3
New Orleans, LA 70112
Phone: 504-988-6203 Fax: 504-585-6428
bishop.tulane.edu
-----------------------------------------
Alternate Contact:
Center for Computational Sciences
Voice: (504) 862-8391
-----------------------------------------
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Wed Apr 07 2004 - 16:53:00 PDT
