X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc-HCD-5D%2Flagrangian_lesyng.F;h=4230e10f84db3b7564511a55ce88a1d41a403eff;hb=b6b1679cfe492bd5b73e03410f82917bd9c3aa7b;hp=f57a432b4b4cb59ed4b8a5fe01c34346675d5f15;hpb=fb31d9cf50bce06dbfae62c92c11f3e9c836f405;p=unres.git

diff --git a/source/unres/src-HCD-5D/lagrangian_lesyng.F b/source/unres/src-HCD-5D/lagrangian_lesyng.F
index f57a432..4230e10 100644
--- a/source/unres/src-HCD-5D/lagrangian_lesyng.F
+++ b/source/unres/src-HCD-5D/lagrangian_lesyng.F
@@ -895,6 +895,7 @@ c---------------------------------------------------------------------------
       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
@@ -919,7 +920,7 @@ Compute accelerations in Calpha and SC
           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
@@ -944,6 +945,7 @@ C Conevert d_a to virtual-bon-vector basis
         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)