爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 7623|回复: 6

[作图] 对海温合成做显著性检验,能不能帮忙看看程序出了什么问题

[复制链接]

新浪微博达人勋

发表于 2015-4-4 20:44:07 | 显示全部楼层 |阅读模式

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

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

x
做10月-2月的海温合成分析的显著性检验,这样出来的结果就没有通过显著性检验的,能不能帮忙看看哪里出问题了,谢谢!
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
filename="data/sst.mnmean.nc"
f=addfile(filename,"float")
year_high=(/1972,1976,1981,1985,1990,2004,2009/)
year_low =(/1971,1975,1980,1984,1989,2003,2008/)
s  = f->sst(0:6,:,:)
s1 = f->sst(0,:,:)
dimx = dimsizes(s1)
nlat = dimx(0)
nlon = dimx(1)
s_high = new((/7,5,nlat,nlon/),float)
s_low  = new((/7,5,nlat,nlon/),float)
ss_high = new((/7,nlat,nlon/),float)
ss_low  = new((/7,nlat,nlon/),float)
do i=0,6
  s_high(i,:,:,:)=f->sst(((year_high(i)-1954)*12+9):((year_high(i)-1954)*12+13),:,:)
  s_low(i,:,:,:) =f->sst(((year_low(i)-1954)*12+9):((year_low(i)-1954)*12+13),:,:)
  ss_high(i,:,:) =dim_avg_n(s_high(i,:,:,:),0)
  ss_low (i,:,:) =dim_avg_n(s_low (i,:,:,:),0)
end do
;-----------------------------------------------------------------------------------
;t-test
;-----------------------------------------------------------------------------------
copy_VarMeta(s,ss_high)
copy_VarMeta(s,ss_low)
sst_high_ave=dim_avg_n_Wrap(ss_high,0)
sst_low_ave =dim_avg_n_Wrap(ss_low ,0)
sst_high_var=dim_variance_n(ss_high,0)
sst_low_var =dim_variance_n(ss_low ,0)
iflag = False
pro = ttest(sst_high_ave,sst_high_var,t1,sst_low_ave,sst_low_var,t2,iflag,False)
copy_VarMeta(s1,pro)
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-4-4 20:46:32 | 显示全部楼层
漏了一个语句:t1=7,t2=7
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-4-4 20:46:49 | 显示全部楼层
漏了一个语句:t1=7,t2=7
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-4-6 16:06:33 | 显示全部楼层
没有人吗?
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2015-4-7 15:31:10 | 显示全部楼层
漏了一个语句:t1=7,t2=7
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-4-8 10:06:56 | 显示全部楼层
不能更赞了~~~~~~~~~~~~~~~~~~~~~~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-4-8 10:07:17 | 显示全部楼层
不能更赞了~~~~~~~~~~~~~~~~~~~~~~
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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