New valence-torsionals completed

[unres.git] / source / unres / src_MD-M / checkder_p.F

diff --git a/source/unres/src_MD-M/checkder_p.F b/source/unres/src_MD-M/checkder_p.F
index 270f4cc..80dd213 100644 (file)
--- a/source/unres/src_MD-M/checkder_p.F
+++ b/source/unres/src_MD-M/checkder_p.F
@@ -2,6 +2,7 @@
 C Check the gradient of Cartesian coordinates in internal coordinates.
       implicit real*8 (a-h,o-z)
       include 'DIMENSIONS'
+      include 'COMMON.CONTROL'
       include 'COMMON.IOUNITS'
       include 'COMMON.VAR'
       include 'COMMON.CHAIN'
@@ -14,8 +15,9 @@ C Check the gradient of Cartesian coordinates in internal coordinates.
 * Check the gradient of the virtual-bond and SC vectors in the internal
 * coordinates.
 *    
-      aincr=1.5d-7  
-      aincr2=2.5d-8   
+      print '("Calling CHECK_ECART",1pd12.3)',aincr
+      write (iout,'("Calling CHECK_ECART",1pd12.3)') aincr
+      aincr2=0.5d0*aincr
       call cartder
       write (iout,'(a)') '**************** dx/dalpha'
       write (iout,'(a)')
@@ -175,6 +177,7 @@ C----------------------------------------------------------------------------
 C Check the gradient of the energy in Cartesian coordinates. 
       implicit real*8 (a-h,o-z)
       include 'DIMENSIONS'
+      include 'COMMON.CONTROL'
       include 'COMMON.CHAIN'
       include 'COMMON.DERIV'
       include 'COMMON.IOUNITS'
@@ -191,8 +194,8 @@ C Check the gradient of the energy in Cartesian coordinates.
       nf=0
       nfl=0                
       call zerograd
-      aincr=1.0D-7
-      print '(a)','CG processor',me,' calling CHECK_CART.'
+      print '("Calling CHECK_ECART",1pd12.3)',aincr
+      write (iout,'("Calling CHECK_ECART",1pd12.3)') aincr
       nf=0
       icall=0
       call geom_to_var(nvar,x)
@@ -282,8 +285,10 @@ c      rlambd=0.3d0
 c      call intcartderiv
 c      call checkintcartgrad
       call zerograd
-      aincr=8.0D-7
-      write(iout,*) 'Calling CHECK_ECARTINT.'
+c      aincr=8.0D-7
+c      aincr=1.0D-7
+      print '("Calling CHECK_ECARTINT",1pd12.3)',aincr
+      write (iout,'("Calling CHECK_ECARTINT",1pd12.3)') aincr
       nf=0
       icall=0
       call geom_to_var(nvar,x)
@@ -298,8 +303,10 @@ c      call checkintcartgrad
         write (iout,*) "exit cartgrad"
         call flush(iout)
         icall =1
+        write (iout,*) "gcard and gxcart"
         do i=1,nres
-          write (iout,'(i5,3f10.5)') i,(gradxorr(j,i),j=1,3)
+          write (iout,'(i5,3f10.5,5x,3f10.5)') i,(gcart(j,i),j=1,3),
+     &      (gxcart(j,i),j=1,3)
         enddo
         do j=1,3
           grad_s(j,0)=gcart(j,0)
@@ -628,6 +635,7 @@ c----------------------------------------------------------------------------
 C Check the gradient of energy in internal coordinates.
       implicit real*8 (a-h,o-z)
       include 'DIMENSIONS'
+      include 'COMMON.CONTROL'
       include 'COMMON.CHAIN'
       include 'COMMON.DERIV'
       include 'COMMON.IOUNITS'
@@ -642,8 +650,9 @@ C Check the gradient of energy in internal coordinates.
       character*6 key
       external fdum
       call zerograd
-      aincr=1.0D-7
-      print '(a)','Calling CHECK_INT.'
+c      aincr=1.0D-7
+      print '("Calling CHECK_INT",1pd12.3)',aincr
+      write (iout,'("Calling CHECK_INT",1pd12.3)') aincr
       nf=0
       nfl=0
       icg=1
@@ -674,12 +683,12 @@ cd    write (iout,'(i3,1pe14.4)') (i,gana(i),i=1,nvar)
         xi=x(i)
         x(i)=xi-0.5D0*aincr
         call var_to_geom(nvar,x)
-        call chainbuild
+        call chainbuild_extconf
         call etotal(energia1(0))
         etot1=energia1(0)
         x(i)=xi+0.5D0*aincr
         call var_to_geom(nvar,x)
-        call chainbuild
+        call chainbuild_extconf
         call etotal(energia2(0))
         etot2=energia2(0)
         gg(i)=(etot2-etot1)/aincr