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