anchor.diff:  
   end subroutine k2coll_scamcs
   
  -subroutine k2coll_iterat(a,b,dh,s)
  +subroutine k2coll_iterat(a, b, dh, s)
   
     real(kind=fPrec), intent(in)    :: a
     real(kind=fPrec), intent(in)    :: b
  
positive.code:   use mathlib_bouncer
    use coll_materials
  
    implicit none
  
    real(kind=fPrec) t,cnorm,cnform
    parameter(cnorm=2.607e-5_fPrec,cnform=0.8561e3_fPrec) ! DM: changed 2.607d-4 to 2.607d-5 to fix Rutherford bug
  
    k2coll_ruth=(cnorm*exp_mb(((-one*real(t,fPrec))*cnform)*emr(mcurr)**2))*(zatom(mcurr)/real(t,fPrec))**2
  end function k2coll_ruth
  
  !>
  
negative.code: !! Select a scattering type (elastic, sd, inelastic, ...)
  !<
  integer function k2coll_ichoix(ma)
  
    use mod_ranlux
    use coll_materials
  
    integer, intent(in) :: ma
    integer i
  
    i = 1
  10 continue
    if(rndm4() > cprob(i,ma)) then
      i = i+1
      goto 10
    end if
  
    k2coll_ichoix = i
  
  end function k2coll_ichoix
  
  !>
  