求助Grads绘图：wrf输出资料如何作垂直剖面（非东-西或南-北向）

gh18350146

wrf输出资料，用grads作（111.5°E,23°N）至（114°E,22°N）垂直剖面图，要素包括雷达回波及风场。
雷达回波能画出，风场无法画出。所得图及对应gs如下：

'reinit'
'open **.ctl'
'enable print **.gmf'
'set grads off'
'set grid off'
'set x 1'
'set y 1'
'set z 1 20'
lon1=111.5
lon2=114
lat1=23
lat2=22
'define alfa=atan2('lat2-lat1','lon2-lon1')'
lon=lon1                             
'collect 1 free'
'collect 2 free'
'collect 3 free'
'collect 4 free'
while(lon<=lon2)
lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
'collect 1 gr2stn(dbz,'lon','lat')'
'collect 2 gr2stn((U),'lon','lat')'
'collect 3 gr2stn((V),'lon','lat')'
'collect 4 gr2stn(W,'lon','lat')'
lon=lon+0.02
endwhile
'set lon 111.5 114'
'set z 1 20'

'set gxout shaded'
'run C:\GrADS20\lib\radarcolor.gs'
'set clevs 5 10 15 20 25 30 35 40 45 50 55 60 65 70'
'set ccols 0 21 22 23 24 25 26 27 28 29 30 31 32 33 34'
'd coll2gr(1,20)'

'set gxout vector'
'set ccolor 1'
'set arrscl 0.35 15'
'd coll2gr(2,20)*cos(alfa)+coll2gr(3,20)*sin(alfa);coll2gr(4,20)'
'cbarn 1 1 9.8 4.39'
'print'
'disable print'

grads错误如下：


求高人指点~

mofangbao

根据提示信息看，是你两个格点的维度不匹配，请尝试设置同一个维度之后画图

gh18350146

mofangbao 发表于 2013-4-24 10:20

                               
登录/注册后可看大图

根据提示信息看，是你两个格点的维度不匹配，请尝试设置同一个维度之后画图

谢谢，弱弱的问一句，如何设置两个格点同一个维度，请高手明示

kongfeng0824

可以看出图还是很漂亮的。

gh18350146

kongfeng0824 发表于 2013-4-24 21:52

                               
登录/注册后可看大图

可以看出图还是很漂亮的。

。。。。么得风场出现，各种伤悲。

river

set z 1 20  你要画的是从1米到20米的高度？试着改改红色的数字d coll2gr(2,20)*cos(alfa)+coll2gr(3,20)*sin(alfa);coll2gr(4,20)

kongfeng0824

gh18350146 发表于 2013-4-25 15:14

                               
登录/注册后可看大图

。。。。么得风场出现，各种伤悲。

你指的是风场矢量的那种有箭头的图吗？

gh18350146

kongfeng0824 发表于 2013-4-25 21:31

                               
登录/注册后可看大图

你指的是风场矢量的那种有箭头的图吗？

是的呢。wrf输出资料，雷达回波能画出来，风场没有画出来~不知道插值哪里出问题。

gh18350146

river 发表于 2013-4-25 15:37

                               
登录/注册后可看大图

set z 1 20  你要画的是从1米到20米的高度？试着改改红色的数字d coll2gr(2,20)*cos(alfa)+coll2gr(3,20)*s ...

set z 1 20 这个是wrf输出资料里面的层次，20层对应17km。好的我试试，谢谢啦
  0.25000
   0.50000
   0.75000
   1.00000
   2.00000
   3.00000
   4.00000
   5.00000
   6.00000
   7.00000
   8.00000
   9.00000
  10.00000
  11.00000
  12.00000
  13.00000
  14.00000
  15.00000
  16.00000
  17.00000

river

gh18350146 发表于 2013-4-26 11:15

                               
登录/注册后可看大图

set z 1 20 这个是wrf输出资料里面的层次，20层对应17km。好的我试试，谢谢啦
  0.25000
   0.50000

coll2gr(2,-u)红色的参数是在你的资料是等压面的情况下才能用的。如果你的实际资料就是实际高度的，那后面就要写成你实际的那个高度吧，不应该用层次了就，你试下