请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 9446|回复: 15

比湿经度-高度剖面图和地形图怎么叠加

[复制链接]

新浪微博达人勋

发表于 2015-4-10 15:39:17 | 显示全部楼层 |阅读模式
GrADS
系统平台: win7
问题截图:
问题概况: 比湿的剖面图和地形叠加,实现不了。求教!
我看过提问的智慧: 看过
自己思考时长(天): 3

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

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

x
本帖最后由 梦纯 于 2015-4-10 15:39 编辑


                               
登录/注册后可看大图
各位大虾们好,求教了!怎么用grads画出来如上图所示的比湿剖面图和地形叠加图啊????

我自己将30度的比湿剖面数据和30度的地形数据分别存成了两个文件:
diff_vapor.ctl'topo_30.ctl'
画图的gs文件如下:'reinit'
'open D:\PhD\run_topo_0.5\vapor\diff_vapor.ctl'
'open D:\PhD\run_topo_1.0\topo\topo_30.ctl'
'set grads off'
'set grid off'
'set gxout shaded'
'set xlopts -10'
*'set ylopts -10'
'set lon 75 135'
'set clevs -0.5 -0.3 -0.1 0 0.1 0.3 0.5'
'set z 1 11'
*'set zlog on'
'd vapor.1'
'cbar_interp'
'set gxout bar'
'set bargap 80'
'set lon 75 135'
'd topo.2'
'run axis.gs  -type b -position o -label on  -interval 10 -suffix `3.`'
*'run axis.gs  -type r  -v  topo.1  -position o -label on -position o -hoffset -0.35 -start 0 -end 10000 -interval 1000 '
grads报错如下: 无标题.png
请各位大侠们帮帮忙啊!!谢谢啦!拜托啦!






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

新浪微博达人勋

发表于 2015-5-6 09:51:03 | 显示全部楼层
19921219 发表于 2015-5-5 16:39
如果不加地形,我的图就可以出来,但是有地形后,就不行了!!


按照下面的GS,把你的资料的路径换了就能出图了,例图也在后面给你贴上去了,你的'set xlopts ' 'set ylopts '在画物理量和地形时设置的不一样,其它的就是要注意一个覆盖的问题就可以了


'reinit'
'sdfopen D:\2013Grads\nc/omega.2013.nc'
'open d:\2013Grads\gs\surface_elevation\dz.ctl'
'set grads off'
'set grid off'
'set csmooth on'
'set map 15 1 9'
'set parea 0.8 10.8 0.8 7.9'
'set xlopts 1 6 0.14'
'set ylopts 1 6 0.14'
'set font 0'
'set lev 1000 300'



'set dfile 2'
'set lat 30 '
'set lon 70 140'
'set t 1'
'set z 1 8'
*'set zlog on'
'set gxout shaded'
'set cint 8000'
'set rgb 18 125 125 125'
'set csmooth on'
'set rbcols 1'
'set cmax 0'
'define zz=dz.2'
'd zz'

'set dfile 1'
'set lat 30 '
'set lon 70 140'
'set gxout contour'
'set cint 0.05'
'set cthick 11'
'set clab forced'
'd omega'
'printim d:\123.png white'
;
123.png
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2015-4-11 15:09:10 | 显示全部楼层
梦纯 发表于 2015-4-11 15:04
这样的图也可以哎。能不能共享一下你的gs文件啊?谢谢了哈!

'reinit'
'open D:\laoshu\fnl2014.ctl'
'open d:\2013Grads\gs\surface_elevation\dz.ctl'

'set grads off'
'set grid off'
'set csmooth on'
'set map 15 1 9'
'set parea 0.8 10.8 0.8 7.9'
'set xlopts 1 6 0.14'
'set ylopts 1 6 0.14'
'set font 0'
'set lev 1000 100'
'set time 20Z10MAY2014'

