求grads批量画图方法

lljly777

要画出122个站点的等值线图，对应有122个dat文件，但是画不出来啊，帮忙看看下面的文件，有什么问题？
gs文件
'reinit'
'open d:\lw\yngrd.ctl'
'open d:\lw\zd.ctl'
t=1
year=2001
while(year<=2003)
'set grads off'
'set grid off'
'set cthick 5'
'set xlopts 1 5 0.16'
'set ylopts 1 5 0.fi16'
'set t 't
'define ay=oacres(g,r.2,12,8,6,4)'
'define by=maskout(ay,g-0.5)'
'define cy=smth9(by)'
'set gxout contour'
'd cy'
'enable print d:\lw\'year'.gmf'
'run d:\lw\YN1.gs'
'print'
'disable print'
'c'
t=t+1
year=year+1
endwhile
'reinit'
ctl文件1
dset  d:\lw\bxgx\%y4.dat
options template
undef 999.9
Title R
dtype station
stnmap d:\lw\st.map
tdef 122 linear jun2001 1yr
vars 1
r 0 99 station data
endvars
ctl文件2
dset d:\lw\yngrd.dat
undef -999.0
title Sample GRIB Data on Yunnan Province by Zhaodi
xdef 101 linear 97 0.1
ydef 91 linear 21 0.1
zdef 1 linear 500 1
tdef 122 linear jun2001 1yr
vars  1
g   0  99  grid data prepared for oacres function
endvars

善人/jw

有什么错误提示呢

river

您这个不是批量不批量的问题。grads的站点作图是有一套规则的，不是您随意想一下就能画出来的。论坛里很多站点作图的例子，先学习一下吧

传说中的谁

你的t是怎么循环的呢？开始时1，结束呢？建议如果你想把t也循环的话，再写一重while（）.......endwhile吧
还有，按照你写的方式输出图，不管循环对不对，最后出来的只可能有2001.gmf、2002.gmf、2003.gmf

nuister

为什么看不到t的终止在哪。。

松下风

期待解决结果

liutaoz21

为啥感觉怪怪的。每个站点文件都是包含三年的数据？

wc619lf

你是说吧122个站点资料格点化吧？那也不叫批处理啊。。。。。

lljly777

就是有122个dat文件，每个文件画一张等值线图，所有dat文件用一个ctl文件来描述，dat文件的命名是按年份来排列的，主要是为了实现批处理

iij

lljly777 发表于 2013-7-7 23:06
就是有122个dat文件，每个文件画一张等值线图，所有dat文件用一个ctl文件来描述，dat文件的命名是按年份来 ...

楼主，请问你有解决么？我也想要做一个这样的批量绘图。