anchor.diff:        real,allocatable,dimension(:)::arrayTmp
         complex(kind=8),allocatable,dimension(:)::complexTmp
         character(len=256)::my_filename,errorMsg
  -      logical::DEBUG=.true.,ok,found
  +      logical::DEBUG=.false.,ok,found
   !
   !
   !     Format statements.
  
positive.code:           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) )
  
negative.code: !           'core hamiltonian'   write the core hamiltonian.
  !           'fock'               write the fock matrix.
  !           'density'            write the density matrix.
  !           'scf density'        write the density matrix with the SCF label
  !           'overlap'            write the overlap matrix.
  !           'wavefunction'       export the wavefunction object.
  !
  