Re: AMBER: Harmonic restraint along z-axis

From: David A. Case <>
Date: Thu, 27 Jul 2006 07:32:52 -0700

On Wed, Jul 26, 2006, Myunggi Yi wrote:
> I'd like to restrain some atoms at a certain value (z=5.0) along
> the z-axis with a harmonic potential during simulations.
> The source code is so complicated for me.
> Would you give me some draft or orientation for adding
> my own potential to the sander code?

Look at subroutine xconst in $AMBERHOME/src/sander/ene.f. That shows you
how a general harmonic constraint is implemented and called. You could
modify this to have forces only along the z-direction.

...good luck...dac

The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" to
Received on Sun Jul 30 2006 - 06:07:09 PDT
Custom Search