AMBER Archive (2002)

Subject: Re: SETVEL, incorrect number of arguments?

From: David A. Case (case_at_scripps.edu)
Date: Tue Sep 17 2002 - 14:18:38 CDT

On Mon, Sep 16, 2002, Giovanni Fossati wrote:
>
> while compiling on IBM-AIX, I have just noticed a possible error in
> sander_classic/sander.f (amber6). There are two calls to SETVEL, and the
> second calls it with one too many argument.
> These are the definition and the calls.
>
> dynlib.f
> SUBROUTINE SETVEL(NRP,NR,NTU,NTX,X,V,F,WINV,
> + TEMPI,HEAT,DT,INIT,IG,iscale,scalm)
>
> sander.f
> CALL SETVEL( NRP,NR,NTU,NTX,X(L30),X(L40),X(L35),X(L20),
> + TEMPI,HEAT,DT,INIT,IG,iscale,scalm)
>
> * CALL SETVEL(NPM,NRP,NR,NTU,NTX,X(L30),X(L40),X(L35),X(L20),
> * TEMPI,HEAT,DT,INIT,IG,iscale,scalm)
>
>

This does indeed seem to be an error, and could lead to unpredictable
results if you tried to set random velocities at the beginning of an MD
run. I've posted a bug-fix to the Amber web site.

Thanks for pointing this out.....dac 

-- 

==================================================================
David A. Case                     |  e-mail:      case_at_scripps.edu
Dept. of Molecular Biology, TPC15 |  fax:          +1-858-784-8896
The Scripps Research Institute    |  phone:        +1-858-784-9768
10550 N. Torrey Pines Rd.         |  home page:                   
La Jolla CA 92037  USA            |    http://www.scripps.edu/case
==================================================================