anchor.diff:          integer,allocatable::natoms,nbasis,nbasisUse,icharge,multiplicity,nelectrons,icgu, &
             NFC,NFV,ITran,IDum9,NShlAO,NPrmAO,NShlDB,NPrmDB,NBTot
           integer,dimension(:),allocatable::atomicNumbers,atomTypes,basisFunction2Atom, &
  -          IBasisFunctionType
  -        real,dimension(:),allocatable::atomicCharges,atomicWeights,cartesians
  +          IBasisFunctionType,IgaussianScalars
  +        real,dimension(:),allocatable::atomicCharges,atomicWeights,cartesians   !hph,gaussianScalarsR
         Contains
           procedure,pass::OpenFile       => MQC_Gaussian_Unformatted_Matrix_Open
           procedure,pass::CloseFile      => MQC_Gaussian_Unformatted_Matrix_Close
  
positive.code:         fileinfo%Multiplicity,fileinfo%nelectrons,Len12L,Len4L,  &
          IOpCl,fileinfo%ICGU
  !
        return
        end subroutine MQC_Gaussian_Unformatted_Matrix_Read_Header
  
  
negative.code: !             This is only an issue for nonsymmetric matrices stored in symmetric form.
                matrixOut = transpose(matrixOut)
                deallocate(complexTmp)
              case('MIXED')
                write(*,1020)
                write(*,1040)' Hrant - LR   = ',LR
  