  A procedure for generating random variates from a relativistic Maxwellian
distribution with arbitrary temperature and drift velocity is presented. The
algorithm is based on the rejection method and can be used to initialize
particle velocities in kinetic simulations of plasmas and gases.
