Re: AMBER: recovering total force array

From: David A. Case <>
Date: Tue, 19 Oct 2004 14:26:31 -0700

On Tue, Oct 19, 2004, Dave S Walker wrote:

> In runmd.f, I had been trying to replace the velocity array with the force
> array when making the call to CORPAC(). The forces (and perhaps the
> velocities and coordinates) in runmd.f I suspect are distributed among the
> nodes and I had been trying to use a reduce statement to recover all the
> forces before they're written over in Step 3. Doing this has failed to
> reproduce the array I get with 1 cpu.

Did you try just forcing the code to go through the "if( init==3 )" block of
code in subroutine fdist()? That should put all of the forces onto all nodes,
including the master node, where they can be printed out.


The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" to
Received on Tue Oct 19 2004 - 22:53:00 PDT
Custom Search