多时次探空资料处理问题

hello-lyn

按照清风的帖子，想做一个多时次探空资料的处理，先用一个站点，两个时次，每个时次17层的数据进行尝试。fortran程序如下：
program stn2grd
parameter(lev=17)
character*8 stid
real lat,lon,alt,year,mon,tt(2),p(2,lev),h(2,lev),t(2,lev),d(2,lev),u(2,lev)
open(3,file='d:\1007\stn01.txt')
open(4,file='d:\1007\stn01.grd',form='binary')
do i=1,2
do k=1,lev
read(3,*) stid,lat,lon,alt,year,mon,tt(i),p(i,k),h(i,k),t(i,k),d(i,k),u(i,k)
enddo;enddo
close(3)

do i=1,2
j=1
tim=0.0
nlev=17
flag=1
write(4) stid,lat,lon,tim,nlev,flag,p(i,j),h(i,j),t(i,j),d(i,j),u(i,j) !j=1时，表示写入地面变量
print*,stid,lat,lon,tim,nlev,flag,p(i,j),h(i,j),t(i,j),d(i,j),u(i,j)
do k=2,lev                                                             !写入高空变量
write(4) p(i,k),h(i,k),t(i,k),d(i,k),u(i,k)
enddo
nlev=0
write(4) stid,lat,lon,tim,nlev,flag
enddo
close(4)
end

生成了grd又做站点映射文件
站点ctl如下：
dset d:\1007\stn01.grd
dtype station
stnmap d:\1007\stn01.map
undef 99999
title tankong
tdef 2 linear 01JAN2010 12hr
VARS 5
p 17 99 pressure
z 17 99 height
t 17 99 temperature
d 17 99 dewpoint temperature
u 17 99 u-wind
endvars


运行grads ！ stnmap命令时出现了这样的错误。 不知是不是fortran程序中写入grd时那一步出错，导致跟ctl不匹配。希望大家指教！