gs文件打开多个nc文件出了问题

hgsndg

求助大神急！我用gs脚本文件打开多个nc文件时报错
'reinit'
'sdfopen e:\aa\uwnd.sig995.1987.nc'
'sdfopen e:\aa\vwnd.sig995.1987.nc'
'sdfopen e:\aa\uwnd.sig995.day.1981-2010.ltm.nc'
'set grads off'
'set lon 105 130'
'set lat 0 20'
'set lev 0'
'set t 1'
'define uadn=ave(uwnd.1,t=130,t=135)-ave(uwnd.3,t=130,t=135)'
'define vadn=ave(vwnd.2,t=130,t=135)'
'set gxout vector'
'd uadn;vadn'
'printim E:\aa\picture\uwnd.png white'
;

报错内容为：“Data Request Warning:  Request is completely outside file limits
如果只打开前两个文件就没有问题，前两个都是vars 1,第三个是vars 2的。以下为ctl文件部分截图：
vars 2
uwnd=>uwnd  0  t,y,x  mean Long Term Mean Daily u-wind at sigma level 995
valid_yr_count=>valid_yr_count  0  t,y,x  count of non-missing values used in mean

请问该怎么解决这个问题，我听说可以用”set nfile n"来解决，是真的吗？到底该怎么用呢？

ldwloveyyz

你最好是发个第三个文件的CTL的截图，问题疑似出现在第三个文件上

ldwloveyyz

有没有查看3个文件的维度是否一样啊

hgsndg

ldwloveyyz 发表于 2019-4-19 06:47
有没有查看3个文件的维度是否一样啊

dset E:\aa\uwnd.sig995.day.1981-2010.ltm.nc
title daily ltm uwnd.sig995 from the NCEP Reanalysis
undef -9.96921e+36
dtype netcdf
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 1 linear 0 1
tdef 365 linear 00Z01JAN0001 1440mn
vars 2
uwnd=>uwnd  0  t,y,x  mean Long Term Mean Daily u-wind at sigma level 995
valid_yr_count=>valid_yr_count  0  t,y,x  count of non-missing values used in mean
endvars

hgsndg

ldwloveyyz 发表于 2019-4-19 06:47
有没有查看3个文件的维度是否一样啊

dset E:\aa\uwnd.sig995.day.1981-2010.ltm.nc
title daily ltm uwnd.sig995 from the NCEP Reanalysis
undef -9.96921e+36
dtype netcdf
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 1 linear 0 1
tdef 365 linear 00Z01JAN0001 1440mn
vars 2
uwnd=>uwnd  0  t,y,x  mean Long Term Mean Daily u-wind at sigma level 995
valid_yr_count=>valid_yr_count  0  t,y,x  count of non-missing values used in mean
endvars

hgsndg

第一个，第二个也差不多
ga-> q ctlinfo
dset E:\aa\uwnd.sig995.1987.nc
title mean daily NMC reanalysis
undef -9.96921e+36
dtype netcdf
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 1 linear 0 1
tdef 365 linear 00Z01JAN1987 1440mn
vars 1
uwnd=>uwnd  0  t,y,x  mean Daily u-wind
endvars
第三个文件
ga-> q ctlinfo
dset E:\aa\uwnd.sig995.day.1981-2010.ltm.nc
title daily ltm uwnd.sig995 from the NCEP Reanalysis
undef -9.96921e+36
dtype netcdf
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 1 linear 0 1
tdef 365 linear 00Z01JAN0001 1440mn
vars 2
uwnd=>uwnd  0  t,y,x  mean Long Term Mean Daily u-wind at sigma level 995
valid_yr_count=>valid_yr_count  0  t,y,x  count of non-missing values used in mean
endvars

ldwloveyyz

第一个，第二个也差不多
ga-> q ctlinfo
dset E:\aa\uwnd.sig995.1987.nc
title mean daily NMC reanalysis
undef -9.96921e+36
dtype netcdf
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 1 linear 0 1
tdef 365 linear 00Z01JAN1987 1440mn
vars 1
uwnd=>uwnd  0  t,y,x  mean Daily u-wind
endvars
第三个文件
ga-> q ctlinfo
dset E:\aa\uwnd.sig995.day.1981-2010.ltm.nc
title daily ltm uwnd.sig995 from the NCEP Reanalysis
undef -9.96921e+36
dtype netcdf
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 1 linear 0 1
tdef 365 linear 00Z01JAN0001 1440mn
vars 2
uwnd=>uwnd  0  t,y,x  mean Long Term Mean Daily u-wind at sigma level 995
valid_yr_count=>valid_yr_count  0  t,y,x  count of non-missing values used in mean
endvars
标红的地方时间维数明显不一样，grads中貌似维数不一样就不能参加运算

hgsndg

ldwloveyyz 发表于 2019-4-19 12:15
第一个，第二个也差不多
ga-> q ctlinfo
dset E:\aa%uwnd.sig995.1987.nc

我发现了，请问怎么修改nc文件呢

hgsndg

ldwloveyyz 发表于 2019-4-19 12:15
第一个，第二个也差不多
ga-> q ctlinfo
dset E:\aa%uwnd.sig995.1987.nc

好的，谢谢

yuhuiqing6

怎么处理的，我遇到了一样的问题