anchor.diff:      ! correct TOT-CSec for energy dependence of qel
       ! TOT CS is here without a Coulomb contribution
       csect(0,ma) = csref(0,ma) + freep(ma) * (pptot - pptref)
  -    bn(ma) = (bnref(ma) * csect(0,ma)) / csref(0,ma)
  +    bn(ma)      = (bnref(ma) * csect(0,ma)) / csref(0,ma)
   
       ! also correct inel-CS
       csect(1,ma) = (csref(1,ma) * csect(0,ma)) / csref(0,ma)
  
positive.code:   ! Do a step for a point-like interaction.
    ! Get monte-carlo interaction length.
  10 continue
    zlm1     = (-one*xintl(mat))*log_mb(real(rndm4(),fPrec))
    nabs_tmp = 0  ! type of interaction reset before following scattering process
    xpBef    = xp ! save angles and momentum before scattering
    zpBef    = zp
  
negative.code:         x  = x  + tiltangle*c_length
          xp = xp + tiltangle
        else if(tiltangle < zero) then
          x  = x  + tiltangle*c_length
          xp = xp + tiltangle
          x  = x  - sin_mb(tiltangle) * c_length
        end if
  
        ! Transform back to particle coordinates with opening and offset
  