求相关做t检验寻找预测因子。。一定是出错了。。。

Sunshining

我用850、500、200三层1月份高度场和sst场与华北地区夏季降水求相关进行t检验，希望找到通过检验的地区做为预测因子。。可是画出来相关系数的图真的是太诡异。。。找了一天不懂自己哪里做的不对啊。。。
这是我求相关系数的程序：
do it=1,yr
      avey=avey+hr(it)
enddo
      avey=avey/yr
do j=1,n
do i=1,m
      do it=1,yr
      avex(i,j)=avex(i,j)+hgt200(i,j,it)
      enddo
      avex(i,j)=avex(i,j)/yr
      enddo;enddo
do it=1,yr
      sy=sy+(hr(it)-avey)*(hr(it)-avey)
enddo
do j=1,n
do i=1,m
      do it=1,yr
      sxy(i,j)=sxy(i,j)+(hgt200(i,j,it)-avex(i,j))*(hr(it)-avey)
      sx(i,j)=sx(i,j)
     *+(hgt200(i,j,it)-avex(i,j))*(hgt200(i,j,it)-avex(i,j))
      enddo;enddo;enddo
do j=1,n
do i=1,m
      r200(i,j)=sxy(i,j)/sqrt(sx(i,j)*sy)
enddo;enddo
do j=1,n
do i=1,m
avex(i,j)=0.0
sxy(i,j)=0.0
sx(i,j)=0.0
      enddo;enddo

这是t检验：
do j=1,n
do i=1,m
t200(i,j)=r200(i,j)*sqrt((yr-2)/(1-r200(i,j)**2))
enddo;enddo

ctl文件：
dset   d:\climate\climate\data\r.grd
undef  -9.96921e+36
xdef   144 linear 0  2.5
ydef   73 linear -90 2.5
zdef   1 linear 500 1
tdef   65 linear  00Z01Jan1951 1mo
vars   6
r850   1   99   r850
r500   1   99   r500
r200   1   99   r200
t850   1   99   t850
t500   1   99   t500
t200   1   99   t200
Endvars
gs文件：
'set grads off'
'set grid off'
'set lat -90 90'
'set lon 0 360'
'set gxout shaded'
'set clevs -1.669 -1.295 1.295 1.669'
'set ccols 9 4 0 8 2'
'd t200'
'cbar'
'set gxout contour'
'd r200'
'printim d:\climate\climate\data\r200.png white'


画出来的图惨不忍睹啊。。。