求问 GRADS怎么画大于某个风速的风向杆图

yanglei_nuist

最近学习使用GRADS处理NCEP资料，想问大家下，我想画一个风场，其中风速和风向用风向杆表示，但是我只想画大于12ms-1的风场，想问下大家怎么设置语句，谢谢了。下面有我编过的GS文件。

'reinit'
'open E:/20130818/fnl201308.ctl'
'set vpage 0 11 0 8.5'
'set parea 1 10.5 1 7'


m=1
while(m<17)
  'set grads off'
  'set grid off'
  'set timelab off'

  'set lat 10 55'
  'set lon 60 150'
  'set mpdset cnworld'
  'set map 1 3 3'
  'draw map'
  'set t ' m

  'set clevs 55 60 65 70'
  'set rbcols 0 4 11 5 7'

  'set gxout shaded'
  'd PWATclm'
  'cbarn 1 0'

  'set z 6'
  'set gxout barb'
  'set digsiz 0.1'
  'd UGRDprs*2.5;VGRDprs*2.5'

  'printim E:\20130818\PWATclm\PWATclm_'m'.jpg white'
   m=m+1
   'c'
endwhile
'print'  
'disable print'
;

yanglei_nuist

fyy123 发表于 2017-3-29 20:50
想问一下，为什么'd UGRDprs*2.5;VGRDprs*2.5'这里要乘以2.5啊？

grads里风标的单位是英寸，咱们是米每秒，乘以2.5正好

夜未央

请问为什么这样画出来的风场是彩色的？
设置ccolor以后也不管用，请问怎么改成黑色？
'set cmin 6'
'set ccolor 1'
'define a=mag(u,v)'
'd u;v;a'

njzqxt

出图前加句
set cmax 12

river

使用set cmax 和set cmin 命令配合display就可以了。比如你的这个这么写就行
set gxout brab
set cmin 12
d UGRDprs*2.5;VGRDprs*2.5;mag(UGRDprs,VGRDprs)



注意红色部分，必须有，否则set cmin 12就不起作用了

yanglei_nuist

river 发表于 2013-11-1 08:18 使用set cmax 和set cmin 命令配合display就可以了。比如你的这个这么写就行 set gxout brab set cmin 12 ...

你好，谢谢你的回复。我试过了还是不行啊，我查了手册mag是算风的矢量和，请问mag的这句话怎么写啊，按照你的回复的写法，不好使啊

river

yanglei_nuist 发表于 2013-11-2 06:59

                               
登录/注册后可看大图

你好，谢谢你的回复。我试过了还是不行啊，我查了手册mag是算风的矢量和，请问mag的这句话怎么写啊，按照 ...

怎么个不好使，有错误提示吗？这是我试验过的，如果是我的猜测或者不确定的一个回复我会说明的

yanglei_nuist

river 发表于 2013-11-2 07:28 怎么个不好使，有错误提示吗？这是我试验过的，如果是我的猜测或者不确定的一个回复我会说明的

‘set gxout brab’                                    ‘set cmin 12’                                               ‘d UGRDprs*2.5;VGRDprs*2.5’             ‘mag(UGRDprs,VGRDprs)’我这样写的话，软件提示‘mag(UGRDprs,VGRDprs)’ unknown。                                 ‘set gxout brab’                                    ‘set cmin 12’                                               ‘d UGRDprs*2.5;VGRDprs*2.5’             ‘d mag(UGRDprs,VGRDprs)’这样写的话就可以出图了，但是相当于画了两遍风场。                                            'set gxout brab’                                    ‘set cmin 12’                                               ‘d mag(UGRDprs,VGRDprs)’，这样写的话，就根本没画不出来风场。不知道我理解您第一条回复是否正确，可能语句哪里写错了，麻烦您指正。                                                                           

yanglei_nuist

river 发表于 2013-11-2 07:28 怎么个不好使，有错误提示吗？这是我试验过的，如果是我的猜测或者不确定的一个回复我会说明的

‘set gxout brab’                                    ‘set cmin 12’                                               ‘d UGRDprs*2.5;VGRDprs*2.5’             ‘mag(UGRDprs,VGRDprs)’我这样写的话，软件提示‘mag(UGRDprs,VGRDprs)’ unknown。                                 ‘set gxout brab’                                    ‘set cmin 12’                                               ‘d UGRDprs*2.5;VGRDprs*2.5’             ‘d mag(UGRDprs,VGRDprs)’这样写的话就可以出图了，但是相当于画了两遍风场。                                            'set gxout brab’                                    ‘set cmin 12’                                               ‘d mag(UGRDprs,VGRDprs)’，这样写的话，就根本没画不出来风场。不知道我理解您第一条回复是否正确，可能语句哪里写错了，麻烦您指正。                                                                           

river

yanglei_nuist 发表于 2013-11-2 08:39

                               
登录/注册后可看大图

‘set gxout brab’                                    ‘set cmin 12’                              ...

你看清楚我的回复，为什么要把好端端的一句分开写呢。我再写一遍
set gxout brab
set cmin 12
d UGRDprs*2.5;VGRDprs*2.5;mag(UGRDprs,VGRDprs)

yanglei_nuist

river 发表于 2013-11-2 12:50 你看清楚我的回复，为什么要把好端端的一句分开写呢。我再写一遍 set gxout brab set cmin 12

谢谢您，我之前看错了，试过啦，可以出图了，万分感谢。

river

yanglei_nuist 发表于 2013-11-2 22:14

                               
登录/注册后可看大图

谢谢您，我之前看错了，试过啦，可以出图了，万分感谢。

小杨同学，没必要那么客气。我还年轻呢，咱俩同辈···