请写一个python程序，计算太阳高度角

以下是一个计算太阳高度角的Python程序：

import math
import datetime

# 经度和纬度（北京）
lat = 39.9
lng = 116.4

# 当前时间
now = datetime.datetime.now()

# 计算儒略日
d = now.toordinal() - datetime.datetime(2000, 1, 1).toordinal() + 366 + (now - datetime.datetime(now.year, 1, 1)).total_seconds() / (24 * 60 * 60)

# 计算真太阳时角
t = (d - 1) / 36525
l0 = math.radians(280.46645 + 36000.76983 * t + 0.0003032 * t ** 2)
m = math.radians(357.52910 + 35999.05030 * t - 0.0001559 * t ** 2 - 0.00000048 * t ** 3)
e = 0.016708617 - 0.000042037 * t - 0.0000001236 * t ** 2
c = (1.914600 - 0.004817 * t - 0.000014 * t ** 2) * math.sin(m) + (0.019993 - 0.000101 * t) * math.sin(2 * m) + 0.000290 * math.sin(3 * m)
theta = l0 + c
nu = m + c
alpha = math.atan2(math.sin(theta) * math.cos(0.00005 * math.pi / 180), math.cos(theta))
delta = math.asin(math.sin(0.00005 * math.pi / 180) * math.sin(alpha))
H = math.radians((now.hour + now.minute / 60 + now.second / 3600 + now.microsecond / 3600000) * 15 - lng - math.degrees(alpha))

# 计算太阳高度角
sin_h = math.sin(math.radians(lat)) * math.sin(delta) + math.cos(math.radians(lat)) * math.cos(delta) * math.cos(H)
h = math.degrees(math.asin(sin_h))

print("太阳高度角：", h)

这个程序首先获取了经度和纬度（这里以北京为例），并获取当前时间。然后计算了儒略日，以及一系列与日地运动相关的参数。最后计算了太阳的高度角，并将其输出。