Library for sampling a poisson distribution in Fortran

I’m looking for a reliable library/module to sample a poisson distribution in Fortran. Thanks for any suggestions!

The interface for the function might be

function sample_poisson(mu) result(k)
  real(dp), intent(in) :: mu
  integer :: k
  ! ...
end function

https://github.com/Jim-215-Fisher/statislib

Alan Miller wrote

random.f90 A module for generating random numbers from a range of distributions [including the Poisson].

ignpoi.f90 Generate random Poisson deviates.