画的剖面风场，没有值的地方是一个个小圆圈，怎么能把它去掉

qxjy009 · 发表于 2016-11-2 15:18:14

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

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

x

qxjy009 · 发表于 2016-11-2 15:20:10

我用的gs就是仿着river的任意剖面
'reinit'
'open E:\wrfout-0712-d02\1010-gfs-20-20\p-axis\gfs-2020.ctl'

'set grid off'
'set grads off'
*'set zlog on'
lon1=114
lon2=116
lat1=30
lat2=32
alfa=math_atan2(lat2-lat1,lon2-lon1)
*say alfa
'set x 1'
'set y 1'
'set lev 1000 300'
'set t 1'
ln=lon1

'collect 1 free'
'collect 2 free'

while(ln<=lon2)
lt=lat1+(lat2-lat1)*(ln-lon1)/(lon2-lon1)
'collect 1 gr2stn(-w*100,'ln','lt')'
'collect 2 gr2stn(u*cos('alfa')+v*sin('alfa'),'ln','lt')'
ln=ln+0.06
endwhile
say 'ok...'

'set lon 114 116'
'set ylevs 1000 850 700 600 500 400 300'
'set csmooth on'
'set clab off'
'set gxout barb'
'd coll2gr(2,-u)*2.5;coll2gr(1,-u)*2.5'

qxjy009 · 发表于 2016-11-2 15:20:49

然后其实我叠加了地形能盖住大部分小圆圈，但是最上面的还是有些盖不住

fengaiyu · 发表于 2016-11-2 19:57:56

问题解决了吗？

qxjy009 · 发表于 2016-11-3 08:57:07

没有呢，没头绪

river · 发表于 2016-11-5 08:05:33

qxjy009 发表于 2016-11-2 15:20
我用的gs就是仿着river的任意剖面
'reinit'
'open E:\wrfout-0712-d02\1010-gfs-20-20\p-axis\gfs-2020.c ...

没有值是说那里都是缺测值吗？如果是缺测的话，grads应该不会画出来的

占戈示申 · 发表于 2016-11-5 12:31:07

画到950呢？不画1000的。

qxjy009 · 发表于 2016-11-9 15:41:36

river 发表于 2016-11-5 08:05
没有值是说那里都是缺测值吗？如果是缺测的话，grads应该不会画出来的

我也觉得奇怪，这到底是个啥，就是wrf转出来的dat画的

signor9 · 发表于 2017-5-7 17:50:08

我遇到了类似的问题，降水图上是一个个的小圆圈，也不知道该怎么解决

		自动登录	找回密码
密码			立即注册

[图形美化] 画的剖面风场，没有值的地方是一个个小圆圈，怎么能把它去掉

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