Algorithm 1. Set time step (increment) 2. set max number of time steps 3. start with ground values of temperature at x=0, y=0, z=0, theta=0, and a particular value of fi 4. co
1 view (last 30 days)
Show older comments
% Corrected Code
% Initialize variables height = F2007120317S7627.Height; wind_speed = F2007120317S7627.WS; wind_direction = F2007120317S7627.WD; temperature = F2007120317S7627.Temp0; humidity = F2007120317S7627.HumiS; min_humid = max(humidity); theta_limit = 45; time_limit = 100; phi_limit = 90; delta_vx = zeros(101, 1); delta_vy = zeros(101, 1); delta_x = zeros(theta_limit, phi_limit); delta_y = zeros(theta_limit, phi_limit); delta_z = zeros(theta_limit, phi_limit); delta_theta = zeros(theta_limit, phi_limit); x = zeros(theta_limit, phi_limit); y = zeros(theta_limit, phi_limit); z = zeros(time_limit, theta_limit); % Fixed dimension delta_t = 0.1; theta_end = zeros(theta_limit + 1, phi_limit); phi = linspace(0, 2 * pi, phi_limit); theta = linspace(0, pi / 2, theta_limit); temp1 = 26.3; delta_t = 1; vx1 = 0; vy1 = 0; alti_matrix = []; vx = zeros(time_limit, theta_limit); vy = zeros(time_limit, theta_limit); cs = zeros(time_limit, theta_limit); z = zeros(time_limit, theta_limit); % Fixed dimension cs(1, :) = 331;
for i = 1:theta_limit delta_z(1, i) = z(1, i) + (cs(1, i) * cos(theta(1, i))) * delta_t; z(1, i) = delta_z(1, i); year = z(1, i); alti_matrix = [alti_matrix, year]; yr = linspace(0, 180, time_limit); altitude = 0:z(1, i):z(1, i) * 4; % alt(k,i) = altitude(1,k); temp1 = spline(height, temperature, altitude); temp_inv = temp1.'; speed1 = interp1(height, wind_speed, altitude); speed_inv = speed1.'; dir1 = interp1(height, wind_direction, altitude); dir_inv = dir1.'; vx(1, i) = speed_inv(1, 1) * cos(deg2rad(dir_inv(1, 1))); vy(1, i) = speed_inv(1, 1) * sin(deg2rad(dir_inv(1, 1))); % delta_vx(1,i) = vx(1,i) - vx(1,i); % delta_vy(1,i) = vy(1,1) - vy(1,i); humd1 = interp1(height, humidity, altitude); humd_inv = humd1.'; q(1, i) = humd_inv(1, 1) / min_humid; temp_kelvin(1, i) = (273.15 + temp_inv(1, 1)) * (1 + 0.608 * q(1, i)); cs(1, i) = 20.047 * sqrt(temp_kelvin(1, i)); % delta_cs(1,i) = cs(k+1,i) - cs(k,i); end
0 Comments
Answers (1)
Walter Roberson
on 25 Sep 2023
Edited: Walter Roberson
on 25 Sep 2023
Remember that comments extend to the end of the physical line, not to the next semi-colon, so your comment about "Initialize variables height" is also commenting out the setting of the wind speed and height and so on that is on the same line.
We recommend against packing so much code onto each physical line: it is difficult for people to read and debug when you do that. You should break it up into separate lines such as
% Initialize variables
height = F2007120317S7627.Height;
wind_speed = F2007120317S7627.WS;
0 Comments
See Also
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!