NCL小波画图遇到的问题，求解决。

hxyj

最近一直用NCL小波分析数据，并画图。但遇到一个问题，就是画出的图的左上角与右上角总是不对称。需要说明的是，如果把 plot2= gsn_csm_contour(wks,power({0:25},:),resL) 里面的{0:25}改为{0:22}就不会出现左右不对称的问题，但这样显示的纵坐标太短。
下面附上我的脚本，数据和图，能否麻烦老师给予解答。不胜感激。祝各位端午节快乐！

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/skewt_func.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/wind_rose.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/run_cor.ncl"

begin

   f3                  =asciiread("D:\paper\20140525\EOF\time\pc1.txt",63,"float")
   nin3                =f3

   ninox               =nin3
   ninox!0             ="time"
   ninox&time          =ispan(1951,2013,1)
;...................................................................wavelet
  y                                        = ninox
  N                                        = dimsizes(y)  ;60
  m                                        = 0
  dt                                       = 1
  param                                    = 6
  s0                                       = 1.5*dt
  dj                                       = 0.125
  jtot                                     = 1 + floattointeger(((log10(N*dt/s0))/dj)/log10(2.)) ;20
  npad                                     = 1024
  noise                                    = 1
  isigtest                                 = 1
  siglvl                                   = 0.1
  nadof                                    = 0
  wave                                     = wavelet(y,m,dt,param,s0,dj,jtot,npad,noise,isigtest,siglvl,nadof)
  power                                    = onedtond (wave@power, (/jtot,N/) )
  phase                                    = onedtond (wave@phase, (/jtot,N/) )
  power!0                                  = "period"                        ; Y axis
  power&period                             = wave@period
  power!1                                  = "time"                          ; X axis
  power&time                               = y&time
  power@long_name                          = " "
  gws                                      = wave@gws
  SIG                                      = power           ; transfer metadata
  SIG                                      = power/conform(power,wave@signif,0)
  SIG@long_name                            = " "
;************************************************create plot
  wks= gsn_open_wks("x11","wavelet")  
  gsn_define_colormap(wks,"MPL_GnBu")

  resL                                 = True
  resL@gsnDraw                         = False
  resL@gsnFrame                        = False   
  resL@trYReverse                      = True                  ; reverse y-axis
  resL@tmXBMode                        = "Explicit"
  resL@tmXBValues                      = ispan(1960,2010,10)
  resL@tmXBLabels                      = (/"1960","1970","1980","1990","2000","2010"/)

  plot2                                = gsn_csm_contour(wks,power({0:25},:),resL)
  res1                                = True
  plot2                                = ShadeCOI(wks,plot2,wave,y&time,res1)
  draw(plot2)
  frame(wks)
end

longlivehj

从shadecoi的ncl源码来来看，的确是会出现左右不对称的情况。但为了达到你提出的要求，我修改了源码。至于这样修改对与否，因为我不太懂小波这块，还请楼主自己把握。防止破坏ncl原始代码，使用的时候将load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"中的路径换为下面这个我修改过的文件实际存放的位置。

shea_util.ncl (156.69 KB, 下载次数: 201)

2014-5-31 12:42 上传

点击文件名下载附件

novelbean

好东西，学习下

hxyj

longlivehj 发表于 2014-5-31 12:45
从shadecoi的ncl源码来来看，的确是会出现左右不对称的情况。但为了达到你提出的要求，我修改了源码。至于 ...

非常感谢您的回复与解决。

Cash

很漂亮，还有显著性检验。{:eb502:}{:eb502:}{:eb502:}{:eb502:}

盈星

你好，我想问一下你这个边界影响区域是怎么画的？对这个函数不太懂

ㄣ韩2723邹ㄜ

一起学习

lilyare_2006

学习了学习了！

cenhui99

好好学习{:eb513:}

yuqianhao0315

简单易行的程序

niuda

好东西，学习下