关于站点数据stnmap命令的使用

特别提示：!stnmap -i d:/eof/td.ctl中路径用反斜杠'/'，而不是'/'!

1、首先使用fortran程序，将数据按照规定格式读入，包含站点名称、经纬度、观测数据值

program bintd
parameter (st=548) ! st means nums of station
character stid(st)*8

integer nflag,nlev
real lat(st),lon(st),std(st),tim,height
open(30,file='d:/eof/sn.txt')
open(40,file='d:/eof/pre_SPACE.txt')
open(50,file='d:/eof/td.dat',form='binary')
do i=1,st
read(30,*)stid(i),lat(i),lon(i),height
read(40,*)std(i)
write(*,*)stid(i),lat(i),lon(i),std(i)
end do
tim=0.0
nflag=1
nlev=1
do i=1,st
write(50)stid(i),lat(i),lon(i),tim,nlev,nflag,std(i)
write(*,*)stid(i),lat(i),lon(i),tim,nlev,nflag,std(i)
end do
nlev=0
write(50)stid(548),lat(548),lon(548),tim,nlev,nflag
write(*,*)stid(548),lat(548),lon(548),tim,nlev,nflag
close(50)
end

2、建立描述文件

dset d:/eof/td.dat
dtype station
stnmap d:/eof/map548.map
undef 999.9
title daily td
TDEF 1 LINEAR Jan1948 1mon
vars 1
std 0 99 tdpoint data
endvars

3、运行grads，运行：!stnmap -i d:/eof/td.ctl生成站点映射文件map548.map

特别提示：!stnmap -i d:/eof/td.ctl中路径用反斜杠'/'，而不是'/'!