纬向分场u的的纬度-高度的EOF的第一模态

一路向南

想要做出文献中的纬向分场u的的纬度-高度的EOF的第一模态，但做出来差别较大（左图是文献中的图），不知道为什么，各位大神能帮我看看哪的错误吗？EOF代码是网上下载的写好的。
代码：path='E:\uwnd.nc';
ncinf=netcdf.open(path,'NOWRITE');
lon=ncread(path,'longitude');
lat=ncread(path,'latitude');
lev=ncread(path,'level');
time=ncread(path,'time');
scale=netcdf.getAtt(ncinf,4,'scale_factor');
offset=netcdf.getAtt(ncinf,4,'add_offset');
n1=length(lon);
latnum1=find(lat==10);
latnum2=find(lat==-80);
n2=latnum2-latnum1+1;
n3=length(lev);
n4=length(time);
j=1;
num=n4/12;
for i=1:num
    x=[];
    x=double(netcdf.getVar(ncinf,4,[0 latnum1-1 0 (i-1)*12],[n1 n2 n3 12])).*scale+offset;
    x=squeeze(mean(x,4));
    var(:,:,:,j)=x;
    j=j+1;
end
var=squeeze(mean(var,1));
varmean=squeeze(mean(var,3));
for i=1:n2
    for j=1:n3
        for k=1:num
            var(i,j,k)=var(i,j,k)-varmean(i,j);
        end
    end
end
[V_eof,PC,per,lamd]=EOF(var,1);
latnew=ncread(path,'latitude',[latnum1],[n2]);
h=contourf(latnew,lev,V_eof(:,:,1)');
colorbar;
set(gca,'YDir','reverse');

文献中的

自己画的

一路向南

没人

蔚蓝情怀

可能是原始数据有问题了

更大华

我也不明白，望懂的人帮忙看看

一路向南

蔚蓝情怀 发表于 2015-2-11 08:29
可能是原始数据有问题了

应该没有，欧洲中心下的

蔚蓝情怀

你换个别的资料分析一下 看看结果合理吗 从而判断程序是否有问题

一水天下

把100hpa上的数据去掉试试呢

一路向南

一水天下 发表于 2015-2-17 11:57
把100hpa上的数据去掉试试呢

不行，还是和文献上的图不一样

沙颖凯

"[V_eof,PC,per,lamd]=EOF(var,1);"

你的EOF.m之前用时有问题吗？

一路向南

沙颖凯 发表于 2015-2-20 06:24
"[V_eof,PC,per,lamd]=EOF(var,1);"

你的EOF.m之前用时有问题吗？

没有，用了其它的EOF程序画出来也是这样，不知道拿出了问题