anchor.diff:  
   contains
   
  -    subroutine read_in_fcidump
  +    subroutine read_in_fcidump(store_in)
   
           ! Read in a FCIDUMP file, which contains the kinetic and coulomb
           ! integrals, one-particle eigenvalues and other system information.
           ! 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 basis, only: write_basis_fn
  +        use basis, only: nbasis, basis_fns, init_basis_fn, end_basis_fns, write_basis_fn
           use molecular_integrals, only: init_molecular_integrals, store_one_e_int_mol, &
                                          store_two_e_int_mol
           use point_group_symmetry, only: init_pg_symmetry
  -        use system, only: fcidump, uhf
  +        use system, only: fcidump, uhf, ecore
   
           use utils, only: get_free_unit
           use checking, only: check_allocate
           use errors, only: stop_all
   
           use, intrinsic :: iso_fortran_env
  +
  +        logical, intent(in), optional :: store_in
  +        logical :: store
   
           ! System data
           ! We don't know how many orbitals we have until we read in the FCI
  
positive.code:         end if
  
          ! Set up symmetry information.
          call init_pg_symmetry()
  
          ! Initialise integral stores.
          call init_molecular_integrals()
  
          ! read integrals and eigenvalues
          ios = 0
  
negative.code:             call write_basis_fn(basis_fns(i), new_line=.true.)
          end do
  
          if (.not.store) then
              ! Should tidy up and deallocate everything we allocated.
              call end_basis_fns()
          end if
  
      end subroutine read_in_fcidump
  
  end module read_in_system
  