点对场相关

mozhiting

用函数tcorr做点对场的相关，出现了问题，现已解决，特总结如下。点的时间序列存放在一个数据文件中，空间场的时间序列存放在另一个数据文件中，描述文件分别如下：1、
DSET  E:\GrADS\homework\17\wy.dat
TITLE monthly  ERA-40 air temperature, from 7/1958 to 7/1997
UNDEF  9.999E+20
XDEF   1   LINEAR  50   5
YDEF   1   LINEAR  -20  2
ZDEF   1  levels 1000
TDEF   40  linear jun1958 1yr
VARS  1
wy      1  99 ** W-Y index
ENDVARS
2、
dset E:\GrADS\homework\17\sst.dat
title sub-surface sea water temperature anomaly data 1958-1997
undef -999
xdef 41 linear 50 5
ydef 34 linear -20 2
zdef 1 levels 1000
TDEF   40  linear jun1958 1yr
vars 1
sst 1 99 sst
endvars

分别打开均无错误。

脚本文件如下：
'reinit'
'open e:\GrADS\homework\17\sst.ctl'
'open e:\GrADS\homework\17\wy.ctl'
'set t 1'
'set x 1 41'
'set y 1 34'
'set z 1'
'set gxout contour'
'define ts=tcorr(wy.2,sst,t=1,t=40)'
'd ts'
say 'end'


错误提示如下：

后进行修改，内容如下：

'reinit'
'open E:\GrADS\homework\17\wy.ctl'
'open E:\GrADS\homework\17\sst.ctl'
'set grid off'
'set grads off'

'set t 1 40'
'define a=wy'
'set dfile 2'
'set lon 50 250'
'set lat -20 45'
'set t 1'
'define b=tcorr(a,sst,t=1,t=40)'
*******************************************draw it
'set gxout shaded'

'set clevs -0.312 0.312'
'set ccols 4 0 2'
'd b'
'cbarn'
'set gxout contour'
'set cstyle 2'
'set clevs -0.4 -0.3 0  0.3 0.4'
'd b'

'set gxout contour'
'set clevs 0 0.3 0.4'
'set ccolor 2'
'd b'
'draw title W-Y'

画出图如下：

river

看看官网的例子是怎么设置的吧
This example calculates the temporal correlation between sea level pressure and the defined variable elnino, an areal average of surface temperature in the equatorial Pacific.
set x 1
set y 1
set z 1
set t 1 100
define elnino = aave(ts,lon=-160,lon=-80,lat=-10,lat=10)
set lon -180 180
set lat -90 90
set z 1
set t 1
d tcorr(elnino, slp, t=1, t=100)

mozhiting

river 发表于 2017-6-17 22:47
看看官网的例子是怎么设置的吧
This example calculates the temporal correlation between sea level pre ...

我的第一个变量已经是区域平均了，所以前半部分不用写了，而后半部分是一样的啊！

river

mozhiting 发表于 2017-6-18 09:13
我的第一个变量已经是区域平均了，所以前半部分不用写了，而后半部分是一样的啊！

谁告诉你不用写了，错误提示很明白的告诉你了，第一个变量必须是 0维的。

mozhiting

river 发表于 2017-6-18 21:29
谁告诉你不用写了，错误提示很明白的告诉你了，第一个变量必须是 0维的。

嗯嗯，已经搞明白了，谢谢~