基于Python的日照时数转太阳辐射计算

#!/usr/bin/env python # -*- encoding: utf-8 -*- ''' @File : 1_Solar_rad_conversion.py @Time : 2022/03/21 18:27:14 @Author : Qi Cheng Yi Shao @Version : 1.0 @Contact : 1441210996@qq.com @License : (C)Copyright 2021-2024, Mutual Aid Commune @Desc : None ''' # here put the import lib from cmath import tan import numpy as np import math Lat=34.299722 #该地区的纬度 ARC=np.pi*Lat/180 TANF=math.tan(ARC) AS=0.25 BS=0.5 def Srconversion(n,julian_day): ''' :param n: 日照时数,单位为小时 :param julian_day: 当前日期在该年度的第几天,如:20210101,则julian_day为1 :return: 返回值为太阳辐射,单位: MJ/m^2/day ''' delta=0.4209*math.sin(0.0172*julian_day-1.39) dr=1+0.033*math.cos(0.0172*julian_day) tanDelta=math.tan(delta) J=-TANF*tanDelta Ws=math.acos(J) L=Ws*math.sin(ARC)*math.sin(delta) M=math.cos(ARC)*math.cos(delta)*math.sin(Ws) Ra=37.6*dr*(L+M) N=24/np.pi*Ws Rs=(AS+BS*(n/N))*Ra return round(Rs,1) if __name__=="__main__": print(Srconversion(7,1))