anchor.diff:          nspx = zero
           nspy = zero
         end if
  -      sampl(ie) = totals
       end do
   
     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:   real(kind=fPrec), allocatable, save :: sqsum_ax(:) !(nblz)
    real(kind=fPrec), allocatable, save :: sum_ay(:) !(nblz)
    real(kind=fPrec), allocatable, save :: sqsum_ay(:) !(nblz)
  
    real(kind=fPrec), allocatable, save :: neffx(:) !(numeff)
    real(kind=fPrec), allocatable, save :: neffy(:) !(numeff)
  