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