I can´t find anything other than the standard expression valid BEFORE recombination which does not depend on the ionization fraction, or expressions for AFTER the precipitous fall in the sound speed after last scattering. Has anyone modelled that "precipitous" fall (in either sound speed or pressure) in detail as a function of the (small and falling) ionisation fraction?