3 char* startab = "/lib/sky/estartab";
8 double lomoon, himoon, sd;
11 static char saoa[100];
14 lomoon = omoon.point[0].ra - sd;
17 himoon = omoon.point[NPTS+1].ra + sd;
26 f = open(startab, OREAD);
28 fprint(2, "%s?\n", startab);
32 epoch = (epoch-1900.0) * 365.24220 + 0.313;
36 * read mean places of stars at epoch of star table
48 alpha = rah + ram/60 + ras/3600;
50 if(alpha < lomoon || alpha > himoon)
53 if(alpha < lomoon && alpha > himoon)
57 sprint(saop, "%ld", sao);
67 * convert rt ascension and declination to internal format
70 delta = fabs(dday) + dmin/60 + dsec/3600;
76 * if(fabs(beta) > 6.55*radian)
79 sd = .0896833e0*cos(beta)*sin(lambda-1.3820+.00092422117*eday)
84 for(i=0; i<=NPTS+1; i++)
85 setobj(&ostar.point[i]);
87 occult(&omoon, &ostar, 0);
88 if(occ.t1 >= 0 || occ.t5 >= 0) {
94 if(occ.t1 >= 0 && occ.e1 >= 0)
95 event("Occultation of SAO %s begins at ",
97 if(occ.t5 >= 0 && occ.e5 >= 0)
98 event("Occultation of SAO %s ends at ",