用Fortran处理站点数据插值成格点数据后用grads作图，grads报错

myx0919

第一次发帖~请大家多多指教！

这个是很初级的问题，但是我已经钻研了特别久也不知道为什么出不了图，所以就想请大家帮忙看一下！谢谢大家了！

我参考了http://bbs.06climate.com/forum.p ... &extra=page%3D1 这个帖子的做法~
我是想做出eof的图！现在是在用傻办法一步一步地走了
几乎都是一模一样地按照这个上面帖子的流程走了，而且在我进行数据处理之前我是成功地出了图的！
（PS:数据处理指平滑啊，距平啊，筛选某些站点之类的）
但是我数据处理之后，grads就开始报错了！（报错附图在下面，还有跟着这个报错一起出现的图）

给大家看看我的各个程序，希望大家可以帮帮我看看哪里出错了
站点处理成格点（Fortran）：
program egvt
      
real,dimension(120)::lat,lon,ind
character*8 stid(120)
open(1,file='regvt4.txt',status='old')
do i=1,120
read(1,*)lat(i),lon(i),ind(i)
enddo
close(1)
do i=1,120
stid(i)=char(i)
enddo
open(3,file='regvt4.grd',form='binary')
tim=0.0
nlev=1
nflag=1
do i=1,120
write(3) stid(i),lat(i),lon(i),tim,nlev,nflag,ind(i)
enddo
nlev=0
write(3) stid(i-1),lat(i-1),lon(i-1),tim,nlev,nflag
close(3)
write(*,*)ind
end


对应的CTL：
dset            E:\MSDEV98\MyProjects\raineof-picture\regvt4.grd
dtype          station
stnmap       E:\MSDEV98\MyProjects\raineof-picture\regvt.map
undef           -999.0
title              rain
tdef          1     linear       Jan1982    1yr
vars          1
re            0     99       rain
endvars



grid的编写（Fortran）：
Program  main
parameter(nx=31, ny=34)
real     lat(ny), lon(nx)
real     s(nx,ny)
open(1, file= 'E:\MSDEV98\MyProjects\raineof-grid\grid1.grd', form='binary')
lat(1)=20.0
lon(1)=105.0
do j=1,ny-1
   lat(j+1)=lat(j)+1.0
enddo
do i=1,nx-1
   lon(i+1)=lon(i)+1.0
enddo
do i=1,nx
     do j=1,ny
           s(i,j)= 1
     enddo
enddo
write(1) s
end


对应的CTL：
dset E:\MSDEV98\MyProjects\raineof-grid\grid1.grd
title grid data
undef -999.9
xdef 31 linear 105 1
ydef 34 linear 20 1
zdef 1 levels 1
tdef 1 linear Jan1982 1yr
vars 1
g     0  99  grid data
endvars


gs文件：
'reinit'
'open    E:\MSDEV98\MyProjects\raineof-picture\grid1.ctl'
'open    E:\MSDEV98\MyProjects\raineof-picture\binary.ctl'
'set lon 105 135'
'set lat 20 53'
'set t 1'
'set mpdset hires cn'
'define   a=oacres(g,re.2,10,7,4,2,1)'
'set gxout contour'
'd a'
'draw title reof lamnda1'
'printim E:\MSDEV98\MyProjects\raineof-picture\regvt1.png white'
;

希望大家可以帮我解答一下为什么我出不了图，我自己已经好努力好努力地去想！也看了好多次原帖，看教材，反复做了好几次，还是不明白哪里出错了！
先谢谢大家了！