地形图叠加出现错误

云一林

以下为中国地形的ctl：
dset D:\ncep\dx.nc
title
undef -9.99e+33
dtype netcdf
xdef 2011 levels 70 70.0333 70.0667 70.1 70.1333等

ydef 1201 levels 15 15.0333 15.0667 15.1等
zdef 1 linear 0 1
tdef 1 linear 00Z01JAN0001 1mn
vars 1
z=>z  0  y,x  z
endvars


地形的gs：
'reinit'
'reset'
'sdfopen d:\ncep\dx.nc'
'set grads off '
'set grid off'
'set lon 73 137'
'set lat 15 55'
define_colors
'set gxout shaded'
'set clevs 3000'
'set ccols 0 16 '
'd z'
'draw title (a)'
'printim d:\ncep\dx.png x1200 y1000 white'
'disable print'

能出图

热源的ctl：
dset d:\ncep\q1pj.dat
undef  -9.99e+33
xdef 144 linear   0 2.5
ydef 69 linear  -85 2.5
zdef 1 linear 1000   1000
tdef 12 linear 0Z1mar1979 1mo
vars 1
q1n  0  99 *Monthly Mean of Q1
endvars


热源的gs：
'reinit'
'reset'
'open d:\ncep\q1pj.ctl'
'set grads off '
'set grid off'
'set lon 73 137'
'set lat 15 55'
'set z 1'
'set gxout contour'
'set clevs  -60  -40  -20  0  20  40  60  80  100  120'
'set ccolor 1 '
'set cthick 4'
'd ave(q1n,t=1,t=3)'
'draw title (a)'
'printim d:\ncep\q1.png x1200 y1000 white'
'disable print'

能出图

两个数据合到一起后的gs：
'reinit'
'reset'
'set grads off '
'set grid off'
'sdfopen d:\ncep\dx.nc'
'open d:\ncep\q1pj.ctl'
'set lon 73 137'
'set lat 15 55'
'set z 1'
define_colors
'set gxout shaded'
'set clevs 3000'
'set ccols 0 16 '
'd z.1'
'set gxout contour'
'set clevs  -60  -40  -20  0  20  40  60  80  100  120'
'set ccolor 1 '
'set cthick 4'
'd ave(q1n.2,t=1,t=3)'
'draw title (a)'
'printim d:\ncep\z.png x1200 y1000 white'
'disable print'

grads报错：
Data Request Error: Invalid grid coordinates
  World coordinates convert to non-integer  grid coordinates
    Variable = q1n.2  Dimension = 3
Operation Error:  Error from ave function
  Error ocurred at column 1
DISPLAY error:  Invalid expression
  Expression = ave(q1n.2,t=1,t=3)

请问这个gs应该怎么改？因为地形本来就好像只有两维，没有高度层和时间层，而热源有一层高度和12个时间层？有没有什么办法将等值线直接叠加到阴影图上？

balfulosa

是是这样可以不'd z.1(t=1)'