主题:站点txt转换站点二进制文件
文件执行是没有问题,但是在grads中map文件不能生成,求前辈指导啊,
创建二进制站点程序:.grd文件
program main
real vec(33,10)
real lat(33),lon(33)
character*8 stid(33)
open(11,file='outtevt.txt')
do i=1,33
read(11,'(10f5.2)') (vec(i,j),j=1,10)
!write(*,'(10f5.2)') (vec(i,j),j=1,10)
enddo
close(11)
open(12,file='station.txt')
do k=1,33
read(12,'(a5,2x,f5.2,2x,f6.2)')stid(k),lat(k),lon(k)
enddo
close(12)
open (13,file='outtevt.grd',form='binary')
do j=1,10
tim=0.0
nlev=1
nflag=1
do i=1,33
write(13) stid(i),lat(i),lon(i),tim,nlev,nflag,vec(i,j)
enddo
nlev=0
write(13) stid(i-1),lat(i-1),lon(i-1),tim,nlev,nflag
enddo
close(13)
end
ctl文件:
dset c:\msdev\projects\zhandianerjinzhi\outtevt.grd
dtype station
stnmap c:\msdev\projects\zhandianerjinzhi\outtevt.map
undef 32767
title daily v
tdef 1 linear Jan1961 12mo
vars 10
v1 0 99 tdpoint data
v2 0 99 tdpoint data
v3 0 99 tdpoint data
v4 0 99 tdpoint data
v5 0 99 tdpoint data
v6 0 99 tdpoint data
v7 0 99 tdpoint data
v8 0 99 tdpoint data
v9 0 99 tdpoint data
v10 0 99 tdpoint data
endvars
然后执行!stnmap
c:\msdev\projects\zhandianerjinzhi\outtevt.ctl
然后提示有错误,invalid station dhr in............