maskout画阴影出现这种情况，求大神指点

猫猫monty

简单直接粗暴的只上gs和图啦，求大神指点！！！
gs：
'reinit'
'open E:\data\ncep_winter_ok\ave_climate\all_level_climate\ave_air_2015.12_2016.01.ctl'
'open E:\data\ncep_winter_ok\2016\2015.12_2016.01\ano\ano_vwnd_2015.12_2016.01.ctl'
'open E:\data\NCEP2\orog.ctl'

'set parea 0.6 8 6.2 10.2'
'set grads off'
'set grid off'
'set xlint 10'
'set font 1'
'set xlopts 1 4 0.16'
'set ylopts 1 4 0.16'
'set clopts 1 1 0.15'
'set mpdset cnworld'

*地形
'set dfile 3'
'set t 1'
'set lon -80'
'set lat 0 90'
'set ylevs 850 500 400 200 100 50'
'set z 1 14'
'define a=(orog.3)*(-1.0)'

'set t 1'
'set lon -80'
'set lat 0 90'
'set ylevs 850 500 400 200 100 50'
'set z 1 14'
'set gxout shaded'
'set cmin 0'
'set csmooth on'
'd orog.3'

'set dfile 1'
'set lat 0 90'
'set lon -180 180'
'set z 1 14'
'set t 1 62'
'define dtx = cdiff(air.1,x)'
'define dty = cdiff(air.1,y)'
'define dx = cdiff(lon,x)*3.1416/180'
'define dy = cdiff(lat,y)*3.1416/180'

'set dfile 2'
'set lat 0 90'
'set lon -180 180'
'set z 1 14'
'set t 1 62'

'define v=vwnd.2'
'define adv=(-1*(v*dty/dy ))/6.37e6'

'set lat 0 90'
'set lon -90'
'set z 1 14'
'set t 1 31'
'define adv1=ave(adv,lon=-100,lon=-80)'

*西经80——100度剖面12月平均

'set lat 0 90'
'set z 1 14'
'set t 1'
'define adv12=ave(adv1,t=1,t=31)'

'set gxout shaded'
'run D:\stratospheric\color\30colors.gs'
'set clevs -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7'
'set ccols 21 22 24 26 27 28 30 32 35 36 39 41 43 45 47 48'
'set csmooth on'
'd maskout(adv12*10000,a)'
'gxprint E:\data\tu\\poumian.500_ano.V_ave.T_tem_advection.eps white'
;


图：

river

既然地形可以直接display，那你就直接先画其他的变量，最后再画地形的阴影不就覆盖了么，还用什么maskout啊

river

既然地形可以直接display，那你就直接先画其他的变量，最后再画地形的阴影不就覆盖了么，还用什么maskout啊

猫猫monty

river 发表于 2016-3-11 21:21
既然地形可以直接display，那你就直接先画其他的变量，最后再画地形的阴影不就覆盖了么，还用什么maskout啊

直接画shaded之后再画阴影，阴影就把之前shaded图全覆盖了啊

river

猫猫monty 发表于 2016-3-11 22:44
直接画shaded之后再画阴影，阴影就把之前shaded图全覆盖了啊

哦，我给忘记了。
出现那种类似马赛克的东西可能是资料分辨率的问题，估计是你两种资料分辨率不同的原因吧

曦尘

http://ncl.ucar.edu/Applications/mask.shtml
官网里masking版块的第14个例子你可以参考一下，但是地形没有那么平滑。

balfulosa

那个地性文件好像只有2.5*2.5的分辨率，所以出现锯齿也是很正常的。

星雨

我觉着你可以尝试着把粗一点的地形资料叠加上去。