walker环流 速度势 比湿想叠加在一张图上

nuistajun

跪求各位大神给小弟看看脚本该怎么修改，因为我的shum只有8层，而速度势和walker环流都有12层，三种要素想叠加在一张图上。。。
'reinit'
'open C:\Users\Lenovo\Desktop\wind\positive\shum1.ctl'
'open C:\Users\Lenovo\Desktop\wind\positive\shum3.ctl'
'set x 1 144'
'set y 1 72 '
'set z 1 8'
'set annot 8 8'
'set mpdraw on'
'set map auto'
'set mpdset cnworld'
'set grads off'
'set grid off'
'set gxout shaded'
'set lon 50 300'
'set lat 0'
'set z 1 8'
'set zlog on'
'define oo=shum.1-shum.2'
'd ave(oo,lat=-10,lat=20)'
'cbarn'
pull dummy


'open C:\Users\Lenovo\Desktop\wind\positive\u1.ctl'
'open C:\Users\Lenovo\Desktop\wind\positive\v1.ctl'
'open C:\Users\Lenovo\Desktop\wind\positive\w1.ctl'
'open C:\Users\Lenovo\Desktop\wind\positive\u3.ctl'
'open C:\Users\Lenovo\Desktop\wind\positive\v3.ctl'
'open C:\Users\Lenovo\Desktop\wind\positive\w3.ctl'
'set x 1 144'
'set y 1 72 '
'set z 1 12'
*'set annot 8 8'
*'set mpdraw on'
*'set map auto'
*'set mpdset cnworld'
*'set grads off'
*'set grid off'
'define aa=uwnd.3-uwnd.6'
'define bb=vwnd.4-vwnd.7'
'define chi=fish_chi(aa,bb)'
*'define chi=fish_chi(uwnd.4,vwnd.5)'
*pull dummy
'chi=chi-aave(chi,global)'
'set gxout contour'
'set cthick 9'
'set lon 50 300'
'set lat 0'
'set z 1 12'
'set zlog on'
'd ave(chi/1e6,lat=-10,lat=20)'
*pull dummy
'set gxout vector'
'set lon 50 300'
'set lat 0'
'set z 1 12'
'set zlog on'
'define cc=uwnd.3-uwnd.6'
'define dd=omega.5-omega.8'
'd ave(cc,lat=-10,lat=20); ave(dd,lat=-10,lat=20)'
*'d ave(uwnd.1,lat=-10,lat=20);ave(omega.3,lat=-10,lat=20)'
*pull dummy
'printim C:\Users\Lenovo\Desktop\wind\positive\grid('b').png white'
;

Ｏ2囡

顶~各位大神帮帮看看

小傅之追梦

可以把CTL贴出来看看吗？

nuistajun

小傅之追梦 发表于 2014-3-30 18:13
可以把CTL贴出来看看吗？

dset C:\Users\Lenovo\Desktop\wind\positive\shum1.grd
title u2  sample
undef -9.99e+33
xdef 144 linear 0 2.5
ydef 73  linear -90 2.5
zdef 8 levels 1000 925 850 700 600 500 400 300
tdef 1 linear jun1971 1yr
vars 1
shum 8 99 grid data
endvars

dset C:\Users\Lenovo\Desktop\wind\positive\u1.grd
title u2  sample
undef -9.99e+33
xdef 144 linear 0 2.5
ydef 73  linear -90 2.5
zdef 12 levels 1000 925 850 700 600 500 400 300 250 200 150 100
tdef 1 linear jun1971 1yr
vars 1
uwnd 12 99 grid data
endvars

nuistajun

Ｏ2囡 发表于 2014-3-30 17:50
顶~各位大神帮帮看看

C:\Users\Lenovo\Desktop\上传

小傅之追梦

nuistajun 发表于 2014-3-30 18:19
dset C:\Users\Lenovo\Desktop\wind\positive\shum1.grd
title u2  sample
undef -9.99e+33

1000 925 850 700 600 500 400 300 这八个层次是所有资料都有的，你可以在这几个层次上进行叠加啊

lqouc

你把画图两部分的顺序颠倒一下试试，只要改动变量后面的文件号，其余不变。

nuistajun

lqouc 发表于 2014-3-30 19:27
你把画图两部分的顺序颠倒一下试试，只要改动变量后面的文件号，其余不变。

好的，我试试。。

nuistajun

lqouc 发表于 2014-3-30 19:27
你把画图两部分的顺序颠倒一下试试，只要改动变量后面的文件号，其余不变。

不行呢大神，这样shaded把前面都给覆盖了

lqouc

nuistajun 发表于 2014-3-30 20:32
不行呢大神，这样shaded把前面都给覆盖了

额，没注意你是填色，话说可以试试画图顺序不变，只是改变open文件的顺序，先open那个12层的。