【求助】Grads批处理AWX文件问题

Crazy_Merlin

大家好，我想用Grads批量处理FY-2E晴空大气可降水量数据。目前使用了11个时次的数据。方法参考了坛子里的一些帖子。
http://bbs.06climate.com/forum.php?mod=viewthread&tid=2528&extra=&page=1
http://bbs.06climate.com/forum.php?mod=viewthread&tid=4978
程序在作图时不知道为什么，只能画第一个时次的图。程序报错如附件所示。望高手帮忙解决！谢谢！
还有一个小小的问题，怎么给图片加上标题（标题名和数据文件一样，如FY2E_TPW_MLT_OTG_20120817_0600）


下面是ctl和gs文件
dset I:\TDDOWNLOAD\TPW\20120826\FY2E_TPW_MLT_OTG_%y4%m2%d2_%h200.AWX
title FY-2E_TPW data
options template
options yrev
undef 99
*99为无效值
xdef 1001 linear  55 0.1
ydef 1001 linear -50 0.1
zdef 1 levels 850
tdef 11 linear 00Z25aug2012 3hr
vars 1
tpw  0 -1,40,1 FY-2E_TPW (mm)
endvars
________________________________________
'reinit'
'open I:\TDDOWNLOAD\TPW\TPW.ctl'
*'enable print I:\TDDOWNLOAD\TPW\TPW.gmf'
*'set mpdraw on'
'set poli on'
'set mpdset  mres hun anh jiangx hub sic guiz'
tt=1
while(tt<=11)
'set t' tt
'set gxout shaded'
'set cint 10'
'set cmin 0.0'
*'set cmax 75'
'set xlopts 1 2 0.14'
'set ylopts 1 2 0.14'
'set lon 75 130'
'set lat 15 50'
'set grads off'
'set grid on'
'set xlint 5'
'set ylint 5'
'I:\TDDOWNLOAD\TPW\rb_matlab.gs'
'd smth9(tpw)'
'I:\TDDOWNLOAD\TPW\cbar_matlab.gs 5 1 1'
*'draw title FY-2E_TPW'
*'draw string 2.3 0.3 time='tt'''
if(rc='0')break;endif;
'printim I:\TDDOWNLOAD\TPW\TPW.png white'
*'print'
*'clear'
tt=tt+1
endwhile
pull dummy
*'disable print'
return
;

跪求高手解决！

mofangbao

楼主把你的代码换行搞一下哈，你可以先贴到word里面，然后再贴过来

瑜之姬

{:eb517:｝太深奥

mofangbao

晕，又是粗心大意的一个
'set t' tt
改成
'set t 'tt
这个问题论坛不知道多少人问过了

Crazy_Merlin

mofangbao 发表于 2012-8-27 12:13

                               
登录/注册后可看大图

楼主把你的代码换行搞一下哈，你可以先贴到word里面，然后再贴过来

不好意思，忘了，已经更改！

Crazy_Merlin

mofangbao 发表于 2012-8-27 17:34

                               
登录/注册后可看大图

晕，又是粗心大意的一个
'set t' tt
改成

谢谢管理员~太爱你了！这个一直没注意过……

njzqxt

楼主，编辑下，把正确的代码放上吧。

Crazy_Merlin

njzqxt 发表于 2012-8-29 07:46

                               
登录/注册后可看大图

楼主，编辑下，把正确的代码放上吧。

正确代码：'set t' tt改成'set t 'tt 改变空格位置

'reinit'
'open I:\TDDOWNLOAD\TPW\TPW.ctl'
'set mpdraw on'
'set poli on'
'set mpdset  mres hun anh jiangx hub sic guiz'
*'draw map'
tt=1
while(tt<=11)
'set t 'tt
'set gxout shaded'
'set cint 10'
'set cmin 0.0'
*'set cmax 75'
'set xlopts 1 2 0.15'
'set ylopts 1 2 0.15'
'set lon 70 135'
'set lat 15 50'
'set grads off'
'set grid on'
'set xlint 5'
'set ylint 5'
'I:\TDDOWNLOAD\TPW\rb_matlab.gs'
'd smth9(tpw)'
'I:\TDDOWNLOAD\TPW\cbar_matlab.gs 5 1 0.8'
'draw title FY-2E_TPW_'tt''
*'draw rec 1 1 2 2'
*'draw string 2.3 0.3 time='tt''
if(rc='0')break;endif;
*'enable print I:\TDDOWNLOAD\TPW\TPW'tt'.gmf'
'printim I:\TDDOWNLOAD\TPW\TPW'tt'.png white'
*'print'
'clear'
tt=tt+1
endwhile
pull dummy
*'disable print'
return
;

Crazy_Merlin

mofangbao 发表于 2012-8-27 17:34

                               
登录/注册后可看大图

晕，又是粗心大意的一个
'set t' tt
改成

管理员同志，有个问题请教一下，怎么才能把AWX文件的文件名批量的赋给Grads生成的图片呢？比如说：title为FY-2E_TPW_MLT_OTG_20120826_0300 .8楼给了我修改后的gs文件。grads不太精通，自学的~谢谢啦！

mofangbao

Crazy_Merlin 发表于 2012-8-29 08:55

                               
登录/注册后可看大图

管理员同志，有个问题请教一下，怎么才能把AWX文件的文件名批量的赋给Grads生成的图片呢？比如说：title为 ...

如果你的ctl里面是批量描述的，这个有点小复杂，不过应该是可以实现的
使用 q ctlinfo获得你当前打开的文件信息，用subline结合subwrd取出文件名
当然批量描述的我没试过，也许取出来的可能是 %d这种新式，如果是这样，那就要再用q dims获得当前的时间信息，还是用上面的方法取出相应的字符串，再拼合成完整的文件名，如果获得的直接是把%d替换好的字符串那就省事了。
如果你不是批量描述的，那通过上面的第一步就能实现你的要求了
我只说了思路，具体实现要你自己摸索了，可能会有点小折腾，应该可以实现，实现了别忘了来分享代码