checkgrad_inc added to input for checkgrad cart

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

diff --git a/source/unres/src_MD/checkder_p.F b/source/unres/src_MD/checkder_p.F
index 4d0379e..05ea860 100644 (file)
--- a/source/unres/src_MD/checkder_p.F
+++ b/source/unres/src_MD/checkder_p.F
@@ -284,8 +284,9 @@ C Check the gradient of the energy in Cartesian coordinates.
 c      call intcartderiv
 c      call checkintcartgrad
       call zerograd
-      aincr=1.0D-5
-      write(iout,*) 'Calling CHECK_ECARTINT.'
+c      aincr=1.0D-4
+      aincr=checkgrad_inc
+      write(iout,*) 'Calling CHECK_ECARTINT.',aincr
       nf=0
       icall=0
       call geom_to_var(nvar,x)
@@ -395,6 +396,7 @@ c          call int_from_cart1(.false.)
           if (.not.split_ene) then
             call etotal(energia1(0))
             etot1=energia1(0)
+c            write (iout,*) "i",i," etot",etot," etot1",etot1
           else
 !- split gradient
             call etotal_long(energia1(0))
@@ -441,6 +443,7 @@ c          write (iout,*)
           if (.not.split_ene) then
             call etotal(energia1(0))
             etot1=energia1(0)
+c            write (iout,*) "i",i," etot",etot," etot1",etot1
           else
 !- split gradient
             call etotal_long(energia1(0))
@@ -686,7 +689,7 @@ cd    write (iout,'(i3,1pe14.4)') (i,gana(i),i=1,nvar)
         call etotal(energia2(0))
         etot2=energia2(0)
         gg(i)=(etot2-etot1)/aincr
-        write (iout,*) i,etot1,etot2
+c        write (iout,*) i,etot1,etot2
         x(i)=xi
       enddo
       write (iout,'(/2a)')' Variable        Numerical       Analytical',