X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD-M%2Fgradient_p.F;h=7295d2ee648ef9c135a9ded784dd59a42ded0192;hb=f94ff29ee2f2c001a0c55100f253cd70a93c33cd;hp=b5ba82afdb8e02533ba03d05723e4f85bc2c913e;hpb=73c361f9b6f2b88612675f975588ee88dbba3a9c;p=unres.git

diff --git a/source/unres/src_MD-M/gradient_p.F b/source/unres/src_MD-M/gradient_p.F
index b5ba82a..7295d2e 100644
--- a/source/unres/src_MD-M/gradient_p.F
+++ b/source/unres/src_MD-M/gradient_p.F
@@ -292,11 +292,9 @@ c If performing constraint dynamics, add the gradients of the constraint energy
          do i=1,nres-3
            gloc(i,icg)=gloc(i,icg)+dugamma(i)
          enddo
-      write(iout,*) "TU JESTEM?"
          do i=1,nres-2
            gloc(nphi+i,icg)=gloc(nphi+i,icg)+dutheta(i)
          enddo
-        write(iout,*) "TU JESTEM?"
       endif 
 #ifdef TIMING
       time01=MPI_Wtime()
@@ -310,7 +308,7 @@ cd      write(iout,*) 'calling int_to_cart'
 #ifdef DEBUG
       write (iout,*) "gcart, gxcart, gloc before int_to_cart"
 #endif
-      do i=1,nct
+      do i=0,nct
         do j=1,3
           gcart(j,i)=gradc(j,i,icg)
           gxcart(j,i)=gradx(j,i,icg)
@@ -351,7 +349,7 @@ C-------------------------------------------------------------------------
 C
 C Initialize Cartesian-coordinate gradient
 C
-      do i=1,nres
+      do i=-1,nres
 	do j=1,3
 	  gvdwx(j,i)=0.0D0
           gradx_scp(j,i)=0.0D0
@@ -383,11 +381,25 @@ C
           gradx(j,i,icg)=0.0d0
           gscloc(j,i)=0.0d0
           gsclocx(j,i)=0.0d0
+          gliptranc(j,i)=0.0d0
+          gliptranx(j,i)=0.0d0
+          gradafm(j,i)=0.0d0
           do intertyp=1,3
            gloc_sc(intertyp,i,icg)=0.0d0
           enddo
         enddo
       enddo
+c
+c Initialize the gradients of local restraints
+c
+      do i=1,nres
+        dutheta(i)=0.0d0
+        dugamma(i)=0.0d0
+        do j=1,3
+          duscdiff(j,i)=0.0d0
+          duscdiffx(j,i)=0.0d0
+        enddo
+      enddo
 C
 C Initialize the gradient of local energy terms.
 C