anchor.diff:    call recuin(scatter_seed1,scatter_seed2)
   
     call alloc(rndVals,   napx*3,             zero,   "rndVals")
  -  call alloc(pLost,     napx,               .false.,"pLost")
     call alloc(pScattered,napx,               .false.,"pScattered")
     call alloc(hasProc,   nGen,scatter_nProc, .false.,"hasProc")
     call alloc(nLost,     nGen,scatter_nProc, 0,      "nLost")
  
positive.code:     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
      else
        ! Calculate new energy, and scattering angle assuming energy >> mass
        if(abs(dPP) >= pieni) then
  
negative.code: 
    real(kind=fPrec), allocatable :: brThreshold(:)
    real(kind=fPrec), allocatable :: rndVals(:)
    logical,          allocatable :: pScattered(:)
    logical,          allocatable :: hasProc(:,:)
    integer,          allocatable :: nLost(:,:)
  