./toolbox
chdoppler.pro

Last modification date:

Thu Jun 2 23:56:08 2005

chdoppler

function chdoppler(ra, dec, julday, [obspos=obspos], [light=light])

computes the projected velocity of the telescope wrt six coordinate systems: geo, helio, bary, lsrk, lsrd, gal negative velocities mean approach The standard LSR is defined as follows: the sun moves at 20.0 km/s toward ra=18.0h, dec=30.0 deg in 1900 epoch coords

CHDOPPLER IS NOT VECTORIZED; ONLY SINGLE NUMBERS

This code came via e-mail from Carl Heiles via Tom Bania on 11/04/2004. Local changes:

• modify this documentation for use by idldoc.
• removed path argument and related code, replaced by obspos argument.
• default position is the GBT.
• Observatory longitude was not being passed to juldaytolmst.
• LSRD added
• Galactocentric added
• Checked against aips++ Measures. Differences were < 20 m/s in one test case (<10m/s for geo, bary, and lsrk).
• Double precision throughout.
• Relativistic addition of velocities.

Returns

The velocity in km/s, or as a faction of c if the keyword /light is specified. the result is a 6-element vector whose elements are [geo, helio, bary, lsrk, lsrd, gal].

Uses

baryvel precess juldaytolmst

Version

$Id: chdoppler.pro,v 1.3 2005/04/08 20:36:06 bgarwood Exp $

Parameters
ra in, required	The source ra in decimal hours, equinox 2000
dec in, required	The source dec in decimal hours, equinox 2000
julday in, required	The julian day

Keywords
obspos in, optional	observatory position [East longitude, latitude] in degrees. Uses the GBT position if not specified.
light in, optional	When set, returns the velocity as a fraction of c