FUNCTION IPS_velocity, rr_earth, R, F, w3darg,  $
        ut_earth    = ut_earth  , $
        pa_earth    = pa_earth  , $
        elo_earth   = elo_earth , $
        elo_sun     = elo_sun   , $
        degrees     = degrees
    @compile_opt.pro        ; On error, return to caller

RETURN, F*sin( (elo_earth+elo_sun)*ToRadians(degrees=degrees) )  &  END