set prnlev 3 prior to the donor/acceptor specification...
> trajin file.mdcrd
> donor mask
> :47,76,78-84,94,192,193,216,260,264,268,336,361,468.NZ,N,NE,NH1,NH2,ND2,OG1
> acceptor mask :1-467.O* .H*
also, this is not a 1-1 mapping as there is one O per 2 H; you need 1-1
prnlev 3
acceptor mask WAT O H1
acceptor mask WAT O H2
prnlev 0
...however, this will likely blow memory if you have lot's of solvent...
> hbond distance 3.2 angle 60.0 series hbonds.dat
...to get around that, use the "solventacceptor" command (which keeps
track of waters interacting generally, rather than by specific ID hence is
less memory intensive)
hbond distance 3.2 angle 60.0 series h1 out hbonds.dat \
solventacceptor WAT O H1 \
solventacceptor WAT O H2
--tom
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Wed Oct 04 2006 - 06:07:19 PDT
