ncread读取NCEP再分析资料疑问

junzhi23

下载了NCEP的再分析资料，但是ncread读取时出了问题。。。不知道是什么缘故。。烦请有经验的小伙伴不吝指导。

南国的佳人

最后一句，从a(1)也就是300开始读，然后多读一层
但是NCEP的资料好像是17层的多，300层的应该没有
楼主可以查一下那套资料里300hPa是第几层，然后对应改一下a(1)那里哦

南国的佳人

b(1),c(1),d(1)那里也有相同的问题

junzhi23

南国的佳人 发表于 2017-2-6 18:23
最后一句，从a(1)也就是300开始读，然后多读一层
但是NCEP的资料好像是17层的多，300层的应该没有
楼主可 ...

C:\Users\bzqxt\Desktop\ncread3.png

level=ncread('rhum.mon.mean.nc','level');
lat=ncread('rhum.mon.mean.nc','lat');
lon=ncread('rhum.mon.mean.nc','lon');
time=ncread('rhum.mon.mean.nc','time');
a=find(level>=300&level<=500);
b=find(lat>=32.5&lat<=35);
c=find(lon>=115&lon<=117.5);
d=find(time>=1569816&time<=1570488);
level_num=length(a);
lat_num=length(b);
lon_num=length(c);
time_num=length(d);
rhum=ncread('rhum.mon.mean.nc','rhum',[a(1) b(1) c(1) d(1)],[level_num lat_num lon_num time_num]);
修改后还是有问题，a(1)的值为6，也就是第六层，应该没问题的啊。。求教
如果把a(1),b(1),c(1),d(1)都修改为6,23,47,2还是同样有问题。。

junzhi23

南国的佳人 发表于 2017-2-6 18:24
b(1),c(1),d(1)那里也有相同的问题

我明白了，ncread中定义数据的起始位置及终止位置时必须同资料的Dimensions顺序，感谢你的建议。

气象万千2011

junzhi23 发表于 2017-2-7 10:00
我明白了，ncread中定义数据的起始位置及终止位置时必须同资料的Dimensions顺序，感谢你的建议。

亲，遇到了相同的问题，但是没懂你是怎么改的啊？

junzhi23

气象万千2011 发表于 2017-3-12 22:28
亲，遇到了相同的问题，但是没懂你是怎么改的啊？

就是如果nc的维度顺序是先纬度，再经度，之后level，最后时间的话，你在用ncread时也要按此顺序。

一夜荒原

a(1) b(1) c(1) d(1)到底什么意思