  We provide a hands-on introduction to Monte Carlo simulations in
nonperturbative lattice quantum gravity, formulated in terms of Causal
Dynamical Triangulations (CDT). We describe explicitly the implementation of
Monte Carlo moves and the associated detailed-balance equations in two and
three spacetime dimensions. We discuss how to optimize data storage and
retrieval, which are nontrivial due to the dynamical nature of the lattices,
and how to reconstruct the full geometry from selected stored data. Various
aspects of the simulation, including tuning, thermalization and the measurement
of observables are also treated. An associated open-source C++ implementation
code is freely available online.
