读取EC风速
ds=ncdataset(filename);
ds.variables;
lat = ds.data(ds.variables{2});
% lat_find = find ( lat <=45 & lat >= 30 );
% lat = lat ( lat_find );
lon = ds.data(ds.variables{3});
% lon_find = find ( lon <=135 & lon >= 125 );
% lon = lon ( lon_find );
for i = 1:size(lon,1)
Lat (:, i) = lat;
end
clear i;
for i = 1:size(lat,1)
Lon (i, :) = lon';
end
clear i lat lon;
gh = ds.data(ds.variables{1});
gh = shiftdim ( double ( gh ) );
clear ds filelist filename i;
cd I:\EC\grib_dir\700u_24h\;
filelist = dir ( strcat ('I:\EC\grib_dir\700u_24h\','*.grib') );
filelist_order = filelist(60:length(filelist));
filelist_order ( (length(filelist_order)+1) : length(filelist) ) = filelist(1:59) ;
filename = filelist_order(k).name;