1 REM Ins{nd av John Stagg <5741>    1985-08-11 18.43.35
10 ! M]NE 850810 
20 ! Ursprungligen EME av Jinge Fl~cht (ABC 4716). F|r 80-tecken ABC800-serien
30 EXTEND : FLOAT 
40 ; CHR$(12)
50 Bl$=STRING$(80,32) : Sp$=SPACE$(25)
60 ; STRING$(80,ASCII("="))
70 ; "           Program f|r utr{kning av elevation och riktning till m}nen"
80 ; STRING$(80,ASCII("=")) : ; 
90 IF S9 THEN GOTO 210
100 ; "Dagens datum (RETN ger datorns interntid): "; : INPUT Date$
110 IF Date$='' THEN GOTO 190
120 ; 'Klockan :'; : INPUT Kl$ : IF LEN(Kl$)=3 THEN Kl$='0'+Kl$
130 T=VAL(LEFT$(Kl$,2))
140 N=VAL(RIGHT$(Kl$,3))
150 D=VAL(MID$(Date$,5,2))
160 M=VAL(MID$(Date$,3,2))
170 Y=VAL(LEFT$(Date$,2))
180 POKE 65519,Y,M,D,T,N
190 IF PEEK(65519)=0 OR PEEK(65520)=255 THEN GOTO 40
200 ; '[r det sommartid '; : INPUT T$ : T$=T$+' ' : S=(LEFT$(T$,1)='j') OR (LEFT$(T$,1)='J')
210 T$=TIME$ : Y=PEEK(65519)+1900 : M=PEEK(65520) : D=PEEK(65521) : S9=-1
220 ; CUR(5,0);Bl$ : ; Bl$ : ; Bl$
230 ; CUR(6,0);'M}nens position (grader) i himlen ';LEFT$(TIME$,16);' :' : ; : ; : ; 
240 A3$=MID$(T$,12,2)
250 B3$=MID$(T$,15,2)
260 B=VAL(A3$+B3$)
270 DEF FNA(X)=INT(X*D5*10+.5)/10 
280 DEF FNB(X)=(X-INT(X))*P5 
290 P5=2*PI
300 D5=360/P5
310 R5=P5/360
320 L5=59
330 U5=10
340 L6=-18
350 U6=4
360 IF L6<0 THEN U6=-U6
370 L5=(L5+U5/60)*R5
380 L6=(L6+U6/60)*R5
390 I6=100
400 Y1=Y-(INT(Y/100)*100)
410 I1=2
420 IF M>3 THEN 500
430 IF INT((Y-1853)/4)<11 THEN 460
440 C1=-1
450 GOTO 470
460 C1=0
470 J1=365*(Y-1853)+D+30*(M+9)+INT((M+10)/2)
480 J2=INT((Y-1853)/4)+1+C1
490 GOTO 610
500 IF INT((Y-1852)/4)<11 THEN 530
510 C1=-1
520 GOTO 540
530 C1=0
540 IF M=9 THEN 580
550 IF M=11 THEN 580
560 C2=0
570 GOTO 590
580 C2=1
590 J1=365*(Y-1852)+D+30*(M-3)+INT((M-2)/2)
600 J2=INT((Y-1852)/4)+C1+C2
610 J=J1+J2
620 T1=J-17472.5
630 T=(B-INT(B/100)*100)/1440+INT(B/100)/24
640 T=T+(S-1)/24 ! Konv till GMT
650 T5=T1+T
660 K1=FNB(.751213+3.6601102E-02*T5)
670 K2=FNB(.822513+3.62916457E-02*T5)
680 K3=FNB(.995766+2.73777852E-03*T5)
690 K4=FNB(.974271+3.38631922E-02*T5)
700 K5=FNB(3.12525E-02+3.67481957E-02*T5)
710 L8=K1+.658*R5*SIN(2*K4)+6.289*R5*SIN(K2)
720 L8=L8-1.274*R5*SIN(K2-2*K4)-.186*R5*SIN(K3)
730 L8=L8+.214*R5*SIN(2*K2)-.114*R5*SIN(2*K5)
740 L8=L8-.059*R5*SIN(2*K2-2*K4)-.057*R5*SIN(K2+K3-2*K4)
750 K6=K5+.6593*R5*SIN(2*K4)+6.2303*R5*SIN(K2)-1.272*R5*SIN(K2-2*K4)
760 L7=5.144*R5*SIN(K6)-.146*R5*SIN(K5-2*K4)
770 D1=COS(L7)*SIN(L8)*.397821+SIN(L7)*.917463
780 D1=ATN(D1/(SQR(1-D1^2))) : ! KOLLA EXPONENTIERINGEN
790 A2=COS(L7)*COS(L8)/COS(D1)
800 A1=(COS(L7)*SIN(L8)*.917463-SIN(L7)*.397821)/COS(D1)
810 A=ATN(A1/A2)
820 GOSUB 1060
830 R1=A
840 L1=6.5709822E-02*T1
850 L=T*24*1.002738+6.646055+(L1-INT(L1/24)*24)
860 L=(L-INT(L/24)*24)
870 G=(L/24)*2*PI-R1
880 IF G<P5 THEN 910
890 G=G-P5
900 GOTO 940
910 IF G<0 THEN 930
920 GOTO 940
930 G=G+P5
940 H=L6-G
950 E3=COS(L5)*COS(H)*COS(D1)+SIN(D1)*SIN(L5)
960 E2=SQR(1-(E3*E3))
970 E=ATN(E3/E2)
980 IF E<-.5 THEN 1220
990 A2=SIN(D1)/(COS(L5)*COS(E))
1000 A2=A2-(SIN(L5)/COS(L5))*(SIN(E)/COS(E))
1010 A1=SIN(L5)*SIN(D1)+COS(L5)*COS(D1)*COS(H)
1020 A1=(SIN(H)*COS(D1))/(SQR(1-A1^2)) : ! KOLLA EXP.
1030 A=ATN(A1/A2)
1040 GOSUB 1060
1050 GOTO 1220
1060 IF A=0 THEN 1080
1070 GOTO 1120
1080 IF A2<0 THEN 1100
1090 GOTO 1210
1100 A=P5/2
1110 GOTO 1210
1120 IF A>0 THEN 1180
1130 IF A2<0 THEN 1160
1140 A=P5+A
1150 GOTO 1210
1160 A=P5+(A-P5/2)
1170 GOTO 1210
1180 IF A2>0 THEN 1210
1190 IF A2=0 THEN 1210
1200 A=A+P5/2
1210 RETURN 
1220 ; 
1230 Nere=FNA(E)<-.5
1240 ; Sp$;"Azimut = ";FNA(A);'  '
1250 ; Sp$;"Elevation = ";FNA(E);'  '
1260 ; Sp$;"GH = ";FNA(G);'  '
1270 ; Sp$;"DC = ";FNA(D1);'  '
1280 IF Nere THEN ; CUR(21,0);'(M}nen {r nu under horisonten)';
1290 ; CUR(21,60);'Avbryt med RETN';
1300 A3$=MID$(TIME$,12,2)
1310 B3$=MID$(TIME$,15,2)
1320 D3$=LEFT$(TIME$,2)
1330 Y$=A3$+B3$
1340 O=VAL(Y$)
1350 IF SYS(5) THEN GET T$ : IF T$=CHR$(13) THEN ; CHR$(12); : GOTO 1380
1360 IF O=B THEN GOTO 1300
1370 GOTO 210
1380 END