matlab數據處理-批量處理nc文件

clc;clear;
datadir = 'F:\GLEAM_E\'; %指定批量數據所在的文件夾
filelist = dir([datadir,'*.nc']);       %列出所有滿足指定類型的文件
k=length(filelist);  %% 統計文件個數
res=[];


for i = 1:k %依次讀取並處理
    [i k]
    %% 批量讀取NC文件
    ncFilePath = ['F:\GLEAM_E\',filelist(i).name]; %設定NC路徑    
    year=str2num(ncFilePath(14:17));  %% 根據path提取年份

    %% 讀取變量值
%     ncdisp('F:\GLEAM_EP\Ep_1980_GLEAM_v3.3a.nc')  %% 顯示數據組成
    lon=ncread(ncFilePath,'lon');          %讀取經度信息（范圍、精度）
    lat=ncread(ncFilePath,'lat');          %讀取維度信息
    dt=ncread(ncFilePath,'E');     %獲取蒸發變量數據

    %% 提取JJA數據
    st=sum(eomday(year, 1:6-1))+1;
    en=sum(eomday(year, 1:8));
    data_JJA=dt(:,:,st:en);
    
    %% 多年份組合
    res=cat(3,res,data_JJA);
   
end