use astropy cosmology, use quantities, but allow dimensionless argument with default units