X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD%2FMD_A-MTS.F;h=30d72a7dd64d8c7039ea7b064648af280b5d78ee;hb=a20d7ba73b9caa74140350c4c97fa62d83711703;hp=29e4fcb1f2761e5cf611e91830d1515e964a8ebc;hpb=a3ef63eda7d771657cc2df1d275e0a03a5cc1966;p=unres.git

diff --git a/source/unres/src_MD/MD_A-MTS.F b/source/unres/src_MD/MD_A-MTS.F
index 29e4fcb..30d72a7 100644
--- a/source/unres/src_MD/MD_A-MTS.F
+++ b/source/unres/src_MD/MD_A-MTS.F
@@ -1681,7 +1681,8 @@ c if the friction coefficients do not depend on surface area
           stdforcp(i)=stdfp*dsqrt(gamp)
         enddo
         do i=nnt,nct
-          stdforcsc(i)=stdfsc(itype(i))*dsqrt(gamsc(itype(i)))
+          stdforcsc(i)=stdfsc(iabs(itype(i)))
+     &                *dsqrt(gamsc(iabs(itype(i))))
         enddo
       endif
 c Open the pdb file for snapshotshots
@@ -1806,6 +1807,7 @@ c      rest2name = prefix(:ilen(prefix))//'.rst'
          write (iout,'(i3,3f10.5,3x,3f10.5)') i,(d_t(j,i),j=1,3),
      &   (d_t(j,i+nres),j=1,3)
        enddo 		 	 
+       call flush(iout)
 c  Zeroing the total angular momentum of the system
        write(iout,*) "Calling the zero-angular 
      & momentum subroutine"
@@ -1823,6 +1825,7 @@ c Removing the velocity of the center of mass
       if(me.eq.king.or..not.out1file)then
        write (iout,*) "vcm right after adjustment:"
        write (iout,*) (vcm(j),j=1,3) 
+       call flush(iout)
       endif
       if (.not.rest) then	  	
          call chainbuild
@@ -2078,6 +2081,7 @@ c          write (iout,*) "i",i," ii",ii," geigen",geigen(i),
 c     &      " d_t_work_new",d_t_work_new(ii)
         enddo
       enddo
+      call flush(iout)
 c diagnostics
 c      Ek1=0.0d0
 c      ii=0
@@ -2116,6 +2120,7 @@ c Transfer to the d_t vector
 c      do i=0,nres-1
 c        write (iout,*) "d_t",i,(d_t(j,i),j=1,3)
 c      enddo
+c      call flush(iout)
       do i=nnt,nct
         if (itype(i).ne.10) then
           do j=1,3