Error reading uo file
1 view (last 30 days)
Show older comments
Error in untitled (line 52)
u1 = ncread(fullfile(filename,'uo'),[1 1 nz 1],[length(lon) length(lat) 1 ntime]);
for yi=1:nyear
ny=ystart+yi-1
filename=fullfile(fp,['June' num2str(ny) '.nc']);
u1 = ncread(fullfile(filename,'uo'),[1 1 nz 1],[length(lon) length(lat) 1 ntime]);
v1 = ncread(fullfile(filename,'uo'),[1 1 nz 1],[length(lon) length(lat) 1 ntime])
% if you want monthly averages, do this
% note this will end up a different dimension if you're extracting
% multiple layers
% might want to predefine these matrices depending on what you want
% https://uk.mathworks.com/help/matlab/ref/mean.html
u2(:,:,yi)=mean(u1,[4]); whos u2
v2(:,:,yi)=mean(v1,[4]);
end
% then take the mean to get the mean over all the years, for one depth
% layer only, remember
u3=mean(u2,[3]); %taking the mean over the third dimension
v3=mean(u2,[3]);
% calculate the magnitude of the vel and plot
uv4=sqrt(u3.^2+v3.^2);
close all, pcolor(LON,LAT,uv4.'), shading flat, colorbar
0 Comments
Answers (1)
VBBV
on 19 Jul 2023
u1 = ncread(filename,'uo',[1 1 nz 1],[length(lon) length(lat) 1 ntime]);
v1 = ncread(filename,'uo',[1 1 nz 1],[length(lon) length(lat) 1 ntime]);
4 Comments
VBBV
on 19 Jul 2023
See the above question / answer which discusses similar issue, I think its related with Version and /or file limit contents
See Also
Categories
Find more on Web Services in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!