#!/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))
- 1
- 2
- 3
- 4
- 5
- 6
前往页