1 REM Ins{nt av 2314
10 REM SAVEOS
20 NOTRACE
30 REM COPYRIGHT DR THOMAS A CLARK
40 REM PHASE3BR
50 REM ELLIPTICAL ORBIT SATELLITE AZIMUTH AND EL OSCAR 10
60 REM SATELLITE TRACKING SOFTWARE FOR THE RADIO AMATEUR
70 REM OMARBETAT F\R ABC 80 AV SM6EFW GUNNAR OLSSON FLODA
80 REM KEPLERIAN PARAMETERS FR]N NASA ENL. AMSAT-SM DEN 1983-11-06
90 ; INP(3) : GOTO 770
100 LET D1=ABS(W5-L0)
110 IF D1>180 THEN LET D1=360-D1
120 LET D9=(SIN(A1*P0))*(SIN(L5*P0))+(COS(A1*P0))
130 LET D0=(COS(A1*P0))*(COS(L5*P0))*(COS(D1*P0))
140 LET D2=(D9)*(D0)
150 LET D2=57.3*(-ATN(D2/SQR((.999999-(D2*D2)+.000001)))+PI/2)
160 IF D3>D2 THEN PRINT TAB(36);"IN"
170 IF D3.1 THEN LET A0=26105.8
1230 IF N0<=.1 THEN LET N0=SQR(G0/(A0*A0*A0))
1240 LET E2=1-(E0*E0)
1250 LET E1=SQR(E2)
1260 LET F1=0
1270 LET Q0=M0/360+K0
1280 PRINT TAB(7);"OSCAR 10 ON";D4;"-";M;"-";Y
1290 PRINT " UTC";TAB(8);"AZ";TAB(13);"EL"
1300 PRINT "----------------"
1310 LET K2=9.95*((R0/A0)^3.5)/(E2*E2)
1320 LET S1=SIN(I0*P0)
1330 LET C1=COS(I0*P0)
1340 LET O=O0-(D-D0)*K2*C1-(T-T0)*K2*C1
1350 LET S0=SIN(O*P0)
1360 LET C0=COS(O*P0)
1370 LET W=W0+(T-T0)*K2*(2.5*(C1*C1)-.5)+(D-D0)*K2*(2.5*(C1*C1)-.5)
1380 LET S2=SIN(W*P0)
1390 LET C2=COS(W*P0)
1400 DIM C(3,2)
1410 LET C(1,1)=(C2*C0)-(S2*S0*C1)
1420 LET C(1,2)=-(S2*C0)-(C2*S0*C1)
1430 LET C(2,1)=(C2*S0)+(S2*C0*C1)
1440 LET C(2,2)=-(S2*S0)+(C2*C0*C1)
1450 LET C(3,1)=(S2*S1)
1460 LET C(3,2)=(C2*S1)
1470 LET Q=N0*(D-D0)+N0*(T-T0)+Q0
1480 LET K=INT(Q)
1490 LET M9=INT((Q-K)*256)
1500 LET M=(Q-K)*P2
1510 LET E=M+E0*SIN(M)+.5*(E0*E0)*SIN(2*M)
1520 LET S3=SIN(E)
1530 LET C3=COS(E)
1540 LET R3=1-E0*C3
1550 LET M1=E-E0*S3
1560 LET M5=M1-M
1570 IF ABS(M5)<.000011 THEN GOTO 1600
1580 LET E=E-M5/R3
1590 GOTO 1520
1600 LET X0=A0*(C3-E0)
1610 LET Y0=A0*E1*S3
1620 LET R=A0*R3
1630 LET X1=X0*C(1,1)+Y0*C(1,2)
1640 LET Y1=X0*C(2,1)+Y0*C(2,2)
1650 LET Z1=X0*C(3,1)+Y0*C(3,2)
1660 LET G7=T*G1+D*G1+G2
1670 LET G7=(G7-(INT(G7)))*P2
1680 LET S7=-SIN(G7)
1690 LET C7=COS(G7)
1700 LET X=(X1*C7)-(Y1*S7)
1710 LET Y=(X1*S7)+(Y1*C7)
1720 LET Z=Z1
1730 LET L8=L9*P0
1740 LET S9=SIN(L8)
1750 LET C9=COS(L8)
1760 LET S8=SIN(-W9*P0)
1770 LET C8=COS(W9*P0)
1780 LET R9=R0*(1-(F/2)+(F/2)*COS(2*L8))+H9/1000
1790 LET L8=ATN((1-F)*(1-F)*S9/C9)
1800 LET Z9=R9*SIN(L8)
1810 LET X9=R9*COS(L8)*C8
1820 LET Y9=R9*COS(L8)*S8
1830 LET X5=(X-X9)
1840 LET Y5=(Y-Y9)
1850 LET Z5=(Z-Z9)
1860 LET R5=SQR(X5*X5+Y5*Y5+Z5*Z5)
1870 LET D3=6378/R
1880 LET D3=57.3*(-ATN(D3/SQR(1-D3*D3))+PI/2)
1890 LET Z8=(X5*C8*C9)+(Y5*S8*C9)+(Z5*S9)
1900 LET X8=-(X5*C8*S9)-(Y5*S8*S9)+(Z5*C9)
1910 LET Y8=(Y5*C8)-(X5*S8)
1920 LET S5=Z8/R5
1930 LET C5=SQR(1-S5*S5)
1940 LET E9=(ATN(S5/C5))/P0
1950 IF E9<-10 THEN GOTO 2350
1960 LET A9=(ATN(Y8/X8))/P0
1970 LET B5=Z/R
1980 LET L5=(ATN(B5/(SQR(1-B5*B5))))*57.3
1990 LET W5=(ATN(Y/X))*57.3
2000 IF X<0 THEN LET W5=180-W5
2010 IF X>0 AND Y<0 LET W5=-W5
2020 IF X>0 AND Y>0 THEN LET W5=360-W5
2030 IF X=0 AND Y>=0 THEN LET W5=270
2040 IF X=0 AND Y<0 THEN LET W5=90
2050 IF X8<0 LET A9=A9+180
2060 IF X8>0 AND Y8<0 THEN LET A9=360+A9
2070 IF X8=0 AND Y8>=0 THEN LET A9=90
2080 IF X8=0 AND Y8<0 LET A9=270
2090 LET T4=INT(T*24)*100
2100 LET T5=INT(T*1440)-T4*.6
2110 LET T6=T4+T5
2120 GOSUB 750
2130 IF F1=1 THEN GOTO 2220
2140 PRINT T$;TAB(6);INT(A9);TAB(11);INT(E9*10)/10
2150 LET T=T+1/96
2160 IF (T-T7)>.17 THEN GOTO 2180
2170 GOTO 1340
2180 LET F1=1
2190 LET T=T7
2200 REM PRINT CUR(21,0);"TO STOP H "A"TILL ""FR""PRINTS.TO UNFRE PR"C" "
2210 GOTO 1340
2220 PRINT CUR(2,19);"DX STATE AT ";T$
2230 PRINT CUR(3,19);"---------------------"
2240 GOSUB 190
2250 GET X$ : IF X$="A" THEN GOTO 2280
2260 IF X$="C" GOTO 2340
2270 IF X$<>"A" GOTO 2350
2280 PRINT CUR(18,25);"FROZEN" : PRINT CUR(21,0)" "
2290 PRINT CUR(21,0);"** SLUT P] DATA **"; : REM GOTO 3610
2300 PRINT CUR(22,0);" ----MERA ? (J/N)----";
2310 GET \$ : IF \$="J" GOTO 10
2320 IF \$<>"J" GOTO 2380
2330 GOTO 2260
2340 PRINT CUR(18,25);'WAIT'
2350 LET T=T+1/96
2360 IF (T-T7)>.17 THEN GOTO 2180
2370 GOTO 1340
2380 REM SLUT