anchor.diff:            end do find_bin
   
             ! Print filter information
  -          write(UNIT=UNIT_TALLY, FMT='(1X,2A,1X,A)') space(1:indent), &
  +          write(UNIT=UNIT_TALLY, FMT='(1X,2A,1X,A)') repeat(" ", indent), &
                  trim(filter_name(j)), trim(get_label(t, j, bins(j)))
   
             ! Determine scoring index for this bin combination -- note that unlike
  
positive.code:           indent = indent + 2
            do k = 1, t % n_macro_bins
               write(UNIT=UNIT_TALLY, FMT='(1X,2A,1X,A,"+/- ",A)') & 
                    space(1:indent), macro_name(abs(t % macro_bins(k) % scalar)), &
                    real_to_str(t % scores(score_index,k) % val), &
                    trim(real_to_str(t % scores(score_index,k) % val_sq))
            end do
  
negative.code:                 else
                     if (has_filter(j)) then
                        ! Print current filter information
                        write(UNIT=UNIT_TALLY, FMT='(1X,2A,1X,A)') repeat(" ", indent), &
                             trim(filter_name(j)), trim(get_label(t, j, bins(j)))
                        indent = indent + 2
                     end if
  