grads画uwnd十年平均

33521

我想用nc资料画uwnd这个变量的十年资料平均。为什么这种方法不行呢
'set lon 50 150'
'set lat 10 80'
'set t 1'
'set lev 500'
'set mpdset hires '

'define a=ave(uwnd.1,t=1,t=366)'
*'define b=ave(uwnd.2,t=1,t=365)'
*'define c=ave(uwnd.3,t=1,t=365)'
*'define d=ave(uwnd.4,t=1,t=365)'
*'define e=ave(uwnd.5,t=1,t=366)'
*'define f=ave(uwnd.6,t=1,t=365)'
*'define g=ave(uwnd.7,t=1,t=365)'
*'define h=ave(uwnd.8,t=1,t=365)'
*'define i=ave(uwnd.9,t=1,t=366)'
*'define g=ave(uwnd.10,t=1,t=365)'
*'define k=ave(uwnd.11,t=1,t=365)'


*'d (a+b+c+d+e+f+g+h+i+g+k)/11'
'd a'

'print'
'disable print'
;

错误提示如图，这是为什么呢

river

贴个完整的gs上来

lxlxllx89

再把q ctlinfo的信息贴上来

斥鷃

看这个错是你读取的超过你打开的文件的维数范围了~同楼上所述，贴个ctl吧~

33521

river 发表于 2013-9-19 13:17

                               
登录/注册后可看大图

贴个完整的gs上来

这个就是完整的gs，下面的工作还没做这里就卡住了，感觉uwnd。2，uwnd。3后面那些都没有值呢

33521

lxlxllx89 发表于 2013-9-19 15:19

                               
登录/注册后可看大图

再把q ctlinfo的信息贴上来

这个就是ctlinfo，感觉后面的那些uwnd。2，uwnd。3那些都没有值，这是咋回事呢

33521

斥鷃 发表于 2013-9-19 18:16

                               
登录/注册后可看大图

看这个错是你读取的超过你打开的文件的维数范围了~同楼上所述，贴个ctl吧~

。nc的资料我用的sdfopen，没有写ctl，是不是这样打开不对呢，感觉第一个能打开，后面的10个都没有值

river

33521 发表于 2013-9-20 10:16

                               
登录/注册后可看大图

这个就是完整的gs，下面的工作还没做这里就卡住了，感觉uwnd。2，uwnd。3后面那些都没有值呢

这个怎么就是完整的了，都没有打开的命令啊

33521

river 发表于 2013-9-20 18:38

                               
登录/注册后可看大图

这个怎么就是完整的了，都没有打开的命令啊

'reinit'

'enable print E:\CH\ncep_2000_2013\u_2000.2000.gmf'
'sdfopen E:\CH\ncep_2000_2013\uwnd.2000.nc'
'sdfopen E:\CH\ncep_2000_2013\uwnd.2001.nc'
'sdfopen E:\CH\ncep_2000_2013\uwnd.2002.nc'
'sdfopen E:\CH\ncep_2000_2013\uwnd.2003.nc'
'sdfopen E:\CH\ncep_2000_2013\uwnd.2004.nc'
'sdfopen E:\CH\ncep_2000_2013\uwnd.2005.nc'
'sdfopen E:\CH\ncep_2000_2013\uwnd.2006.nc'
'sdfopen E:\CH\ncep_2000_2013\uwnd.2007.nc'
'sdfopen E:\CH\ncep_2000_2013\uwnd.2008.nc'
'sdfopen E:\CH\ncep_2000_2013\uwnd.2009.nc'
'sdfopen E:\CH\ncep_2000_2013\uwnd.2010.nc'

'set lon 50 150'
'set lat 10 80'
'set t 1'
'set lev 500'
'set mpdset hires '

'define a=ave(uwnd.1,t=1,t=366)'
'define b=ave(uwnd.2,t=1,t=365)'
'define c=ave(uwnd.3,t=1,t=365)'
'define d=ave(uwnd.4,t=1,t=365)'
'define e=ave(uwnd.5,t=1,t=366)'
'define f=ave(uwnd.6,t=1,t=365)'
'define g=ave(uwnd.7,t=1,t=365)'
'define h=ave(uwnd.8,t=1,t=365)'
'define i=ave(uwnd.9,t=1,t=366)'
'define g=ave(uwnd.10,t=1,t=365)'
'define k=ave(uwnd.11,t=1,t=365)'


'd (a+b+c+d+e+f+g+h+i+g+k)/11'
*'d a'

'print'
'disable print'
;
不好意思，是这个

CHch

斥鷃 发表于 2013-9-19 18:16

                               
登录/注册后可看大图

看这个错是你读取的超过你打开的文件的维数范围了~同楼上所述，贴个ctl吧~

'reinit'

'enable print E:\CH\ncep_2000_2013\u_2000.2000.gmf'
'sdfopen E:\CH\ncep_2000_2013\uwnd.2000.nc'
'sdfopen E:\CH\ncep_2000_2013\uwnd.2001.nc'
'sdfopen E:\CH\ncep_2000_2013\uwnd.2002.nc'
'sdfopen E:\CH\ncep_2000_2013\uwnd.2003.nc'
'sdfopen E:\CH\ncep_2000_2013\uwnd.2004.nc'
'sdfopen E:\CH\ncep_2000_2013\uwnd.2005.nc'
'sdfopen E:\CH\ncep_2000_2013\uwnd.2006.nc'
'sdfopen E:\CH\ncep_2000_2013\uwnd.2007.nc'
'sdfopen E:\CH\ncep_2000_2013\uwnd.2008.nc'
'sdfopen E:\CH\ncep_2000_2013\uwnd.2009.nc'
'sdfopen E:\CH\ncep_2000_2013\uwnd.2010.nc'

'set lon 50 150'
'set lat 10 80'
'set t 1'
'set lev 500'
'set mpdset hires '

'define a=ave(uwnd.1,t=1,t=366)'
'define b=ave(uwnd.2,t=1,t=365)'
'define c=ave(uwnd.3,t=1,t=365)'
'define d=ave(uwnd.4,t=1,t=365)'
'define e=ave(uwnd.5,t=1,t=366)'
'define f=ave(uwnd.6,t=1,t=365)'
'define g=ave(uwnd.7,t=1,t=365)'
'define h=ave(uwnd.8,t=1,t=365)'
'define i=ave(uwnd.9,t=1,t=366)'
'define g=ave(uwnd.10,t=1,t=365)'
'define k=ave(uwnd.11,t=1,t=365)'


'd (a+b+c+d+e+f+g+h+i+g+k)/11'
*'d a'

'print'
'disable print'
;
完整的是这个