grads用tcorr函数出错，不知道怎么回事，以前就好好的。

笨北极熊笨

DSET   h:\data2\newsummerpmi.dat
*options template big_endian 365_day_calendar
UNDEF  -9.96921e+36
TITLE  k         
XDEF 1 LINEAR 1 1
YDEF 1 LINEAR 1 1
ZDEF 1 LEVELS 0 0
TDEF 36 LINEAR 00z01jan1979 1yr   
VARS 1           
pmi  0 -999 plateau monsoon index     
ENDVARS

********求相关********
'reinit'
'set gxout fwrite'
'set fwrite h:\1\sm1\1\1\678corr.dat'
'open h:\data2\newsummerpmi.ctl'
'sdfopen h:\1\sm1\1\36mean.nc'
'set x 1'
'set y 1'
'set z 1'
'set t 1 36'
'define pmi=pmi.1'
'set dfile 2'
'set lon 60 120'
'set lat 0 60'
'set z 1'
'set t 1'
'd tcorr(pmi,swvl1,t=1,t=36)'
'disable fwrite'
提示的错误在图片中，其中36mean.nc的ctl也在图片中，找了好久的错没找出来。赶脚自己的好像似乎没错啊？
各位有没遇到这种情况，求大神啊。。。。

katrina

你看一下是不是你两个ctl的时间维不一样！

carpenter3

'set t 1'
'd tcorr(pmi,swvl1,t=1,t=36)'

这个写法貌似不对吧

笨北极熊笨

katrina 发表于 2014-11-14 11:14
你看一下是不是你两个ctl的时间维不一样！

一样啊，都是36个

笨北极熊笨

carpenter3 发表于 2014-11-14 11:31
'set t 1'
'd tcorr(pmi,swvl1,t=1,t=36)'

那应该怎么写呢？

katrina

你的ctl文件一个是36年，一个36个月

liuzf406

katrina 发表于 2014-11-14 20:38
你的ctl文件一个是36年，一个36个月

我 也 是 用 这 个 命 令 出 错 了 ,你 用 的 是 CMIP5数 据 吗 ?
CMIP5数 据 :
lon_bnds
lat_bnds  
这 咋 都 成 变 量 了 呢  

katrina

你把那份描述文件改一下，时间维一定要相同！

pqman3

36mean.nc是从1979年8月1日开始，而newsummerpmi.dat是从1979年1月1日开始，时间维不一致。

笨北极熊笨

pqman3 发表于 2014-11-14 23:45
36mean.nc是从1979年8月1日开始，而newsummerpmi.dat是从1979年1月1日开始，时间维不一致。

真的，改了就可以了，以前都不知道。万分感谢。。。。。。