求助一个eof分解的问题

Amyhiganbana

求助各位大神
这是我写的eof分解求NAO指数的程序，可是求出的时间序列的值在1到-1之间，但我看大多数文章都在-2到2之间，大神们帮我看看是为什么吗

begin
        fili = "test1.nc"
    f = addfile (fili, "r")
    psl = f->psl(:,{29.68586:81.51833},{-100:40})
    year = ispan(1850, 2005, 1)
    dims = dimsizes(psl)

  optEOF = True      
  optEOF@jopt = 0   ; This is the default; most commonly used; no need to specify.
;;optEOF@jopt = 1   ; **only** if the correlation EOF is desired

  optETS = False

    w = sqrt(cos(0.01745329*psl&lat))
    wp = psl * conform(psl, w, 1)
    copy_VarMeta(psl, wp)

    x=wp(lat|:,lon|:,year|:)
    neof = 3
    eof = eofunc_Wrap(x, neof, optEOF)
    printVarSummary(eof)


    ntim = dims(0)
    sig_pcv = eofunc_north(eof@pcvar, ntim, True)

    eof_ts = eofunc_ts_Wrap(x, eof, optETS)
    print(eof_ts)
    eof_ts = dim_standardize_n(eof_ts, 1, 1)

    nedf = addfile("pcl_winter-2.nc", "c")
    eof_ts!1 = "year"
    eof_ts&year = year

    nedf->NAO = eof_ts(0,:)

end
;;;;;;;;;next ncl;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
begin
        fili = "pcl_winter-2.nc"
    f = addfile (fili, "r")
    NAO = f->NAO
    NAO_mean = runave_Wrap(NAO, 9, 0)

    wks = gsn_open_wks("png","plot-xy-NAO_EOF-2")
  gsn_define_colormap(wks,"rainbow")  

  res           = True              


  res@vpHeightF = 0.35  ; 改变viewport的大小
  res@vpWidthF  = 0.6   ;

  res@trXMinF     = 1850  ; 设定X轴坐标范围
  res@trXMaxF     = 2005
  res@trYMaxF   = 2 ; 设定Y轴坐标范围
  res@trYMinF   = -2
  res@tiXAxisString   = "year"   ; 添加X轴名称   
  res@tiYAxisString = "NAO"    ; 添加Y轴名称

  plot = gsn_csm_xy(wks,NAO&year,NAO_mean,res)   

end

Amyhiganbana

{:eb302:}

Amyhiganbana

{:lxm_26:}

Amyhiganbana

weiyuntao

你做了9点平滑，是否别人的文章中也做了9点平滑，是否你的处理手段和别人文章中一样，得重视！为了验证你的结果是否正确，你可以画出空间pattern来和标准NAO模态做对比，若一致，说明是正确的，指数是否在-2~2，又或是-1~1之间无所谓的吧，主要是看位相以及相对大小！祝好！

Amyhiganbana

谢谢！我的NAO模态空间分布和标准NAO模态是一致的，只是等值线数值不同，那是否说明，我的结果是正确，只是方法和所看文章的方法不一致呢？

风之牧语

dim_standardize_n你标准化了啊

Amyhiganbana

风之牧语 发表于 2017-10-19 14:23
dim_standardize_n你标准化了啊

时间序列不都要标准化吗，我试了没有标准化，那值很大。。

weiyuntao

对的  只是方法不一致而已   只要提取出来了  就可以了   

wawatou83

占位…………