想利用站点资料画图，但不知为什么grads报错？

飞行＊想念

我按照马红云的书练习使用站点资料画图，已经将站点资料处理为了二进制文件，但建立站点映射文件时grads总是报错，我不知是什么原因。还请各位指教。我现在把所有内容都传上来，大家帮忙看看是哪的问题。以下是站点资料处理为二进制文件的fortran程序。
real vec(16)
open(1,file='d:\2\r16.txt',status='old')
read(1,*)(vec(i),i=1,16)
close(1)
call stntogrd(vec)
end
subroutine stntogrd(x)
real lat(16),lon(16),x(16)
character * 8 stid(16)
open(2,file='d:\2\china.dat')
do 20 k=1,16
20 read(2,*) lat(k),lon(k)
close(2)
do 2 i=1,16
2 stid(i)=char(i)
open(3,file='d:\2\201006.grd',form='binary')
tim=0.0
nlev=1
nflag=1
do 40 i=1,16
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显示的截屏

然后身为菜鸟的我就不知道该怎么处理了，求大神们指教！

飞行＊想念

补充：grads提示说二进制文件有问题，但不知是哪有问题，怎样解决。本人才学grads不久，还十分不熟悉，还望大家指教，我一定好好学习！

飞行＊想念

再补充一下我的ctl文件

做个霸气的木头

你的程序没有问题，我用你的程序就正常映射地图了，是不是你的数据文件格式或者你文件没有放对位置。好好检查

mofangbao

建议你在我们论坛找一下 grads站点作图详解这篇帖子 照着做一遍就会了

le_soleil

会不会的你的CHINA.DAT文件已经包含了站点的经纬度资料，我记得我老师给我的站点程序打开CHINA.DAT时顺带把经纬度也输入了，你试一试呢

nuister

映射文件出来了吗？

river

程序没问题，ctl也是对的。

飞行＊想念

做个霸气的木头 发表于 2013-5-25 17:04

                               
登录/注册后可看大图

你的程序没有问题，我用你的程序就正常映射地图了，是不是你的数据文件格式或者你文件没有放对位置。好好检 ...

数据文件格式应该是对的呀，文件也应该放对位置了的。

飞行＊想念

le_soleil 发表于 2013-5-25 18:46

                               
登录/注册后可看大图

会不会的你的CHINA.DAT文件已经包含了站点的经纬度资料，我记得我老师给我的站点程序打开CHINA.DAT时顺带把 ...

china.dat是我自己建的文件，不会已经包含了的