- 积分
- 614
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-7-27
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
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) |
|