********> bugfix.10 Author: Dave Case Date: Mar 16, 2000 Programs: sander_classic Description: solvent temperature is computed incorrectly if ntt=5 and MPI is used Fix: Make the following change to amber6/src/sander_classic/runmd.f *** OLD runmd.f. Tue Aug 3 11:13:18 1999 --- NEW runmd.f Thu Mar 16 15:01:00 2000 *************** *** 774,783 **** c #ifdef MPI i3 = 3*(istart-1) ! ilimit = nsolut ! if (ilimit .gt. iend) ilimit = iend ! do j=istart,ilimit ! #else I3 = 0 DO J = 1,nsolut --- 774,781 ---- c #ifdef MPI i3 = 3*(istart-1) ! ibiggest = min( iend, nsolut ) ! do j=istart,ibiggest #else I3 = 0 DO J = 1,nsolut *************** *** 799,805 **** c --- solvent c #ifdef MPI ! do j=nsolut+1,iend #else do j=nsolut+1,nrp #endif --- 797,804 ---- c --- solvent c #ifdef MPI ! ismallest = max( nsolut+1, istart ) ! do j=ismallest,iend #else do j=nsolut+1,nrp #endif *************** *** 822,830 **** c #ifdef MPI i3 = 3*(istart-1) ! ilimit = nsolut ! if (ilimit .gt. iend) ilimit = iend ! do j=istart,ilimit #else I3 = 0 DO J = 1,nsolut --- 821,828 ---- c #ifdef MPI i3 = 3*(istart-1) ! ibiggest = min( iend, nsolut ) ! do j=istart,ibiggest #else I3 = 0 DO J = 1,nsolut *************** *** 842,848 **** c --- solvent c #ifdef MPI ! do j=nsolut+1,iend #else do j=nsolut+1,nrp #endif --- 840,847 ---- c --- solvent c #ifdef MPI ! ismallest = max( nsolut+1, istart ) ! do j=ismallest,iend #else do j=nsolut+1,nrp #endif *************** *** 912,920 **** c #ifdef MPI i3 = 3*(istart-1) ! ilimit = nsolut ! if (ilimit .gt. iend) ilimit = iend ! do j=istart,ilimit #else I3 = 0 DO J=1,nsolut --- 911,918 ---- c #ifdef MPI i3 = 3*(istart-1) ! ibiggest = min( iend, nsolut ) ! do j=istart,ibiggest #else I3 = 0 DO J=1,nsolut *************** *** 929,935 **** c --- solvent c #ifdef MPI ! do j=nsolut+1,iend #else DO j=nsolut+1,NRP #endif --- 927,934 ---- c --- solvent c #ifdef MPI ! ismallest = max( nsolut+1, istart ) ! do j=ismallest,iend #else DO j=nsolut+1,NRP #endif *************** *** 1006,1014 **** c #ifdef MPI i3 = 3*(istart-1) ! ilimit = nsolut ! if (ilimit .gt. iend) ilimit = iend ! do j=istart,ilimit #else I3 = 0 DO J = 1,nsolut --- 1005,1012 ---- c #ifdef MPI i3 = 3*(istart-1) ! ibiggest = min( iend, nsolut ) ! do j=istart,ibiggest #else I3 = 0 DO J = 1,nsolut *************** *** 1025,1031 **** c --- solvent c #ifdef MPI ! do j=nsolut+1,iend #else DO j=nsolut+1,NRP #endif --- 1023,1030 ---- c --- solvent c #ifdef MPI ! ismallest = max( nsolut+1, istart ) ! do j=ismallest,iend #else DO j=nsolut+1,NRP #endif