grads画图提示Data Request Error: Invalid grid coordinates

樱花阿喵

文件pr_Amon_BCC-CSM2-MR_historical_1x1regrid_1961-2014.DJF.mean.nc和pr_monthly_CN05.1_1x1regrid_1961-2016.DJF.mean.nc分别是青藏高原冬季降水的模式模拟值和观测值，我想计算1995-2014年的降水距平百分率。但出现如下错误：Data Request Error: Invalid grid coordinates
  World coordinates convert to non-integer  grid coordinates
    Variable = pr.2  Dimension = 3
Operation Error:  Error from ave function
  Error ocurred at column 1
DEFINE error:  Invalid expression.

请各位帮帮忙，已经困在这好几天了。

第一个nc文件需要用ctl打开，它的ctl为：
DSET  E:\prr\pr\historical\pr_Amon_BCC-CSM2-MR_historical_1x1regrid_1961-2014.DJF.mean.nc
dtype netcdf
TITLE MEAN PRECIPITATION
undef 1e+20
OPTIONS 365_day_calendar
xdef 51 linear 60 1
ydef 31 linear 20 1
zdef 1 linear 0 1
tdef 54 linear 12Z16JUL1961 12mo
vars 1
pr=>pr  0  t,y,x  Precipitation
endvars

计算距平百分率的gs为：
'reinit'
'open E:\prr\pr\historical\pr_Amon_BCC-CSM2-MR_historical_1x1regrid_1961-2014.DJF.mean.ctl'
'set lon 70 110'
'set lat 25 40'
'define aa=ave(pr,t=34,t=53)'
'sdfopen E:\prr\pr\observe\pr_monthly_CN05.1_1x1regrid_1961-2016.DJF.mean.nc'
'set lon 70 110'
'set lat 25 40'
'define bb=ave(pr.2,t=34,t=53)'
'define c=(aa-bb)/bb'
'set gxout shaded'
'd c'
;

苏拉苏拉

pr.2有问题，可以看一下第二个文件的描述信息，是不是三维的数据？是不是与第一个文件维度不匹配？之类的

樱花阿喵

苏拉苏拉 发表于 2021-3-8 19:18
pr.2有问题，可以看一下第二个文件的描述信息，是不是三维的数据？是不是与第一个文件维度不匹配？之类的

这是第二个文件的文件信息，也是三维的数据，就是时间维多了一年，维度都是一样的，不知道哪里出了问题？

苏拉苏拉

樱花阿喵 发表于 2021-3-9 15:23
这是第二个文件的文件信息，也是三维的数据，就是时间维多了一年，维度都是一样的，不知道哪里出了问题？

第一个文件起始时间1961.12.16，第二个文件1962.01.01，有可能是这里没对上。你计算bb之前，'set dfile 2'试试。

樱花阿喵

苏拉苏拉 发表于 2021-3-10 15:12
第一个文件起始时间1961.12.16，第二个文件1962.01.01，有可能是这里没对上。你计算bb之前，'set dfile 2 ...

set dfile 2之后又出现了新错误，还是提示维度不对：gancgrid error: nc_get_vara_double failed; NetCDF: Index exceeds dimension bound

                               
登录/注册后可看大图

如果我单独打开第二个文件，define bb的话就没有错误；
但先打开第二个文件，再打开第一个的话还是会有Data Requst Error的错误提示，这两个文件的数据维度是一样的呀，真的不知道错误出在哪儿了

苏拉苏拉

樱花阿喵 发表于 2021-3-10 16:27
set dfile 2之后又出现了新错误，还是提示维度不对：gancgrid error: nc_get_vara_double failed; NetCDF ...

分别单独打开可以，调换位置还是一样的错误提示，说明就是两个文件的维度信息不匹配的问题。经纬度格点一个个都是对应的，时间网格起始点不对应。把第一个文件的ctl里的12Z16JUL1961改成和第二个文件的00Z01JAN1962一样即可。

瑾珂

你好，你的问题解决了吗？？怎么解决的