1 cd5bae78 2004-04-21 devnull #include "astro.h"
3 cd5bae78 2004-04-21 devnull char* solstr[] =
5 cd5bae78 2004-04-21 devnull "Fall equinox",
6 cd5bae78 2004-04-21 devnull "Winter solstice",
7 cd5bae78 2004-04-21 devnull "Spring equinox",
8 cbeb0b26 2006-04-01 devnull "Summer solstice"
13 cd5bae78 2004-04-21 devnull double beta;
16 cd5bae78 2004-04-21 devnull char *betstr;
17 cd5bae78 2004-04-21 devnull } bettab[] =
19 cd5bae78 2004-04-21 devnull -1.3572, 231, 50, "Quadrantid",
20 cd5bae78 2004-04-21 devnull 0.7620, 336, 0, "Eta aquarid",
21 cd5bae78 2004-04-21 devnull 1.5497, 260, -20, "Ophiuchid",
22 cd5bae78 2004-04-21 devnull 2.1324, 315, -15, "Capricornid",
23 cd5bae78 2004-04-21 devnull 2.1991, 339, -17, "Delta aquarid",
24 cd5bae78 2004-04-21 devnull 2.2158, 340, -30, "Pisces australid",
25 cd5bae78 2004-04-21 devnull 2.4331, 46, 58, "Perseid",
26 cd5bae78 2004-04-21 devnull -2.6578, 95, 15, "Orionid",
27 cd5bae78 2004-04-21 devnull -1.8678, 15, -55, "Phoenicid",
28 cd5bae78 2004-04-21 devnull -1.7260, 113, 32, "Geminid",
33 cd5bae78 2004-04-21 devnull search(void)
35 cd5bae78 2004-04-21 devnull Obj2 *p, *q;
36 cd5bae78 2004-04-21 devnull int i, j;
37 cd5bae78 2004-04-21 devnull double t;
39 cd5bae78 2004-04-21 devnull for(i=0; objlst[i]; i++) {
40 cd5bae78 2004-04-21 devnull p = objlst[i];
41 cd5bae78 2004-04-21 devnull if(p == &oshad)
42 cd5bae78 2004-04-21 devnull continue;
43 cd5bae78 2004-04-21 devnull t = rise(p, -.833);
44 cd5bae78 2004-04-21 devnull if(t >= 0.)
45 cd5bae78 2004-04-21 devnull event("%s rises at ", p->name, "", t,
46 cd5bae78 2004-04-21 devnull i==0? PTIME: PTIME|DARK);
47 cd5bae78 2004-04-21 devnull t = set(p, -.833);
48 cd5bae78 2004-04-21 devnull if(t >= 0.)
49 cd5bae78 2004-04-21 devnull event("%s sets at ", p->name, "", t,
50 cd5bae78 2004-04-21 devnull i==0? PTIME: PTIME|DARK);
51 cd5bae78 2004-04-21 devnull if(p == &osun) {
52 cd5bae78 2004-04-21 devnull for(j=0; j<4; j++) {
53 cd5bae78 2004-04-21 devnull t = solstice(j);
54 cd5bae78 2004-04-21 devnull if(t >= 0)
55 cd5bae78 2004-04-21 devnull event("%s at ", solstr[j], "", t,
56 cd5bae78 2004-04-21 devnull SIGNIF|PTIME);
58 cd5bae78 2004-04-21 devnull for(j=0; bettab[j].beta!=0; j++) {
59 cd5bae78 2004-04-21 devnull t = betcross(bettab[j].beta);
60 cd5bae78 2004-04-21 devnull if(t >= 0)
61 cd5bae78 2004-04-21 devnull event("%s meeteeor shouwer",
62 cd5bae78 2004-04-21 devnull bettab[j].betstr, "", t, SIGNIF);
64 cd5bae78 2004-04-21 devnull t = rise(p, -18);
65 cd5bae78 2004-04-21 devnull if(t >= 0)
66 cd5bae78 2004-04-21 devnull event("Twilight starts at ", "", "", t, PTIME);
67 cd5bae78 2004-04-21 devnull t = set(p, -18);
68 cd5bae78 2004-04-21 devnull if(t >= 0)
69 cd5bae78 2004-04-21 devnull event("Twilight ends at ", "", "", t, PTIME);
71 cd5bae78 2004-04-21 devnull if(p == &omoon)
72 cd5bae78 2004-04-21 devnull for(j=0; j<NPTS; j++) {
73 cd5bae78 2004-04-21 devnull if(p->point[j].mag > .75 && p->point[j+1].mag < .25)
74 cd5bae78 2004-04-21 devnull event("New moon", "", "", 0, 0);
75 cd5bae78 2004-04-21 devnull if(p->point[j].mag <= .25 && p->point[j+1].mag > .25)
76 cd5bae78 2004-04-21 devnull event("First quarter moon", "", "", 0, 0);
77 cd5bae78 2004-04-21 devnull if(p->point[j].mag <= .50 && p->point[j+1].mag > .50)
78 cd5bae78 2004-04-21 devnull event("Full moon", "", "", 0, 0);
79 cd5bae78 2004-04-21 devnull if(p->point[j].mag <= .75 && p->point[j+1].mag > .75)
80 cd5bae78 2004-04-21 devnull event("Last quarter moon", "", "", 0, 0);
82 cd5bae78 2004-04-21 devnull if(p == &omerc || p == &ovenus) {
83 cd5bae78 2004-04-21 devnull t = melong(p);
84 cd5bae78 2004-04-21 devnull if(t >= 0) {
85 cd5bae78 2004-04-21 devnull t = rise(p, 0) - rise(&osun, 0);
86 cd5bae78 2004-04-21 devnull if(t < 0)
87 cd5bae78 2004-04-21 devnull t += NPTS;
88 cd5bae78 2004-04-21 devnull if(t > NPTS)
89 cd5bae78 2004-04-21 devnull t -= NPTS;
90 cd5bae78 2004-04-21 devnull if(t > NPTS/2)
91 cd5bae78 2004-04-21 devnull event("Morning elongation of %s", p->name,
92 cd5bae78 2004-04-21 devnull "", 0, SIGNIF);
94 cd5bae78 2004-04-21 devnull event("Evening elongation of %s", p->name,
95 cd5bae78 2004-04-21 devnull "", 0, SIGNIF);
98 cd5bae78 2004-04-21 devnull for(j=i; objlst[j]; j++) {
99 cd5bae78 2004-04-21 devnull if(i == j)
100 cd5bae78 2004-04-21 devnull continue;
101 cd5bae78 2004-04-21 devnull q = objlst[j];
102 cd5bae78 2004-04-21 devnull if(p == &omoon || q == &omoon) {
103 cd5bae78 2004-04-21 devnull occult(p, q, 0);
104 cd5bae78 2004-04-21 devnull if(occ.t3 < 0)
105 cd5bae78 2004-04-21 devnull continue;
106 cd5bae78 2004-04-21 devnull if(p == &osun || q == &oshad) {
107 cd5bae78 2004-04-21 devnull if(occ.t1 >= 0)
108 cd5bae78 2004-04-21 devnull event("Partial eclipse of %s begins at ", p->name, "",
109 cd5bae78 2004-04-21 devnull occ.t1, SIGNIF|PTIME);
110 cd5bae78 2004-04-21 devnull if(occ.t2 >= 0)
111 cd5bae78 2004-04-21 devnull event("Total eclipse of %s begins at ", p->name, "",
112 cd5bae78 2004-04-21 devnull occ.t2, SIGNIF|PTIME);
113 cd5bae78 2004-04-21 devnull if(occ.t4 >= 0)
114 cd5bae78 2004-04-21 devnull event("Total eclipse of %s ends at ", p->name, "",
115 cd5bae78 2004-04-21 devnull occ.t4, SIGNIF|PTIME);
116 cd5bae78 2004-04-21 devnull if(occ.t5 >= 0)
117 cd5bae78 2004-04-21 devnull event("Partial eclipse of %s ends at ", p->name, "",
118 cd5bae78 2004-04-21 devnull occ.t5, SIGNIF|PTIME);
119 cd5bae78 2004-04-21 devnull } else {
120 cd5bae78 2004-04-21 devnull if(occ.t1 >= 0)
121 cd5bae78 2004-04-21 devnull event("Occultation of %s begins at ", q->name, "",
122 cd5bae78 2004-04-21 devnull occ.t1, SIGNIF|PTIME);
123 cd5bae78 2004-04-21 devnull if(occ.t5 >= 0)
124 cd5bae78 2004-04-21 devnull event("Occultation of %s ends at ", q->name, "",
125 cd5bae78 2004-04-21 devnull occ.t5, SIGNIF|PTIME);
127 cd5bae78 2004-04-21 devnull continue;
129 cd5bae78 2004-04-21 devnull if(p == &osun) {
130 cd5bae78 2004-04-21 devnull if(q != &omerc && q != &ovenus)
131 cd5bae78 2004-04-21 devnull continue;
132 cd5bae78 2004-04-21 devnull occult(p, q, -1);
133 cd5bae78 2004-04-21 devnull if(occ.t3 >= 0.) {
134 cd5bae78 2004-04-21 devnull if(occ.t1 >= 0)
135 cd5bae78 2004-04-21 devnull event("Transit of %s begins at ", q->name, "",
136 cd5bae78 2004-04-21 devnull occ.t1, SIGNIF|LIGHT|PTIME);
137 cd5bae78 2004-04-21 devnull if(occ.t5 >= 0)
138 cd5bae78 2004-04-21 devnull event("Transit of %s ends at ", q->name, "",
139 cd5bae78 2004-04-21 devnull occ.t5, SIGNIF|LIGHT|PTIME);
141 cd5bae78 2004-04-21 devnull continue;
143 cd5bae78 2004-04-21 devnull t = dist(&p->point[0], &q->point[0]);
144 cd5bae78 2004-04-21 devnull if(t > 5000)
145 cd5bae78 2004-04-21 devnull continue;
146 cd5bae78 2004-04-21 devnull event("%s is in the house of %s",
147 cd5bae78 2004-04-21 devnull p->name, q->name, 0, 0);
150 cd5bae78 2004-04-21 devnull if(flags['o'])
151 cd5bae78 2004-04-21 devnull stars();
152 cd5bae78 2004-04-21 devnull if(flags['a'])
153 cd5bae78 2004-04-21 devnull satels();
154 cd5bae78 2004-04-21 devnull evflush();
