I think this was written for an older version?  Right now, ``wavelength = wavelength * (1. + heliocentric/c)`` *does* work.  So this discussion need to be updated to reflect that.  @hamogu, do you want to re-write this, or should @migueldvb do it?