爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 2747|回复: 3

[脚本编辑] 求助!用grads转化为二进制后画图出错

[复制链接]

新浪微博达人勋

发表于 2015-4-25 10:15:51 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
各位好,我下载的逐日数据格式有点奇怪,但是网站配有ctl可以打开,每一个逐日数据是一个文件,数据是V1.0_0.50deg.lnx.19790501,V1.0_0.50deg.lnx.19790502,V1.0_0.50deg.lnx.19790503。。。我现在想用让它们共用一个crl文件把它们合并到一个grd文件中,(将19790501—19790930这153个数据合并)。但是合并后的数据画某一天的图和原始数据不一样,只有第一天是一样的,后面都不一样,不知道是什么问题??请教各位!!!!谢谢!!!
ctl 如下:
dset e:\1979\V1.0_0.50deg.lnx.1979%m2%d2
options template
title  Pre
undef -999.0
xdef 720 linear   0.25 0.50
ydef 360 linear -89.75 0.50
zdef 1   linear 1 1
tdef 153 linear 01May1979 1dy
vars 2
r    1  00 the grid analysis (0.1mm/day)
gnum     1  00 the number of stn
ENDVARS

gs 如下:
'reinit'
'open E:\1979\Pre1979.ctl'
'set gxout fwrite'
'set fwrite E:\1979\1979_56789.grd'
'set grads off'
tn=1
while(tn<=153)
'set x 1 720'
'set y 1 360'
'set z 1'
'set t 'tn''
'd r'
tn=tn+1
endwhile
'disable fwrite'
''
;
合并后的数据的ctl:
dset e:\1979\1979_56789.grd
title Pre
undef -999.0
xdef 720 linear 0.25 0.50
ydef 360 linear -89.75 0.50
zdef 1   linear 1 1
tdef 153 linear 01May1979 1dy
vars 2
r    1  00 the grid analysis (0.1mm/day)
gnum     1  00 the number of stn  
ENDVARS

用合并后数据的ctl画图会有这样错误信息 QQ截图20150425101507.png
请教各位是什么原因???
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-4-25 11:00:01 | 显示全部楼层
你只输出了一个变量,为什么描述文件要有两个变量
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-4-25 14:50:48 | 显示全部楼层
传说中的谁 发表于 2015-4-25 11:00
你只输出了一个变量,为什么描述文件要有两个变量

e:\11a.png如果ctl里不写另一个变量只写输出变量,ctl就打不开!!我也觉得好奇怪,不知道为什么!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-4-25 14:56:02 | 显示全部楼层
这个是错误信息
11a.png
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表