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'
SITENAME= 'Mt. Hopkins 60-inch'
SITELAT = '31:40:51.4'
The results are summarized to standard output as follows:
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
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.