关于绘制站点文件的一点问题

mrxiaoan

通过学习用grads绘制站点文件等值线图的，收获很大，但现在一直还存在一个问题。教程gs原文如下：
'reinit'
'open e:\program\example\sta2grd\grid.ctl'
'open e:\program\example\sta2grd\sta.ctl'
'set lon 70 140'
'set lat 15 55'
'set mpdset hires'
'd  oacres(g,rain.2)'
'printim e:\program\example\sta2grd\sta_grid.png white'

这里grid.ctl是格点背景场，全部都是1的值；sta.ctl就是站点要素值。但按照教程gs这样的顺序，一直出不了图，显示错误如下：
但如果将grid.ctl和sta.ctl打开顺序调换下，先打开sta.ctl，再打开grid.ctl，再把'd  oacres(g,rain.2)'改成'd  oacres(g.2,rain.1)'就可以画出图来了。
以前画单张图还没什么问题，顺序不一样就不一样。但现在要连续绘制多张图，就有点麻烦了，因为grads只能关闭最后一个文件。我就无法对sta.ctl进行循环读取了，这会是什么问题呢？

做个霸气的木头

因为原始站点数据要插值到格点上面，所以必须先打开格点文件，再打开站点文件，这个顺序不能变。绘制多张图也可以，只是你是哪种多张图呢，如果是同一个变量不同时间，你可以全部插值完成后再循环绘制图形，就像下面这样的（1-600个时间点的全部插值完再循环画）。
'set lat 15 55'
'set lon 70 140'
'set t 1 600'
'define a=oacres(g(t=1),pre.2)'
'define aa=smth9(a)'
'set t 1 12'
'save=ave(aa,t+0,t=600,12)'
'modify save seasonal'
'set t 1 600'
'have=aa-save'
i=1
while(i<=12)
'set grads off'
'set t 'i''
'set gxout shaded'
'cnbasemap save'
'cbarn'
'draw title pre_month-'%i%''
'printim c:\pre_yue_ave_'%i%'.png white'
'c'
i=i+1
endwhile

river

木头的回答很详细了······

mrxiaoan

好啊，谢谢大神的回复。我赶紧试下看

香草拿铁

顺便学习了
多谢