- 积分
- 774
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2020-11-5
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
请问一下,各位老师 我将2020年7-9月份台风发生时刻的垂直风切变做一个合成分析图,想要进行显著性检验,下面是我根据气象家园的例子写的代码,
其中a是所有台风发生时刻的垂直风切变,b为2020年7-9月的垂直风切变,输出的结果如图1所示,蓝色是通过显著性检验的部分,黄色为未通过显著性的部分,输出结果是错误的。
请问一下,是我的方法选择有问题,还是其它的地方有问题,请各位老师帮忙看一下
do i=0,200
do j=0,320
avga(i,j) = avg (a(:,i,j))
avgb(i,j) = avg (b(:,i,j))
vara(i,j) = variance (a(:,i,j))
varb(i,j) = variance (b(:,i,j))
na(i,j) = dimsizes (a(:,i,j))
nb(i,j) = dimsizes (b(:,i,j))
df(i,j) = na(i,j)+nb(i,j)-2
varavarb_ts(i,j) = (na(i,j)-1)*vara(i,j)+(nb(i,j)-1)*varb(i,j)
t_ts(i,j) = (avga(i,j)-avgb(i,j))/(sqrt(varavarb_ts(i,j)/df(i,j))*sqrt(1.0/na(i,j)+1.0/nb(i,j)))
s(i,j) = student_t(t_ts(i,j),df(i,j))
end do
end do
|
-
图1
|