bcvcorr Example 1

Compute Barycentric Velocity Correction Using Header Parameters

bcvcorr will compute the barycentric velocity correction and heliocentric Julian date for an observation based on position and observatory keywords in the image header.

bcvcorr Task Parameters

The names of the keywords to be used are set in the parameter file:

  (keyra = "RA")            Right ascension header keyword
  (keydec = "DEC")          Declination header keyword
  (keyeqnx = "EPOCH")       Coordinate equinox header keyword
...
  (keydate = "DATE-OBS")    Header keyword for date of observation (dd/mm/yy)
  (keystart = "UTOPEN")     Header keyword for UT at start of observation
  (keymid = "UTMID")        Header keyword for UT in middle of observation
  (keyend = "UT")           Header keyword for UT at end of observation
  (keyexp = "EXPOSURE")     Header keyword for duration of observation in sec
...
  (obsname = "file")        Observatory name--file to read from image heade
  (keyobs = "SITENAME")     Observatory name header keyword
  (keylat = "SITELAT ")     Observatory latitude header keyword
  (keylong = "SITELONG")    Observatory longitude header keyword
  (keyalt = "SITEELEV")     Observatory altitude header keyword

Spectrum Image Header

If the parameter obsname is file, the position of the observatory is read from the specified header keywords. Here is the relevant portion of the image header:

...
    RA      = '09:36:48.0'
    HA      = '00:53:18.9'
    DEC     = '16:59:00.0'
    EPOCH   =                1950.
    UT      = '04:32:26'
    ST      = '10:32:24.6'
    DATE-OBS= '13/04/91'
...
    EXPOSURE=                3600.
...
    SITENAME= 'Mt. Hopkins 60-inch'
    SITELONG= '110:52:39.0'
    SITELAT = '31:40:51.4'
    SITEELEV=                2320.
...

Results

The results are summarized to standard output as follows:

ar> bcvcorr
List of input spectra (22408): 
22408: 4096 x 1 x 1 1-D image
UT: 4:32:26.0, UTOPEN: INDEF, UTMID: INDEF, EXP: 3600
4/13/1991 4:02:26.0 UT
RA: 9:36:48.0, Dec: 16:59:00.0  1950.0
Mt. Hopkins 60-inch lat 31:40:51.4 , long 110:52:39.0, alt 2320.0
Julian date is 2448359.68919 at 4/13/1991 4:02:26.0 UT
Object at ra 9:36:48.000 dec 16:59:00.00 eq 1950
Heliocentric Julian date: 2448359.69197
gbcvel = -25.7858  ghcvel = -25.7864  geovel = -0.0875
bcv = -25.8732  hcv = -25.8739 computed
bcv = INDEF  hcv = -25.8721 from file
ar>

As only the exposure time and the UT of the end of the observation are given, the midtime of the event is computed by subtracting half of the exposure time from the end time. This time is given on the Julian Date line. The heliocentric Julian date is the time which the measured light from the observed object passes the sun. A correction from the center of the earth to the observatory site (geovel) is added to the geocentric barycentric (gbcvel) and heliocentric (ghcvel) velocity corrections to give the topocentric barycentric (bcv) and heliocentric (hcv corrections. The computed corrections can be compared to those already in the header (keywords BCV and HCV) which are given on the last line.

Last updated October 30, 2003 by Doug Mink dmink@cfa.harvard.edu

[rvsao] [bcvcorr]