基于Python的可照时数的计算

from math import *

#计算太阳赤纬
t=int(input("请输入t："))
sita=2*pi*t/365#单位弧度
chiwei=0.006918-0.399912*cos(sita)+0.070257*sin(sita)-0.006758*cos(2*sita)+0.000907*sin(2*sita)-0.002697*cos(3*sita)+0.00148*sin(3*sita)
print(chiwei)

b=degrees(chiwei)#b为度数
print(b)

#计算日可照时数
weidu=float(input("请输入纬度："))
g1=radians(45+(weidu-b+0.566667)/2)
s1=sin(g1)
g2=radians(45-(weidu-b-0.566667)/2)
s2=sin(g2)
shang=s2*s1
weiduhu=radians(weidu) #转弧度
c1=cos(weiduhu)
c2=cos(chiwei)
xia=c1*c2
Ta=degrees(asin(sqrt(shang/xia))) #转为度数 单位度
T=4*Ta/15 #单位hour
print(T)

 

仅在最后两行做了修改，感谢原作者和他底下的评论者，代码内容归原作者所有仅供学习参考，欢迎交流提供改进意见。