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

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3760|回复: 5

[图形美化] 如何使南风为红色箭头,北风为蓝色箭头??

[复制链接]

新浪微博达人勋

发表于 2018-8-8 17:27:17 | 显示全部楼层 |阅读模式

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

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

x
!!!求助,我想令南风为红色箭头,北风为蓝色箭头,用的是nc资料,grads该咋么做呢?我想到一种办法,就是用atan2函数判断,但是我画出来都是一个颜色,红色。算出来的风向值都远远大于360,自己实在没想出哪错了?(这个NC资料单位是m/s)。还有办法给南北风分别上色吗?
参考过论坛里的一篇(下面链接)
http://bbs.06climate.com/forum.p ... 1&fromuid=72572
(出处: 气象家园)

'reinit'
'set grads off'
'set grid off'
'sdfopen D:/ncep/daily/uwnd.2018.nc'
'sdfopen D:/ncep/daily/vwnd.2018.nc'
'set x 1 144'
'set y 1 73'
'set lev 500'
'set t 12'
***!!!用u,v分量计算风向(值域为0-360)
'define radian=180./3.1415926'     
'define wind=atan2(vwnd.2,uwnd)*radian+180'
***!!!判断风向
if(wind>180)  
   'set gxout vector'
   'set arrscl 0.3'
   'set ccolor 2'   ;**红色
   'set arrowhead -0.25'
   'd skip(uwnd,4,4);skip(vwnd.2,4,4)'
else
   'set gxout vector'
   'set ccolor 11'   ;**蓝色
   'set arrscl 0.3'
   'set arrowhead -0.25'
   'd skip(uwnd,4,4);skip(vwnd.2,4,4)'
endif
'printim d:\coldwave\h/AAA.png white'
;




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

新浪微博达人勋

 楼主| 发表于 2018-8-13 09:22:18 | 显示全部楼层
本帖最后由 董伟是我 于 2018-8-13 09:33 编辑

H:\coldwave\IPV\E/20180124.PNG   
已解决,如图。我现在所知道的最简单的方法是可利用abs函数
如:   'set gxout vector'
  'run D:\GRADS\redbule.gs'
  'set rbcols 60 50'
  'set arrscl 0.7 50'
  'set arrowhead -0.6'
  'set cthick 8'
  'd  maskout(skip(ue,4,4),mag(ue,ve)-15);maskout(skip(ve,4,4),mag(ue,ve)-15);maskout(skip(ve,4,4),mag(ue,ve)-15)/(abs(maskout(skip(ve,4,4),mag(ue,ve)-15)))'


简单点就是   d u;v;v/abs(v)
201824.png
20180124.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-8-13 09:24:15 | 显示全部楼层
H:\coldwave\IPV\E/20180124.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-8-14 22:24:04 | 显示全部楼层
感谢楼主无私分享
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-8-15 09:32:44 | 显示全部楼层
厉害~感谢楼主无私分享!!!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2019-3-6 17:12:55 | 显示全部楼层
谢谢楼主分享
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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