[AIT logo]

Institut für Astronomie und Astrophysik

Abteilung Astronomie

Sand 1, D-72076 Tübingen, Germany
[Uni logo]


binradvel Source code in binradvel.pro

binradvel

Name
         binradvel
Purpose
         compute radial velocity of a binary system as a function of
         orbital phase or time (if t0 or t90 are given)
Category
         astronomy
Calling Sequence
         vel=binradvel(phase,k=k,period=period,asini=asini,eccen=eccen, $
                   omega=omega,gamma=gamma,degrees=degrees)
Input Parameters
          phase: array with the orbital phase points (from 0 to 1)
                 for which the radial velocity is to be computed;
                 if t0 or t90 are given: time in days (same time
                 system as t0 or t90)
          k    : velocity amplitude of the system
          asini: value of a sin i in kilometers
          period: orbital period, in days
          NOTE: either k or  both, asini and period, must be given
Optional Input Parameters
          t0   : epoch of periastron (see above, description of "phase" argument)
          t90  : epoch when mean longitude = 90 degrees
          omega: longitude of periastron (little omega; radians,
                 default: 0.)
          eccen: eccentricity (default: 0.)
          gamma: systemic velocity (same units as k or in km/s; default:0.)
Keyword Parameters
          degrees: if set, angular arguments are in degrees instead
                 of radian
Output Parameters
      the function returns the radial velocity as a function of the
      phase in the units of k or of  asini/period
Restrictions
      * the eccentricity must be 0<=e<1
      * only one of the arguments t0 and t90 is allowed
Procedure
      See chapter 3 of
         R.W. Hilditch, An Introduction to Close Binary Stars,
         Cambridge Univ. Press, 2001
Example
 Radial Velocity of HR7000
 (Hilditch, Fig. 3.27; Griffin et al., 1997, Fig. 2)
 npt=100
 phase=findgen(npt)/(npt-1)
 vel=binradvel(phase,k=21.68,gamma=-22.52,eccen=0.372,omega=306.5,/degrees)
 ; Show two orbits for clarity
 pp=[phase,1.+phase]
 vv=[vel,vel]
 plot,pp,vv,xtitle='Orbital Phase',ytitle='Radial Velocity [km/s]'
 xrange=!x.crange
 oplot,xrange,[0.,0.],linestyle=2
Revision History
 $Log: binradvel.pro,v $
 Revision 1.2  2005/02/18 21:41:23  wilms
 * added t0 and t90 arguments
 Revision 1.1  2002/09/09 14:52:04  wilms
 initial release

Last modified by pro2html on 2005 February 19 at 04:10 UTC

[Home Page] [Software, Documentation] [IDL Documentation] [Quick Reference] [Feedback]

Jörn Wilms (wilms@astro.uni-tuebingen.de)
Updated automatically