anchor.diff:    double precision          :: Ecctol = 1.0d-7
     double precision          :: Eold   = 0.0d0
     double precision          :: deltaE = 99.0d0
  -  double precision          :: in1
  +  double precision          :: in1,in2
   
     Elccd  = 0.0d0
   
  
positive.code:               enddo
              enddo
  
  !            write(*,*) "t2   ",p3,p4,h1,h2,AMO(p3,p4,h1,h2),AMO(h1,h2,p3,p4)
  !            Dijab = (Eps_SO(p3) + Eps_SO(p4) - Eps_SO(h1) - Eps_SO(h2))
              t2(a,b,i,j) =  1.0d0*in1/Dijab
  
negative.code:   use module_ints,      only : Index2e
    use module_wavefun,   only : Fock_to_MO,Fock_to_AO
    implicit none
    logical                   :: CCD = .true.
    integer                   :: iSpin,MOZero,iT2
    integer                   :: p3,p4,h1,h2,h5,h6,p5,p6
    integer                   :: i,j,a,b 
  