Use new nH function in recombinationXe
Make zmin,zmax options to recombinationXe, with support for zmax -> Automatic
Push default zmin down to zero, since this seems to work fine
Add option to extend Xe above zmax using XeEq