| 16 |
|
YEAR = Calendar year AD (1950-2050) |
| 17 |
|
|
| 18 |
|
Outputs: |
| 19 |
+ |
STIME = Solar time (hours past midnight) |
| 20 |
|
SAZI = Solar altitude (degrees above horizon) |
| 21 |
|
SALT = Solar azimuth (degrees west of south) |
| 22 |
|
} |
| 52 |
|
{ Greenwich mean sidereal time } |
| 53 |
|
gmst = normcir(6.697375 + 0.0657098242*n + UTIME, 24); |
| 54 |
|
lmst = normcir(gmst - LON/15, 24); |
| 55 |
< |
ha = 15*(normcir(lmst - ra/15 + 12, 24) - 12); |
| 55 |
> |
STIME = normcir(lmst - ra/15 + 12, 24); |
| 56 |
> |
ha = 15*(STIME - 12); |
| 57 |
|
{ Final solar angles } |
| 58 |
|
sin_SALT = sin_dec*sin_LAT + cos_dec*cos(LAT*DEG)*cos(ha*DEG); |
| 59 |
|
SALT = asin(sin_SALT) / DEG; |
