Merge branch 'adasko' of mmka.chem.univ.gda.pl:unres into adasko
[unres.git] / source / unres / src_MD-M / kinetic_lesyng.f
index 8535f5d..db959b3 100644 (file)
@@ -44,7 +44,7 @@ c to the velocities of the first Calpha.
         incr(j)=d_t(j,0)
       enddo
       do i=nnt,nct
-        iti=itype(i)
+        iti=iabs(itype(i))
         if (itype(i).eq.10) then
           do j=1,3
             v(j)=incr(j)
@@ -81,7 +81,7 @@ c       write(iout,*) 'KEr_p', KEr_p
 c  The rotational part of the side chain virtual bond
        KEr_sc=0.0D0
        do i=nnt,nct
-        iti=itype(i)
+        iti=iabs(itype(i))
         if (itype(i).ne.10) then
         do j=1,3
          incr(j)=d_t(j,nres+i)