OF THE ROGER N. CLARK
OF THE ROGER N. CLARK
OF THE ROGER N. CLARK
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
VISUAL ASTRONOMY <strong>OF</strong> <strong>THE</strong> DEEP SKY<br />
Table C.2 (cont. )<br />
--<br />
APPENDIX C: AIR MASS, ATMOSPHERIC EXTINCTION<br />
Table C.2 (cont. )<br />
1605 XSIGN = 1: Al = ABS(VAL(A$ »<br />
1610 IF LEFT$(A$ ,l)="-" <strong>THE</strong>N XSIGN= -1<br />
1615 A# = XSIGN* (A1+A2/60#+A3/3600# )<br />
1620 RETURN<br />
1650 REM ** ******** Ddeg to DMS<br />
1651 REM subroutine to convert decimal degrees (+ or -) to Deg , Min, Sec.<br />
1652 REM<br />
1655 XSIGN =SGN(DDEG#): X#=ABS (DDEG#): IF XSIGN=O <strong>THE</strong>N XSIGN=l<br />
1660 DEG INT(X#)*XSIGN: X#=(X#-INT(X#» *60#<br />
1665 MIN = INT(X#): X#=(X#-INT(X#» *60#<br />
1670 SEC = X#<br />
1675 IF SEC > 59.99 <strong>THE</strong>N SEC SEC-60 : MIN MIN + XSIGN<br />
1680 IF MIN > 59.99 <strong>THE</strong>N MIN MIN=60: DEG DEG + XSIGN<br />
1685 RETURN<br />
1700 REM ** ******** compute Alt , Az, Airm<br />
1705 REM subroutine to compute altitude , azimuth, airmass, hour angle<br />
1710 REM<br />
1715 REM input: ST, RA, DEC , LAIR , SO<br />
1720 REM output is corrected for refraction<br />
1725 REM<br />
1730 GOSUB 1200 REM compHA<br />
1735 GOSUB 1250 REM Zenithangle<br />
1736 GOSUB 2500 REM refcorr<br />
1740 GOSUB 2550 REM Altitude<br />
1745 GOSUB 2600 REM Azimuth<br />
1750 GOSUB 2650 REM Airmass<br />
1755 RETURN<br />
1799 REM ********** Cal to JD<br />
1800 REM SUBROUTINE to convert Calendar day to Julian Day<br />
1805 REM adapted from Sky and Telescope May 1984 , page 454-455.<br />
1810 REM input: year=y month=m day=d<br />
1812 REM g=flag for gregorian calendar (begin 1752 in england , colonies )<br />
1813 REM (begin 1873 in japan , 1927 in turkey )<br />
1814 REM g=l gregorian calendar , =0 for Julian calendar<br />
1815 G=l<br />
1816 IF Y=O <strong>THE</strong>N RETURN<br />
1865 F=F+1 : J=J- 1<br />
1870 RETURN<br />
1899 REM ** *** ***** JD to Cal<br />
1900 REM SUBROUTINE to convert Julian Day to Ca lendar Day .<br />
1901 REM<br />
1902 REM Adapted from Sky and Telescope May , 1984 page 454-455 .<br />
1905 REM input : J= Julian Day (integer part ), F= fractional part of day<br />
1906 REM ou tput : Y= year , M= month , D=day (including fractional part )<br />
-<br />
-<br />
19 10 G=l :FTMP=F : JTMP=J<br />
19 1 5 IF J 59.995 <strong>THE</strong>N MIN MIN-60 : HOURS = HOURS + XSIGN<br />
2090 RETURN<br />
2100 REM ********** compSO<br />
2110 REM subroutine to compute side real time at Greenwich at o hours UT<br />
2120 REM on the date DD MM YY (day , month, year )<br />
2130 REM J= Julian Centuries elapsed from DD MM YY<br />
2140 Y= YY<br />
2150 M= MM<br />
2160 D= DD<br />
2170 GOSUB 1800 REM Cal to JD<br />
2173 JJ#<br />
2176 JF#<br />
2180 Jl#<br />
2185 REM<br />
2190 JO#<br />
2200 SOH<br />
220 1 SOU<br />
220 5 SOil<br />
2210 SOli<br />
2220 RETURN<br />
J<br />
F<br />
JJ/I + JFII<br />
Julian day 2451545.0 is January 1, 2000 at 12 hours UT<br />
(Jll1 - 2451 54511)/3652511<br />
24 110. 5484111+8640184.81286611*JOII+ .09310411*JOII*JOII<br />
SOli - .000006211*JOII*JO#*JOII<br />
SOil /8640011<br />
SOil - INT( SOil)<br />
276<br />
277