请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 831|回复: 4

[脚本编辑] 水汽散度数据的提取

[复制链接]

新浪微博达人勋

发表于 2016-10-31 17:34:58 | 显示全部楼层 |阅读模式

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

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

x
我想把水汽散度的数据提取出来,但是出的图有问题,有没有人可以帮忙看下啊
gs文件:
'reinit'
'open C:\Users\Administrator\Desktop\js\qua.ctl'
'open C:\Users\Administrator\Desktop\js\qva.ctl'
'set x 1 144'
'set y 1 73'
'set gxout fwrite'
'set fwrite C:\Users\Administrator\Desktop\js\sd.grd'
i=1
while (i<=639)
'set t 'i''
'define qdiv=hdivg(qua.1,qva.2)*100000'
'd qdiv'
i=i+1
endwhile
'disable fwrite'
;
ctl文件:
DSET C:\Users\Administrator\Desktop\js\sd.grd
TITLE sd
UNDEF -9.96921e+36
XDEF 144 LINEAR 0 2.5
YDEF 73 LINEAR -90 2.5
ZDEF 1  LINEAR 1 1
TDEF 639 LINEAR DEC1961 1mo
vars 1
qdiv 1 99 sd
endvars
qua.ctl和qva.ctl是经纬向分量,可以正常出图,但是计算散度出的图就这样了
         

任一时次的散度

任一时次的散度

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-10-31 18:00:06 | 显示全部楼层
改你的缺测值
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-10-31 19:18:51 | 显示全部楼层

缺测值怎么改啊?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-10-31 20:05:44 | 显示全部楼层
谢谢river,这个问题已经解决了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-10-31 20:06:39 | 显示全部楼层

谢谢river,问题解决了
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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