利用UTC时间和经纬度计算太阳高度角的python代码

以下是使用Python计算太阳高度角的代码示例：

import datetime
import math

# 输入UTC时间和经纬度
utc_time = datetime.datetime.utcnow()
longitude = 120.1536  # 南京市的经度
latitude = 31.9818   # 南京市的纬度

# 计算儒略日
J2000 = datetime.datetime(2000, 1, 1, 12, 0, 0)
delta_JD = utc_time - J2000
JD = 2451545.0 + delta_JD.days + delta_JD.seconds / 86400.0

# 计算世界时
T = (JD - 2451545.0) / 36525.0
UT = utc_time.hour + utc_time.minute / 60.0 + utc_time.second / 3600.0
TU = (UT - longitude / 15.0) / 24.0

# 计算格林威治平恒星时
GMST = 6.697374558 + 0.06570982441908 * JD + 1.00273790935 * UT + 0.000026 * T * T
GMST = GMST % 24

# 计算黄赤交角
epsilon = 23.439 - 0.0000004 * JD

# 计算太阳真黄经
M = (357.5291 + 0.98560028 * (JD - 2451545)) % 360
C = 1.9148 * math.sin(math.radians(M)) + 0.02 * math.sin(math.radians(2*M)) + 0.0003 * math.sin(math.radians(3*M))
lambda_sun = (M + 102.9372 + C + 180) % 360

# 计算太阳直角坐标
alpha_sun = math.atan2(math.cos(math.radians(epsilon))*math.sin(math.radians(lambda_sun)), math.cos(math.radians(lambda_sun)))
delta_sun = math.asin(math.sin(math.radians(epsilon))*math.sin(math.radians(lambda_sun)))

# 计算太阳高度角和方位角
h = math.asin(math.sin(math.radians(latitude))*math.sin(delta_sun) + math.cos(math.radians(latitude))*math.cos(delta_sun)*math.cos(math.radians(GMST*15)))
A = math.atan2(-math.cos(delta_sun)*math.sin(math.radians(GMST*15)), math.sin(delta_sun)*math.cos(math.radians(latitude)) - math.cos(delta_sun)*math.sin(math.radians(latitude))*math.cos(math.radians(GMST*15)))

# 将太阳高度角转换为角度制
h_deg = math.degrees(h)

# 输出结果
print("太阳高度角：", h_deg)

这个代码在计算太阳高度角时，会先计算出当前的世界时和格林威治平恒星时，然后计算黄赤交角和太阳真黄经，最后根据经纬度和时间计算太阳直角坐标和太阳高度角、方位角。这个代码可能存在一些舍入误差，但是可以用于一般的应用。