This example is generated in ``examples/cooling_rate`` with the Grackle cooling.
The plot are the same than in [Smith2016]_.
The code has two different modes 1D or 2D.
In the first mode, the code generates particles at different temperatures but same density and computes the cooling rate associated to theses conditions.