anchor.diff:    use cross_section, only: get_macro_xs
     use error,         only: fatal_error
     use global
  -  use mesh,          only: get_mesh_bin
  +  use mesh,          only: get_mesh_bin, bin_to_mesh_indices
     use output,        only: message
     use search,        only: binary_search
     use string,        only: int_to_str, real_to_str
  
positive.code:                    if (has_filter(j)) then
                        ! Print current filter information
                        write(UNIT=UNIT_TALLY, FMT='(1X,2A,1X,A)') space(1:indent), &
                             trim(filter_name(j)), trim(get_uid(t, j, bins(j)))
                        indent = indent + 2
                     end if
                     j = j + 1
  
negative.code: 
            ! Print filter information
            write(UNIT=UNIT_TALLY, FMT='(1X,2A,1X,A)') space(1:indent), &
                 trim(filter_name(j)), trim(get_label(t, j, bins(j)))
  
            ! Determine scoring index for this bin combination -- note that unlike
            ! in the score_tally subroutine, we have to use max(bins,1) since all
  