计算赤纬、时角与理论太阳可照时数（附代码）

标签：赤纬 可照 sita rad 时角 lat deg

基于Python的可照时数的计算_可照时数的微机计算程序-CSDN博客

这篇博客帮我初期找到了计算方向，但他用于计算赤纬（θ）和时角（w）的公式来历比较奇怪，没有沿用。

（一）代码

from math import *

def possibleSunshineDuration(N, lat_deg):
    """
    :param N: 日数（N）
    :param lat_deg: 纬度（°）
    :return: T:
    """
    # 计算太阳赤纬
    sita_deg = -23.45 * cos(2 * pi / 365 * (N + 10))  # degrees
    # print(sita_deg)
    sita_rad = radians(sita_deg)  # radians
    # print(sita_rad)

    # 计算可照日数
    lat_rad = radians(lat_deg)
    w_rad = acos(-tan(lat_rad) * tan(sita_rad))
    T = w_rad / pi * 24  # hour
    return T

（二）公式
① 太阳赤纬：

其中N为一年的日数，1月1日的日数为1。
参考：赤纬

② 日出方程式

其中是日出（当数值为负数时）或日落（当数值为正值时）时的时角，以度为单位；
是在地球上观测者的纬度；
 是太阳的赤纬。
参考：日出方程式

 

标签：赤纬,可照,sita,rad,时角,lat,deg
From： https://www.cnblogs.com/RumuRyl/p/18610583