anchor.diff:        complex(kind=real64),allocatable,dimension(:)::compVectorTmp
         type(MQC_Matrix)::matrixInUse 
         character(len=256)::my_filename,my_storage
  -      logical::DEBUG=.true.,ok
  +      logical::DEBUG=.false.,ok
         Parameter(LenBuf=4000)
   !
   !     Format statements.
  
positive.code:           call mqc_error_a('Unknown wavefunction type in writeESTObj', 6, &
                 'my_integral_type', my_integral_type )
          endIf
        case('density')
          if(.not.(Present(est_integral))) call mqc_error_L('wrong EST type in writeESTOBj', 6, &
               'Present(est_integral)', Present(est_integral) )
          if(my_integral_type.eq.'space') then
  
negative.code:           call mqc_error_a('Unknown wavefunction type in writeESTObj', 6, &
                 'my_integral_type', my_integral_type )
          endIf
        case('density')
           if(.not.(Present(est_integral))) call mqc_error_L('wrong EST type in writeESTOBj', 6, &
                'Present(est_integral)', Present(est_integral) )
           if(my_integral_type.eq.'space') then
             call fileInfo%writeArray('ALPHA DENSITY MATRIX', &
               matrixIn=est_integral%getBlock('alpha'))
           elseIf(my_integral_type.eq.'spin') then
             call fileInfo%writeArray('ALPHA DENSITY MATRIX', &
               matrixIn=est_integral%getBlock('alpha'))
             call fileInfo%writeArray('BETA DENSITY MATRIX', &
               matrixIn=est_integral%getBlock('beta'))
           elseIf(my_integral_type.eq.'general') then
             call mqc_matrix_undoSpinBlockGHF(est_integral,tmpMatrix)
             if(.not.mqc_matrix_haveComplex(tmpMatrix)) call MQC_Matrix_Copy_Real2Complex(tmpMatrix)
             call fileInfo%writeArray('ALPHA DENSITY MATRIX',matrixIn=tmpMatrix)
           else
             call mqc_error_a('Unknown wavefunction type in writeESTObj', 6, &
                  'my_integral_type', my_integral_type )
           endIf
        case('overlap')
          if(.not.(Present(est_integral))) call mqc_error_L('wrong EST type in writeESTOBj', 6, &
               'Present(est_integral)', Present(est_integral) )
  