NCL处理nc就是会更方便，但是避免不了一些用不了NCL的情况

写了一个简单用matlab提取nc文件站点数据的脚本，通过判断经纬度与格点经纬度距离，选取最近格点值。

后面也附上了提取wrfout中站点数据的NCL脚本，如果是再分析资料的话，wrf系列函数不知道是否可以，可能还是需要通过距离判断。

clc;         %清屏

clear all;   %清空

%%

%nc批量读取的数据

datadir1='C:\Users\sunxy\Desktop\';                                                    %指定批量数据所在的文件夹

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

%%

%读取数据

lat_station=32.5;lon_station=90.8;

  filename1=[datadir1,filelist1.name];

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

  

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

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

  [lat_2D lon_2D]=meshgrid(lat,lon);

  distance=sqrt((lat_2D-lat_station).^2+(lon_2D-lon_station).^2);          %算出站点与每个格点距离

  [ilon ilat]=find(distance==min(min(distance)));                          %寻找距离最近格点

  lat_back=lat(ilat);                                                      %输出经纬度检查是否为希望读取站点

  lon_back=lon(ilon,:);

  blh= ncread(filename1,'blh');

  blh_station=blh(ilon,ilat,:);

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

%   

% % ncdisp('C:\Users\sunxy\Desktop\200901.nc','/','full')