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: !! Select a scattering type (elastic, sd, inelastic, ...)
  !<
  integer function k2coll_ichoix(ma)
    use mod_ranlux
    use coll_materials
    integer ma,i
    real(kind=fPrec) aran
    aran=real(rndm4(),fPrec)
    i=1
  10 if( aran.gt.cprob(i,ma) ) then
      i=i+1
      goto 10
    end if
    k2coll_ichoix=i
  end function k2coll_ichoix
  
  !>
  
negative.code:   use mathlib_bouncer
    use coll_materials
  
    real(kind=fPrec), intent(in) :: t
  
    ! DM: changed 2.607d-4 to 2.607d-5 to fix Rutherford bug
    real(kind=fPrec), parameter :: cnorm=2.607e-5_fPrec
    real(kind=fPrec), parameter :: cnform=0.8561e3_fPrec
  
    k2coll_ruth = (cnorm*exp_mb(((-one*t)*cnform)*emr(mcurr)**2)) * (zatom(mcurr)/t)**2
  
  end function k2coll_ruth
  
  !>
  