X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD-M%2FMD_A-MTS.F;fp=source%2Funres%2Fsrc_MD-M%2FMD_A-MTS.F;h=42b3af40c1d307c6b697be93a23f8171f06e8f18;hb=f038962fbb96e0c2c6f1ccb910373ceddd5b387b;hp=bb8709057f33dc8319a8a020e6e51467a1e4bfc6;hpb=0ef28a0babbfafd06d3977c622ecbe98a5f41e86;p=unres.git

diff --git a/source/unres/src_MD-M/MD_A-MTS.F b/source/unres/src_MD-M/MD_A-MTS.F
index bb87090..42b3af4 100644
--- a/source/unres/src_MD-M/MD_A-MTS.F
+++ b/source/unres/src_MD-M/MD_A-MTS.F
@@ -1798,7 +1798,7 @@ c-----------------------------------------------------------
       include 'COMMON.IOUNITS'
       include 'COMMON.NAMES'
       include 'COMMON.TIME1'
-      double precision xv,sigv,lowb,highb
+      double precision xv,sigv,lowb,highb,vec_afm(3)
 c Generate random velocities from Gaussian distribution of mean 0 and std of KT/m 
 c First generate velocities in the eigenspace of the G matrix
 c      write (iout,*) "Calling random_vel dimen dimen3",dimen,dimen3
@@ -1812,10 +1812,27 @@ c      call flush(iout)
           lowb=-5*sigv
           highb=5*sigv
           d_t_work_new(ii)=anorm_distr(xv,sigv,lowb,highb)
+
 c          write (iout,*) "i",i," ii",ii," geigen",geigen(i),
 c     &      " d_t_work_new",d_t_work_new(ii)
         enddo
       enddo
+       if (SELFGUIDE.gt.0) then
+       distance=0.0
+       do j=1,3
+       vec_afm(j)=c(j,afmend)-c(j,afmbeg)  
+       distance=distance+vec_afm(j)**2
+       enddo
+       distance=dsqrt(distance)
+       do j=1,3
+         d_t_work_new(j+(afmbeg-1)*3)=-velAFMconst*vec_afm(j)/distance
+         d_t_work_new(j+(afmend-1)*3)=velAFMconst*vec_afm(j)/distance
+         write(iout,*) "myvel",d_t_work_new(j+(afmbeg-1)*3),
+     &    d_t_work_new(j+(afmend-1)*3)
+       enddo
+
+       endif
+
 c diagnostics
 c      Ek1=0.0d0
 c      ii=0