关于province-basemap出图色标出错求教

子龙

   我利用站点资料绘制降水量图，利用province-basemap 将不需要区域白化，可色标出现错误，我的白色色标设置的是0以下，可这个色标绘制20以下，我的gs文件色标设置没有错，求教，是不是province-basemap这个gs文件出错呢？
gs文件如下：
'reinit'
'open D:\20160724\grid.ctl'
'open D:\20160724\2016072400.ctl'
'set grads off'
'set grid off'
palette()
'set xlopts 1 4 0.14'
'set ylopts 1 4 0.14'
'set clopts 1 4 0.12'
'set lat 31 40'
'set lon 103 113'
'set z 1'
'set t 1'
*突出显示陕西地图
'set mpdset shannx '
'set poli on'
'set map 0 1 13'
*计算散度
'define a=oacres(g.1,u.2,10,7,4,2,1)'
'define b=oacres(g.1,v.2,10,7,4,2,1)'
'define d1=hdivg(a,b)'
*差值物理量
'define ct=oacres(g.1,t.2,10,7,4,2,1)'
'define ctd=oacres(g.1,td.2,10,7,4,2,1)'
'define cetd=oacres(g.1,etd.2,10,7,4,2,1)'
'define cprs=oacres(g.1,prs.2,10,7,4,2,1)'
'define cp=oacres(g.1,rr.2,10,7,4,2,1)'
'define cctd=ct-ctd'   
*水汽通量散度
'define  qcetd=1000*0.622*cetd/(cprs-0.378*cetd)'
'define q1=qcetd*a/9.8'
'define q2=qcetd*b/9.8'
'define qfdiv=hdivg(q1,q2)*1e6'
*辐合上升和湿区和净辐合区计算
'define p1=const(maskout(cp,5-cctd),0,-u)'      
'define p2=const(maskout(p1,0.0-d1),0,-u)'
'define p3=const(maskout(p2,0.0-qfdiv),0,-u)'
'define p4=smth9(p3)'
'set xlint 2'
'set ylint 2'
*绘图设置
'set gxout shaded'
'set clevs 0.1 10 25 50 100 250'
'set rbcols 40 41 42 43 44 45 46'
*'d smth9(qfdiv)'
'd p4'
'province-basemap shannx p4 0'
*'set gxout stream'
*'set ccolor 24'
*'set arrscl 0.1'
*'set arrowhead 0.02'
*'set strmden 5'
*'d u.2*2.5;v.2*2.5'
'cbarn  1 1'
'printim   D:\20160724\2016072408p.png  white'
'reinit'
function palette()
'set rgb 40 255 255 255'
'set rgb 41 166 242 142'
'set rgb 42 61 185 61'
'set rgb 43 97 184 255'
'set rgb 44 1 0 224'
'set rgb 45 250 0 250'
'set rgb 46 129 0 64'

四叶草

'd p4'
'province-basemap shannx p4 0'
把‘d  p4’注释看看

river

直接用'province-basemap shannx p4 0'，province-basemap里面已经包含了display