anchor.diff:    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, sampl(i), torbx(i), torby(i)
  +    write(coll_orbitCheckUnit,"(i8,3(1x,1pe15.7))") i, dcum(i), torbx(i), torby(i)
     end do
     call f_close(coll_orbitCheckUnit)
   
  
positive.code:           nspx = zero
            nspy = zero
          end if
  
            sampl(ie)    = totals
        end do
      end if
    end if
  
negative.code:   call alloc(sqsum_ax, nblz_new, zero, "sqsum_ax") !(nblz_new)
    call alloc(sum_ay,   nblz_new, zero, "sum_ay") !(nblz_new)
    call alloc(sqsum_ay, nblz_new, zero, "sqsum_ay") !(nblz_new)
  
    call alloc(secondary,            npart_new, 0, "secondary") !(npart_new)
    call alloc(tertiary,             npart_new, 0, "tertiary") !(npart_new)
  