求助：nc文件批量读取

突破的渴望

请问各位大神 我求的是ncep数据的nc文件中特定经纬度（lon>=120 & lon<=290；lat>=-30 & lat<=30）下 高度为第一层 （level=1）的多年平均，请问这个代码错在哪里？ 平均的应该是怎么写呢？写在哪里呢？谢谢大家~~~

1. clc;  %清屏
   
2. clear; %清空
   
3. datadir='C:\Users\xuchu\Desktop\air\'; %指定批量数据所在的文件夹
   
4. filelist=dir([datadir,'*.nc']); %指定批量数据的类型
   
5. a=filelist(1).name; %查看你要读取的文件的编号。filelist(1).name在window下为第一个标号数据
   
6. b=filelist(2).name; %查看你要读取的文件的编号。filelist(2).name在window下为第二个标号数据
   
7. k=length(filelist);
   
8. for s=1:k
   
9.    filename=[datadir,filelist(s).name];
   
10.    lon=ncread(filename,'lon');
    
11.    lat=ncread(filename,'lat');
    
12.    c=find(lon>=120 & lon<=290);
    
13.    d=find(lat>=-30 & lat<=30);
    
14.    lon_num=length(c);
    
15.    lat_num=length(d);
    
16.    level_lim=ncread(filename,'level',[1], [1]);
    
17.    level_num=length(level_lim);
    
18.    ai  = ncread(filename,'air',[level_lim c d],[level_num lon_num lat_num]); %读入变量air
    
19. end

复制代码

突破的渴望

问题已解决~~~