slalink - functions available with the slalink extensions to tawk.

SYNOPSYS

preces(sys, ep0, ep1, ra, dec)
convert from equinox0 to equinox1.
fk45(r1950, d1950, bepoch)
convert from FK4 to FK5.
fk54(r2000, d2000, bepoch)
convert from FK5 to FK4.
r2d(radians)
convert radians to degrees.
d2r(degrees)
convert degrees to radians.
r2h(radians)
convert radians to hours.
h2r(hours)
convert hours to radians.

DESCRIPTION

preces(sys, ep0, ep1, ra, dec)
convert from equinox0 to equinox1.

Precession - either FK4 (Bessel-Newcomb, pre-IAU1976) or FK5 (Fricke, post-IAU1976) as required.

sys
precession to be applied: "FK4" or "FK5"
ep0,ep1
starting and ending epoch
ra,dec
RA,Dec, mean equator & equinox of epoch ep0
RETURNS
[ ra,dec ]
RA,Dec, mean equator & equinox of epoch ep1
Notes
1)
The epochs are Besselian if sys='FK4' and Julian if 'FK5'. For example, to precess coordinates in the old system from equinox 1900.0 to 1950.0 the call would be.
	[ ra, dec ] = preces ("FK4", 1900.0, 1950.0, ra, dec)
2)
This routine will not correctly convert between the old and the new systems - for example conversion from B1950 to J2000. For these purposes see fk425, fk524, fk45z and fk54.
3)
If an invalid sys is supplied, values of -99.0,-99.0 will be returned for both ra and dc.
P.T.Wallace
Starlink 22 December 1993
fk45(r1950, d1950, bepoch)
convert from FK4 to FK5.

Convert B1950.0 FK4 star data to J2000.0 FK5 assuming zero proper motion in an inertial frame (double precision)

This routine converts stars from the old, Bessel-Newcomb, FK4 system to the new, IAU 1976, FK5, Fricke system, in such a way that the FK5 proper motion is zero. Because such a star has, in general, a non-zero proper motion in the FK4 system, the routine requires the epoch at which the position in the FK4 system was determined.

r1950,d1950
B1950.0 FK4 RA,Dec at epoch (rad)
bepoch
Besselian epoch (e.g. 1979.3)
RETURNS
[ r2000, d2000 ]
J2000.0 FK5 RA,Dec (rad)
Notes
1)
The epoch BEPOCH is strictly speaking Besselian, but if a Julian epoch is supplied the result will be affected only to a negligible extent.
2)
Conversion from Besselian epoch 1950.0 to Julian epoch 2000.0 only is provided for. Conversions involving other epochs will require use of the appropriate precession, proper motion, and E-terms routines before and/or after FK425 is called.
3)
In the FK4 catalogue the proper motions of stars within 10 degrees of the poles do not embody the differential E-term effect and should, strictly speaking, be handled in a different manner from stars outside these regions. However, given the general lack of homogeneity of the star data available for routine astrometry, the difficulties of handling positions that may have been determined from astrometric fields spanning the polar and non-polar regions, the likelihood that the differential E-terms effect was not taken into account when allowing for proper motion in past astrometry, and the undesirability of a discontinuity in the algorithm, the decision has been made in this routine to include the effect of differential E-terms on the proper motions for all stars, whether polar or not. At epoch 2000, and measuring on the sky rather than in terms of dRA, the errors resulting from this simplification are less than 1 milliarcsecond in position and 1 milliarcsecond per century in proper motion.
Example
	[ r2000, d2000 ] = fk45(r1950, d1950, 1950)
fk54(r2000, d2000, bepoch)
convert from FK5 to FK4.

Convert a J2000.0 FK5 star position to B1950.0 FK4 assuming zero proper motion and parallax.

This routine converts star positions from the new, IAU 1976, FK5, Fricke system to the old, Bessel-Newcomb, FK4 system.

r2000,d2000
J2000.0 FK5 RA,Dec (rad)
bepoch
Besselian epoch (e.g. 1950)
RETURNS
[ r1950, d1950 ]
B1950.0 FK4 RA,Dec (rad) at epoch BEPOCH
Notes
1)
The proper motion in RA is dRA/dt rather than cos(Dec)*dRA/dt.
2)
Conversion from Julian epoch 2000.0 to Besselian epoch 1950.0 only is provided for. Conversions involving other epochs will require use of the appropriate precession routines before and after this routine is called.
3)
Unlike in the slaFK524 routine, the FK5 proper motions, the parallax and the radial velocity are presumed zero.
4)
It is the intention that FK5 should be a close approximation to an inertial frame, so that distant objects have zero proper motion; such objects have (in general) non-zero proper motion in FK4, and this routine returns those fictitious proper motions.
5)
The position returned by this routine is in the B1950 reference frame but at Besselian epoch bepoch. For comparison with catalogues the bepoch argument will frequently be 1950.0.
P.T.Wallace
Starlink 30 October 1993
Example
	[ r1950, d1950 ] = fk54(r2000, d2000, 1950)

SEE ALSO

slalink.1
programs available with the slalink extensions to tawk.
slalib.3
The Star Link Astrometry library.
tabtable.1
starbase data tables introduction (this page).
tablefun.3
starbase data tables introduction (this page).
tawk
awk language extensions implimented at SAO.
mawk
an implimentation of the awk programming language by Mike Brennan (brennanboeing.com).