matlab打开hgt文件,科学网—MATLAB读取ERA和NCEP数据以及写入新的nc文件 - 张乐乐的博文...

function[ value pressure latitude longitude] = get_NCEP_NCAR(year,month,field)

%Field is the type of variable you want : geopotential_height,temperature,u,v,w%Dimension of output value is (time,pressure_level,latitude,longitude)

folder='/where_the_netcdf_are/';

%List the files in selected folder

air_list=dir([folder,'air.',num2str(year,'%04i'),'*']);

hgt_list=dir([folder,'hgt.',num2str(year,'%04i'),'*']);

omega_list=dir([folder,'omega.',num2str(year,'%04i'),'*']);

uwnd_list=dir([folder,'uwnd.',num2str(year,'%04i'),'*']);

vwnd_list=dir([folder,'vwnd.',num2str(year,'%04i'),'*']);

if(strcmp(field,'temperature'))

filename=air_list.name;

filepath=[folder,filename];

ncid=netcdf.