If you use the COOLING in your simulation, there will be a field ElectronAbundance, and you can use it
as described. In case you don’t, you need to choose a mean atomic/molecular weight \mu, which depends on assumed chemical composition and ionization state.
