NCL绘图示例（三）小波图.PDFVIP

NCL 绘图示例（三）：小波图 施宁 （南京信息工程大学大气科学学院） begin f = addfile(./ENSO-index.nc, r) ensoi = f-ensoi time = ensoiyear N = dimsizes(ensoi) ;; 小波计算 mother = 0 ; 母小波类型，通常为 0，即Morlet小波。其余两中被分 别为1，Paul小波和2 ，DOG (derivative of Gaussian)小波 dt = 1 ; 数组中数值之间的时间间隔，通常为 1。本例中表示间 隔1 年。 param = -1 ; 母小波参数。 如果param 0 ，则使用默认数值，即采 用Morlet小波时为6；Paul小波为4;DOG小波为2 s0 = dt ; Morlet小波s0 = dt ; Paul小波s0 = dt/4 dj = 0.25 ; 常用设定 jtot = 1+floattointeger(((log10(N*dt/s0))/dj)/log10(2.)) ; 常用设定 npad = N ; 常用设定 nadof = 0 ; 常用设定 noise = 1 ; 常用设定,h 红噪声检验 siglvl = .05 ; 置信度水平 isigtest= 0 ; 采用chi-square 检验；若为 1 则是对全部波谱进行时间平 均检验 w = wavelet(ensoi,mother,dt,param,s0,dj,jtot,npad,noise,isigtest,siglvl,nadof) ;************************************ power = onedtond(w@power,(/jtot,N/)) ; 功率谱 power!0 = period ; Y axis powerperiod = w@period power!1 = time ; X axis powertime = time power@long_name = Power Spectrum power@units = 1/unit-freq ;计算显著性 ( = 1 则显著) SIG = power ; 复制元数据 SIG = power/conform (power,w@signif,0) SIG@long_name = Significance SIG@units = ;************************************************* wks = gsn_open_wks(eps,plot-enso-wavelet) gsn_define_colormap(wks,BlAqGrYeOrReVi200) YLValues = (/1,2,4,8,16/) YLLabels = (/1,2,4,8,16/) res = True res@gsnDraw = False res@gsnFrame = False res@gsnRightString = res@gsnLeftString = res@trYReverse = True ; 倒置 y-axis res@tmYL