anchor.diff:            call mqc_integral_allocate(est_integral,'core hamiltonian','space',tmpMatrixAlpha)
           elseIf(fileinfo%isUnrestricted()) then
             call fileInfo%getArray('CORE HAMILTONIAN ALPHA',tmpMatrixAlpha,foundOut=found)
  -          call mqc_error_l('CORE HAMILTONIAN ALPHA not present on file - aborting',6,'found',found)
  +          if(.not.found) call mqc_error_l('CORE HAMILTONIAN ALPHA not present on file - aborting',6,'found',found)
             call fileInfo%getArray('CORE HAMILTONIAN BETA',tmpMatrixBeta,foundOut=found)
  -          call mqc_error_l('CORE HAMILTONIAN BETA not present on file - aborting',6,'found',found)
  +          if(.not.found) call mqc_error_l('CORE HAMILTONIAN BETA not present on file - aborting',6,'found',found)
             if(MQC_Matrix_HaveComplex(tmpMatrixAlpha)) then
               call mqc_matrix_symm2full(tmpMatrixAlpha,'hermitian')
               tmpMatrixAlpha = transpose(tmpMatrixAlpha)
  
positive.code:       case('mo coefficients')
          if(fileinfo%isRestricted()) then
            call fileInfo%getArray('ALPHA MO COEFFICIENTS',tmpMatrixAlpha,foundOut=found)
            call mqc_error_l('ALPHA MO COEFFICIENTS not present on file - aborting',6,'found',found)
            call mqc_integral_allocate(est_integral,'mo coefficients','space',tmpMatrixAlpha)
          elseIf(fileinfo%isUnrestricted()) then
            call fileInfo%getArray('ALPHA MO COEFFICIENTS',tmpMatrixAlpha,foundOut=found)
            call mqc_error_l('ALPHA MO COEFFICIENTS not present on file - aborting',6,'found',found)
            call fileInfo%getArray('BETA MO COEFFICIENTS',tmpMatrixBeta,foundOut=found)
            call mqc_error_l('BETA MO COEFFICIENTS not present on file - aborting',6,'found',found)
            call mqc_integral_allocate(est_integral,'mo coefficients','spin',tmpMatrixAlpha, &
              tmpMatrixBeta)
          elseIf(fileinfo%isGeneral()) then
            call fileInfo%getArray('ALPHA MO COEFFICIENTS',tmpMatrixAlpha,foundOut=found)
            call mqc_error_l('ALPHA MO COEFFICIENTS not present on file - aborting',6,'found',found)
            nBasis = fileInfo%getVal('nBasis')
            call mqc_matrix_spinBlockGHF(tmpMatrixAlpha,fileInfo%getVal('nElectrons'), &
              fileInfo%getVal('multiplicity'),elist)
  
negative.code:           call mqc_integral_allocate(est_integral,'overlap','space',tmpMatrixAlpha)
          elseIf(fileinfo%isUnrestricted()) then
            call fileInfo%getArray('OVERLAP',tmpMatrixAlpha,foundOut=found)
            if(.not.found) call mqc_error_l('OVERLAP not present on file - aborting',6,'found',found)
            if(MQC_Matrix_HaveComplex(tmpMatrixAlpha)) then
              call mqc_matrix_symm2full(tmpMatrixAlpha,'hermitian')
              tmpMatrixAlpha = transpose(tmpMatrixAlpha)
  