Dense suspensions are widespread in nature, manufacturing and process engineering. Particle-based simulations have proven to be an invaluable complement to experimental rheological characterisation, serving as a virtual rheometer that enables rapid exploration of parameter space and detailed scrutiny of microscopic dynamics. To maximise the utility of such simulations, it can be advantageous to exploit pre-existing, well-optimised, well-documented codes. Here we provide a simple description of how to use LAMMPS to study the rheology of dense, granular suspensions.