On Sun, Sep 16, 2007, David Cerutti wrote:
>
> Fixing this, however, does not do much. Like I said earlier, however,
> SHAKE does seem to work on the protein atoms
Sander or pmemd reports a line like this:
Number of triangulated 3-point waters found: 11324
Are you getting this printout, with the proper number of water molecules?
Water SHAKE is done in the routine setlep() [called shake_fastwater() in
pmemd.] You might put some print statements in there to see what parameters
this is being called with.
Can you make a small system that shows this problem? It's hard to debug just
based on the information we have so far.
....dac
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Wed Sep 19 2007 - 06:07:21 PDT
