合成分析的显著性检验问题

推开世界的门

请问一下，各位老师      我将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