ind=ind+3
enddo
do i=nnt,nct
- if (itype(i).ne.10) then
+ if ((itype(i).ne.10).and.(itype(i).ne.ntyp1)) then
do j=1,3
d_t_work(ind+j)=d_t(j,i+nres)
enddo
ind=ind+3
enddo
do i=nnt,nct
- if (itype(i).ne.10) then
+ if ((itype(i).ne.10).and.(itype(i).ne.ntyp1)) then
do j=1,3
friction(j,i+nres)=fric_work(ind+j)
enddo
stochforc(j,0)=ff(j)+force(j,nnt+nres)
enddo
do i=nnt,nct
- if (itype(i).ne.10) then
+ if ((itype(i).ne.10).and.(itype(i).ne.ntyp1)) then
do j=1,3
stochforc(j,i+nres)=force(j,i+nres)
enddo
ind=ind+3
enddo
do i=nnt,nct
- if (itype(i).ne.10) then
+ if ((itype(i).ne.10).and.(itype(i).ne.ntyp1)) then
do j=1,3
stochforcvec(ind+j)=stochforc(j,i+nres)
enddo