grads对两个序列求相关并画出图

hwt5813023

这个问题纠结了我一天了,问题如下
对温度和湿度的异常序列之间求相关系数
相关系数我求出来了 可是不知道如何用这些数据画出一维的图
'sdfopen e:\air.mon.mean.nc'
'sdfopen e:\rhum.mon.mean.nc'
'set lon 120 260'
'set lat 0 60'
'set lev 1000'
'set t 1 12'
'define airave=ave(air,t+0,t=396,12)'
'modify airave seasonal'
'define shuave=tvrh2q(ave(air,t+0,t=396,12),ave(rhum.2,t+0,t=396,12))'
'modify shuave seasonal'
'set t 1 396'
'd scorr(air-airave,tvrh2q(air,rhum.2)-shuave,lon=120,lon=260,lat=0,lat=60)'
我应该在哪加上写什么????????

♂雨已~○

在你画相关系数之前把经度、纬度、时间、高度都设定成定值~~~自然就是一维的图了

hwt5813023

♂雨已~○ 发表于 2012-5-20 10:09

                               
登录/注册后可看大图

在你画相关系数之前把经度、纬度、时间、高度都设定成定值~~~自然就是一维的图了

谢谢 我先试试

hwt5813023

♂雨已~○ 发表于 2012-5-20 10:09

                               
登录/注册后可看大图

在你画相关系数之前把经度、纬度、时间、高度都设定成定值~~~自然就是一维的图了

可是我用区域与区域求的相关系数啊..如果再把这些设置成定值该如何操作啊?

hwt5813023

♂雨已~○ 发表于 2012-5-20 10:09

                               
登录/注册后可看大图

在你画相关系数之前把经度、纬度、时间、高度都设定成定值~~~自然就是一维的图了

哈哈 谢谢 出图了 可以交差了

mofangbao

hwt5813023 发表于 2012-5-20 10:18

                               
登录/注册后可看大图

哈哈 谢谢 出图了 可以交差了

把正确的脚本发一下吧，谢谢啦

hwt5813023

'sdfopen e:\air.mon.mean.nc'
'sdfopen e:\rhum.mon.mean.nc'
'set grads off'
'set lon 120 260'
'set lat 0 60'
'set lev 300'
'set t 1 12'
'define airave=ave(air,t+0,t=396,12)'
'modify airave seasonal'
'define shuave=tvrh2q(ave(air,t+0,t=396,12),ave(rhum.2,t+0,t=396,12))'
'modify shuave seasonal'
'set lon 120'
'set lat 0'
'set lev 300'
'set t 1 396'
'define res=scorr(air-airave,tvrh2q(air,rhum.2)-shuave,lon=120,lon=260,lat=0,lat=60)'
'set x 1'
'set y 1'
'set lev 300'
'set t 1 396'
*'set xaxis  1979 2011 3'
*'set cmark 0'
*'set ccolor 2'
'd tloop(res)'
星号打头无特殊作用

老师居然还让我做两个区域点对点求相关并画出图....不过有了这些基础 那个还挺好求的...
再上张图吧 庆祝一下

善人/jw

最近好多朋友都在问相关系数，可以给大家做个例子了

易小凯

原来还有scorr这个函数，学习了

♂雨已~○

善人/jw 发表于 2012-5-20 13:32

                               
登录/注册后可看大图

最近好多朋友都在问相关系数，可以给大家做个例子了

师兄说得好啊~~~要不师兄带队做一个例子？