;+
; Convert velocity(m/s) to frequency(MHz) using the given rest
; frequency and velocity definition.
;
; @param vel {in}{required} Velocity in m/s
; @param restfreq {in}{required} Rest frequency in MHz
; @keyword veldef {in}{optional}{type=string} The velocity definition
; which must be one of OPTICAL, RADIO, or TRUE.  Defaults to RADIO.
;
; @returns frequency in MHz
;
; @version $Id: veltofreq.pro,v 1.4 2005/04/08 20:36:06 bgarwood Exp $
;-
function veltofreq, vel, restfreq, veldef=veldef
    compile_opt idl2

    if (not keyword_set(veldef)) then veldef = "RADIO"

    result = double(vel) / !gc.light_speed
    case veldef of
        'RADIO': result = (1.0d - result) * restfreq
        'OPTICAL': result = restfreq / (result + 1.0d)
        'TRUE': result = restfreq * sqrt((1.d - result)/(1.d + result))
        else: message, 'unrecognized velocity definition'
    endcase

    return, result
end