grads错误求帮忙解答

qinyi0521

上面的错误不理解，我想知道错误在哪里。我的fortran转换程序如下：
program main
implicit none
integer,parameter::m=160,n=60
real lon(160),lat(160),rain(m,n),i,j,aver(m),jup(m,n),time,iflag,ilev
integer stid(160)
character*8::std(160)
open(1,file='d:\data\fujia\id160.txt',status='old')
do i=1,m
read(1,*) stid(i),lon(i),lat(i)
enddo
close(1)
open(2,file='d:\data\fujia\r1607.txt')
do j=1,n
read(2,*) (rain(i,j),i=1,m)
enddo
close(2)
!求六十年的平均值
do i=1,m
aver(i)=0
do j=1,n
aver(i)=aver(i)+rain(i,j)
enddo
aver(i)=aver(i)/60
enddo
!求距平值
do i=1,m
do j=1,n
jup(i,j)=rain(i,j)-aver(i)
enddo
enddo
do i=1,m
std(i)=char(stid(i))
enddo
open(4,file='d:\grads\rain.grd',form='binary')
do i=1,n
iflag=1
ilev=1
time=0.0
do j=1,m
write(4) std(j),lon(j),lat(j),time,ilev,iflag,jup(j,i)
enddo
ilev=0
write(4) std(m),lon(m),lat(m),time,ilev,iflag
enddo
close(4)
end

做个霸气的木头

ctl，gs一并贴出来吧

qinyi0521

ctl文件：
dset d:\grads\rain.grd
dtype station
stnmap d:\grads\rain.map
undef -999.0
title rain
tdef 60 linear jul1951 1yr
vars 1
rain 0 99 rainfall data
endvars

gs文件没写，直接执行就出现我截屏出现的问题了