请问有人绘制过IPO或PDO的图吗

旺仔女孩 · 发表于 2020-10-26 21:16:40

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

求助各位大佬，我想对海温数据进行EOF分解，文献里是按照画IPO的图来处理数据的，然后我就对海温数据求了距平，标准化，面积权重，然后求年平均消除年际变化，但是我的数据处理好像存在问题，画出的图总是和文献中的不一样，空间模态的量级差了大概十倍，解释方差我的是5.5%，文献中是9.5%
小白第一次发帖，如果发帖内容不符合规范，请见谅

begin
  f = addfile("b.e11.BRCP85C5CNBDRD.f09_g16.033.cam.h0.TS.200601-208012.nc","r")
  TS= f->TS(:,:,:)
  printVarSummary(TS)
  ;距平
  TS_avg = clmMonTLL(TS)
  printVarSummary(TS_avg)
  TS_anom= calcMonAnomTLL(TS, TS_avg)
  printVarSummary(TS_anom)

  ;annual mean
  TS_annual=month_to_annual(TS_anom,1)
  printVarSummary(TS_annual)

  ;标准化
  TS_annual= dim_standardize_n_Wrap(TS_annual,1,0)

  ;权重
  w = sqrt(cos(0.01745329*TS_annual&lat))
  wu = TS_annual*conform(TS_annual, w, 1)
  copy_VarCoords(TS_annual, wu)
  printVarSummary(wu)

  x = wu(lat|:,lon|:,year|:)
  neof = 2
  eof = eofunc_Wrap( x,neof, False)
  eof_ts = eofunc_ts_Wrap( x, eof, False)
  eof_ts = dim_standardize_n( eof_ts, 1, 1)