14.01.2015 Views

OF THE ROGER N. CLARK

OF THE ROGER N. CLARK

OF THE ROGER N. CLARK

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!