X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2FMD.F90;h=52d7ec3e7190cf52ab5e8802eca0386ee4847201;hb=bc23440fbe68672d430f71f22f46b11265f003db;hp=b60dd2726c1f43369a22175503161a2cb0e02c0d;hpb=3d3780d89d2885ee8700324ef832ef07ee35c45b;p=unres4.git

diff --git a/source/unres/MD.F90 b/source/unres/MD.F90
index b60dd27..52d7ec3 100644
--- a/source/unres/MD.F90
+++ b/source/unres/MD.F90
@@ -2058,6 +2058,7 @@
           do j=1,3    
             adt=(d_a_old(j,inres)+d_af_work(ind+j))*d_time
             adt2=0.5d0*adt+sqrt13*d_as_work(ind+j)*d_time
+!            write(iout,*) "adt",adt,"ads",adt2
             dc(j,inres)=dc_old(j,inres)+(d_t_old(j,inres)+adt2)*d_time
             d_t_new(j,inres)=d_t_old(j,inres)+0.5d0*adt
             d_t(j,inres)=d_t_old(j,inres)+adt
@@ -2531,9 +2532,14 @@
        write (iout,*) "vcm right after adjustment:"
        write (iout,*) (vcm(j),j=1,3) 
       endif
-      if (.not.rest) then	  
-  122   continue        	
+
+
+  
 !         call chainbuild
+
+      if ((.not.rest).or.(forceminim)) then	  	
+         if (forceminim) call chainbuild_cart
+  122   continue                
          if(iranconf.ne.0 .or.indpdb.gt.0.and..not.unres_pdb .or.preminim) then
           if (overlapsc) then 
            print *, 'Calling OVERLAP_SC'
@@ -5168,11 +5174,14 @@
 !el       ginvfric(2*nres,2*nres)	!maxres2=2*maxres
 !el      common /przechowalnia/ ginvfric
       
-      logical :: lprn = .false., checkmode = .false.
+      logical :: lprn, checkmode
       integer :: i,j,ind,k,nres2,nres6,mnum
       nres2=2*nres
       nres6=6*nres
-
+      lprn=.false.
+      checkmode=.false.
+!      if (large) lprn=.true.
+!      if (large) checkmode=.true.
       if(.not.allocated(gamvec)) allocate(gamvec(nres6)) !(MAXRES6)
       if(.not.allocated(ginvfric)) allocate(ginvfric(nres2,nres2)) !maxres2=2*maxres
       do i=0,nres2