（提问）站点数据作图出不来图

littleqi

我想画160站各站回归系数等值线图，但grads出不来图.....
这是程序：
program sx6
implicit none
integer,parameter::m=160,n=25
integer i,j
integer::sta(m),t(25)=(/(i,i=1,25)/)
real::f(m,n),jd(m),wd(m),s1(m)=0,s2(m)=0,avef(m)=0,avet=0,b(m)=0
real tim
integer nlev,nflag
!读数据
open(1,file='e:\tongji\shixi6\160zhan-rainfall-summer.txt')
read(1,*)
do i=1,m
read(1,*)sta(i),wd(i),jd(i),(f(i,j),j=1,n)
enddo
!计算各站回归系数b
do i=1,160
  do j=1,25
  avef(i)=avef(i)+f(i,j)/n
  enddo
enddo
do i=1,25
avet=avet+t(i)*1.0/n
enddo
do i=1,160
  do j=1,25
  s1(i)=s1(i)+f(i,j)*t(j)
  s2(i)=s2(i)+t(j)**2
  enddo
enddo
do i=1,m
b(i)=(s1(i)-n*avef(i)*avet)/(s2(i)-n*(avet**2))
write(*,*)b(i)
enddo
!将各站回归系数b写入
open(3,file='e:\tongji\shixi6\rt.grd',form='binary')
tim=0.0
nlev=1
nflag=1
do i=1,160
write(3) sta(i),wd(i),jd(i),tim,nlev,nflag,b(i)
enddo
nlev=0
write(3) sta(i-1),wd(i-1),jd(i-1),tim,nlev,nflag

end

可以运行出结果并生成rt.grd
[img]file:///C:/Users/Administrator.USER-20141106MT/AppData/Roaming/Tencent/Users/1425213625/QQ/WinTemp/RichOle/(IUGQUETYLUZ(9GAXR0Y2[S.png[/img]
但是grads会显示：
[img]file:///C:/Users/Administrator.USER-20141106MT/AppData/Roaming/Tencent/Users/1425213625/QQ/WinTemp/RichOle/XZCR5WM4[CFWFE[$INNSW.png[/img]
这是我的ctl：
dset  e:/tongji/shixi6/rt.grd
dtype station
stnmap e:/tongji/shixi6/rt.map
undef 99999.9
title the 160 station Jan t and ind rc   
tdef 1  linear jan1982 1mo
vars 1
b 0 99 grid data1
endvars
我觉得应该是程序哪里出错了，但实在想不出来，求帮助！！

littleqi

C:\Users\Administrator.USER-20141106MT\Desktop

littleqi

不会发图片