生成站点映射文件时出错，有人提过类似问题，但按其修改无果。多个时次在一个文件里

晚晴山

数据是46年，每年9月份的，因此时次应该是46，每个时次有56个站点。每个时次有5列：站点号，经度，纬度，年份，月份，降水。（见附件）

fortran程序是这样的：
program main
implicit none
integer,parameter::nt=46,nn=56
real lat(nt,nn),lon(nt,nn),rain(nt,nn),tim
integer stan(nt,nn),yr(nt,nn),mon(nt,nn)
integer nlev,nflag,i,k
open(11,file='F:\cy\rain99.txt')
do i=1,nt
do k=1,nn
read(11,*)stan(i,k),lon(i,k),lat(i,k),yr(i,k),mon(i,k),rain(i,k)
enddo

enddo
close(11)

open(12,file='f:\cy\obs.grd',form='binary')
do i=1,nt
tim=0.0
nlev=1
nflag=1
do k=1,nn

write(12) stan(i,k),lat(i,k),lon(i,k),tim,nlev,nflag,rain(i,k)
enddo
nlev=0
write(12)stan(i,k-1),lat(i,k-1),lon(i,k-1),tim,nlev,nflag
enddo
close(12)
end

对应的ctl是这样的：
dset F:\cy\obs.grd
dtype station
stnmap F:\cy\obs.map
undef -999.0
title ob
tdef 46 linear 00Z01JAN1961 1mo
vars 1
rain 0 99
endvars



GrADS错误提示是这样的（见图）：

mofangbao

你参考资料看的不仔细吧，一般都会强调 第一个序号（站号）必须是character*8的才行

river

这是你没有仔细看别人的帖子，按照别人的教程修改没看到点子上去，站号必须是字符型*8

晚晴山

mofangbao 发表于 2013-9-22 08:19

                               
登录/注册后可看大图

你参考资料看的不仔细吧，一般都会强调 第一个序号（站号）必须是character*8的才行

教训得极是

晚晴山

river 发表于 2013-9-22 08:30

                               
登录/注册后可看大图

这是你没有仔细看别人的帖子，按照别人的教程修改没看到点子上去，站号必须是字符型*8

教训得即是