tcorr用法

叩容

个人觉得其实这函数实际意义不大，因为同期相关毕竟。。。不是大多数。但一些时候大概估计下还是可以用到的。

虽然LASG手册已经很详细，但是一开始照着例子用tcorr仍然常常会出错，调试很久。下面记录我自己用此函数时的要点，希望给各位一些帮助。

1.
X 、Y、Z层必须统一。即 linear 对应 linear  ，leves 对应 leves 。
X 、Y、Z层初始位置不一致倒是没关系。如 xdef  360 linear 0.500000 1.000000 对应 XDEF 1 linear 1 1 是可以通过的。

2.
简单计算两个时间序列的相关系数用这个函数是很简单的。
'set x 1'
'set y 1'
'set z 1'
'set t 1'
'd tcorr(a,b,t=1,t=365)'
就可以得到结果。

3.
如果需要计算点对场相关，就必须严格按照例子。特别注意，
'define a=时间序列'
这步不可以省略。

4.
'run define_colors.gs'
'set clevs   -0.5 -0.4 -0.3 -0.2 -0.1 -0.05 0.05 0.1 0.2 0.3 0.4 0.5'
'set ccols 48 47 46 44 43 41 0 21 23 24 26 27 28'
'd tcorr(a,b,t=1,t=73)'
直接出图也很漂啦！

mofangbao

是不是从你空间拷过来的啊，要重新上传一下图片的
ps：好久没看到你了哈

叩容

mofangbao 发表于 2013-4-10 10:25

                               
登录/注册后可看大图

是不是从你空间拷过来的啊，要重新上传一下图片的
ps：好久没看到你了哈

收到，已经重新上传。。。
俺不热爱学术，整天都在忙别的事，消失很正常。。。。

river

谢谢分享经验

蔚蓝的天

表示没看明白啊

开始起飞1

'reinit'
*'set gxout fwrite'
*'set fwrite c:\GrADS19\data\bbbb.grd'
'open c:\GrADS19\data\2000data.ctl'
'open c:\GrADS19\data\1988data.ctl'
*'set lon 73 135.5'
*'set lat 17 55'
'set x 1'
'set y 1'
'set z 1'
'set t 1'
'define res=tcorr(a,a.2,t=1,t=100)'
'd res'
'enable print c:\GrADS19\data\bbb.gmf'
'disable print'
;
看来我还是没弄明白啊，求指点啊

叩容

开始起飞1 发表于 2013-4-30 10:11

                               
登录/注册后可看大图

'reinit'
*'set gxout fwrite'
*'set fwrite c:\GrADS19\data\bbbb.grd'

你出现的问题是啥呢？

hikingbaby

一直没弄明白tcorr中时间系列是一个格点的系列还是这个空间范围内的系列，空间场是一个时间的场还是多个时间的场。

chuanhao

如果想把求出来的这个相关系数写入grd，怎么办呢，我写进去老出错

lucy_9208

hikingbaby 发表于 2014-4-28 16:27
一直没弄明白tcorr中时间系列是一个格点的系列还是这个空间范围内的系列，空间场是一个时间的场还是多个时 ...

多个时间的空间场