Calculating the circulation given specific parameters

Question

Jonathon Klepatzki 2023-11-14

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2047347-calculating-the-circulation-given-specific-parameters

评论： Jonathon Klepatzki 2023-11-15

I am trying to calculate the circulation using the following code. Although the code works, the calculated data doesn't make sense as I get a linear plot. I would also like to create a image plot, if possible. Any help would be greatful!

close all;

clear all;

clc;

h = 200; % height of stack in meters

z = (0:10); % average wind in the vertical

mr = 2.5; % mean mixing ratio

sigmaz = z*mr; %vertical mixing

dt = (24:72)*60*60; %change in time

for i = 1:length(z)

for j = 1:length(dt)

c = exp(-(z+h).^2.)/((2*sigmaz.^2))./dt;

end

plot((1:numel(dt)),(1:numel(c)))

4 个评论
显示 2更早的评论隐藏 2更早的评论

Jonathon Klepatzki 2023-11-15

You were correct, I made an error.

However, perhaps you can help me. I changed the equation to include either:

C = (2./pi).^1/2.*(Q./sigmaz.*meanu.*W.*dt) or C = Q./meanu.*h.*W.*sigmaz.*dt

and the surf plot provided by Matt J doesn't work. Any ideas?

Jonathon Klepatzki 2023-11-15

Never mind, I figured out what I did wrong.

请先登录，再进行评论。

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

Answer 1

Matt J 2023-11-14

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2047347-calculating-the-circulation-given-specific-parameters#answer_1352602

编辑：Matt J 2023-11-14

在 MATLAB Online 中打开

Perhaps this is what you meant. The expression for c looks like it's supposed to be Gaussian, and the double for-loop suggests this is supposed to be a surface plot.

h = 200; % height of stack in meters

z = (0:10); % average wind in the vertical

mr = 2.5; % mean mixing ratio

sigmaz = z*mr; %vertical mixing

dt = (24:72)'*60*60; %change in time

c = exp( -0.5 * ( (z+h)./sigmaz).^2 )./dt;

surf(z,dt,c)

3 个评论
显示 1更早的评论隐藏 1更早的评论

Jonathon Klepatzki 2023-11-14

编辑：Jonathon Klepatzki 2023-11-14

@Matt J it is indeed Gaussian. I narrowed down the equation that spoke of dispersion to only account for the upper-levels insteand of the surface. This is an attempt for me to conduct a study of aerosols in the upper-levels at a given depth in the arctic (Master's degree) when significant snow fall events are taking place.

I like your plots btw! And thank you.

Jonathon Klepatzki 2023-11-15

@Matt J, quick question. I had to change the equation to include the full dispersion (mathematicl error) and noticed that when I changed the variable "mr" to either 1 or 3, it gives me a linear plot once again. Could this be an issue with the dimensions?

请先登录，再进行评论。