5 POKE55,.:POKE56,56:CLR
10 DV=PEEK(186):IFDV<8THENDV=8
15 POKE53280,.:POKE53281,.:PRINT"[147]"
16 POKE53371,0
20 AD=49152
22 SYSAD:SYSAD+12
25 POKE53272,29
30 GOSUB400
40 POKE53272,31:SYSAD+9,12
45 BX$="_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_"
50 PRINT"[147][150]^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^"
55 PRINTBX$;""TAB(38)BX$
60 PRINT"[150]^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^"
80 PRINT"":PRINTTAB(6)"[150] [198] [193] [210] [211] [201] [199] [200] [212] [206] [197] [211] [211] "
85 PRINTTAB(7)"[153] [211][213][206]-[205][207][207][206] [211][197][212][213][208] "
90 PRINTTAB(2)"[156][205]OON [208]HASE [159]([198]=0,[209]TS=90,[206]=180):";:L9%=3:GOSUB1355:AM=Q9
100 PRINTTAB(2)"[156][205]OON FROM [218]ENITH [159]([219]'S): ";:L9%=3:GOSUB1355:ZM=Q9
110 PRINTTAB(2)"[156][205]OON'S [193]NGULAR DISTANCE [159]([219]'S)";:L9%=3:GOSUB1355:RM=Q9
120 PRINTTAB(2)"[156][211]UN FROM [218]ENITH [159]([219]'S): ";:L9%=3:GOSUB1355:ZS=Q9
130 PRINTTAB(2)"[156][211]UN'S ANGULAR DISTANCE [159]([219]'S): ";:GOSUB1355:RS=Q9
140 PRINTTAB(7)"[153] [211][201][212][197] [197][206][214][201][210][207][206][205][197][206][212] "
150 PRINTTAB(2)"[156][210]ELATIVE [200]UMIDITY: ";:L9%=6:GOSUB1355:RH=Q9
160 PRINTTAB(2)"[156][212]EMPERATURE IN [198]AHRENHEIT: ";:L9%=3:GOSUB1355:FH=Q9
165 TE=(FH-32)*5/9
170 PRINTTAB(2)"[156][217]OUR [204]ATITUDE: ";:L9%=6:GOSUB1355:LA=Q9
180 PRINTTAB(2)"[156][193]LTITUDE IN FEET: ";:L9%=5:GOSUB1355:AL=Q9
185 AL=AL/3.3
188 PRINTTAB(7)"[153] [196][193][212][197] "
190 PRINTTAB(2)"[156][212]HE [205]ONTH- [159](1-12):";:L9%=2:GOSUB1355:M=Q9
200 PRINTTAB(2)"[156][212]HE [217]EAR: ";:L9%=4:GOSUB1355:Y=Q9
210 SN=1:REM SNELLEN RATIO (20/20=1, GOOD 20/10=2
215 PRINTTAB(7)"[153] [207][194][202][197][195][212]'[211] [204][207][195][193][212][201][207][206] "
220 PRINTTAB(2)"[156][218]ENITH DISTANCE [159]([219]'S):";:L9%=4:GOSUB1355:Z=Q9
222 POKE214,20:PRINT:PRINTTAB(9)"[158][201]S THIS CORRECT? [154][217][150]/[154][206]":POKE198,0
223 GETA$:IFA$<>"Y"ANDA$<>"N"THEN223
225 IFA$="N"THEN40
229 PRINT"[147]":SYSAD+9,14
230 PRINT"[159]^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^"
232 PRINTBX$:PRINT""TAB(38)BX$
235 PRINT"[159]^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^"
248 GOSUB1000:REM EXTINCTION
249 GOSUB2000:REM SKY
250 SYSAD+9,14:REM VISUAL LIMITING MAGNITUDE
260 BL=B(3)/1.11E-15:REM IN NANOLAMBERTS
270 IFBL>1500THEN300
280 C1=10^-9.8:C2=10^-1.9
290 GOTO310
300 C1=10^-8.350001:C2=10^-5.9
310 TH=C1*((1+SQR(C2*BL))^2):REM IN FOOT-CANDLES
320 MN=-16.57-2.5*(LOG(TH)/LOG(10))-DM(3)+5*(LOG(SN)/LOG(10))
340 PRINT:PRINTTAB(2)"[150]-[153][211]KY [194]RIGHTNESS ([155]N[204][153])[158]:";BL
350 PRINTTAB(2)"[150]-[153][204]IMITING [205]AGNITUDE[158]:";MN
360 GOSUB13000
365 RESTORE:PRINT"[147]":GOTO30
400 FORI=1TO5:READWA(I):NEXTI
410 DATA.365,.44,.55,.7,.9
420 FORI=1TO5:READMO(O):NEXTI
425 DATA-10.93,-10.45,-11.05,-11.9,-12.7
430 RD=(null)/180
435 RETURN
1000 REM EXTINCTION SUBROUTINE
1010 FORI=1TO5:READOZ(I):NEXTI
1020 DATA0,0,.031,.008,0
1030 FORI=1TO5:READWT(I):NEXTI
1040 DATA.074,.045,.031,.02,.015
1050 LT=LA*RD
1060 RA=(M-3)*30*RD
1070 SL=LA/ABS(LA)
1080 REM AIRMASS FOR EACH COMPONENT
1090 ZZ=Z*RD
1100 XG=1/(COS(ZZ)+.0286*EXP(-10.5*COS(ZZ)))
1110 XA=1/(COS(ZZ)+.0123*EXP(-24.5*COS(ZZ)))
1120 XO=1/SQR(1-(SIN(ZZ)/(1+(20/6378)))^2)
1130 REM UBVRI EXTINCTION FOR EACH COMPONENT
1140 FORI=1TO5
1150 KR=.1066*EXP(-1*AL/8200)*((WA(I)/.55)^-4)
1160 KA=.1*((WA(I)/.55)^-1.3)*EXP(-1*AL/1500)
1170 KA=KA*((1-.32/LOG(RH/100))^1.33)*(1+SL*SIN(RA))
1180 KO=OZ(I)*(3+.4*(LT*COS(RA)-COS(3*LT)))/3
1190 KW=WT(I)*.94*(RH/100)*EXP(TE/15)*EXP(-1*AL/8200)
1200 K(I)=KR+KA+KO+KW
1210 DM(I)=KR*XG+KA*XA+KO*XO+KW*XG
1220 NEXTI
1230 REM WRITE RESULTS AND RETURN
1240 PRINT"":PRINTTAB(2)"[153] [193][201][210][205][193][211][211] [146][159]-[199]AS-  [158]-[193]EROSOL-  [154]-[207]ZONE-:"
1250 PRINTTAB(1)"[159]"XG,"[158]"XA
1255 PRINTTAB(1)"[154]"XO
1260 PRINTTAB(2)"[153] [197]XTINCTION [195]OEFFICIENTS [146] [154][213]-[159][194]-[214]-[150][210]-[129][201]:"
1270 PRINTTAB(1)"[154]"K(1),"[159]"K(2)
1273 PRINTTAB(1)""K(3),"[150]"K(4)
1275 PRINTTAB(1)"[129]"K(5)
1280 PRINTTAB(2)"[153] [212]OTAL [197]XTINCTION [146] [154][213]-[159][194]-[214]-[150][210]-[129][201]:"
1290 PRINTTAB(1)"[154]"DM(1),"[159]"DM(2)
1292 PRINTTAB(1)""DM(3),"[150]"DM(4)
1294 PRINTTAB(1)"[129]"DM(5)
1295 SYSAD+9,4
1300 RETURN
1355 Q9$="":POKE198,.
1360 GETA$
1365 POKE646,RND(1)*15+1:PRINT"*[157]";:IFA$=""THEN1360
1370 IFA$=CHR$(13)THENPRINT" ":Q9=VAL(Q9$):RETURN
1375 IF(A$=CHR$(20)ANDLEN(Q9$))THENQ9$=LEFT$(Q9$,LEN(Q9$)-1):GOTO1397
1380 IFLEN(Q9$)>=L9%THENSYSAD+9,6:GOTO1360
1385 IF(A$>="0"ANDA$<="9")ORA$="."ORA$="-"THEN1390
1387 GOTO1360
1390 Q9$=Q9$+A$
1395 PRINT"[158]"A$;:SYSAD+9,13:GOTO1360
1397 PRINT" [157][157] [157]";:GOTO1360
2000 REM SKY ROUTINE
2010 FORI=1TO5:READBO(I):NEXTI
2020 DATA8E-14,7E-14,1E-13,1E-13,3E-13
2030 FORI=1TO5:READCM(I):NEXTI
2040 DATA1.36,.91,0,-.76,-1.17
2050 FORI=1TO5:READMS(I):NEXTI
2060 DATA-25.96,-26.09,-26.74,-27.26,-27.55
2070 X=1/(COS(ZZ)+.025*EXP(-11*COS(ZZ))):REM AIR MASS
2080 XM=1/(COS(ZM*RD)+.025*EXP(-11*COS(ZM*RD))):REM AIR MASS MOON
2090 IFZM>90THENXM=40
2100 XS=1/(COS(ZS*RD)+.025*EXP(-11*COS(ZS*RD))):REM AIR MASS SUN
2110 IFZS>90THENXS=40
2120 FORI=1TO5
2130 REM DARK NIGHT SKY BRIGHTNESS
2140 BN=BO(I)*(1+.3*COS(6.283*(Y-1992)/11))
2150 BN=BN*(.4+.6/SQR(1-.96*((SIN(ZZ))^2)))
2160 BN=BN*(10^(-.4*K(I)*X))
2170 REM MOONLIGHT BRIGHTNESS
2180 MM=-12.73+.026*ABS(AM)+4E-09*(AM^4):REM MOON MAG IN V
2190 MM=MM+CM(I):REM MOON MAG
2200 C3=10^(-.4*K(I)*XM)
2210 FM=6.2E+07*(RM^-2)+(10^(6.16-RM/40))
2220 FM=FM+(10^5.36)*(1.06+((COS(RM*RD))^2))
2230 BM=10^(-.4*(MM-MO(I)+43.27))
2240 BM=BM*(1-10^(-.4*K(I)*X))
2250 BM=BM*(FM*C3+440000*(1-C3))
2260 REM TWIGHT BRIGHTNESS
2270 HS=90-ZS:REM HEIGHT OF SUN
2280 BT=10^(-.4*(MS(I)-MO(I)+32.5-HS-(Z/(360*K(I)))))
2290 BT=BT*(100/RS)*(1-10^(-.4*K(I)))
2300 REM DAYLIGHT BRIGHTNESS
2310 C4=10^(-.4*K(I)*XS)
2320 FS=6.2E+07*(RS^-2)+(10^(6.15-RS/40))
2330 FS=FS+(10^5.36)*(1.06+((COS(RS*RD))^2))
2340 BD=10^(-.4*(MS(I)-MO(I)+43.27))
2350 BD=BD*(1-10^(-.4*K(I)*X))
2360 BD=BD*(FS*C4+440000*(1-C4))
2370 REM TOTAL SKY BRIGHTNESS
2380 IFBD>BTTHEN2410
2390 B(I)=BN+BD
2400 GOTO2420
2410 B(I)=BN+BT
2420 IFZM<90THENB(I)=B(I)+BM
2430 NEXTI
2440 PRINTTAB(2)"[153] [211]KY [194]RIGHTNESS [146] [154][213]-[159][194]-[214]-[150][210]-[129][201]:"
2450 PRINTTAB(1)"[154]"B(1),"[159]"B(2)
2452 PRINTTAB(1)""B(3),"[150]"B(4)
2455 PRINTTAB(1)"[129]"B(5)
2460 RETURN
10000 D=PEEK(186):N$="0:FARSIGHTED":OPEN15,D,15,"S"+N$:CLOSE15:SAVEN$,D:END
13000 POKE214,20:PRINT:PRINTTAB(8)"(1) [212]RY ANOTHER ONE
13010 [153][163]8)"(2) (null)O (null)(null)RIGHT$(null)VAL(null)(null)VAL (null)ENU
13020 POKE198,0
13030 GETA$:IFA$<"1"ORA$>"2"THEN13030
13040 IFA$="1"THENRETURN
13050 SYSAD+15
13060 PRINT"[147][144]LOAD"CHR$(34)"B.UNIVERSE III"CHR$(34)","DV
13070 PRINT"RUN28"
13080 POKE631,13:POKE632,13:POKE198,2:END