anchor.diff:  !           '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.
   !
  
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:       real,allocatable,dimension(:)::arrayTmp
        complex(kind=8),allocatable,dimension(:)::complexTmp
        character(len=256)::my_filename,errorMsg
        logical::DEBUG=.false.,ok,found
  !
  !
  !     Format statements.
  