4000 '  Lunar position for the MC-10 and CoCo

4005 ' Paul J. Burke  P.O. Box 613 Palmdale CA 93550  CIS 74656,2333 

4010 '  PJB on Delphi

4015 ' Computes the position of the Moon for any recent date

4020 ' Lists the position in RA and Dec. Merge this program with

4025 ' ASTRO5.BAS and add prompts to use ASTRO5 features

4030 '

4035 INPUT"Turn on your printer ";A$

4040 GOTO 4150

4045 '

4050 'Julian day computation/day of year

4055 Y=INT(A):M=INT(100*(A-Y)):D=100*(100*(A-Y)-M)

4060 IF M>2 THEN 4070

4065 Y=Y-1:M=M+12

4070 A1=INT(Y/100):B=2-A1+INT(A1/4)

4075 D1=INT(365.25*Y)+INT(30.6001*(M+1))+D+1720994.5+B

4080 RETURN

4085 '

4090 ' HMS conversion

4095 H=INT(G1):M=!(ABS((G1-INT(G1)*60)):S=ABS((G1-INT(G1)*60-M)*60

4100 RETURN

4105 '

4110 ' HMS to Hours

4115 M=ABS(H1-INT(H1)*100)

4120 '

4125 INPUT"Your time zone ";TZ

4130 INPUT"Desired time of the event ";LT$:LT=VAL(LT$)

4135 LT=(LT+TZ)/24

4140 RETURN

4145 '

4150 INPUT"Date:YYYY.MMDD ";A$:A=VAL(A$):GOSUB 4125

4155 A=A+LT/10000

4160 GOSUB 4050

4165 Y$=STR$(Y):M$=STR$(M):D$=STR$(INT(100*D+.5)/100)

4170 D3=Y*365+(M-1)*31+D:IF M>2 THEN 4180

4175 D3=D3+INT((Y-1)/4)-INT(.75*INT((Y-1)/100)):GOTO4185

4180 D3=D3-INT(2.3+.4*M)+INT(Y/4)-INT(.75*INT((Y/100)+1))

4185 N2=D3-715875.5 :' days of Epoch 1960

4190 MZ=311.1687:ME=178.699:MP=255.7433

4195 P2=.111404*N2+MP:P2=ABS(P2-INT(P2/360)*360)

4200 RH=MZ+13.17636*N2:P2=RH-P2

4205 D4=6.283185*SIN(.017455328*P2):RH=RH+D4:RH=ABS(RH-INT(RH/360)*360)

4210 R1=RH/15:R1=(R1-INT(R1/24)*24):' Lunar RA

4215 A2=ME-N2*.052954:A2=ABS(A2-INT(A2/360)*360)

4220 A2=RH-A2:IF A2<0 THEN A2=A2+360

4225 IF A2>360 THEN A2=A2-360

4230 A3=5.1454*SIN(A2*.01745329)

4235 D2=A3+23.1444*SIN(RH*.01745329):' Lunar dec

4240 ' REM lines 4255 thru 4275 if no printer is attached

4245 PRINT"   Event date  ";Y$;"/";M$;"/";D$

4250 PRINT"      Lunar Position"

4255 PRINT#-2

4260 PRINT#-2,"    Event date :"Y$"/"M$"/"D$

4265 PRINT#-2,"     Lunar Position"

4270 PRINT#-2,"    Lunar R.A.     ";INT(R1*100+.5)/100

4275 PRINT#-2,"    Lunar dec      ";INT(D2*100+.5)/100

4280 PRINT"  Lunar R.A. ";INT(R1*100+.5)/100

4285 PRINT"  Lunar dec  ";TAB(15);INT(D2*100+.5)/100

4290 PRINT

4295 A=A+.0001:PRINTA:GOTO4160