Extra features: allow constant alphaS; add guards for massless quarks; add speed of sound computation (to be double checked)
