- 积分
- 120
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-4-1
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
看了清风大侠的帖子之后,综合了一些其他帖子内容,画近半个月时间,来实现自动读取micaps第三类站点资料文件(本帖...\r24-8-p),把读取的数据(本贴数据是2000年1月1日08时到31日08时,共计31天的降水资料)存到一个数据格点文件里,文中附有多时次站点资料转格点资料的Fortran程序(【终版】批量读取+txt转dat)。在程序里用一个文档文件来检验数据读写是否正常((jianyan)sta2000-01-end)。
这是多时次的sta2000-01-end.ctl:
DSET f:\sta2000-01-end.grd
DTYPE station
STNMAP f:\sta2000-01-end.map ********注意.用grads生成.map文件 语句:!stnmap -i D:/000/6/sta.ctl (注意/的方向,更改你使用的地址
UNDEF 9.999e+20
TITLE 2000-01
TDEF 31 linear 01jan2000 1dy
VARS 1
rain 0 99 Rain Data
ENDVARS
根据文本文档(jianyan)sta2000-01-end,检验了数据是按31个时次,每个时次后都有结束标志(nlev=0)。
但是,一到grads运行 !stnmap -i f:/sta2000-01-end.ctl,就提示错误(如图grads提示的错误)
不知道错在哪儿,请大神指点迷津。
|
|