anchor.diff:          isym = 0
           syml = 0
           symlz = 0
  +
  +        if (present(store_in)) then
  +            store = store_in
  +        else
  +            store = .false.
  +        end if
   
           ! FCIDUMP file format is as follows:
   
  
positive.code: 
  contains
  
      subroutine read_in_fcidump
  
          ! 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 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 utils, only: get_free_unit
          use checking, only: check_allocate
          use errors, only: stop_all
  
          use, intrinsic :: iso_fortran_env
  
          ! System data
          ! We don't know how many orbitals we have until we read in the FCI
  
negative.code:         end if
  
          ! Set up symmetry information.
          if (store) call init_pg_symmetry()
  
          ! Initialise integral stores.
          if (store) call init_molecular_integrals()
  
          ! read integrals and eigenvalues
          ios = 0
  