grads用tcorr算出来的相关系数数量级不对，求大神指导

笨北极熊笨

第一个数据是ccismmon45.nc，是土壤湿度数据，d sm 显示数据是正常的。

第二个数据是一个36年的时间序列
数据也是正常的，我就写了一个gs算相关：
'reinit'
'set gxout fwrite'
'set fwrite E:\zj\201501\GPCP\ccism\smcor.dat'
'sdfopen E:\zj\201501\GPCP\ccism\spmizhou.nc'
'sdfopen E:\zj\201501\GPCP\ccism\ccismmon45.nc'
'set x 1'
'set y 1'
'set z 1'
'set t 1 36'
'define pmi=pmi.1'
'set dfile 2'
'set lon 0 360'
'set lat -89.875 89.875'
'set z 1'
'set t 1'
'd tcorr(pmi,0.0001*sm,t=1,t=36)'
'disable fwrite'
成功生成一个dat数据，根据数据写了一个ctl文件
DSET  E:\zj\201501\GPCP\ccism\smcor.dat
*options template big_endian 365_day_calendar
UNDEF  -9999
TITLE  k         
XDEF 1440 LINEAR -179.875 0.25
YDEF 720 LINEAR -89.875 0.25
ZDEF 1 LEVELS 0 1
TDEF 1 LINEAR 00z01MAY1979 1mo   
VARS 1         
r  0 -999 plateau monsoon index   
ENDVARS
可以打开文件，但生成的数据数量级不对，相关系数应该在正负1之间

不知道哪儿出错了。。。。。。。。。。。。。。。
PS： 用NCL 的esccr函数也算了，但是相关系数都到正负2.4了
求大神指导我到底是哪儿有问题啊



来自群组: 国科大

笨北极熊笨

'd tcorr(pmi,0.0001*sm,t=1,t=36)'
这个乘以或者不乘以0.0001结果都一样，数量级不对，0.0001是scale factor，add_offset 为0