- 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
+C if (SELFGUIDE.gt.0) then
+C distance=0.0
+C do j=1,3
+C vec_afm(j)=c(j,afmend)-c(j,afmbeg)
+C distance=distance+vec_afm(j)**2
+C enddo
+C distance=dsqrt(distance)
+C do j=1,3
+C d_t_work_new(j+(afmbeg-1)*3)=-velAFMconst*vec_afm(j)/distance
+C d_t_work_new(j+(afmend-1)*3)=velAFMconst*vec_afm(j)/distance
+C write(iout,*) "myvel",d_t_work_new(j+(afmbeg-1)*3),
+C & d_t_work_new(j+(afmend-1)*3)
+C enddo