anchor.diff:      if(scatter_allowLosses .and. iLost == 1) then
         ! If lost, flag it, count it, but no need to update energy and angle
         nLost(iGen,procID) = nLost(iGen,procID) + 1
  -      pLost(j) = .true.
  -      phi      = zero
  +      llostp(j) = .true.
  +      phi       = zero
       else
         ! Calculate new energy, and scattering angle assuming energy >> mass
         if(abs(dPP) >= pieni) then
  
positive.code: 
    real(kind=fPrec), allocatable :: brThreshold(:)
    real(kind=fPrec), allocatable :: rndVals(:)
    logical,          allocatable :: pLost(:)
    logical,          allocatable :: pScattered(:)
    logical,          allocatable :: hasProc(:,:)
    integer,          allocatable :: nLost(:,:)
  
negative.code: 
    call dealloc(rndVals,   "rndVals")
    call dealloc(hasProc,   "hasProc")
    call dealloc(pScattered,"pScattered")
    call dealloc(nLost,     "nLost")
    call dealloc(nScattered,"nScattered")
  