NCEP:

clc;         %清屏

clear all;   %清空

%nc批量读取的数据

datadir1='D:\data\ncep_monthly\pressure\';                                                    %指定批量数据所在的文件夹

filelist1=dir([datadir1,'air.mon.mean.nc']);                      %指定批量数据的类型

minlat=-90;maxlat=90;minlon=0;maxlon=357.5;

filename1=[datadir1,filelist1.name];

ncid1=netcdf.open(filename1,'NC_NOWRITE');                                 %打开nc文件                                 

lat=ncread(filename1,'lat');                                         %读入变量lat

lon=ncread(filename1,'lon');                                        %读入变量lon

la=find((lat<=maxlat)&(lat>=minlat));

lo=find((lon<=maxlon)&(lon>=minlon));

latitude=lat(la,:);

longitude=lon(lo,:);

t_1= ncread(filename1,'air');

t(:,:,:,:)=t_1(lo,la,1:12,:);

netcdf.close(ncid1);                                                      %关闭nc文件

 

fileID1=fopen('airtemp.bin','w');

fwrite(fileID1,t,'single');

fclose(fileID1);

 

 

ECMWF:

clc;         %清屏

clear all;   %清空datadir1='D:\sunxy\data\1979-201805_monthly_nc\1\';        %指定批量数据所在的文件夹

filelist1=dir([datadir1,'sounding*.nc']);                      %指定批量数据的类型

k1=length(filelist1);

% 读取数据

minlat=-90;maxlat=90;minlon=0;maxlon=360;

n=1;

for s=1:k1

  filename=[datadir1,filelist1(s).name];

  ncid=netcdf.open(filename,'NC_NOWRITE');                                 %打开nc文件         

 

  lat=ncread(filename,'latitude');                                         %读入变量lat

  lon=ncread(filename,'longitude');                                        %读入变量lon

  time=ncread(filename,'time');

  l=length(time);

  la=find((lat<=maxlat)&(lat>=minlat));

  lo=find((lon<=maxlon)&(lon>=minlon));

  latitude=lat(la,:);

  longitude=lon(lo,:);

 

w_1= ncread(filename,'w');

w(:,:,:,n:n+l-1)=w_1(lo,la,37:-1:11,:);

  netcdf.close(ncid);                                                      %关闭nc文件

  n=n+l

  t_1=[];

end;

 

fileID5=fopen('D:\sunxy\study\Q1Q2\Fortran热源\glo\1\omega.bin','w');

fwrite(fileID5,w,'single');

fclose(fileID5)

 

注意：用MATLAB读取ECMWF文件，气压层会倒置。