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: !           This CASE block uses NI, NR, N1-N5, and NRI to determine the data
  !           type (integer, real, etc.) and data structure (scalar, vector,
  !           matrix, etc.).
              select case(MQC_Gaussian_Unformatted_Matrix_Array_Type(NI,NR,N1,N2,N3,N4,N5,NRI,ASym))
              case('INTEGER-VECTOR')
                allocate(integerTmp(LR))
  
negative.code:         allocate(fileinfo%NBTot)
          fileinfo%NBTot = 0
        endIf
        nAt3 = fileinfo%natoms*3
        if(.not.allocated(fileinfo%atomicNumbers)) then
          allocate(fileinfo%atomicNumbers(fileinfo%natoms))
          fileinfo%atomicNumbers = 0
  