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

gh18350146

river 发表于 2013-5-2 13:42

                               
登录/注册后可看大图

coll2gr(2,-u)红色的参数是在你的资料是等压面的情况下才能用的。如果你的实际资料就是实际高度的，那后面 ...

好的谢谢呢，我试试看。可是coll2gr(1,20)能画出雷达回波，风场应该也能画出来的。（雷达回波和风场是同等格式资料）

qazgsy90

来学习下，也想画一张。

banna

漂亮。。。。学习了

晓艾

你好，请问你的问题解决了么？？我也需要画这个图，可是改了之后还是画不出来。

gh18350146

晓艾 发表于 2014-4-15 19:18
你好，请问你的问题解决了么？？我也需要画这个图，可是改了之后还是画不出来。

解决了。红色为改动的部分，改后风场能画出来。你试试
'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*cos(alfa)),'lon','lat')'
'collect 3 gr2stn((V*sin(alfa)),'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)+coll2gr(3,20);coll2gr(4,20)'
'cbarn 1 1 9.8 4.39'
'print'
'disable print'

晓艾

谢谢了！！

大气太阳茄

请问一下，gs中
lon=lon+0.02
这一句的0.02是怎么确定的呢？我也在做剖面图，发现lon=lon+1和lon=lon+2画出的图不一样啊，不知道应该取多少？

gh18350146

大气太阳茄 发表于 2014-6-22 21:09
请问一下，gs中
lon=lon+0.02
这一句的0.02是怎么确定的呢？我也在做剖面图，发现lon=lon+1和lon=lon+2画 ...

这一句表示剖面时候经纬度分辨率。结合自己需要把。在原资料分辨率允许的情况下选取。

偏爱路人乙

初学者进来学习

smilehfr

请问楼主  垂直方向上的风矢量疏密程度对么 我现在等间距的 不知道是什么原因