读取ERA5 grib数据

clear;
clc; 
path =  'I:\ERA5\';
matlist = dir(strcat(path,'*.mat'));
for i = 1:size(matlist,1)
    matname = matlist(i).name;
    load(matname);
end
value_all = cat(1,value16,value17,value18,value19',value20);
plot(1:length(value_all), value_all(:,1));

%%
% cd('D:\代码储备\nctoolbox-1.1.3');
% setup_nctoolbox;
% cd ('F:\ERA5\');

% %%
% clear;
% clc;
% path =  'F:\ERA5\';
% filelist = dir(strcat(path,'2019*.grib'));
% 
% for i = 1:size(filelist,1)
%     filename = filelist(i).name;
%     ds=ncdataset(filename);
%     ds.variables;
%     data19(i).ssrd = ds.data(ds.variables{8});
%     data19(i).lat = ds.data(ds.variables{12});
%     data19(i).lon = ds.data(ds.variables{13});
% end
% clear i;
% 
% for m = 1:size(data19,2)
%     data19(m).value = squeeze ( nanmean ( nanmean ( data19(m).ssrd,2 ), 3 ) );
% end
% clear m;
% value19 = data19(1).value;
% for m = 2:size(data19,2)
%     value19 = cat (1,value19,data19(m).value);
% end
% value19 = value19./3600;
% value19 (isnan(value19)==1)= [];
% plot(1:length(value19), value19(:,1));
% save('value19.mat','value19');