anchor.diff:  
     integer j
     real(kind=fPrec) average,sigma,zpj,sum,sqsum
  -
  -#ifdef HDF5
  -  ! For tracks2
  -  integer hdfturn,hdfpid,hdftyp
  -  real(kind=fPrec) hdfx,hdfxp,hdfy,hdfyp,hdfdee,hdfs
  -#endif
  -
   
   !++  Output information:
   !++
  
positive.code:       end if
      end if
    end do
  
  !GRD THIS LOOP MUST NOT BE WRITTEN INTO THE "IF(FIRSTRUN)" LOOP !!!!!
    if(dowrite_tracks) then
      call coll_writeTracks2(1)
    end if
  
    do j=1,napx
      if(part_hit_pos(j) == ie .and. part_hit_turn(j) == iturn) then
  
negative.code: !++  End of check for hit this turn and element
      end if
    end do ! end do j = 1, napx
  
    if(dowrite_tracks) then
      call coll_writeTracks2(1)
    end if
  
  !++  Calculate statistical observables and save into files...
    if (n_impact.gt.0) then
  