A model for the large-scale-circulation (LSC) of turbulent Rayleigh-Bénard convection in cylindrical samples is presented. It consists of two physically motivated stochastic ordinary differential equations, one each for the strength and the azimuthal orientation of the LSC. Stochastic forces represent phenomenologically the influence of turbulent fluctuations. Consistent with measurements, the model yields an azimuthally meandering LSC with occasional rotations, and with more rare cessations. As in experiment, cessations have a uniform distribution of LSC orientation changes.