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

diff --git a/source/unres/src_MD/MD_A-MTS.F b/source/unres/src_MD/MD_A-MTS.F
index 30d72a7..995b558 100644
--- a/source/unres/src_MD/MD_A-MTS.F
+++ b/source/unres/src_MD/MD_A-MTS.F
@@ -1096,7 +1096,7 @@ c Backup the coordinates, velocities, and accelerations
 
        if(tnp .or. tnp1) then
 #ifndef G77
-        write (iout,'(a3,7f)') "TTT",EK,s_np,potE,pi_np,Csplit,
+        write (iout,'(a3,7f20.10)') "TTT",EK,s_np,potE,pi_np,Csplit,
      &                          E_long,energia_short(0)
 #else
         write (iout,'(a3,7f20.10)') "TTT",EK,s_np,potE,pi_np,Csplit,
@@ -1681,8 +1681,7 @@ c if the friction coefficients do not depend on surface area
           stdforcp(i)=stdfp*dsqrt(gamp)
         enddo
         do i=nnt,nct
-          stdforcsc(i)=stdfsc(iabs(itype(i)))
-     &                *dsqrt(gamsc(iabs(itype(i))))
+          stdforcsc(i)=stdfsc(itype(i))*dsqrt(gamsc(itype(i)))
         enddo
       endif
 c Open the pdb file for snapshotshots
@@ -1829,7 +1828,8 @@ c Removing the velocity of the center of mass
       endif
       if (.not.rest) then	  	
          call chainbuild
-         if(iranconf.ne.0) then
+         write (iout,*) "PREMINIM ",preminim
+         if(iranconf.ne.0 .or. preminim) then
           if (overlapsc) then 
            print *, 'Calling OVERLAP_SC'
            call overlap_sc(fail)
@@ -1851,14 +1851,18 @@ c Removing the velocity of the center of mass
            call minimize(etot,varia,iretcode,nfun)
            call var_to_geom(nvar,varia)
           endif
-          if(me.eq.king.or..not.out1file)
-     &       write(iout,*) 'SUMSL return code is',iretcode,' eval ',nfun
+          if(me.eq.king.or..not.out1file) then
+             write(iout,*) "Minimized energy is",etot
+             write(iout,*) 'SUMSL return code is',iretcode,' eval ',nfun
+             call etotal(potEcomp)
+             call enerprint(potEcomp)
+          endif
          endif
       endif 	  
       call chainbuild_cart
       call kinetic(EK)
       if (tbf) then
-        call verlet_bath(EK)
+        call verlet_bath
       endif	  
       kinetic_T=2.0d0/(dimen3*Rb)*EK
       if(me.eq.king.or..not.out1file)then