*---------计算假相当位温------------
'define br=(0.1158*log(rhprs)-1.5332787)*tmpprs-789.92816'                              
'define cr=789.92816*tmpprs'                                                         
'define tc2=(-br-sqrt(br*br-4*cr))/2.0'                                               
'define pc2=lev*pow(tc2/tmpprs,3.5005574)'                                                        
'define etc2=6.1078*pow(273.16/tc2,5.1347779)*exp(3148.0973*(tc2-273.16)/(126.06334*tc2))'   
'define wc2=0.622*etc2/(pc2-etc2)'                                                        
'define ed2=tc2*pow(1000.0/(pc2-etc2),0.28765924)'                                          
'define lc2=2500.79-2.3697*(tc2-273.16)'                                                  
'define eqt=ed2*exp(wc2*lc2/1.0048/tc2)'


*---------输出地形------------
'set strsiz 0.2'
'set xlopts 1 4 0.15'
'set ylopts 1 4 0.15'
'set grads off'
'set grid off'
'set dfile 2'
'set lat 30 '
'set lon 70 140'
'set t 1'
'set z 1 12'
*'set zlog on'
'set gxout shaded'
'set cint 8000'
'set rgb 18 125 125 125'
'set csmooth on'
'set rbcols 1'
'set cmax 0'
'define zz=dz.2'
'd zz'


*---------输出图形------------
'set lat 30'
'set lon 70 140 '

'set gxout contour'
'set ccolor 1'
'set cthick 6'
'set clskip 1'
'set cint 3'
'd eqt'

'printim D:\0.gif white'




不知道你能不能用得上?
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2015-4-11 10:35:25 | 显示全部楼层
画了一张假相当位温的   不过只画了线条和地形  你看看是不是你要的那样的形式?
0.gif
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-4-11 15:04:22 | 显示全部楼层
youyou 发表于 2015-4-11 10:35
画了一张假相当位温的   不过只画了线条和地形  你看看是不是你要的那样的形式?

这样的图也可以哎。能不能共享一下你的gs文件啊?谢谢了哈!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-4-11 15:13:42 | 显示全部楼层
我试试哈,多谢啦!!!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-4-11 15:23:49 | 显示全部楼层
youyou 发表于 2015-4-11 15:09
'reinit'
'open D:\laoshu\fnl2014.ctl'
'open d:\2013Grads\gs\surface_elevation\dz.ctl'

youyou,你的地形资料是12层的么?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-4-11 23:26:05 | 显示全部楼层
梦纯 发表于 2015-4-11 15:23
youyou,你的地形资料是12层的么?

多少层都可以 ,只要和你想画的那个的那个东西的层次配合好就行……要不竖坐标就乱了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-4-16 08:42:18 | 显示全部楼层
努力学习
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2015-5-5 15:39:29 | 显示全部楼层
youyou 发表于 2015-4-11 15:09
'reinit'
'open D:\laoshu\fnl2014.ctl'
'open d:\2013Grads\gs\surface_elevation\dz.ctl'

想问一下,我叠加画omege,但是只能出地形图,其他物理量出不来!是为什么?
'reinit'
'open h:\data\dz.ctl'
'sdfopen h:/data/omega.2014.nc'
'set fwrite h:/data/omega.grd'
*'enable print h:/data/omega/0720.gmf'

'set grads off'
'set grid off'
'set csmooth on'
'set map 15 1 9'
'set parea 0.8 10.8 0.8 7.9'
'set xlopts 1 6 0.14'
'set ylopts 1 6 0.14'
'set font 0'
'set lev 1000 300'

'set dfile 2'
'set lat 30 '
'set lon 70 140'
'set gxout contour'
'set cint 0.05'
'set cthick 11'
'set clab forced'

'd omega.2'

'set strsiz 0.2'
'set xlopts 1 4 0.15'
'set ylopts 1 4 0.15'
'set grads off'
'set grid off'
'set dfile 1'
'set lat 30 '
'set lon 70 140'
'set t 1'
'set z 1 8'
*'set zlog on'
'set gxout shaded'
'set cint 8000'
'set rgb 18 125 125 125'
'set csmooth on'
'set rbcols 1'
'set cmax 0'
'define zz=dz.1'
'd zz'


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

新浪微博达人勋

发表于 2015-5-5 15:40:04 | 显示全部楼层
19921219 发表于 2015-5-5 15:39
想问一下,我叠加画omege,但是只能出地形图,其他物理量出不来!是为什么?
'reinit'
'open h:\data\d ...

不仅仅是omege其他物理量也是的
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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