NCL做SSTA的EOF分解

Hurricane_Hu

数据：1961-2000年DJF的SSTA数据（1996的DJF平均场，1997DJF平均场....以此类推）
           数据原始来源：SODA各年DJF的平均场
不知道为什么出现这个情况？ 北纬30度左右看起来是叠加了一副图上去的样子，但是实际上并没有plot...
把分析区域改成20N-60N也是同样的情况..
是所选的区域不对还是我的数据有问题？谢谢！！
+++++++++++++++++++++++++++++++++++++++++++++
NCL代码：
latS=-40
latN=40
lonL=120
lonR=280

neval=3
optEOF=True
optEOF@jopt=0
optETS = False

; transfer the surface ssta data to "data"
data=temp({lat|latS:latN},{lon|lonL:lonR},time|:)

data!0="time"
data!1="lat"
data!2="lon"

ev=eofunc_Wrap(data,neval,optEOF)
ts=eofunc_ts_Wrap(data,ev,optETS)

printVarSummary(ev)
printVarSummary(ts)

wks=gsn_open_wks("eps","EOF")
res=True
res@cnFillOn=True
res@cnLinesOn=False
res@gsnSpreadColors = True
res@mpCenterLonF=180
res@mpMinLatF=latS
res@mpMaxLatF=latN
res@mpMinLonF=lonL
res@mpMaxLonF=lonR
res@pmTickMarkDisplayMode = "Always"
; important for select regions
res@gsnAddCyclic         = False

symMinMaxPlt(ev, 16, False, res)

plot=gsn_csm_contour_map_ce(wks,ev(0,:,:),res)


year=ispan(1961,2000,1)
wks_ts=gsn_open_wks("eps","ts")
plot2=gsn_csm_xy(wks_ts,year,ts(0,:),False)

niuda

Hurricane_Hu 发表于 2015-6-7 22:56
单独画非常正常啊。。。哭

单独画正常的话，是不是EOF过程及前处理的问题，你可以换一个数据用相同的代码处理后看看出图效果看看

niuda

感觉是不是数据处理的问题，没有EOF前单独画数据看看有没有问题？

Hurricane_Hu

niuda 发表于 2015-6-7 22:55
感觉是不是数据处理的问题，没有EOF前单独画数据看看有没有问题？

单独画非常正常啊。。。哭

当时年少春衫薄

楼主请问你那个time维度的数值是年份吗

zongqing123

解决了，楼主

zongqing123

谢谢分享，找到了

zongqing123

请问楼主是如何解决的？能不能分享一下，谢谢

大头袁

谢谢楼主！

sxtylh

很可能是数据预处理的问题，请楼主详查