anchor.diff:              tmpMatrixAlpha)
           elseIf(fileinfo%isGeneral()) then
             call fileInfo%getArray('OVERLAP',tmpMatrixAlpha,foundOut=found)
  -          call mqc_error_l('OVERLAP not present on file - aborting',6,'found',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)
  
positive.code:           call mqc_integral_allocate(est_integral,'density','space',tmpMatrixAlpha)
          elseIf(fileinfo%isUnrestricted()) then
            call fileInfo%getArray('ALPHA SCF DENSITY MATRIX',tmpMatrixAlpha,foundOut=found)
            call mqc_error_l('ALPHA SCF DENSITY MATRIX not present on file - aborting',6,'found',found)
            call fileInfo%getArray('BETA SCF DENSITY MATRIX',tmpMatrixBeta,foundOut=found)
            call mqc_error_l('BETA SCF DENSITY MATRIX not present on file - aborting',6,'found',found)
            if(MQC_Matrix_HaveComplex(tmpMatrixAlpha)) then
              call mqc_matrix_symm2full(tmpMatrixAlpha,'hermitian')
              tmpMatrixAlpha = transpose(tmpMatrixAlpha)
  
negative.code:       case('fock')
          if(fileinfo%isRestricted()) then
            call fileInfo%getArray('ALPHA FOCK MATRIX',tmpMatrixAlpha,foundOut=found)
            if(.not.found) call mqc_error_l('ALPHA FOCK MATRIX not present on file - aborting',6,'found',found)
            if(MQC_Matrix_HaveComplex(tmpMatrixAlpha)) then
              call mqc_matrix_symm2full(tmpMatrixAlpha,'hermitian')
              tmpMatrixAlpha = transpose(tmpMatrixAlpha)
  