1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
|
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<!--Converted with LaTeX2HTML 97.1 (release) (July 13th, 1997)
by Nikos Drakos (nikos@cbl.leeds.ac.uk), CBLU, University of Leeds
* revised and updated by: Marcus Hennecke, Ross Moore, Herb Swan
* with significant contributions from:
Jens Lippman, Marek Rouchal, Martin Wilck and others -->
<HTML>
<HEAD>
<TITLE>SLA_PREC - Precession Matrix (FK5)</TITLE>
<META NAME="description" CONTENT="SLA_PREC - Precession Matrix (FK5)">
<META NAME="keywords" CONTENT="sun67">
<META NAME="resource-type" CONTENT="document">
<META NAME="distribution" CONTENT="global">
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso_8859_1">
<LINK REL="STYLESHEET" HREF="sun67.css">
<LINK REL="next" HREF="node153.html">
<LINK REL="previous" HREF="node151.html">
<LINK REL="up" HREF="node13.html">
<LINK REL="next" HREF="node153.html">
</HEAD>
<BODY >
<BR> <HR>
<A NAME="tex2html1946" HREF="node153.html">
<IMG WIDTH="37" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="next" SRC="next_motif.gif"></A>
<A NAME="tex2html1944" HREF="node13.html">
<IMG WIDTH="26" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="up" SRC="up_motif.gif"></A>
<A NAME="tex2html1938" HREF="node151.html">
<IMG WIDTH="63" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="previous" SRC="previous_motif.gif"></A> <A HREF="sun67.html#stardoccontents"><IMG ALIGN="BOTTOM" BORDER="0"
SRC="contents_motif.gif"></A>
<BR>
<B> Next:</B> <A NAME="tex2html1947" HREF="node153.html">SLA_PRECES - Precession</A>
<BR>
<B>Up:</B> <A NAME="tex2html1945" HREF="node13.html">SUBPROGRAM SPECIFICATIONS</A>
<BR>
<B> Previous:</B> <A NAME="tex2html1939" HREF="node151.html">SLA_PREBN - Precession Matrix (FK4)</A>
<BR> <HR> <P>
<P><!--End of Navigation Panel-->
<H2><A NAME="SECTION0004139000000000000000">SLA_PREC - Precession Matrix (FK5)</A>
<A NAME="xref_SLA_PREC"> </A><A NAME="SLA_PREC"> </A>
</H2>
<DL>
<DT><STRONG>ACTION:</STRONG>
<DD>Form the matrix of precession between two epochs (IAU 1976, FK5).
<DT><STRONG>CALL:</STRONG>
<DD><TT>CALL sla_PREC (EP0, EP1, RMATP)</TT>
<P> </DL>
<P> <DL>
<DT><STRONG>GIVEN:</STRONG>
<DD>
<BR>
<TABLE CELLPADDING=3>
<TR VALIGN="TOP"><TD ALIGN="LEFT"><EM>EP0</EM></TD>
<TH ALIGN="LEFT"><B>D</B></TH>
<TD ALIGN="LEFT" NOWRAP>beginning epoch</TD>
</TR>
<TR VALIGN="TOP"><TD ALIGN="LEFT"><EM>EP1</EM></TD>
<TD ALIGN="LEFT"><B>D</B></TD>
<TD ALIGN="LEFT" NOWRAP>ending epoch</TD>
</TR>
</TABLE></DL>
<P> <DL>
<DT><STRONG>RETURNED:</STRONG>
<DD>
<BR>
<TABLE CELLPADDING=3>
<TR VALIGN="TOP"><TD ALIGN="LEFT"><EM>RMATP</EM></TD>
<TH ALIGN="LEFT"><B>D(3,3)</B></TH>
<TD ALIGN="LEFT" NOWRAP>precession matrix</TD>
</TR>
</TABLE></DL>
<P> <DL>
<DT><STRONG>NOTES:</STRONG>
<DD><DL COMPACT>
<DT>1.
<DD>The epochs are TDB Julian epochs.
<DT>2.
<DD>The matrix is in the sense:
<BLOCKQUOTE><B>v</B><SUB>1</SUB> = <B>M</B><IMG WIDTH="7" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
SRC="img72.gif"
ALT="$\cdot$"><B>v</B><SUB>0</SUB></BLOCKQUOTE>
where <B>v</B><SUB>1</SUB> is the star vector relative to the
mean equator and equinox of epoch EP1, <B>M</B> is the
<IMG WIDTH="39" HEIGHT="25" ALIGN="MIDDLE" BORDER="0"
SRC="img18.gif"
ALT="$3\times3$"> matrix RMATP and
<B>v</B><SUB>0</SUB> is the star vector relative to the
mean equator and equinox of epoch EP0.
<DT>3.
<DD>Though the matrix method itself is rigorous, the precession
angles are expressed through canonical polynomials which are
valid only for a limited time span. There are also known
errors in the IAU precession rate. The absolute accuracy
of the present formulation is better than
<IMG WIDTH="23" HEIGHT="18" ALIGN="BOTTOM" BORDER="0"
SRC="img25.gif"
ALT="$0\hspace{-0.05em}^{'\hspace{-0.1em}'}\hspace{-0.4em}.1$"> from
1960AD to 2040AD, better than <IMG WIDTH="17" HEIGHT="17" ALIGN="BOTTOM" BORDER="0"
SRC="img140.gif"
ALT="$1\hspace{-0.05em}^{'\hspace{-0.1em}'}$"> from 1640AD to 2360AD,
and remains below <IMG WIDTH="17" HEIGHT="18" ALIGN="BOTTOM" BORDER="0"
SRC="img75.gif"
ALT="$3\hspace{-0.05em}^{'\hspace{-0.1em}'}$"> for the whole of the period
500BC to 3000AD. The errors exceed <IMG WIDTH="25" HEIGHT="18" ALIGN="BOTTOM" BORDER="0"
SRC="img74.gif"
ALT="$10\hspace{-0.05em}^{'\hspace{-0.1em}'}$"> outside the
range 1200BC to 3900AD, exceed <IMG WIDTH="33" HEIGHT="18" ALIGN="BOTTOM" BORDER="0"
SRC="img141.gif"
ALT="$100\hspace{-0.05em}^{'\hspace{-0.1em}'}$"> outside 4200BC to
5600AD and exceed <IMG WIDTH="41" HEIGHT="18" ALIGN="BOTTOM" BORDER="0"
SRC="img142.gif"
ALT="$1000\hspace{-0.05em}^{'\hspace{-0.1em}'}$"> outside 6800BC to 8200AD.
The SLALIB routine sla_PRECL implements a more elaborate
model which is suitable for problems spanning several
thousand years.
</DL></DL>
<P> <DL>
<DT><STRONG>REFERENCES:</STRONG>
<DD><DL COMPACT>
<DT>1.
<DD>Lieske, J.H., 1979. <I>Astr.Astrophys.</I> <B>73</B>, 282;
equations 6 & 7, p283.
<DT>2.
<DD>Kaplan, G.H., 1981. <I>USNO circular no. 163</I>, pA2.
</DL></DL>
<BR> <HR>
<A NAME="tex2html1946" HREF="node153.html">
<IMG WIDTH="37" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="next" SRC="next_motif.gif"></A>
<A NAME="tex2html1944" HREF="node13.html">
<IMG WIDTH="26" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="up" SRC="up_motif.gif"></A>
<A NAME="tex2html1938" HREF="node151.html">
<IMG WIDTH="63" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="previous" SRC="previous_motif.gif"></A> <A HREF="sun67.html#stardoccontents"><IMG ALIGN="BOTTOM" BORDER="0"
SRC="contents_motif.gif"></A>
<BR>
<B> Next:</B> <A NAME="tex2html1947" HREF="node153.html">SLA_PRECES - Precession</A>
<BR>
<B>Up:</B> <A NAME="tex2html1945" HREF="node13.html">SUBPROGRAM SPECIFICATIONS</A>
<BR>
<B> Previous:</B> <A NAME="tex2html1939" HREF="node151.html">SLA_PREBN - Precession Matrix (FK4)</A>
<BR> <HR> <P>
<P><!--End of Navigation Panel-->
<ADDRESS>
<I>SLALIB --- Positional Astronomy Library<BR>Starlink User Note 67<BR>P. T. Wallace<BR>12 October 1999<BR>E-mail:ptw@star.rl.ac.uk</I>
</ADDRESS>
</BODY>
</HTML>
|