从tcorr函数的用法到grads对时间维的控制

小草

早上用grads画一个时间序列和另外一个三维变化的场的变量求相关时，查阅了以前发的帖子，很是受益。但是在查阅了相关讨论后依然还是会有问题出现，遂自己想办法予以解决，现和大家分享一下，共同学习。
gs如下：
'c'
'reinit'
say ''
'open F:\aa\plot\aa\bb.ctl'
'sdfopen E:\data\sst\111\sst.mon.anom.nc'
'set x 1'
'set y 1'
'set t 1 636'
'define funny=aa'
'set dfile 2'
'set lon 1 360'
'set lat -30 80'
'set t 1'
'set vpage 0 11 0 8.5'
'set grid off'
'set grads off'
'd tcorr(funny,sst,time=JAN1960,time=DEC2012)'   *重点在这一句
'draw title the relationship between aa and sst'
'enable print F:\aa\plot\aar\aa.gmf'
'print'
'disable print'
'reset'
;

    由于两个数据的时间起止范围是不一样的，所以当将'd tcorr(funny,sst,time=JAN1960,time=DEC2012)' 改成'd tcorr(funny,sst,t=1,t=636)' ,则将会出现全场缺测的现象，但是将时间控制写成绝对的年月时间时，则可以顺利出图。
这样就不用再把sst的相应时间段的资料再额外切出来，重新配ctl文件了。

善人/jw

多谢分享经验~~

river

好经验，谢谢分享

kongfeng0824

感谢分享！！！！！

斥鷃

小问题初学者很容易犯，谢谢楼主分享啦~

liutaoz21

感谢楼主分享经验！

麻花

谢谢楼主啦

叩容

楼主这个重点句……

“
'd tcorr(funny,sst,time=JAN1960,time=DEC2012)'
“

sst不写成sst.2 ，你的grads找得到吗？

小草

叩容 发表于 2013-12-18 15:01

                               
登录/注册后可看大图

楼主这个重点句……

“

'set dfile 2' 了，呵呵

chuanhao

小草 发表于 2013-12-25 20:56
'set dfile 2' 了，呵呵

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