********>Bugfix 25:
Author: Dave Mathews
Date: 09/29/2004
Programs: sander
Description: If weight changes are used to vary the value of TEMP0 during
the run, Langevin dynamics will not work properly, since the
random force properties are not updated correctly.
Fix: apply the following patch to amber8/src/sander/runmd.f
------------------------------------------------------------------------------
*** runmd.f 2004/08/08 17:32:08 7.258
--- runmd.f 2004/09/29 17:37:03
***************
*** 1012,1017 ****
--- 1012,1027 ----
! Loncharich, Brooks and Pastor, Biopolymers 32:523-535 (1992),
! Eq. 11. (Note that the first term on the rhs of Eq. 11b
! should not be there.)
+
+ ! Update the Langevin parameters, since temp0 may have changed:
+ if (nmropt>0) then
+ sdfac = sqrt( 4.d0*gammai*boltz2*temp0/dtx )
+ #ifdef LES
+ sdfacles = sqrt( 4.d0*gammai*boltz2*temp0les/dtx )
+ #endif
+ endif
+
+
#ifdef PSANDER
do idcmp=0,nneedatm-1
------------------------------------------------------------------------------
Temporary workarounds: none.