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

梦纯

                               
登录/注册后可看大图

各位大虾们好，求教了！怎么用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报错如下：
请各位大侠们帮帮忙啊！！谢谢啦！拜托啦！

youyou

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'
;

youyou

梦纯 发表于 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'




不知道你能不能用得上？

youyou

画了一张假相当位温的   不过只画了线条和地形  你看看是不是你要的那样的形式？

梦纯

youyou 发表于 2015-4-11 10:35
画了一张假相当位温的   不过只画了线条和地形  你看看是不是你要的那样的形式？

这样的图也可以哎。能不能共享一下你的gs文件啊？谢谢了哈！

梦纯

我试试哈，多谢啦！！！

梦纯

youyou 发表于 2015-4-11 15:09
'reinit'
'open D:\laoshu\fnl2014.ctl'
'open d:\2013Grads\gs\surface_elevation\dz.ctl'

youyou,你的地形资料是12层的么？

youyou

梦纯 发表于 2015-4-11 15:23
youyou,你的地形资料是12层的么？

多少层都可以 ，只要和你想画的那个的那个东西的层次配合好就行……要不竖坐标就乱了

笑哈乐呵

努力学习

19921219

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'


;

19921219

19921219 发表于 2015-5-5 15:39
想问一下，我叠加画omege，但是只能出地形图，其他物理量出不来！是为什么？
'reinit'
'open h:\data\d ...

不仅仅是omege其他物理量也是的