grads 如何显示函数计算结果

yingkr

利用grads计算全球平均温度的时间序列，然后想计算两个时间序列的相关系数

程序如下：

'reinit'
'open G:/T639GSI2GRIB2_ORIG_20140701-20141020_00/2014.ctl'
'open G:/ECMWF/pl-20140701-20141020-temp.ctl'

'set z 1'
'set vpage 0 11 0.7 7.7'
'set grads off'
'set grid off'
'set t 1 112'
'set x 1'
'set y 1'
'set ccolor 1'
'set cmark 0'
'define aa1=tloop(aave(TMPprs,lon=0,lon=360,lat=-90,lat=90))'
'd aa1'
'set ccolor 4'
'set cmark 0'
'define aa2=tloop(aave(Tprs.2,lon=0,lon=360,lat=-90,lat=90))'
'd aa2'
'c'
'set t 1'
'set lon -180 180'
'set lat -90 90'
'd tcorr(aa1,aa2,t=1,t=112)'
;

程序没有问题，但不知道最后的结果'd tcorr(aa1,aa2,t=1,t=112)'怎么显示到屏幕上？

如果不编写gs文件的话，单独敲入d  tcorr(aa1,aa2,t=1,t=112)可以出现
Result Value=0.996

那gs文件要怎么写，能在频幕上直接出现Result Value呢？

lhaikun@163.com

帮顶 同问

Figo

两个序列算相关系数直接用excel就可以了啊，何必这么麻烦呢

yingkr

excel读不了grib数据吧  grads可以直接读进去 然后我就想直接用grads算一下相关 因为不知道怎么输出至屏幕 我只能在刚刚的程序后面加了一句'grads2ascii corr G:/test.txt' 把结果写出来看。。

wan_wan

在脚本里面加上‘draw string  x y  r=0.996'不就可以了？ x，y是要写上去的位置。

wan_wan

不过楼主，我想问一下，你看的是哪两段的时间序列？为什么研究相关系数？不太懂哎

yingkr

不是在图里显示 是在屏幕上显示相关系数的计算数值 因为我都不知道相关等于0.996 怎么写在图里呢？

yingkr

因为我花了两条时间序列 后面想算一下时间序列的相关

sunflower

其实很多办法的，要么w2xy，要么干脆在文档里加个文本框