how to draw a section of current at a given latitude according to longitude and depth?

Section of current at a single latitude.i have a NetCDF file with current components as u and v.

[KSSV]

What data you have? Give us more details.

Meanwhile read about slice

[Khalifa ababacar Ndoye]

NetCDF file with current components u and v then the latitudes and longitudes and depth

Answer 1

[Vinay]

Hii Khalifa,

The target latitude where current section to be defined needs to be compared with the latitude in the dataset and closet latitude has to be identified.Extract the U and V components of the current at this latitude and compute the magnitude. Use the “pcolor” function to plot the current magnitude, with longitude represented on the x-axis and depth on the y-axis.

Kindly refer to the following documentations for “pcolor” function:

https://www.mathworks.com/help/matlab/ref/pcolor.html [AO4]

I hope this helps!

%Custom dataset
longitudes = linspace(-180, 180, 50); 
latitudes = linspace(-90, 90, 30);   
depths = linspace(0, 5000, 20);       
% Generate synthetic current data
[lon_grid, lat_grid, depth_grid] = ndgrid(longitudes, latitudes, depths);

U = 0.1 * sin(2 * pi * lon_grid / 360) .* cos(2 * pi * lat_grid / 180) .* exp(-depth_grid / 2000);
V = 0.1 * cos(2 * pi * lon_grid / 360) .* sin(2 * pi * lat_grid / 180) .* exp(-depth_grid / 2000);
% Save the dataset
save('custom_current_data.mat', 'U', 'V', 'longitudes', 'latitudes', 'depths');

Code for plotting the current section

% Load the dataset
load('custom_current_data.mat'); % Replace with your actual data file

% Define the target latitude
target_lat = 30; 

% Find the index of the closest latitude
[~, lat_idx] = min(abs(latitudes - target_lat));

% Extract the current data at the given latitude
U_section = squeeze(U(:, lat_idx, :));
V_section = squeeze(V(:, lat_idx, :));

% current magnitude
current_magnitude = sqrt(U_section.^2 + V_section.^2);

% Plot the section
figure;
pcolor(longitudes, depths, current_magnitude'); 
shading interp;
colorbar;
xlabel('Longitude');
ylabel('Depth');
title(['Current Magnitude at Latitude ', num2str(target_lat)]);
set(gca, 'YDir', 'reverse'); % Reverse the Y-axis to have depth increasing downwards