This is a python implementation of the Rock'n'Roll resuspension model formulated by Reeks and Hall (2001). The implementation uses the parametrization proposed by Biasi (2001). The main features of this implementation are: Discretized distributions for adhesion forces and particle size (all scipy.stats distributions are available for use!) Unsteady friction velocity inputs (both generated in place and read from .npy files.) Several drag force models. Plotting capabilities Although the focus here is on the Rock'n'Roll model, any other kinetic models (models that compute a resuspension rate from a set of parameters) can be integrated into the code quite easily.