爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 6056|回复: 10

[脚本编辑] Grads沿任意两点作倾斜剖面时垂直层数出现问题

[复制链接]

新浪微博达人勋

发表于 2016-9-30 14:09:41 | 显示全部楼层 |阅读模式

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

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 花小菜 于 2016-9-30 22:12 编辑

请教各位大神,用Grads作斜剖时,画出来的剖面图中垂直层数出现问题,不均匀的高度层变成了均匀的高度层,如下图:
17_1.png
我的sigma层是不等间距的35层,画出来应该是下面的风场比较密上面比较疏,但是剖面图中每层风矢量间距相同,变成了均匀的,不知道哪里出现了问题?

35层sigma如下:(单位:km)   
0.01528 0.03056 0.06118 0.10207 0.14307 0.18418 0.22543 0.28765 0.37111 0.45526 0.54012 0.62569 0.71194 0.79888 0.88651 0.97484 1.10891 1.28984 1.47375 1.66076 1.85102 2.04467 2.34324 2.75291 3.17928 3.62424 4.15096 4.76847 5.4264 6.13116 6.89191 7.72143 9.16967 11.50275 14.75942





































gs如下:

'reinit'
'open F:\0531\cntl\cntl.ctl'
'set grads off'
'set grid off'
'set mpdset cnworld'

lon1=109.68
lon2=110.46
lat1=19.91
lat2=19.24
'define alfa=atan2('lat2-lat1','lon2-lon1')'
'set x 1'
'set y 1'
'set z 1 27'
'set t 15'
lon=lon1
'collect 1 free'
'collect 2 free'
while(lon<=lon2)
lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
'collect 1 gr2stn(w,'lon','lat')'
'collect 2 gr2stn(u*cos(alfa)+v*sin(alfa),'lon','lat')'
lon=lon+0.01
endwhile

'set lon 109.68 110.46'
'set z 1 27'
'set xlabs 19.90N,109.68E|19.68N,109.94E|19.46N,110.20E|19.24N,110.46E'

'set gxout shaded'
'set clevs 0.2 0.4 0.6 0.8 1 1.2 1.4'
'set ccols 0 13 10 7 12 8 2 6'
'd smth9(coll2gr(1,27))'

'set gxout vector'
'drawvect skip(coll2gr(2,27),4,1) skip(coll2gr(1,27)*10,1,1) tr 0.4 8 1 5 5 0.18 0 m`31`5s`a-1`n'

'cbar_interp 1 1 1 9.8 4.2'
'print F:\0531\cntl\section\oblique\14.eps'
;

求大神解答~~~多谢!!!

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-9-30 21:11:18 | 显示全部楼层
加上set zlog on试试
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-9-30 21:56:03 | 显示全部楼层
本帖最后由 花小菜 于 2016-9-30 22:02 编辑
river 发表于 2016-9-30 21:11
加上set zlog on试试


谢谢~~刚才试了一下,发现加上之后好像是把纵坐标轴变成不均匀了,还是很奇怪,是不是还缺点什么??
17_2.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-9-30 22:22:37 | 显示全部楼层
花小菜 发表于 2016-9-30 21:56
谢谢~~刚才试了一下,发现加上之后好像是把纵坐标轴变成不均匀了,还是很奇怪,是不是还缺点什么??
...

这是啥情况,我也搞不懂了·····
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-9-30 22:25:07 | 显示全部楼层
river 发表于 2016-9-30 22:22
这是啥情况,我也搞不懂了·····

好吧,其实我是看过你之前发的帖子,程序也是参照你的~~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-9-30 23:22:31 | 显示全部楼层
花小菜 发表于 2016-9-30 22:25
好吧,其实我是看过你之前发的帖子,程序也是参照你的~~

应该是插值以后就插值成均匀的了吧,我觉得
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-10-1 08:19:30 | 显示全部楼层
学习了,谢谢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-10-1 08:20:24 | 显示全部楼层
river 发表于 2016-9-30 23:22
应该是插值以后就插值成均匀的了吧,我觉得

嗯呢,我觉得也是这样,但是不知道怎么改
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-10-1 16:19:49 | 显示全部楼层
学习了斜剖,感谢楼主
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-10-1 16:21:12 | 显示全部楼层
yujie 发表于 2016-10-1 16:19
学习了斜剖,感谢楼主


密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表