Add the cddf from Altay 2012
Correct the units on the cold gas fraction (internal energy was 1e-10 ergs/g) and generally refactor a bit. Now cold_gas returns the density in PHYSICAL atoms /cm^3.
Add script for plotting the gas fraction with density