anchor.diff:    do j=1,iu
   ! this transformation gives the right marker/name to the corresponding
   ! beta-dunctions or vice versa ;)
  -    if(ic(j).le.nblo) then
  -      do jb=1,mel(ic(j))
  -        myix=mtyp(ic(j),jb)
  -      end do
  +    if(ic(j) <= nblo) then
  +      myix = mtyp(ic(j),mel(ic(j)))
       else
  -      myix=ic(j)-nblo
  +      myix = ic(j)-nblo
       end if
   
       if(cdb_elemMap(myix) > 0) then
  
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:     write(coll_ampUnit,"(a1,1x,a6,1x,a16,19(1x,a20))") "#","ielem",chr_rPad("name",20),"s","AX_AV","AX_RMS","AY_AV","AY_RMS",&
        "alphax","alphay","betax","betay","orbitx","orbity","dispx","dispy","xbob","ybob","xpbob","ypbob","mux","muy"
      do i=1,iu
        if(ic(i) <= nblo) then
          ix = mtyp(ic(i),mel(ic(i)))
        else
          ix = ic(i)-nblo
        end if
        write(coll_ampUnit,"(i8,1x,a16,19(1x,1pe20.13))") i, bez(ix)(1:20), sampl(i),                      &
          sum_ax(i)/real(max(nampl(i),1),fPrec),                                                           &
          sqrt(abs((sqsum_ax(i)/real(max(nampl(i),1),fPrec))-(sum_ax(i)/real(max(nampl(i),1),fPrec))**2)), &
          sum_ay(i)/real(max(nampl(i),1),fPrec),                                                           &
  