anchor.diff:           write(fn,'(E22.15)')enrep
            end if
            close(fn)
  -c.. Now that the integrals are in g_lu1e and g_lu2e, fully packed. 
  +c.. Now that the integrals are in KINT1 and KINT2, fully packed. 
   
          ELSE IF (INTIMP .EQ. 5) THEN
   *
  
positive.code:        open(fn,file=filename,form='formatted',status='unknown')
         rewind (fn)
  c.. generate GA to store 1-e ints for LUCIA
           if(.not.ga_create(mt_dbl,1,nint1,'lucia_1e',minchunk,1,g_lu1e))
       &       call errquit('scf: cannot create lucia_1e',0,GA_ERR)
           if(.not.ga_create(mt_dbl,1,nint2,'lucia_2e',minchunk,1,g_lu2e))
       &       call errquit('scf: cannot create lucia_2e',0,GA_ERR)
  c..  1-e integrals
           g_hcore = ga_create_atom_blocked(geom,basis, 'rohf: hcore')
  
negative.code:          call ga_get(g_hcore,i,i,j,j,int1e,1)
  c        write(6,'(E22.15)')int1e
  c        write(fn,'(E22.15)')int1e
           call ga_put(KINT1,1,1,ind1e,ind1e,int1e,1)
            end do
           end do 
  c.. 2e ints
  