anchor.diff:          ! File format (partially) defined in Comp. Phys. Commun. 54 (1989) 75.
           ! See also notes below.
   
  +        use basis, only: nbasis, basis_fns, init_basis_fn
           use system, only: fcidump
   
           use utils, only: get_free_unit
  +        use checking, only: check_allocate
           use errors, only: stop_all
   
           ! System data
  -        integer :: norb, nelec, ms2, orbsym(1000), isym
  +        ! We don't know how many orbitals we have until we read in the FCI
  +        ! namelist, so have to hardcode the array sizes.
  +        ! It's reasonably safe to assume that we'll never use more than 1000
  +        ! orbitals!
  +        integer :: norb, nelec, ms2, orbsym(1000), isym, syml(1000), symlz(1000)
           logical :: uhf = .false.
   
           ! Integrals
  
positive.code:         real(dp) :: x
  
          ! reading in...
          integer :: ir, ios
          logical :: t_exists
  
          namelist /FCI/ norb, nelec, ms2, orbsym, uhf
  
          ! FCIDUMP file format is as follows:
  
  
negative.code: use const
  
  implicit none
  
  real(p) :: Ecore
  
  contains
  
  