grads的问题

虎跃平阳

run gs文件后会出现unknown command 0,但是计算出来的结果是对的，这是什么原因呢？
这是其中的一部分gs

river

大雨的小伙伴了，怎么标题还写成这样······gs直接贴上来吧

虎跃平阳

river 发表于 2014-1-4 12:07

                               
登录/注册后可看大图

大雨的小伙伴了，怎么标题还写成这样······gs直接贴上来吧

gs在上面的

传说中的谁

想看看的，可是要花钱下载附件

river

虎跃平阳 发表于 2014-1-4 12:36

                               
登录/注册后可看大图

gs在上面的

我是说贴上来，像你帖子里那些字一样，不是传附件

虎跃平阳

传说中的谁 发表于 2014-1-4 12:46

                               
登录/注册后可看大图

想看看的，可是要花钱下载附件

莫急，我是嫌占屏幕，都是土豪还在乎这么点银子。
'reinit'
'open G:\data\era\1000u\1979-1989.ctl'
'open G:\data\era\1000u\1990-1999.ctl'
'open G:\data\era\1000u\2000-2011.ctl'
'set fwrite G:\work\trackdata\climate\circulation\data\e-6\1000u.grd'
'set gxout fwrite'
'set lon 70 140'
'set lat 10 50'
i=1
while(i<=31)
    ret=read('G:\work\trackdata\climate\circulation\time\e-6.txt')
rec=sublin(ret,2)
tt=subwrd(rec,3)
yr=subwrd(rec,1)
'set t 'tt''
if(yr<1990)
   'd (u1000mb(t='tt')+u1000mb(t='tt+1')+u1000mb(t='tt+2')+u1000mb(t='tt+3'))/4.0'
endif
if(yr>1989&yr<2000)
    'd (u1000mb.2(t='tt')+u1000mb.2(t='tt+1')+u1000mb.2(t='tt+2')+u1000mb.2(t='tt+3'))/4.0'
endif
if(yr>1999)
    'd (u1000mb.3(t='tt')+u1000mb.3(t='tt+1')+u1000mb.3(t='tt+2')+u1000mb.3(t='tt+3'))/4.0'
endif
i=i+1
endwhile
close('G:\work\trackdata\climate\circulation\time\e-6.txt')
'disable fwrite'
;
say 1000V
'reinit'
'open G:\data\era\1000v\1979-1989.ctl'
'open G:\data\era\1000v\1990-1999.ctl'
'open G:\data\era\1000v\2000-2011.ctl'
'set fwrite G:\work\trackdata\climate\circulation\data\e-6\1000v.grd'
'set gxout fwrite'
'set lon 70 140'
'set lat 10 50'
i=1
while(i<=31)
    ret=read('G:\work\trackdata\climate\circulation\time\e-6.txt')
rec=sublin(ret,2)
tt=subwrd(rec,3)
yr=subwrd(rec,1)
'set t 'tt''
if(yr<1990)
   'd (v1000mb(t='tt')+v1000mb(t='tt+1')+v1000mb(t='tt+2')+v1000mb(t='tt+3'))/4.0'
endif
if(yr>1989&yr<2000)
    'd (v1000mb.2(t='tt')+v1000mb.2(t='tt+1')+v1000mb.2(t='tt+2')+v1000mb.2(t='tt+3'))/4.0'
endif
if(yr>1999)
    'd (v1000mb.3(t='tt')+v1000mb.3(t='tt+1')+v1000mb.3(t='tt+2')+v1000mb.3(t='tt+3'))/4.0'
endif
i=i+1
endwhile
close('G:\work\trackdata\climate\circulation\time\e-6.txt')
'disable fwrite'
;

虎跃平阳

river 发表于 2014-1-4 17:12

                               
登录/注册后可看大图

我是说贴上来，像你帖子里那些字一样，不是传附件

在楼下了，只是其中的一部分，但能说明问题。

虎跃平阳

怎么跑到楼下了，sorry，请看6楼。

虎跃平阳

gs的目的是：由于数据太大，分段下载的，txt里面有挑出来的时次，从数据中选出相应时次的数据。

river

虎跃平阳 发表于 2014-1-4 22:02

                               
登录/注册后可看大图

在楼下了，只是其中的一部分，但能说明问题。

这个有点高级了，我的grads能力不足了·····等高人吧