anchor.diff:      else
         scatterProb = ratioTot*elemScale
       end if
  +
  +    ! Some accounting of probabilities
  +    sumProb = sumProb + scatterProb
  +    if(scatterProb > one) nAbove = nAbove + 1
  +
       if(rndVals(k) > scatterProb) then
         cycle
       else
  
positive.code:   iLost    = 0
    isDiff   = .false.
    updateE  = .false.
  
    pScattered(:) = .false.
  
  
negative.code:   integer, intent(in) :: iElem
    integer, intent(in) :: nTurn
  
    integer i, j, k, idElem, idPro, iGen, nGen, idGen, iError, unitDens, iLost, procID, fixedID, nAbove
    logical updateE, autoRatio, isDiff, isExact, densDump, beamDump, pScattered(napx)
    real(kind=fPrec) t, dEE, dPP, theta, phi, pVec(3), pNew, xRot, yRot, elemScale, sigmaTot,         &
      ratioTot, crossSection, scatterProb, targetDensity, scRatio, brRatio, rndVals(napx*3), sumProb
  
    real(kind=fPrec), allocatable :: brThreshold(:)
    logical,          allocatable :: hasProc(:,:)
  