matlab读取JRA-55数据（已解决）

小莹子

JRA-55的数据格式为文件名.日期，是grib格式，依照论坛贴吧博客等经验，安装MinGW64编译器，安装read_grib插件，但是出现问题：
read_grib(gribname,'invent')
错误使用 fileparts
输出参数太多。
出错 read_grib (line 77)
[fpath,fname,fext,fver]=fileparts(gribname);
论坛上说fileparts只能输出3个返回值，可以删掉多余返回值，可是怎么删除那？


换种思路：将grib文件修改为nc文件，也方便自己后续处理。
1）在windows上直接添加后缀，显示文件可能会被破坏
ncFilePath='F:\MATLAB\JRA\anl_surf125.197901.nc';
ncid=netcdf.open(ncFilePath,'NC_NOWRITE');
ncdisp(ncFilePath,'lat');
lat=ncread(ncFilePath,'lat');

错误使用 netcdflib
无法打开文件。文件可能损坏或文件名可能包含无效字符。如果数据源为 OPeNDAP URL，请参阅 "用户指南" 中的 "OPeNDAP 故障排除" 部分。
出错 netcdf.open (line 59)
        [varargout{:}] = netcdflib ( 'open', filename, varargin{1} );

2）Linux上修改后缀为nc,在panoply工具中可以打开文件，但是matlab中依旧无法处理
错误使用 netcdflib
无法打开文件。文件可能损坏或文件名可能包含无效字符。如果数据源为 OPeNDAP URL，请参阅 "用户指南" 中的 "OPeNDAP 故障排除" 部分。
出错 netcdf.open (line 59)
        [varargout{:}] = netcdflib ( 'open', filename, varargin{1} );



解决方法：
在Linux上将grib数据转为nc数据，转化代码为：ncl_convert2nc

代码如下：
clc;        
clear all;
datadir='/data/JRA/anl_surf.nc/';
file=dir([datadir,'*.nc']);
%%
%插值为0.5*0.5分辨率
for i=1:length(file)
ncid=[datadir,file(i).name];
tem=ncread(ncid,'TMP_GDS4_HTGL_S123')-273.15;  %TMP_GDS4_HTGL_S123;PWAT_GDS4_EATM_S123kg m-2  RH_GDS4_HTGL_S123相对湿度；PRES_GDS4_SFC_S123表面压力
lat=ncread(ncid,'g4_lat_0');
lat=double(lat);
lon=ncread(ncid,'g4_lon_1');
lon=double(lon);
% time=ncread(ncid,'time');
Lat=double(lat(78:116,:));
Lon=double(lon(114:188,:));
T=tem(114:188,78:116);
T=double(T);
%插值为统一分辨率
% [X,Y]=meshgrid(Lat,Lon);
% [x,y]=meshgrid(20:0.5:45,70:0.5:105);
% vi=interp2(X,Y,T,x,y);
B(:,:,i)=T;
end
% b=B(:,:,1:480);
Ttb=reshape(B(:,:,13:480),[75 39 12 39]);
T0=reshape(mean(Ttb(:,:,1:12,:),3),[75 39 39]);%年平均
T2=reshape(mean(Ttb(:,:,6:8,:),3),[75 39 39]);%夏季平均
%冬季 1-12----3月-2月
T=reshape(B(:,:,15:482),[75 39 12 39]);
T4=reshape(mean(T(:,:,10:12,:),3),[75 39 39]);%冬季平均 %3月-2月

smean=mean(T0,3); %气候态平均
% %%
% %斜率
x=1980:1:2018;
k=size(T0,1);
for j =1:k
    y=T0(j,:);
    [a m b]=regression(x,y);
    m=m*10;
    M(:,j)=m;
end
    N=M';