恒星の出没現象の計算を行うパッケージ
[パッケージrs.lzhの構成] rs.exe - DOS/V用の恒星の出没現象計算プログラム 使用法: []は省略化 set RS_TBL=xxxx.tbl (この指定が必須になりました) rs 計算年(西暦) FK5星表番号 観測地の緯度 [恒星の高度 [太陽の伏角]] 子プロセスとしてwhenhv.exeを使用しますので、whenhv.exeに コマンドパスが通っている必要があります。 rs.c - rs.exe のソースリスト rss.bat - 主な恒星に対する計算をまとめて行うバッチファイル neuge.tbl - P.V.Neugebauerの太陽伏角表 aveni.tbl - A.F.Aveniの結果から推定した太陽伏角表 nipe.fk5 - ノーバスさん(GBA02113@nifty.com)作成のFK5星表 (転載をご了解いただいたノーバスさんに感謝します)
[計算原理 - 現象の起こる時の太陽黄経を求める] 1) 現象の起こる恒星の出没時の太陽の伏角 β を知る Tafel 28(m,type -> β) RS.EXE は 環境変数RS_TBLに指定される表を参照して計算する (SET RS_TBL=NEUGE.TBL など) neuge.tbl: P.V.Neugebauerに従う aveni.tbl: A.F.Aveni に近い 2) 歳差を考慮して恒星の赤経 α 赤緯 δ を知る Tafel I and IV (in “Sterntafeln von 4000 vor Chr. bis zur Gegenwart”) 3) 観測地の地理学的経度 λ 緯度 φ を知る 4) 恒星の出没時の天頂と恒星の赤経の差 t を求める 図1-現代風の計算方法, 補足-Neugebauerの方法(下記) Tafel 1(δ,φ -> t) tan^2 (t'/2) = cos(φ-δ)/cos(φ+δ) 赤緯がφで恒星と90度の角距離のある点 t'' = 0.039h/(cosφcosδsint') 大気差による補正(35'8"なら0.039037h) t = t' + t'' 5) 恒星の出没時の天頂の赤道座標 θ = α +/- t 赤経 φ 赤緯 6) 恒星の出没時の天頂の黄道座標(表を使える様に工夫しているが、 現代では6')の方が判りやすい) 6.1) 計算用パラメータ p,r, 垂直引数 s を求める Tafel 24(θ -> p,r,s) p = +2.5sin 2θ tan r = -sinθtanε sin s = +cosθsinε (ε=23.70 <- ちょっと大きすぎる.疑問!) 6.2) 水平引数 Σ を求める Σ = φ + r 6.3) 計算用パラメータ P を求める Tafel 25(s,Σ -> P) tan P = sin s tan Σ 6.4) 黄経 L を求める L = θ + p + P 6.5) 黄緯 B を求める Tafel 26(s,Σ -> B) sin B = cos s sinΣ 7) 現象が起こる時の地平線から太陽までの黄道の長さ L1 を求める(図2) Tafel 27(B,β -> L1) sin L1 = sinβ/cos B 8) 求める太陽黄経は L +/- (90 + L1) 6') 恒星の出没時の天頂の黄道座標(現代風) 6'.1) 赤道座標を直交座標で表す x = cosφcosθ y = cosφsinθ z = sinφ 6'.2) 黄道傾角だけ x 軸の回りに回転し、黄道直交座標に直す X = cosφcosθ Y = cosφsinθcosε + sinφsinε Z = -cosφsinθsinε + sinφcosε 6'.3) 黄道極座標は tan L = Y/X sin B = Z