anchor.diff:              else
                 call mqc_error_l('errorMsg',6,'found',found)
               endIf
  -          endIf
  -          call fileInfo%getArray('BETA FOCK MATRIX',tmpMatrixBeta,foundOut=found)
  -          if(present(foundObj)) foundObj = found
  -          if(.not.found) then
  -            errorMsg = 'BETA FOCK MATRIX not present on file'
  -            if(present(foundObj)) then
  -              write(6,'(A)') errorMsg
  +          else
  +            call fileInfo%getArray('BETA FOCK MATRIX',tmpMatrixBeta,foundOut=found)
  +            if(present(foundObj)) foundObj = found
  +            if(.not.found) then
  +              errorMsg = 'BETA FOCK MATRIX not present on file'
  +              if(present(foundObj)) then
  +                write(6,'(A)') errorMsg
  +              else
  +                call mqc_error_l('errorMsg',6,'found',found)
  +              endIf
               else
  -              call mqc_error_l('errorMsg',6,'found',found)
  -            endIf
  -          endIf
  -          if(MQC_Matrix_HaveComplex(tmpMatrixAlpha)) then
  -            call mqc_matrix_symm2full(tmpMatrixAlpha,'hermitian')
  -            tmpMatrixAlpha = transpose(tmpMatrixAlpha)
  -          endIf
  -          if(MQC_Matrix_HaveComplex(tmpMatrixBeta)) then
  -            call mqc_matrix_symm2full(tmpMatrixBeta,'hermitian')
  -            tmpMatrixBeta = transpose(tmpMatrixBeta)
  -          endIf
  -          call mqc_integral_allocate(est_integral,'fock','spin',tmpMatrixAlpha, &
  -            tmpMatrixBeta)
  +              if(MQC_Matrix_HaveComplex(tmpMatrixAlpha)) then
  +                call mqc_matrix_symm2full(tmpMatrixAlpha,'hermitian')
  +                tmpMatrixAlpha = transpose(tmpMatrixAlpha)
  +              endIf
  +              if(MQC_Matrix_HaveComplex(tmpMatrixBeta)) then
  +                call mqc_matrix_symm2full(tmpMatrixBeta,'hermitian')
  +                tmpMatrixBeta = transpose(tmpMatrixBeta)
  +              endIf
  +              call mqc_integral_allocate(est_integral,'fock','spin',tmpMatrixAlpha, &
  +                tmpMatrixBeta)
  +            endIf
  +          endIf
           elseIf(fileinfo%isGeneral()) then
             call fileInfo%getArray('ALPHA FOCK MATRIX',tmpMatrixAlpha,foundOut=found)
             if(present(foundObj)) foundObj = found
  
positive.code:             else
                call mqc_error_l('errorMsg',6,'found',found)
              endIf
            endIf
            call fileInfo%getArray('CORE HAMILTONIAN BETA',tmpMatrixBeta,foundOut=found)
            if(present(foundObj)) foundObj = found
            if(.not.found) then
              errorMsg = 'CORE HAMILTONIAN BETA not present on file'
              if(present(foundObj)) then
                write(6,'(A)') errorMsg
              else
                call mqc_error_l('errorMsg',6,'found',found)
              endIf
            endIf
            if(MQC_Matrix_HaveComplex(tmpMatrixAlpha)) then
              call mqc_matrix_symm2full(tmpMatrixAlpha,'hermitian')
              tmpMatrixAlpha = transpose(tmpMatrixAlpha)
            endIf
            if(MQC_Matrix_HaveComplex(tmpMatrixBeta)) then
              call mqc_matrix_symm2full(tmpMatrixBeta,'hermitian')
              tmpMatrixBeta = transpose(tmpMatrixBeta)
            endIf
            call mqc_integral_allocate(est_integral,'core hamiltonian','spin',tmpMatrixAlpha, &
              tmpMatrixBeta)
          elseIf(fileinfo%isGeneral()) then
            call fileInfo%getArray('CORE HAMILTONIAN ALPHA',tmpMatrixAlpha,foundOut=found)
            if(present(foundObj)) foundObj = found
  
negative.code:             else
                call mqc_error_l('errorMsg',6,'found',found)
              endIf
            else
              if(MQC_Matrix_HaveComplex(tmpMatrixAlpha)) then
                call mqc_matrix_symm2full(tmpMatrixAlpha,'hermitian')
                tmpMatrixAlpha = transpose(tmpMatrixAlpha)
              endIf
              call mqc_integral_allocate(est_integral,'overlap','space',tmpMatrixAlpha)
            endIf
          elseIf(fileinfo%isUnrestricted()) then
            call fileInfo%getArray('OVERLAP',tmpMatrixAlpha,foundOut=found)
            if(present(foundObj)) foundObj = found
  