anchor.diff:  
       call calc_Energy()
   
  +!    call calc_DynDamp()
  +
       dE = Etot - Eold     
   
       !  iSCF   Etot  dE   Drmsd
  -    write(*,'("SCF  ",I5,"  ",3(F18.10,"  "),F5.2)') iSCF, Etot, dE, Drmsd, Damp
  +    write(*,'("  SCF  ",I5,"  ",3(F18.10,"  "),F5.2)') iSCF, Etot, dE, Drmsd, Damp
   
       call check_Conv()
   
  
positive.code:   double precision              :: Par_ii(8,3)
  
    write(*,*) ""
    write(*,*) "      Extended Hückel..."
    write(*,*) ""
  
  !  allocate(Hii(dim_1e))
  
negative.code: !  write(*,*) "    Initial guess..."
  !  write(*,*) ""
  
  !  write(*,*) "guess = ", guess
  
    if(guess == "core")then
      write(*,*) ""
      write(*,*) "    Core initial guess..."
      write(*,*) ""
      do i=1,Spins
        Fock(:,:,i) = Hcore(:,:)
  