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:   end subroutine write_tallies
  
  !===============================================================================
  ! GET_UID returns the uid for a cell/surface/etc given a tally, filter type, and
  ! corresponding bin
  !===============================================================================
  
    function get_uid(t, filter_type, bin) result(uid)
  
      type(TallyObject), pointer :: t           ! tally object
      integer, intent(in)        :: filter_type ! type of filter
      integer, intent(in)        :: bin         ! bin in filter array
      character(30)              :: uid         ! user-specified identifier
  
      integer :: index ! index in cells/surfaces/etc array
      real(8) :: E0    ! lower bound for energy bin
      real(8) :: E1    ! upper bound for energy bin
  
      select case(filter_type)
      case (T_UNIVERSE)
         index = t % universe_bins(bin) % scalar
         uid = int_to_str(universes(index) % uid)
      case (T_MATERIAL)
         index = t % material_bins(bin) % scalar
         uid = int_to_str(materials(index) % uid)
      case (T_CELL)
         index = t % cell_bins(bin) % scalar
         uid = int_to_str(cells(index) % uid)
      case (T_CELLBORN)
         index = t % cellborn_bins(bin) % scalar
         uid = int_to_str(cells(index) % uid)
      case (T_SURFACE)
         index = t % surface_bins(bin) % scalar
         uid = int_to_str(surfaces(index) % uid)
      case (T_MESH)
         uid = "Bin " // int_to_str(bin)
      case (T_ENERGYIN)
         E0 = t % energy_in(bin)
         E1 = t % energy_in(bin + 1)
         uid = "[" // trim(real_to_str(E0)) // ", " // trim(real_to_str(E1)) // ")"
      case (T_ENERGYOUT)
         E0 = t % energy_out(bin)
         E1 = t % energy_out(bin + 1)
         uid = "[" // trim(real_to_str(E0)) // ", " // trim(real_to_str(E1)) // ")"
      end select
  
    end function get_uid
  
  !===============================================================================
  ! TALLY_STATISTICS computes the mean and standard deviation of the mean of each
  
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
  