| ');
document.write('日出 '+ Rs +'');
document.write(' | ');
document.write('日没 '+ Ss +'');
if (disp > 1) {
document.write(' |
');
document.write('方位 '+ ARs +'
');
document.write(' | ');
document.write('方位 '+ ASs +'
');
}
if (disp > 0) {
document.write(' |
| ');
document.write('南中 '+ Ds +'');
document.write(' | ');
document.write('高度 '+ hs +'');
}
// 月出/月没表示
document.write(' |
| ');
document.write('月出 '+ Rm +'');
document.write(' | ');
document.write('月没 '+ Sm +'');
if (disp > 1) {
document.write(' |
');
document.write('方位 '+ ARm +'
');
document.write(' | ');
document.write('方位 '+ ASm +'
');
}
if (disp > 0) {
document.write(' |
| ');
document.write('南中 '+ Dm +'');
document.write(' | ');
document.write('高度 '+ hm +'');
}
document.write(' |
');
}
//=========================================================================
// 日の出/日の入り計算
// 引数 .... 計算対象地点の経度、緯度 lng lat
// .... 出没フラグ flag(0)=出、flag(1)=没
// 戻り値 .... 出没時刻（0.xxxx日）
//=========================================================================
function calc_SUN(lng,lat,flag) {
var ans = new Array();
var T,rm_sun,r_sun,theta;
var S,W,k,dt;
var dD = 1; //補正値初期値
var D = 0.5; //逐次計算初期時刻(日)
while( Math.abs(dD) > converge ) { // 逐次計算
T = ( K + D + dT ) / 365.25; //Dの経過ユリウス年(日)
rm_sun = lng_SUN(T); //太陽の黄経
r_sun = dist_SUN(T); //太陽の距離
theta = calc_THETA(lng,T,D); //恒星時
Kou2Seki(rm_sun,0,T,ans); //黄道 -> 赤道変換
S = 0.266994 / r_sun; //太陽の視半径
W = 0.0024428 / r_sun; //太陽の視差
k = -S - R - E + W; //太陽の出没高度
dt = calc_dt(ans[0],ans[1],theta,k,lat,flag); //時角差(dt=tk-t)計算
dD = dt / 360; //仮定時刻に対する補正値
D = D + dD;
}
return D;
}
//=========================================================================
// 月の出/月の入り計算
// 引数 .... 計算対象地点の経度、緯度 lng lat
// .... 出没フラグ flag(0)=出、flag(1)=没、flag(2)=南中
// 戻り値 .... 出没/南中時刻（0.xxxx日）
//=========================================================================
function calc_MOON(lng,lat,flag) {
var ans = new Array();
var T,rm_moon,bt_moon,theta;
var W,k,dt;
var dD = 1; //補正値初期値
var D = 0.5; //逐次計算初期時刻(日)
while( Math.abs(dD) > converge ) { // 逐次計算
T = ( K + D + dT ) / 365.25; //Dの経過ユリウス年(日)
rm_moon = lng_MOON(T); //月の黄経
bt_moon = lat_MOON(T); //月の黄緯
theta = calc_THETA(lng,T,D); //恒星時
Kou2Seki(rm_moon,bt_moon,T,ans);
if (flag != 2) { //南中のときは計算しない
W = para_MOON(T); //月の視差
k = - R - E + W; //月の出没高度
}
dt = calc_dt(ans[0],ans[1],theta,k,lat,flag); //時角差(dt=tk-t)計算
dD = dt / 347.8; //仮定時刻に対する補正値
D = D + dD;
}
if (D < 0 || D >= 1) { D = ""; } //出没なし
return D;
}
//=========================================================================
// 時刻(D)における黄経、黄緯(λ(T),β(T))の天体の方位角(A)計算
// 引数 .... 天体の黄経、黄緯(λ(T),β(T))(度) ramda,beta
// .... J2000.0からの経過ユリウス年(x.xxx日) T
// .... 日本時によるその日０時からの経過時刻(0.xxx日) D
// .... 観測地点の経度、緯度 lng,lat
// 戻り値 .... 高度(xx.x度)
//=========================================================================
function calc_AZ(ramda,beta,T,D,lng,lat) {
var ans = new Array();
Kou2Seki(ramda,beta,T,ans); //黄道 -> 赤道変換
return calc_AZe(ans[0],ans[1],T,D,lng,lat) ;
}
//=========================================================================
// 時刻(D)における赤経、赤緯(α(T),δ(T))(度)の天体の方位角(A)計算
//=========================================================================
function calc_AZe(alpha,delta,T,D,lng,lat) {
var theta,t,A,A0,A1;
theta = calc_THETA(lng,T,D); //恒星時
t = theta - alpha; //天体の時角
with(Math) {
A0 =-cos(pai * delta) * sin(pai * t);
A1 = sin(pai * delta) * cos(pai * lat) - cos(pai * delta) * sin(pai * lat) * cos(pai * t);
A = atan(A0 / A1) / pai;
if (A1 > 0 && A < 0) { A = 360 + A; } //分母がプラスのときは-90