anchor.diff:              if (max(ii,jj,aa,bb) <= nbasis) then
   
                   ! Have integrals involving only core or active orbitals.
  -                if (ii == 0 .and. jj == 0 .and. aa == 0 .and. bb == 0) then
  +                if (i == 0 .and. j == 0 .and. a == 0 .and. b == 0) then
   
                       ! Nuclear energy.
                       Ecore = Ecore + x
   
  -                else if (ii > 0 .and. jj == 0 .and. aa == 0 .and. bb == 0) then
  +                else if (ii > 0 .and. j == 0 .and. a == 0 .and. b == 0) then
   
                       ! \epsilon_i
                       if (uhf) then
  
positive.code:         integer :: ir, ios, ierr
          logical :: t_exists
          integer :: active_basis_offset, rhf_fac
          integer, allocatable :: seen_ij(:), seen_iab(:,:)
  
          namelist /FCI/ norb, nelec, ms2, orbsym, uhf, isym, syml, symlz
  
  
negative.code:         call check_allocate('basis_fns', nbasis, ierr)
  
          ! Set up basis functions.
          call init_basis_fns_read_in(norb, uhf, orbsym(active_basis_offset/rhf_fac+1:norb+active_basis_offset/rhf_fac), basis_fns) 
  
          ! Was a symmetry found for all basis functions?  If not, then we must
          ! turn symmetry off.
  