anchor.diff:            sqsum_ay(ie) = sqsum_ay(ie) + nspy**2
             nampl(ie)    = nampl(ie) + 1
             sampl(ie)    = totals
  -          ename(ie)    = bez(myix)
           else
             nspx = zero
             nspy = zero
  
positive.code:   use hdf5_tracks2
  #endif
  
    integer i
  
    ! Just call it here since samples are no longer supported
    call collimate_end_sample(1)
  
negative.code:     end do
    end if
  
    if(iturn == 1 .and. firstrun) then
      write(outlun,*) ' '
      write(outlun,*)   'Collimator information: '
  