求大神帮忙啊！grads 的 tcorr 函数出错！

aiai1992

我把文件都上传了，希望各位大神帮帮忙，小弟在此谢过！压缩包里面有CTL GS,DAT
'open d:\vf\sst\'tt'.ctl'
'open d:\vf\grads\'yy'.ctl'
'set dfile 2'
'set x 1 1'
'set y 1 1'
'set z 1 1'
'set t 2 66'
'define corr=n'
'set dfile 1'
'set x 1 360'
'set y 1 180'
'set z 1'
'set t 1'
'd tcorr(corr,sst,t=1,t=65)'
这是我的程序，我之前用过，没有问题，可是最近两天运行时总是出错，自己实在找不到原因了，求各位大侠帮着看一下，tcorr这个函数有什么特殊的要求吗？还是怎么，为什么总是提示我corr变量是三维的啊？

柠檬夏

两个ctl文件的tdef起始时间不一致，也许得修改成一样的吧？

aiai1992

柠檬夏 发表于 2014-4-27 22:05
两个ctl文件的tdef起始时间不一致，也许得修改成一样的吧？

还是不可以的，我觉得和那个时间没什么关系，我要做超前和滞后相关。你运行的时候提示什么错误啊？

柠檬夏

aiai1992 发表于 2014-4-28 10:25
还是不可以的，我觉得和那个时间没什么关系，我要做超前和滞后相关。你运行的时候提示什么错误啊？

我把d:\sst\1.ctl的tdef项改成“tdef 66 linear 00Z16Jan1948 1yr”就能运行了,但是这违背你做超前和滞后相关的要求。我觉得你还先得是用Fortran处理grd，截取、输出你要研究的数据段，再编写新的ctl文件来运行吧。
grads中一些函数对不同ctl变量中的tdef要求起始时间、时间步长一致，但tcorr函数我就不清楚了。

lqouc

这样不行的，在define前定义了时间范围，那么你定义的变量的时间范围就是那个，后面的相关计算时自然会出错。先要算超前之后相关，你可以用fortran先把前后的部分截断，使得正好是这样，或者就修改两个ctl文件好了让他们错开一个时间点，这样就能计算了，反正物理意义你明白就好了。

aiai1992

柠檬夏 发表于 2014-4-28 12:39
我把d:\sst\1.ctl的tdef项改成“tdef 66 linear 00Z16Jan1948 1yr”就能运行了,但是这违背你做超前和滞后 ...

好的，谢谢了！

aiai1992

lqouc 发表于 2014-4-28 12:46
这样不行的，在define前定义了时间范围，那么你定义的变量的时间范围就是那个，后面的相关计算时自然会出错 ...

哦哦，原来是这样，这样就增加了很多工作量啊！也许这就是我程序做出的图和Matlab不一样的原因了！多谢！

lqouc

aiai1992 发表于 2014-4-28 21:14
哦哦，原来是这样，这样就增加了很多工作量啊！也许这就是我程序做出的图和Matlab不一样的原因了！多谢！

其实没什么工作量，直接把ctl中的一个文件时间错后一个时刻，然后计算时用time设定纬度就行了。
可以自己揣摩下，grads其实是很灵活的。

aiai1992

lqouc 发表于 2014-4-28 21:58
其实没什么工作量，直接把ctl中的一个文件时间错后一个时刻，然后计算时用time设定纬度就行了。
可以自 ...

好的，谢谢学长！

aiai1992

lqouc 发表于 2014-4-28 21:58
其实没什么工作量，直接把ctl中的一个文件时间错后一个时刻，然后计算时用time设定纬度就行了。
可以自 ...

弄好了，十分感谢！