站点资料画图生成映射文件是出错

飞天漫雪

program transfer
!将十进制数据转换为二进制，成功生成二进制文件
parameter n=30
integer i,k
character*8 name(n) ! 数据为字符型，最大长度为8，最小为4
open(1,file='d:\lengkongqi\5qu\sname.txt')
do i=1,n
read(1,*) name(i)
print*,name(i)
enddo
close(1)
  call stntogrd(name)
end

  subroutine stntogrd(x)
  character x(30)
  parameter m=30
character*8 st(m)
real lon(m),lat(m)
open(2,file='d:\lengkongqi\5qu\30stl.txt')
do  j=1,m
read(2,*) st(j),lon(j),lat(j)! 台站号、经度、维度
enddo
close(2)
open(3,file='d:\lengkongqi\5qu\30st.grd',form='binary')
open(4,file='d:\lengkongqi\5qu\30st.txt')
time=0.0
nlev=1
nflag=1
do 40 k=1,m
write(3) st(k),lat(k)/100,lon(k)/100,time,nlev,nflag,trim(x(k))
write(4,*) st(k),lat(k)/100,lon(k)/100,time,nlev,nflag,x(k)
40 continue
nlev=0
write(3) st(k-1),lat(k-1),lon(k-1),time,nlev,nflag
close(3)
return
end
ctl  文件
dset d:\lengkongqi\5qu\30st.grd
dtype station
stnmap d:\lengkongqi\5qu\30st.map
undef -999.0
title 30st
tdef 1 linear jun1961 1mo
vars 1
name 0 999 stationname
endvars
问题详见附件，跪求大神帮忙。
主要图片粘不上来只好放在附件

洪hong

CTL文件 不需要你自己写，GRADS书上有，你直接用书上那个就好了，不容易出错。

熊猫支队007

你好，请问问题解决了吗？遇到同样的错误