Work around API change in Astropy 5
In Astropy 4, z_at_value
returns a float, but in Astropy 5,
returns a quantity with dimensionless redshift units. Work around
this by making sure that its return value is a Quantity before we
convert it to a float.
Note: remove this workaround once we drop support for astropy < 5.