爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 9817|回复: 15

[分享资料] sst海温NC资料用grads挑选年份求平均图

[复制链接]
发表于 2013-5-10 17:38:56 | 显示全部楼层 |阅读模式

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

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

x
下面为绘此图的gs,由于本人对于aave,ave的嵌套始终尝试错误,于是采取下面无脑式编辑计算了九个年份冬季sst的平均状况,画出此图,但是所用的sst.mean.anom.nc文件lat 为87.5南北纬,图中看出有资料缺失,而且图形非常不美观,求助大神们……我也不知道这样不编辑是否正确,gs该怎么修改以及图形怎样美化……
'reinit'
'sdfopen d:\pressure\sst.mean.anom.nc'
'set cint 1'
'set grid off'
'set grads off'
'set gxout shaded'
'define a=ave(sst,t=1284,t=1286)'
'define b=ave(sst,t=1308,t=1310)'
'define c=ave(sst,t=1344,t=1346)'
'define d=ave(sst,t=1536,t=1538)'
'define e=ave(sst,t=1560,t=1562)'
'define f=ave(sst,t=1572,t=1574)'
'define g=ave(sst,t=1680,t=1682)'
'define h=ave(sst,t=1716,t=1718)'
'define i=ave(sst,t=1836,t=1838)'
'd (a+b+c+d+e+f+g+h+i)/9'
'cbarn'
'set gxout contour'
'set clab off'
'enable print d:\sst-1.gmf'
'print'
'disable print'
;
1.png

密码修改失败请联系微信:mofangbao
发表于 2013-5-10 18:15:04 | 显示全部楼层
你这九个年份没有挨着,也没有规律,只能这样求平均。图形不美观,跟你的数据精度有关,试试smth9()平滑看看。
密码修改失败请联系微信:mofangbao
发表于 2013-5-10 20:37:06 | 显示全部楼层
对于没有什么规律的年份,这样求倒也没什么,挺好的方法。不过少了一个
'set gxout contour'
'set clab off'
'd (a+b+c+d+e+f+g+h+i)/9.'
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2013-5-11 13:26:34 | 显示全部楼层
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2013-5-11 13:27:16 | 显示全部楼层
嗯,我试一试
密码修改失败请联系微信:mofangbao
发表于 2013-5-11 14:15:54 | 显示全部楼层
对于grads,我也不是很懂,就说一下看法啦~~感觉你的'set gxout contour''set clab off'应该都位于d的前面吧,还有river把9改成了9.0,可能也和数据类型有关~~我想或许可以先画等值线图,再设置阴影区域~~
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2013-5-11 16:37:14 | 显示全部楼层
逸儿~~ 发表于 2013-5-11 14:15
对于grads,我也不是很懂,就说一下看法啦~~感觉你的'set gxout contour''set clab off'应该都位于d的前面吧 ...

额,其实set gxout contour是应该在d之后,这个set clab off位置倒可变化,关于9应该是整型还是浮点型刚刚已经再次做过对比,发现没有什么变化……不过还是感谢大家为我解疑
密码修改失败请联系微信:mofangbao
发表于 2013-5-11 21:39:39 | 显示全部楼层
月是故乡明 发表于 2013-5-11 16:37
额,其实set gxout contour是应该在d之后,这个set clab off位置倒可变化,关于9应该是整型还是浮点型刚刚 ...

set gxout contour怎么可能放在d之后呢。这个命令和set gxout shaded一样,本来就是要告诉display要画的是阴影图或者等值线图,你都d完了才告诉他画什么样的图有什么用啊
密码修改失败请联系微信:mofangbao
发表于 2013-5-11 21:40:40 | 显示全部楼层
逸儿~~ 发表于 2013-5-11 14:15
对于grads,我也不是很懂,就说一下看法啦~~感觉你的'set gxout contour''set clab off'应该都位于d的前面吧 ...

grads里面必须先画阴影图再画等值线图,反过来的话阴影会覆盖掉等值线
密码修改失败请联系微信:mofangbao
发表于 2013-5-11 22:34:34 | 显示全部楼层
river 发表于 2013-5-11 21:40
grads里面必须先画阴影图再画等值线图,反过来的话阴影会覆盖掉等值线

嗯哈,多谢啦~~以前都是画等值线图,特定部分才加阴影的,嘿嘿~~
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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