请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 6190|回复: 0

[其他] 功率谱图,求大神帮我看看哪里错了,画出来的图很奇怪

[复制链接]

新浪微博达人勋

发表于 2018-2-22 12:18:31 | 显示全部楼层 |阅读模式

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

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 MARZOLINI 于 2018-2-26 10:19 编辑

画功率谱图错误,不知道错在哪里,求大神帮忙看看,谢谢!
原始资料是中国地区1995-2014年共20年的夏季降水日资料,然后我在中国某一个区域上对这20年夏季日降水求区域权重平均,这步我是通过grads求得。脚本如下:'reinit''open c:\data1\ldobs\summer_precip.ctl'

'set lon 72 135.5'
'set lat 18 53.5'
'set t 1 1840'

'define regavg=aave(precip,lon=107,lon=125,lat=18,lat=30)'

'set fwrite c:\data1\ldobs\summer_regavg1.grd'
'set gxout fwrite'

'set t 1 1840'

'd regavg'
'disable fwrite'
;
下面是summer_regavg1.ctl文件:
DSET C:/data1/ldobs/summer_regavg1.grd
UNDEF -9.99e8

TITLE  China Hourly Merged Precipitation Analysis

*
xdef  1 levels 1
*
ydef  1 levels 1
*
ZDEF   1 levels 1

TDEF 1840 LINEAR 00Z01jun1995 1dy
*
VARS 1
regavg      0 99  CH01   combined analysis (mm/Hour)
ENDVARS

接着我用ncl画功率谱的图,但是图画出来特别的奇怪,找不出来自己错在哪了,希望哪位大神能帮我看看,谢谢。
ncl脚本如下:
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"
begin
  preci1=fbindirread("C:/data1/ldobs/summer_regavg1.grd", 0,1840,"float")
  preci1!0="t"
  sdof_1=specx_anal(preci1,0,2,0.10)
  splt_1=specx_ci(sdof_1,0.05,0.95)
  printVarSummary(splt_1)

  wks=gsn_open_wks("eps","spec")
   r=True
   r@gsnDraw                 = False
   r@gsnFrame                = False  
   r@tiYAxisString = "Variance/frq_interval"
   f = sdof_1@frq
   p   = 1/f                                    ; *highly non-linear*
   p!0 = "f"
   p&f = f
   p@long_name = "period"
   p@units     = "day"
   r@tiXAxisString = "Period (day)"          ; xaxis
   r@tiYAxisString = "Variance/freq"
   plot= gsn_csm_xy(wks,p, splt_1,r)
   draw(wks)
   frame(wks)
end

这张图是功率谱图
而且如果选取中国另一区域进行区域权重平均得出的图与这张图一样。

捕获.PNG
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表