double precision forces(3*ndim),accel(3,0:maxres2),rs(ndim),
& xsolv(ndim),d_a_vec(6*nres)
integer i,j,ind,ichain,n,iposc,innt,inct,inct_prev
+ accel=0.0d0
do j=1,3
Compute accelerations in Calpha and SC
do ichain=1,nchain
enddo
enddo
enddo
-C Conevert d_a to virtual-bon-vector basis
+C Convert d_a to virtual-bon-vector basis
#ifdef DEBUG
write (iout,*) "accel in CA-SC basis"
do i=1,nres
end if
enddo
accel(:,nres)=0.0d0
+ accel(:,nct)=0.0d0
accel(:,2*nres)=0.0d0
if (nnt.gt.1) then
accel(:,0)=accel(:,1)