急！用grads处理awx文件遇到的问题

繁星一点

本人用grads处理FY-2E 1hr每次的TBB资料，ctl文件是依照http://bbs.06climate.com/forum.php?mod=viewthread&tid=2528&page=1#pid20008  所编，gs文件是自己编写。在处理单个文件的时候，能画出图来，但是在gs文件中做循环批处理文件时发生错误，错误提示 set error：invalid operand。请教各位高手错误原因。
一下是我在批处理时的gs文件和ctl文件：
***********gs文件**********************
'reinit'
'clear'
'open f:/FY2E/2010tbb.ctl'
'set grads off'
'set lat 15 55'
'set lon 70 140'
'set t 1 48'
'set mpdset cnhimap'
'set map 1 1 6'
'draw map'
t=1
while(t<49)
'set xlint -10'
'set ylint -10'
'set t' t
'set gxout shaded'
'set clevs 180 200 220 240 260 280 '
'set ccols  9 14 4 13 3  12 2 '
'set grads off'
'set grid off'
'd TBB+100'
'cbarn.gs'
'draw string 2.3 0.3 time='t*1''
if(rc='0')break;endif;
'enable print f:/FY2E/2010_tbb.gmf'
'print'
'clear'
t=t+1
endwhile
pull dummy
'disable print'
return

**************ctl文件**********
dset H:\AWX\06\FY2E_TBB_IR1_OTG_2010%m2%d2_%h200.AWX
title 2010 twpice fy2e_tbb
options template
options yrev
* yrev表示y轴反向
undef 0.0
xdef 1201 linear  45 0.1
ydef 1201 linear -60 0.1
zdef 1 linear 1000 1
tdef 288 linear 00z17JUN2010 1hr
vars  1
tbb   1    -1,40,1    Geopotential Heights (m)
endvars

冬日情愫+

你试试将'set t' t 中的t换一个字母代替，因为t在grads中默认代表时间维数。

空谷幽兰

去掉'set t 1 48'试试~

繁星一点

冬日情愫+ 发表于 2011-12-12 08:45

                               
登录/注册后可看大图

你试试将'set t' t 中的t换一个字母代替，因为t在grads中默认代表时间维数。

谢谢，问题已解决。去掉‘set t 1 48’, 将'set t' t 改成  'set t' tt，相应的循环变量也改成tt

繁星一点

空谷幽兰 发表于 2011-12-12 09:20

                               
登录/注册后可看大图

去掉'set t 1 48'试试~

谢谢，问题已解决。去掉‘set t 1 48’, 将'set t' t 改成  'set t' tt，相应的循环变量也改成tt

baggio

请问‘d tbb+100’是为什么？

繁星一点

baggio 发表于 2012-10-7 11:33

                               
登录/注册后可看大图

请问‘d tbb+100’是为什么？

这个是指云顶亮温（单位：K）
tbb+100  是根据论坛里大家的帖子来的
如果是d tbb  则反演出来的温度明显较实际温度低很多  加上100后与实际结果一致
至于为什么加100  我也不是很清楚
你可以参考论坛里其他人发的类似的帖子

baggio

繁星一点 发表于 2012-10-8 10:42

                               
登录/注册后可看大图

这个是指云顶亮温（单位：K）
tbb+100  是根据论坛里大家的帖子来的
如果是d tbb  则反演出来的温度明显 ...

tbb的单位是开尔文，转换为摄氏度的话是‘d tbb-273.15’,但micaps与实际有100k的差值，我昨天刚刚做出来，应该是tbb-273.15+100=tbb-173.15，这样就跟实际差不多吻合了。至于那100k的差值是怎么出来的，我也不知道！跪求高人解释！@mofangbao

繁星一点

baggio 发表于 2012-10-8 13:39

                               
登录/注册后可看大图

tbb的单位是开尔文，转换为摄氏度的话是‘d tbb-273.15’,但micaps与实际有100k的差值，我昨天刚刚做出来 ...

其实你已经解释清楚了呀
如果tbb单位用摄氏度的话  grads语句是 d tbb-173.15
所以如果用开尔文 那grads语句就是 d tbb-173.15+275.15 即 d tbb+100

baggio

呃...这位同志，开氏温度转换为摄氏温度是需要k-273.15，再加上100的差值，就该是k-173.15，而不是d tbb + 100，该是tbb-273.15+100=tbb-173.15