anchor.diff:    maxrmsd(1) = 0.d0
     maxediff(1) = 0.d0
     
  +  write(6,'(A)') 'Starting analysis ...'
  +  
  +  tmpratio=0
  +  write(6,'(A)') '---|---|---|---|---|---|---|---|---|---|'
  +  
     do i=2,num_files
       rmsd = 1111.d0
       cand = 0
       call rmsd_set_anc_s1(struc_rep(i))
       !~write(*,*)
  +	curratio = i*40/num_files
  +	call progressbar_add(tmpratio,curratio)
  +	tmpratio = curratio
       do j=1,size(tl)
         ediff_tmp = abs(struc_get_energy(struc_rep(i)) - struc_get_energy(struc_rep(tl(j)%i(1))))*h2kcm
         if (ediff_tmp .gt. ediff_lim) cycle
  
positive.code:   stop
  666 write(*,*) 'bad input'
  
  end program
  
negative.code:   read(123,*) num_files
    read(123,*) rmsd_lim
    read(123,*) anceps
    
    write(6,'(A)') 'Reading structures ...'
    
    allocate(struc_rep(num_files))
    
  