Calculating the circulation given specific parameters

9 次查看（过去 30 天）
Jonathon Klepatzki
Jonathon Klepatzki 2023-11-14
评论： Jonathon Klepatzki 2023-11-15
Ran in:
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
end
plot((1:numel(dt)),(1:numel(c)))
  4 个评论
显示 2更早的评论
Jonathon Klepatzki
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?

请先登录，再进行评论。

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

采纳的回答

Matt J
Matt J 2023-11-14
编辑：Matt J 2023-11-14
Ran in:
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更早的评论
Jonathon Klepatzki
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
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?

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Surface and Mesh Plots 的更多信息

标签

产品


版本

R2023a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by