爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 8071|回复: 18

[分享资料] 多文件循环~求助

[复制链接]

新浪微博达人勋

发表于 2012-10-16 13:28:33 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 Eee萱 于 2012-10-16 13:30 编辑

想做一个while循环,同时输出12个月数据的叠加图,gs文件如下:
'reinit'
'open /cluster/g/ctl/so4.ctl'
'open /cluster/g/ctl/u_G4.ctl'
'open /cluster/g/ctl/v_G4.ctl'
i=1
while(i<=12)
'set dfile 1'
'set x 1 144'
'set y 1 91'
'set z 14'
'set t'i''
'set gxout shaded'
'd SO4'
'q time'
time=subwrd(result,3)
date=substr(time,4,9)
'run /cluster/g/grads/grads-2.0.1/lib/cbarn_interp.gs 2'
'draw title Sulfate_200hpa_'date'[ug/m3]'
'enable print so4_200hpa_200401-12.gmf'
'set dfile 2'
'set gxout vector'
'set x 1 144''set y 1 91'
'set z 9'
'set t'i''
'set arrscl 0.15 8'
'd Skip(u.2,2);v.3'
'print'
'clear'
i=i+1
endwhile
'disable print'
但是总是出现如下错误:
SET error: Invalid operand  
Operand = t1
SET error: Invalid operand  
Operand = t1
Data Request Error: Invalid grid coordinates  
World coordinates convert to non-integer  grid coordinates   
Variable = u.2  Dimension = 3
Operation Error:  Error from skip function  
Error ocurred at column 1DISPLAY error:  
Invalid expression   Expression = skip(u.2,2)
不用循环是可以画图的,数据应该没什么问题,是我循环设置哪里出错了呢?拜托各位大家了!谢谢谢谢!


密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-10-16 13:38:26 | 显示全部楼层
t和'i'中间不用加空格么?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-10-16 13:52:56 | 显示全部楼层
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-10-16 14:04:36 | 显示全部楼层
我试了下前半部分的循环,中间加空格的话SET error: Invalid operand  
Operand = t1
会消去;
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-10-16 14:10:10 | 显示全部楼层
'set x 1 144''set y 1 91'
把这两句分两行写试试吧
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-10-16 16:47:56 | 显示全部楼层
风飘摇 发表于 2012-10-16 14:10
'set x 1 144''set y 1 91'
把这两句分两行写试试吧

这个文件里是分行写的哦,复制过来不知道怎么就这样了..
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-10-16 17:11:52 | 显示全部楼层
本帖最后由 river 于 2012-10-16 17:13 编辑

一楼说得对,是'set t 'i''或者'set t 'i  t和'i'之间的空格一定要加的。你一共有两个set t 'i' 。记着都要加空格啊
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-10-16 17:13:43 | 显示全部楼层
'set dfile 2' 以后,u v就不是u.2 v.3了……
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-10-16 17:22:53 | 显示全部楼层
改过之后还是不行啊,错误是一样的...
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 19710
发表于 2012-10-16 17:33:12 | 显示全部楼层
本帖最后由 兰溪之水 于 2012-10-16 17:34 编辑

'reinit'
'open /cluster/g/ctl/so4.ctl'
'open /cluster/g/ctl/u_G4.ctl'
'open /cluster/g/ctl/v_G4.ctl'

'enable print so4_200hpa_200401-12.gmf'

i=1
while(i<=12)
'set x 1 144'
'set y 1 91'
'set z 14'
'set t 'i''
'set gxout shaded'
'd SO4'
'q time'
time=subwrd(result,3)
date=substr(time,4,9)
'run /cluster/g/grads/grads-2.0.1/lib/cbarn_interp.gs 2'
'draw title Sulfate_200hpa_'date'[ug/m3]'

'set gxout vector'
'set x 1 144'
'set y 1 91'
'set z 9'
'set t 'i''
'set arrscl 0.15 8'
'd skip(u.2,2);v.3'
'print'
'clear'
i=i+1
endwhile
'disable print'
;
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表