Cache ageOfUniverse for each cosmology
Add rsdrag = rs[zdrag] with caching