anchor.diff:    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 :: sampl(:) !(nblz)
   
     real(kind=fPrec), allocatable, save :: neffx(:) !(numeff)
     real(kind=fPrec), allocatable, save :: neffy(:) !(numeff)
  
positive.code:                 (ejv(j)-myenom)/myenom,secondary(j)+tertiary(j)+other(j)+scatterhit(j)
  
                write(coll_tracksUnit,'(1x,i8,1x,i4,1x,f10.2,4(1x,e12.5),1x,e11.3,1x,i4)') &
                  ipart(j),iturn,sampl(ie)+half*c_length,           &
                  xv1(j)+half*c_length*yv1(j),yv1(j),                             &
                  xv2(j)+half*c_length*yv2(j),yv2(j),(ejv(j)-myenom)/myenom,      &
                  secondary(j)+tertiary(j)+other(j)+scatterhit(j)
  
negative.code:               hdftyp  = secondary(j)+tertiary(j)+other(j)+scatterhit(j)
                call h5tr2_writeLine(hdfpid,hdfturn,hdfs,hdfx,hdfxp,hdfy,hdfyp,hdfdee,hdftyp)
  
                hdfs  = dcum(ie)+half*c_length
                hdfx  = xv1(j) + half*c_length*yv1(j)
                hdfxp = yv1(j)
                hdfy  = xv2(j) + half*c_length*yv2(j)
  