利用maskout突出显示某省份，遇到错误

tiantianyixin

本人想将内蒙古相对湿度突出显示，因而使用了MASKOUT函数。但是运行grads，错误提示为无法打开maskout的ctl文件，当我单独打开maskout的ctl文件时没有问题的，为什么放在这个gs文件下就不能运行了呢？是因为前三条命令是直接打开nc文件的原因吗？
另外，'set mpdset neimeng'，这个语句对吗？


'reinit'
'sdfopen G:\ncepdaily\rhum.2012.nc'
'sdfopen G:\ncepdaily\uwnd.2012.nc'
'sdfopen G:\ncepdaily\vwnd.2012.nc'
'OPEN E:\haze\maskout.ctl'
'enable print E:\haze\20120105_850_rh.gmf white'
st= 17
while(st<=17)
'set t   ' st
'set grid off'
'set grads off'
'set csmooth on'
'set lon 90 130'
'set lat 30 55'
'set xlopts 1 4 0.18'
'set ylopts 1 4 0.18'
'set xlint 5'
'set ylint 5'
'set lev 850'
'set mpdset neimeng'
'define a=oacres(mask.4,rhum)'

'set gxout contour'
'd smth9(maskout(a,mask.4))'
'set lon 90 130'
'set lat 30 55'
'set ccolor 1'
'set cthick 4'
'set gxout barb'
'set digsize 0.06'
'd skip(uwnd.2*2.5,1,1);skip(vwnd.3*2.5,1,1)'

'set mpdset neimeng_q'
'draw map'
'writehz 1.5 7.8 850hPa'
'set strsiz 0.15'
'q time'
x=subwrd(result, 3)
'draw title 'x''
'print'
'c'
st=st+2
endwhile
'disable print'
'reinit'
;

river

既然说是打不开，那就看看路径什么的对不对。错误截图为什么不发上来呢？
set mpdset neimeng' 本身不存在问题，容易有问题的是你的相应文件件下有没有neimeng这个地图文件

tiantianyixin

river 发表于 2014-9-1 12:34
既然说是打不开，那就看看路径什么的对不对。错误截图为什么不发上来呢？
set mpdset neimeng' 本身不存在 ...

好的，我在调整一下，多谢啦