********>Bugfix 31: Author: found by John D. Chodera Date: 01/18/03 Programs: sander Description: A typo in a units conversion factor makes randomized velocities come out too cold. This is especially troublesome if you are using the "vrand" option to periodically re-randomize velocities. Fix: apply the following patch to amber7/src/sander/newvel.f ------------------------------------------------------------------------------ *** newvel.f 2001/11/07 18:07:38 7.7 --- newvel.f 2003/01/18 18:49:40 *************** *** 55,61 **** DIMENSION V(*),WINV(*),ISCR(*),IBEL(*) C ISEED = IG ! BOLTZ = (8.31441D-3/4.814D0)*TEMP0 NTOTAT = NRP C C Clear scratch array; after any velocity is reset, ISCR(I) is set to 1. --- 55,61 ---- DIMENSION V(*),WINV(*),ISCR(*),IBEL(*) C ISEED = IG ! BOLTZ = (8.31441D-3/4.184D0)*TEMP0 NTOTAT = NRP C C Clear scratch array; after any velocity is reset, ISCR(I) is set to 1. ------------------------------------------------------------------------------ Temporary workarounds: none