请教大神，图中黑白外框如何显示出来，该用什么脚本

yulong0631

请教大神，图中黑白外框如何显示出来，该用什么脚本命令，原帖应该是没有将画此部分的脚本列上，请指教

扇形图的边框明显使用GRADS制作的，请教如何完成？

原帖地址：http://bbs.06climate.com/forum.php?mod=viewthread&tid=29853gs文件如下，请指点：
'run rgbset2'
'reinit'
'set grid on'
'set grads off'

'open year.ctl'
'sdfopen NA1901.2011.tmp.djf.nc'
'set grid on'
'set grads off'
'set poli on'


'set poli on'
'set mpdraw on'

'set vpage 0 11 0 8.5'

'set parea 1 10 1 8'
'set mpdset hires'
'set frame circle'
'set mproj nps'

'set lon -140 -60'
'set lat  15 75'
'set xlint 10'
'set ylint 10'

'set ylopts 1 4 0.20'
'set xlopts 1 4 0.20'

'set gxout shaded'
'set csmooth on'

'set rbcols 49 48 47 46 45 44 43 42 22 23 24 25 26 27 28 29'
'set clevs -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5'


'set t 1 110'
'ltrend djft.2 air slope rmse'
'set t 1'
'd slope*10'
'run cbarn'
'run gmt.gs 0.1 10 10'

'set x 1'
'set y 1'
'set z 1'
'set t 1 110'
'define year= aave(var.1,global)'
'set lon -140 -60'
'set lat  15 75'

'set z 1'
'set t 1'
'define rr=tcorr(year,djft.2, t=1, t=110)'
*'d rr'
********************************test the significance of correlation correlation coefficient （rr）
* degree of freedom=110-2=108 (from 1902 to 2011), 95% significance (both sides)=0.195
* degree of freedom=110-2=108 (from 1902 to 2011), 90% significance (both sides)=0.164

'define tt=rr*sqrt(108/(1-pow(rr,2)))'
'define ts=abs(tt)-0.195'

'set gxout shp'
'set shp -pt shppt'
'd skip(maskout(tt,ts),2,2)'
'set shpopts -1 2 0.001'
'draw shp shppt'
'drawbox -95 -75 30 40'
'drawbox -130 -100 45 65'

river

你这基本跟别人的一样，出不来别人这个框，难道不会报错吗？

言深深

你确定下面的不是IDL绘图的吗？

Leslie

这个黑白外框idl可以做出来，但是帖子中的命令不是idl语言。

yulong0631

扇形图的边框明显使用GRADS制作的，请教如何完成？

yulong0631

有同学知道吗？求指教

yulong0631

有同学知道吗？求指教

喵咪O(∩_∩)O~

谁有drawbox.gs脚本~求脚本