X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD%2Fintcartderiv.F;h=c22054007e6da27e700a4de033e76349bb2e2cd7;hb=4866688f3356059b9a933ec2be6da24e6784fa9f;hp=dd4db242526e19cf929bf818510c67a0f2b1d74d;hpb=87c21417a317385a49127c6e1a9b4a74541a6c6e;p=unres.git

diff --git a/source/unres/src_MD/intcartderiv.F b/source/unres/src_MD/intcartderiv.F
index dd4db24..c220540 100644
--- a/source/unres/src_MD/intcartderiv.F
+++ b/source/unres/src_MD/intcartderiv.F
@@ -151,7 +151,7 @@ c   Obtaining the gamma derivatives from cosine derivative
 
 Calculate derivative of Tauangle
 #ifdef PARINTDER
-      do i=iphi1_start-1,iphi1_end
+      do i=itau_start,itau_end
 #else
       do i=3,nres
 #endif
@@ -226,7 +226,7 @@ c        enddo
       enddo
 CC Second case Ca...Ca...Ca...SC
 #ifdef PARINTDER
-      do i=iphi1_start,iphi1_end
+      do i=itau_start,itau_end
 #else
       do i=4,nres
 #endif
@@ -291,7 +291,7 @@ c   Obtaining the gamma derivatives from cosine derivative
      &     dcosomicron(j,1,2,i)-fac0*(dc_norm(j,i-3)-scalp*
      &     dc_norm(j,i-1+nres))/vbld(i-1+nres)
            dtauangle(j,2,3,i)=-1/sing*dcostau(j,2,3,i)
-        write(iout,*) i,j,"else", dtauangle(j,2,3,i) 
+c        write(iout,*) i,j,"else", dtauangle(j,2,3,i) 
          enddo
         endif                                                                                            
       enddo
@@ -300,7 +300,7 @@ c   Obtaining the gamma derivatives from cosine derivative
 CCC third case SC...Ca...Ca...SC
 #ifdef PARINTDER
 
-      do i=iphi1_start-1,iphi1_end
+      do i=itau_start,itau_end
 #else
       do i=3,nres
 #endif
@@ -364,7 +364,7 @@ c   Obtaining the gamma derivatives from cosine derivative
      &     dcosomicron(j,1,2,i)-fac0*(dc_norm2(j,i-2+nres)-scalp*
      &     dc_norm(j,i-1+nres))/vbld(i-1+nres)
            dtauangle(j,3,3,i)=-1/sing*dcostau(j,3,3,i)
-          write(iout,*) "else",i 
+c          write(iout,*) "else",i 
          enddo
         endif                                                                                            
       enddo
@@ -458,12 +458,12 @@ c   obtaining the derivatives of omega from cosines
 #if defined(MPI) && defined(PARINTDER)
       if (nfgtasks.gt.1) then
 #ifdef DEBUG
-cd      write (iout,*) "Gather dtheta"
+       write (iout,*) "Gather dtheta"
 cd      call flush(iout)
-      write (iout,*) "dtheta before gather"
-      do i=1,nres
-        write (iout,'(i3,3(3f8.5,3x))') i,((dtheta(j,k,i),k=1,3),j=1,2)
-      enddo
+c      write (iout,*) "dtheta before gather"
+c      do i=1,nres
+c        write (iout,'(i3,3(3f8.5,3x))') i,((dtheta(j,k,i),k=1,3),j=1,2)
+c      enddo
 #endif
       call MPI_Gatherv(dtheta(1,1,ithet_start),ithet_count(fg_rank),
      &  MPI_THET,dtheta(1,1,1),ithet_count(0),ithet_displ(0),MPI_THET,