用grads的maskout语句画西南地区降水站点数据

云一林

Grads降水站点数据只画西南地区出错，请问谁知道那个maskout语句到底怎么个用法阿？

站点插值成格点已经画出图来，画的是全中国的，以下是所有资料。

降水站点数据是地面层76个站的1961-2017年57年的逐年降水数据。

站点数据grd的ctl：

DSET E:\lunwen\3.grd

DTYPE station

STNMAP E:\lunwen\3.map

UNDEF -32766

TITLE Rain Data Sample

TDEF 57 linear 01jan1961 1yr

VARS 1

sd 0 99 Rain Data

ENDVARS

Fortran生成的格点背景场的ctl：

DSET E:\lunwen\grid.grd

TITLE Grid Data Sample

UNDEF -9.99E33

XDEF 29 LINEAR 70.0 2.5

YDEF 17 LINEAR 15.0 2.5

ZDEF 1 LEVELS 1000

TDEF 1 LINEAR 01jan1961 1yr

VARS 1

g 0 99 Grid Data

ENDVARS

画全中国的gs：

'reinit'

'open e:\lunwen\grid.ctl'

'open e:\lunwen\3.ctl'

'set lon 70 140'

'set lat 15 55'

'set t 12'

'set grads off'

'set mpdset hires'

'define rgrid=oacres(g,sd.2)'

'set gxout shaded'

'd rgrid'

'set gxout contour'

'cnbasemap rgrid'

'cbarn'

'printim e:\lunwen\1972.png white'

图片：

现在需要只画西南地区的降水站点数据，从MeteoInfo根据论坛步骤制作好西南地区地图。

西南地区地图的ctl：

DSET E:\lunwen\map.dat

TITLE Mask data

UNDEF -9999

XDEF 29 LINEAR 70 2.5

YDEF 17 LINEAR 15 2.5

ZDEF 1 LINEAR 1 1

TDEF 57 linear 01jan1961 1yr

VARS 1

mask 0 99  background mask data

ENDVARS

西南地区画图gs：

'reinit'

'open e:\lunwen\3.ctl'

'open e:\lunwen\map.ctl'

'set lon 70 140'

'set lat 15 55'

'set t 12'

'set mpdset hires'

'define rgrid=oacres(mask.2,sd.1)'

'set gxout shaded'

'd maskout(rgrid,mask.2)'

'set gxout contour'

'd maskout(rgrid,mask.2)'

'cbarn'

'printim e:\lunwen\xinan.png white'

Grads运行结果：

Define error:  Define not yet valid for station data

    Default file is a station data file

Syntax Error:  Invalid Operand

  'rgrid' not a variable or function name

Operation Error:  Error from maskout function

  Error ocurred at column 1

DISPLAY error:  Invalid expression

  Expression = maskout(rgrid,mask.2)

Syntax Error:  Invalid Operand

  'rgrid' not a variable or function name

Operation Error:  Error from maskout function

  Error ocurred at column 1

DISPLAY error:  Invalid expression

  Expression = maskout(rgrid,mask.2)

我的问题：我把站点数据3.grd的ctl时间那一行贴到map的ctl里了，还有就是已经有西南地区的地图了，画图时还需要fortran生成的背景场吗？是打开3个ctl画图还是两个ctl就行了？？最后谁能帮我看看那个西南地区的gs怎么改啊？有点看不太懂论坛那个，看着好像maskout只用上了两个ctl，但是打开了3个ctl。