请教grads画相关图的问题

金色牛神

为什么我画南亚高压强度与北半球冬季气温资料做相关图的时候，画出来的是这样的：
资料是这样的：

蒙古高压1951-2009年59年的强度（p.dat）；

北半球气温数据为NCAR/NCEP再分析全球月平均温度（air_win.grd），网格距为2.5*2.5。

请教一下大家我哪里出错了，附上我的gs文件和fortran相关程序
fortran：

program main

   implicit none

   integer i,k,m

   real x(59),y(59,73,144),r(73,144)

   real sum1,sum2(73,144),av1,av2(73,144),f1(73,144),f2,f3(73,144)

open(1,file='d:qihoubianhua\data\p.dat')

open(2,file='d:\qihoubianhua\data\air_win.grd',form='binary')

open(3,file='d:\qihoubianhua\data\r.grd',form='binary')

do i=1,59

read(1,*)x(i)

enddo

do m=1,144

do k=1,73

  do i=1,59

         read(2)y(i,k,m)

  enddo

enddo  

enddo

!计算均值

sum1=0

do i=1,59

sum1=sum1+x(i)

enddo

av1=sum1/59

      do m=1,144

do k=1,73

sum2(k,m)=0

do i=1,59

sum2(k,m)=sum2(k,m)+y(i,k,m)

enddo

av2(k,m)=sum2(k,m)/59

enddo

Enddo

do m=1,144

do k=1,73

       f1(k,m)=0

       do i=1,59

      f1(k,m)=f1(k,m)+(x(i)-av1)*(y(i,k,m)-av2(k,m))

       end do

enddo

enddo

do m=1,144

do k=1,73

f2=0

       f3(k,m)=0

      do i=1,59

       f2=f2+(x(i)-av1)*(x(i)-av1)

      f3(k,m)=f3(k,m)+(y(i,k,m)-av2(k,m))*(y(i,k,m)-av2(k,m))

      end do

      r(k,m)=f1(k,m)/(sqrt(f2*f3(k,m)))

        print *,'R=',r(k,m)

      write(3)r(k,m)

enddo

enddo

      end

ctl程序：

dset D:\qihoubianhua\data\r.grd
title xiangguanxishutu20121368043
undef -9.96921e+36
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 1  linear 0 1
tdef  59 linear  00Z01JAN1951 1yr
vars 1
r 0 99 r
endvars

pearal

是不是缺省值设置错误呀。程序中并没有判断缺省。其实这个时间相关系数可以直接用grads的tcorr函数来算，把每一点的tcor系数写出来，缺省为-9.99*e+08。再画图就方便多啦~

金色牛神

pearal 发表于 2015-6-18 10:54
是不是缺省值设置错误呀。程序中并没有判断缺省。其实这个时间相关系数可以直接用grads的tcorr函数来算，把 ...

好的 谢谢啊 我试一试

金色牛神

pearal 发表于 2015-6-18 10:54
是不是缺省值设置错误呀。程序中并没有判断缺省。其实这个时间相关系数可以直接用grads的tcorr函数来算，把 ...

懂了  谢谢啊