grads,fortran站点资料的转换

LIfighting

正在学习如何用grads和fortran处理站点资料，参考过清风大侠还有书上的程序，但在用Fortran，grads建立站点映射文件时出现错误并还是一直出错找不到问题在哪儿！
其中Fortran程序：
program ex1
implicit none
integer i
real vec(160)
!!!!!!!!!!!!!从文件"r1606.txt"中读入数据!!!!!!!!!
open(1,file='E:\self\5\r1606.txt',status='old')
read(1,*) (vec(i),i=1,160)
close(1)
!!!!!!!!!!!!!将文本记录转换为GrADS所支持的二进制记录!!!!!!!!!!!!
call stntogrd(vec)
end
!!!!!!!!!!!!!stntogrd子程序!!!!!!!!!!!!!!!!!
subroutine stntogrd(x)
real lat(160),lon(160),x(160)
character*8 stid(160)
!!!!!!!!!!!!!从文件"china.dat"中读入各站的站号!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!从文件"china.dat"中读入160站的经纬度!!!!!!!!!!!!!!
open(2,file='E:\self\5\china.dat')
do i=1,160
read(2,*) stid(i),lat(i),lon(i)
enddo
close(2)
!!!!!!!!!!!!!按排序向二进制文件"201006.grd"中写入记录!!!!!!!!!!!!!!!!
open(3,file='E:\self\5\201006.grd',form='binary')
tim=0.0
nlev=1
nflag=1
do  40 i=1,160
write(3) stid(i),lat(i),lon(i)
!# ,tim,nlev,nflag,x(i)
40 continue
!!!!!!!!!!!!!!!在文件最后，给出一个特别记录，表示这个时间组的记录结束!!!!!!!!!!!!!!!
nlev=0
write(3) stid(i-1),lat(i-1),lon(i-1),tim,nlev,nflag
close(3)
return
end


grads编辑的数据描述文件：
dset E:\self\5\201006.grd
dtype station
stnmap E:\self\5\rain.map
undef -999.0
title rain
tdef 1 linear Jun2010 1mo
vars 1
p 0 99 rainfall data
endvars



在用Grads生成映射文件rain.map时出现如下问题：

grads的报错

LIfighting

希望能得到大家的帮忙！！！！{:5_213:}{:5_213:}