解决错误data request warning:request beyond file limit

cjjj

有资料uwnd.2003.nc和uwnd.2005.nc
经查询uwnd.2003.nc的ctl是
dset xx
title xx
undef XX
xdef   144 linear    0  2.5
ydef    73 linear  -90  2.5
zdef    17 levels  1000 950 850 700 600 500 400 300
200 150 100 70 50 30 20 10   
tdef   365    linear   00Z0JAN2003 1440mn   
vars 1
uwnd 17 -999 mean daily U-wind

uwnd.2005.nc的ctl与上者仅一行不同
tdef   365    linear   00Z0JAN2005 1440mn

欲画出这两年平均U风场动能的经向平均
gs如下：
'reinit'
*'enable print F:\uwnd\ku500hPa.gmf'

'sdfopen F:\uwnd.2003.nc'
'sdfopen F:\uwnd.2005.nc'
'set lat -90 90'
'set lev 500'
'define k=(ave(uwnd.1,t=1,t=365)*ave(uwnd.1,t=1,t=365)/2+ave(uwnd.2,t=1,t=365)*ave(uwnd.2,t=1,t=365)/2)/2'

'd ave(k,lon=0,lon=360)'

*'print'
*'disable print'
;

Grads运行后报错如下：
data request warning:request beyond file limit
cannot plot data-all undefined values

请问应该怎么改呢？谢谢

小傅之追梦

我觉得'd ave(k,lon=0,lon=360)'有问题的吧，应该改成'd ave(k,lon=0,lon=360，-b)'。但是你的问题应该主要不是这个问题~~~

小傅之追梦

主要问题还是在于你的两个数据的时间不一致的原因吧?
'define k=(ave(uwnd.1,t=1,t=365)*ave(uwnd.1,t=1,t=365)/2+ave(uwnd.2,t=1,t=365)*ave(uwnd.2,t=1,t=365)/2)/2'

river

在'd ave(k,lon=0,lon=360)' 的前面给经度随便固定一个值试试，比如set lon 1

cjjj

river 发表于 2013-5-20 21:15

                               
登录/注册后可看大图

我好像建议你加上的是set lon 1 而不是setlon1吧。我的空格呢？
还有另外一个建议，你这个是不 ...

谢谢 改成gs如下：
'reinit'
*'enable print F:\uwnd\ku500hPa.gmf'

'sdfopen F:\uwnd.2003.nc'
'sdfopen F:\uwnd.2005.nc'
'set lat -90 90'
'set lev 500'

define aveu1=ave(ave(uwnd.1,t=1,t=365),lon=0,lon=360)
define aveu2=ave(ave(uwnd.2,t=1,t=365),lon=0,lon=360)
define k=(aveu1**2/2+aveu2**2/2)/2
d k

*'print'
*'disable print'
;

报错：
unmatched parents on function call
error occurred processing function arguments
error occurred on line 9（gs中是：define aveu1=ave(ave(uwnd.1,t=1,t=365),lon=0,lon=360)）

cjjj

小傅之追梦 发表于 2013-5-19 16:53

                               
登录/注册后可看大图

我觉得'd ave(k,lon=0,lon=360)'有问题的吧，应该改成'd ave(k,lon=0,lon=360，-b)'。但是你的问题应该主要 ...

确实不是这个问题。。。不过还是个人很感谢你~~

cjjj

river 发表于 2013-5-19 17:57

                               
登录/注册后可看大图

在'd ave(k,lon=0,lon=360)' 的前面给经度随便固定一个值试试，比如set lon 1

还是不行 报错中还加了一条：unkown command:setlon1
这两个NC单独画图，都是可以出图的。。。

river

cjjj 发表于 2013-5-20 18:13

                               
登录/注册后可看大图

还是不行 报错中还加了一条：unkown command:setlon1
这两个NC单独画图，都是可以出图的。。。

我好像建议你加上的是set lon 1 而不是setlon1吧。我的空格呢？
还有另外一个建议，你这个是不是可以先分别求时间平均和精度平均，然后再利用K的那个公式也行吧
define aveu1=ave(ave(uwnd.1,t=1,t=365),lon=0,lon=360)
define aveu2=ave(ave(uwnd.2,t=1,t=365),lon=0,lon=360)
define k=(aveu1**2/2+aveu2**2/2)/2
d k

cjjj

river 发表于 2013-5-20 21:15

                               
登录/注册后可看大图

我好像建议你加上的是set lon 1 而不是setlon1吧。我的空格呢？
还有另外一个建议，你这个是不 ...

写命令的时候是写了空格的~报错的时候就没有空格了
好的 我来试试 快断网了 下次讲~谢谢谢谢啦

river

cjjj 发表于 2013-5-21 15:15

                               
登录/注册后可看大图

谢谢 改成gs如下：
'reinit'
*'enable print F:\uwnd\ku500hPa.gmf'

难道非得我从头到尾写一遍你才能明白啊，我只是给你一个新的思路，你自己放进gs里面是不是应该把应该加的都加上啊，单引号都哪去了啊。