Add optional parameter Omega_cdm to createCosmology that can be used
instead of Omega_b*h^2.