可以在阴影图上叠加阴影图吗？为什么我的出不来？

schliezer

打算先画假相当位温垂直剖面阴影图，再叠加上地形，可输出的GIF图片上只有黑色的地形，
后来试了试颠倒顺序，先画地形再画假相当位温，结果恰好相反，GIF图片上只有假相当位温，
附gs如下：

'reinit'
'open D:\gradsdata\ncepfnl_1_6h\dz.ctl'
'open D:\gradsdata\ncepfnl_1_6h\fnl_20120728_00_00_c.ctl'
'set display greyscale white'
'set grid off'
'set grads off'
'set csmooth on'
'set mpdset cnworld'
'set lat 30 45'
'set lon 105'
tt=1
while(tt<=13)
'c'
'set grads off'
'set dfile 2'
'set t 'tt' '
'set z 1 21'
'define prs=lev.2'
'define es=(6.112*exp(17.67*(tmpprs.2-273.15)/(tmpprs.2-29.65)))'
'define q=RHprs.2*(0.62197*es/(prs-es))/100.'
'define e=prs*q/(0.62197+q)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(tmpprs.2)-log(e)-4.805)'
'define theta=tmpprs.2*pow((1000/prs),(0.2854*(1.0-0.28*q)))'
'define eqt=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))'
'set gxout shaded'
'set cint 4'
'set rbcols auto'
'd eqt'
'cbarn 1 1 0.2 4.5'

'set dfile 1'
'set t 1'
'set z 1 12'
'set gxout shaded'
'set csmooth on'
'set rbcols 1'
'set cmax 0'
'define zz=dz.1'
'd zz'
'printim D:\gradsdata\ncepfnl_1_6h\0728图\不稳定条件\假相当位温纬度-高度'%tt%'.gif  white'
tt=tt+1
endwhile
;

其中dz是地形文件，其ctl：
DSET  D:\gradsdata\ncepfnl_1_6h\dz.dat
TITLE  dz
UNDEF  -9.99E33
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  JAN2006 1yr
VARS   1
dz     12  99  dz
ENDVARS

不知道是不是因为
'set gxout shaded'
'set csmooth on'
'set rbcols 1'
'set cmax 0'
'define zz=dz.1'
'd zz'
这个步骤里，地形大于0时grads的颜色是白色而不是空白，所以盖住了先画的shaded图？
请教各位又该如何处理呢？

river

这个需要利用maskout来画

schliezer

river 发表于 2014-2-17 17:59

                               
登录/注册后可看大图

这个需要利用maskout来画

maskout也试过了，把'd eqt' 改成 'd maskout(eqt,zz)' ，但还是不行，出图只有地形，提示错误来自maskout函数，且网格不匹配（incompatible grid）。我觉得原因有可能是eqt的垂直层次（21层）与zz的层次（12层）不匹配造成的，不知道怎样选出eqt与zz共有的层次？或者还有什么更好的办法？谢谢你了！

river

schliezer 发表于 2014-2-18 10:55

                               
登录/注册后可看大图

maskout也试过了，把'd eqt' 改成 'd maskout(eqt,zz)' ，但还是不行，出图只有地形，提示错误来自maskou ...

还是去看一下官方文档，研究一下maskout的用法吧

idfall0803

可以試地形一個顏色，假相当位温是地形以外的顏色，應該是可以

山_脉

我看好多人在说要把DZ用re()函数重新插值后才能使用.也不知道你现在弄出来没有.