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

小草 · 发表于 2013-6-29 11:14:13

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由小草于 2013-6-29 11:39 编辑

早上用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 · 发表于 2013-6-29 11:20:37

多谢分享经验~~

river · 发表于 2013-6-29 12:19:33

好经验，谢谢分享

kongfeng0824 · 发表于 2013-6-29 22:09:05

感谢分享！！！！！

斥鷃 · 发表于 2013-6-29 22:48:27

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

liutaoz21 · 发表于 2013-6-30 01:27:39

感谢楼主分享经验！

麻花 · 发表于 2013-9-2 19:57:32

谢谢楼主啦

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

楼主这个重点句……

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

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

小草 · 发表于 2013-12-25 20:56:50

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

登录/注册后可看大图

楼主这个重点句……

“

'set dfile 2' 了，呵呵

chuanhao · 发表于 2014-5-20 20:57:07

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

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

		自动登录	找回密码
密码			立即注册

[分享资料] 从tcorr函数的用法到grads对时间维的控制

登录后查看更多精彩内容~

评分

呵呵