用wrf结果画任意两点间斜剖面，画出图来了，求指导

亽小笨蜑■ɑ

WRF模拟结果经过ARWPOST后处理后，想画任意两点间斜剖面的垂直气流，看到论坛的帖子，学习后，画了两张图，贴出脚本和图，图中有一大块蓝色，是地形的影响吗？求高手指点。。。画图过程中一直报错。。说超出文件限制
'reinit'
'open F:\data\test.ctl'
'set grads off'
'set grid off'
'set zlog off'
lon1=106
lon2=126
lat1=21.0
lat2=41.0
'define alfa=atan2('lat2-lat1','lon2-lon1')'
'set x 1'
'set y 1'
'set lev 1000 100'
'set t 40'
lon=lon1
'collect 1 free'
'collect 2 free'
'collect 3 free'
while(lon<=lon2)
lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
'collect 1 gr2stn(w.1*100,'lon','lat')'
'collect 2 gr2stn(u*cos(alfa)+v*sin(alfa),'lon','lat')'
'collect 3 gr2stn(hgt(t=40),'lon','lat')'
lon=lon+1
endwhile
'set lon 106 126'
'set lev 1000 100'
'set ylevs 1000 925 850 700 600 500 400 300 200 100'
'set xlabs 41N,106E|37N,110E|33N,114E|29N,118E|25N,122E|21N,126E'
'set gxout shaded'
'set csmooth on'
'set clab off'
'set clevs 0'
'set ccols 0 5 0'
'd coll2gr(3,-u)'
'set gxout stream'
'set clab on'
'd maskout(coll2gr(2,-u),-coll2gr(3,-u));maskout(coll2gr(1,-u),-coll2gr(3,-u));mag(coll2gr(2,-u),coll2gr(1,-u))'

'draw title along (41N,106E) to (21N,126E)'
'printim F:\data\2013.10.05.18.png x1000 y800 white'
'c'
;

大气太阳茄

最近也画过这样的图。首先根据系统提示，可能是设定的经纬度范围太大，而且你还有计算部分，grads估计不能承受，建议缩小范围；其次，最后的图上出现色块的原因是因为lon=lon+1这一句，1代表资料的经向分辨率，你要改成你的ctl里面所示的分辨率。

guo871015

同求

亽小笨蜑■ɑ

guo871015 发表于 2014-6-26 20:46
同求

我一看有回复，真激动，结果你说同求！！！

guo871015

我看过river大神的这篇帖子，我也没调出来，你的程序里面test.ctl应该是wrfout的输出文件，
'collect 1 gr2stn(-omega.1*100,'lon','lat')'
'collect 2 gr2stn(u.2*cos(alfa)+v.3*sin(alfa),'lon','lat')'
'collect 3 gr2stn(orog.5(t=1),'lon','lat')'
这一段程序中首先-omega.1是指的打开的omega文件，在wrf中应该是w分量来代替omega，不用加"-"号，其次u.2文件在wrf中指的是u分量，v.3在wrf中指的是v分量，后面的".*"不用加，因为你只打开了一个ctl文件，最后orog.5是地形文件，在wrf中应该对应的变量是hgt。我之前也调试过，但是没成功。

river大神的这篇帖子中'collect 3 gr2stn(orog.5(t=1),'lon','lat')'用的地形文件和他的其他文件相对应，好像都是2.5*2.5°分辨率17个高度层的地形文件，与他用的其他文件相对应。

但是同一经度或纬度的地形垂直剖面比较简单，如果楼主需要，我可以把gs文件发给你，如果楼主自己调试出来任意直线的最好。

亽小笨蜑■ɑ

guo871015 发表于 2014-6-26 22:12
我看过river大神的这篇帖子，我也没调出来，你的程序里面test.ctl应该是wrfout的输出文件，
'collect 1 gr ...

同一经度或者纬度的我画过，谢谢你，我再研究研究这个，要不探讨探讨？

guo871015

亽小笨蜑■ɑ 发表于 2014-6-26 22:19
同一经度或者纬度的我画过，谢谢你，我再研究研究这个，要不探讨探讨？

好的，我也在研究，我不是学气象的，所以这些软件用的很不专业，嘿嘿

亽小笨蜑■ɑ

guo871015 发表于 2014-6-26 22:21
好的，我也在研究，我不是学气象的，所以这些软件用的很不专业，嘿嘿

我也不是学气象的。。。。。。你的QQ？

guo871015

亽小笨蜑■ɑ 发表于 2014-6-26 22:24
我也不是学气象的。。。。。。你的QQ？

331984865

亽小笨蜑■ɑ

guo871015 发表于 2014-6-26 22:54
331984865

不知道你叫啥。。。我按照你说的改了 出来了两张图。。。发给你看看 我的674057237