grads作图时出现了问题，希望大家能帮我看一下

irides

160站气温资料，先将1月份的气温资料01.txt转化为grd文件：real vec(160)
open(1,file='d:\160\temp\01.txt',status='old')
read(1,*)(vec(i),i=1,160)
close(1)

call stntogrd(vec)
end

subroutine stntogrd(x)
real lat(160),lon(160),x(160)
character*8 stid(160)

open(2,file='d:\160\china.dat')
do 20 k=1,160
20 read(2,*) lat(k),lon(k)
close(2)

do 2 i=1,160
2 stid(i)=char(i)

open(3,file='d:\160\grd\01.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


生成01.grd文件后编写了ctl文件：
DSET d:\160\grd\01.grd
dtype station
stnmap d:\160\grd\temp.map
undef -999.0
title temp
TDEF 63 LINEAR JUN1950 1yr
vars 1
t 0 0 temps
endvars


然后生成temp.map文件时出了问题，下面是grads提示的错误：


我看了很多相似问题的高人的解答，一一尝试还是没有解决，请问问题出在哪里呢？

f117hqh

貌似我只看到了一个时次的资料，然而ctl里面tdef写的63

irides

是1951—2013这63年里1月份的气温资料，我把63改成1还是有相同的问题。

irides

f117hqh 发表于 2016-4-21 15:10
貌似我只看到了一个时次的资料，然而ctl里面tdef写的63

是1951—2013这63年里1月份的气温资料，我把63改成1还是有相同的问题。

pearal

我看到你的程序里有 2 stid(i)=char(i)  ，可是char(i)在前面没有 被赋值啊~你不妨 把stid(i)输出到屏幕看看？

irides

pearal 发表于 2016-4-21 16:43
我看到你的程序里有 2 stid(i)=char(i)  ，可是char(i)在前面没有 被赋值啊~你不妨 把stid(i)输出到屏幕看 ...

char应该是一个把十进制转换成二进制的函数，stid是存放二进制数据的数组吧好像是

pearal

http://bbs.06climate.com/forum.p ... p;extra=&page=1
这个可以看一下

irides

pearal 发表于 2016-4-21 17:42
http://bbs.06climate.com/forum.php?mod=viewthread&tid=13663&extra=&page=1
这个可以看一下

好的，谢谢大神