Re: AMBER: Problem with Getting SHAKE to work

From: David A. Case <>
Date: Mon, 17 Sep 2007 09:38:26 -0700

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.


The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" to
Received on Wed Sep 19 2007 - 06:07:21 PDT
Custom Search