anchor.diff:        write(*,*) '    '
         exit
       endif
  -    write(*,*) "      LCCD: ",iT2,Elccd,deltaE 
  - 
  +    if(CCD)then
  +      write(*,*) "      CCD: ",iT2,Elccd,deltaE
  +    else
  +      write(*,*) "      LCCD: ",iT2,Elccd,deltaE 
  +    endif
     enddo ! t2 loop
   
     elseif(spins == 2)then
  
positive.code:     Dijab = 1.0d0
      Elccd = 0.0d0
  
    do iT2=1,maxT ! T amplitude loop
      Eold  = Elccd
      t2old  = t2 !R0.5d0*(t2+t2old) ! uses mixing
      !T2 equations
      do i=MOZero,nOcc
        do j=MOZero,nOcc 
          do a=nOcc+1,nmo
  
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 
  