读取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;