incr(j)=d_t(j,0)
enddo
do i=nnt,nct
- iti=iabs(itype(i))
+ iti=itype(i)
if (itype(i).eq.10) then
do j=1,3
v(j)=incr(j)
c The rotational part of the side chain virtual bond
KEr_sc=0.0D0
do i=nnt,nct
- iti=iabs(itype(i))
+ iti=itype(i)
if (itype(i).ne.10) then
do j=1,3
incr(j)=d_t(j,nres+i)