20 |
|
{ Default to Tregenza sky } |
21 |
|
MF = 1; |
22 |
|
{ Default axis orientation (+Y == north) } |
23 |
< |
rNx = 0; rNy = 0; rNz = 1; |
23 |
> |
rNx = 0; rNy = 0; rNz = -1; |
24 |
|
Ux = 0; Uy = 1; Uz = 0; |
25 |
|
{ Compute oriented axis angles } |
26 |
|
inc_dz = -Dx*rNx-Dy*rNy-Dz*rNz; |
27 |
|
inc_rx = -Dx*(Uy*rNz-Uz*rNy) - Dy*(Uz*rNx-Ux*rNz) - Dz*(Ux*rNy-Uy*rNx); |
28 |
< |
inc_ry = -Dx*Ux-Dy*Uy-Dz*Uz - inc_dz*(rNx*Ux+rNy*Uy+rNz*Uz); |
28 |
> |
inc_ry = Dx*Ux+Dy*Uy+Dz*Uz + inc_dz*(rNx*Ux+rNy*Uy+rNz*Uz); |
29 |
|
|
30 |
|
r_alt = Asin(inc_dz); |
31 |
|
r_azi = Atan2(inc_rx,inc_ry); |