1 #include "astro.h"
2
3 char* startab;
4
5 void
stars(void)6 stars(void)
7 {
8 double lomoon, himoon, sd;
9 int wrap, f, i;
10 char *saop;
11 static char saoa[100];
12
13 sd = 1000*radsec;
14 lomoon = omoon.point[0].ra - sd;
15 if(lomoon < 0)
16 lomoon += pipi;
17 himoon = omoon.point[NPTS+1].ra + sd;
18 if(himoon > pipi)
19 himoon -= pipi;
20 lomoon *= 12/pi;
21 himoon *= 12/pi;
22 wrap = 0;
23 if(lomoon > himoon)
24 wrap++;
25
26 f = open(startab, OREAD);
27 if(f < 0) {
28 fprint(2, "%s?\n", startab);
29 return;
30 }
31 epoch = 1950.0;
32 epoch = (epoch-1900.0) * 365.24220 + 0.313;
33 saop = saoa;
34
35 /*
36 * read mean places of stars at epoch of star table
37 */
38
39 loop:
40 if(rline(f)) {
41 close(f);
42 return;
43 }
44 rah = atof(line+17);
45 ram = atof(line+20);
46 ras = atof(line+23);
47
48 alpha = rah + ram/60 + ras/3600;
49 if(wrap == 0) {
50 if(alpha < lomoon || alpha > himoon)
51 goto loop;
52 } else
53 if(alpha < lomoon && alpha > himoon)
54 goto loop;
55
56 sao = atof(line+0);
57 sprint(saop, "%ld", sao);
58 da = atof(line+30);
59 dday = atof(line+37);
60 dmin = atof(line+41);
61 dsec = atof(line+44);
62 dd = atof(line+50);
63 px = atof(line+57);
64 mag = atof(line+61);
65
66 /*
67 * convert rt ascension and declination to internal format
68 */
69
70 delta = fabs(dday) + dmin/60 + dsec/3600;
71 if(dday < 0)
72 delta = -delta;
73
74 star();
75 /*
76 * if(fabs(beta) > 6.55*radian)
77 * goto loop;
78 */
79 sd = .0896833e0*cos(beta)*sin(lambda-1.3820+.00092422117*eday)
80 + 0.99597*sin(beta);
81 if(fabs(sd) > .0183)
82 goto loop;
83
84 for(i=0; i<=NPTS+1; i++)
85 setobj(&ostar.point[i]);
86
87 occult(&omoon, &ostar, 0);
88 if(occ.t1 >= 0 || occ.t5 >= 0) {
89 i = PTIME;
90 if(mag > 2)
91 i |= DARK;
92 if(mag < 5)
93 i |= SIGNIF;
94 if(occ.t1 >= 0 && occ.e1 >= 0)
95 event("Occultation of SAO %s begins at ",
96 saop, "", occ.t1, i);
97 if(occ.t5 >= 0 && occ.e5 >= 0)
98 event("Occultation of SAO %s ends at ",
99 saop, "", occ.t5, i);
100 while(*saop++)
101 ;
102 }
103 goto loop;
104 }
105