anchor.diff:            nspx = zero
             nspy = zero
           end if
  -
  -        sampl(ie) = totals
         end do
       end if
     end if
  
positive.code:           if(h5_writeTracks2) then
              hdfpid=ipart(j)
              hdfturn=iturn
              hdfs=sampl(ie)
              hdfx=xv1(j)
              hdfxp=yv1(j)
              hdfy=xv2(j)
  
negative.code:   call f_open(unit=coll_orbitCheckUnit,file=coll_orbitCheckFile,formatted=.true.,mode="w")
    write(coll_orbitCheckUnit,"(a1,1x,a6,3(1x,a15))") "#","s","torbitx","torbity"
    do i=1,iu
      write(coll_orbitCheckUnit,"(i8,3(1x,1pe15.7))") i, dcum(i), torbx(i), torby(i)
    end do
    call f_close(coll_orbitCheckUnit)
  
  