如何将一维图和二维图叠加？

梳子

我的数据是这样的，在有地形的情况下绘制一个垂直剖面图，物理量是三维的，剖面图就是二维的纬向z轴剖面，同时我的数据里包含有一个二维的地形，取垂直剖面就是一维图，于是我就用set gxout bar来画，但是结果两者的坐标轴不能重合，各位有没有什么好办法？非常感谢！

下面附上我的ctl和gs，以及画出的图：
DSET G:\20090730\mref\dat\0907%d2%h2%n2.dat
options template
TITLE 3DNet_radar
UNDEF -9999
XDEF 1116 linear   99.410 0.01
YDEF 1053 linear   24.879 0.01
ZDEF 28 LEVELS 300 500 800 1000 1200 1500 2000 2500 3000 3500 4000 4500 5000 5500 6000 6500 7000 7500 8000 8500 9000 9500 10000 11000 12000 13000 14000 15000
TDEF 480 LINEAR 00:00Z30Jul2009 6mn
VARS 2
ref  28 99 3DRadarData
hgt  0 99 height
ENDVARS

'reinit'
'c'
'open G:\20090730\mref\dat\mref_cs.ctl'
'enable print G:\20090730\mref\dat\2009'


'set lon 102 108'
'set lat 32'
it=1
while(it<=1)
'set grads off'
'set grid off'

'set t 'it''
'set gxout shaded'
'set z 1 28'
'run radarcolor.gs'
'set rbcols 21 22 23 24 25 26 27 28 29 30 31 32 33 34'
'set clevs -10 0 5 10 15 20 25 30 35 40 45 50 55 60 65'
'd ref'
'run cbarn.gs'

'set z 1'
'set vrange 500 15000'
'set gxout bar'
'set ccolor 1'
'd hgt'

'q time'
tt=subwrd(result,3)
'draw title 'tt' REF'
if(it<10)
'printim G:\20090730\mref\dat\cs09073100'it'.png white x1000 y800 '
else
if(it<100)
'printim G:\20090730\mref\dat\cs0907310'it'.png white x1000 y800'
else
'printim G:\20090730\mref\dat\cs090731'it'.png white x1000 y800'
endif
endif
'c'
it=it+1
endwhile
'disable print'

kakaxiao

你把两个数据分开处理，然后读取两个ctl，画在一张图里成不？

梳子

kakaxiao 发表于 2012-12-31 10:21

                               
登录/注册后可看大图

你把两个数据分开处理，然后读取两个ctl，画在一张图里成不？

试了一下，还是不行，为什么会错开呢？想不通。。

kakaxiao

梳子 发表于 2012-12-31 10:35

                               
登录/注册后可看大图

试了一下，还是不行，为什么会错开呢？想不通。。

你好像设置了两次坐标轴啊，只设置一次

river

set vrange 500 15000去掉试试

jacci_zou

地形高度需要跟位势高度对比一下，才能使地形在图形上显示出来比较平滑，具体操作，可以搜索下论坛，有一个关于任意两点剖面图的帖子有讲到！

梳子

jacci_zou 发表于 2012-12-31 17:37

                               
登录/注册后可看大图

地形高度需要跟位势高度对比一下，才能使地形在图形上显示出来比较平滑，具体操作，可以搜索下论坛，有一个 ...

位势高度？这个纵轴都是海拔高度，任意两点剖面的那个是用collect和gr2stn实现的，我没懂怎么对比这两个高度啊。。。能再说明一下吗？谢谢！

river

梳子 发表于 2013-1-1 22:24

                               
登录/注册后可看大图

位势高度？这个纵轴都是海拔高度，任意两点剖面的那个是用collect和gr2stn实现的，我没懂怎么对比这两个高 ...

如果你的hgt资料也是海拔高度的资料，不需要对比啊，直接画就行。如果不是的话用它减掉500hPa的多年平均值对比一下就行吧

njzqxt

应该是处理成一个数据吧。

梳子

river 发表于 2012-12-31 16:54

                               
登录/注册后可看大图

set vrange 500 15000去掉试试

去掉这个两个坐标轴就不是统一的范围了，海拔高度的轴会只在0-5000范围内。。。这个数据的垂直层数就是按照海拔高度来分的，地形数据也是海拔高度，两个单位一样，不过在ref这个变量中，高度是一个维数，而地形数据中，高度是变量。。。