update 5D

[unres.git] / source / unres / src-HCD-5D / MD_A-MTS.F

diff --git a/source/unres/src-HCD-5D/MD_A-MTS.F b/source/unres/src-HCD-5D/MD_A-MTS.F
index ca52aaa..08852ba 100644 (file)
--- a/source/unres/src-HCD-5D/MD_A-MTS.F
+++ b/source/unres/src-HCD-5D/MD_A-MTS.F
@@ -1256,7 +1256,6 @@ c      write (iout,*) "friction accelerations"
       call fivediaginv_mult(dimen,fric_work, d_af_work)
 c      write (iout,*) "stochastic acceleratios"
       call fivediaginv_mult(dimen,stochforcvec, d_as_work)
-c      write (iout,*) "Leaving sddir_precalc"
 #else
       call ginv_mult(fric_work, d_af_work)
       call ginv_mult(stochforcvec, d_as_work)
@@ -1266,6 +1265,7 @@ c      write (iout,*) "Leaving sddir_precalc"
       write (iout,'(3f10.5)') (d_af_work(i),i=1,dimen3)
       write (iout,*) "d_as_work"
       write (iout,'(3f10.5)') (d_as_work(i),i=1,dimen3)
+      write (iout,*) "Leaving sddir_precalc"
 #endif
       return
       end
@@ -2187,7 +2187,7 @@ c-----------------------------------------------------------
       double precision xv,sigv,lowb,highb,vec_afm(3),Ek1,Ek2,Ek3,aux
       integer i,ii,j,k,l,ind
       double precision anorm_distr
-      logical lprn /.true./
+      logical lprn /.false./
 #ifdef FIVEDIAG
       integer ichain,n,innt,inct,ibeg,ierr
       double precision work(8*maxres6)
@@ -2797,8 +2797,8 @@ c      enddo
         dc(j,0)=dc_work(j)
         d_t(j,0)=d_t_work(j)
       enddo
-      ind=3    
-      do i=nnt,nct-1   
+      ind=3
+      do i=nnt,nct-1
         do j=1,3
           dc(j,i)=dc_work(ind+j)
           d_t(j,i)=d_t_work(ind+j)