Calculate solar zenith angle

Question

Arvind Gauns 2022-5-16

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1720540-calculate-solar-zenith-angle

评论： Lei 2023-12-24

I want to calculate solar zenith angle for inputting a function. i tried with the available codes but didnt work out for me. Any references/guides for estimating solar zenith angle

2 个评论
显示无隐藏无

Image Analyst 2022-5-16

Probably. Did you try a web search?

TUTORIAL: How to ask a question (on Answers) and get a fast answer

Arvind Gauns 2022-5-17

Yes. I had the methodology and also found some codes but it didnt work on my system

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Sam Chak 2022-5-16

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1720540-calculate-solar-zenith-angle#answer_965090

在 MATLAB Online 中打开

Hi @Arvind Gauns

I did a web search and found the formula. It's a little strange that your browser didn't show references/guides to calculate the solar zenith angle.

Anyhow, I also tried the code and it works:

prompt1 = "What is the value of the local latitude? ";
Phi = input(prompt1)
prompt2 = "What is the value of the current declination of the Sun? ";
delta = input(prompt2)
prompt3 = "What is the value of the hour angle, in the local solar time? ";
h = input(prompt3)
theta_s = acos(sin(Phi)*sin(delta) + cos(Phi)*cos(delta)*cos(h))

2 个评论
显示无隐藏无

Arvind Gauns 2022-5-17

在 MATLAB Online 中打开

with a lot of trial and error the whole night, I found this way to get the solar zenith angle.

Lat         = 52.65;
Long        = 5.6;
height      = 47.5;
f1 = A(i,:);
            f = char(f1);
            yea    = f(1:4); year = str2double(yea);
            mont   = f(6:7); month = str2double(mont);
            do     = f(9:10); dom = str2double(do);
            date = datenum(year,month,dom);
            Doy     = datenum(year,month,dom)-datenum(year-1,12,31);
            hou = f(11:12); hour = str2double(hou);    
            mi = f(14:15); min= str2double(mi);
            se = f(17:18); sec=str2double(se);
            UTC = +1;
            time = ((hour - UTC)/24) + (min/(60*24)) + (sec/(60*60*24));
            cos_sza = cos(calczenithangle(Doy,time,0,0,Long,Lat));