Example for testing COLIBRE stellar winds by massive stars (early kinetic feedback). This consists of a uniform box of gas with a star in the center. The amount of momentum injected can be checked by summing the PartType0/Momentum Received variable in the snapshots over all gas particles. 

Use the python script, read_momentum.py to compare the total amount of momentum received by the gas particles in the simulation against the model input. 

Currently (08/2020) the code does not compile without cooling when including the COLIBRE feedback/chemisty, but it works with the full subgrid model:

./configure --with-hydro=sphenix --with-kernel=quartic-spline --disable-hand-vec --with-subgrid=COLIBRE

in this case, the gas of course starts cooling, which is not desired for this example. 
A hack to avoid this is to TEMPORARILY set the net cooling rate to 0. This can be done by adding "return 0.;" at the very beginning of
routine colibre_cooling_rate in cooling_rates.h.