time00=MPI_Wtime()
#endif
icg=1
+#ifdef DEBUG
+ write (iout,*) "Before sum_gradient"
+ do i=1,nres-1
+ write (iout,*) i," gradc ",(gradc(j,i,icg),j=1,3)
+ write (iout,*) i," gradx ",(gradx(j,i,icg),j=1,3)
+ enddo
+#endif
call sum_gradient
#ifdef TIMING
#endif
gel_loc_long(j,i)=0.0d0
ghpbc(j,i)=0.0D0
ghpbx(j,i)=0.0D0
+ gsaxsc(j,i)=0.0D0
+ gsaxsx(j,i)=0.0D0
gcorr3_turn(j,i)=0.0d0
gcorr4_turn(j,i)=0.0d0
gradcorr(j,i)=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