  A Monte Carlo method for the collisional guiding-center Fokker-Planck kinetic
equation is derived to include the effects of background magnetic-field
nonuniformity. It is shown that, in the limit of a homogeneous magnetic field,
the guiding-center Monte Carlo collision operators reduce to the standard
Coulomb operators. The coefficients required for the implementation of the
method are also given.
