关于grads中scorr函数的使用

苜蓿

目的：求两次降水过程某区域的空间相关系数

问题：lz实在不会使用scorr函数，因为是新手手册和上网都找不到解决方法。。。于是只能向论坛里的大牛求助。。。希望不要被喷

菜鸟自己编写的gs：
'reinit'
'open d:\1\ncep\dat\3.ctl'
'open d:\1\ncep\dat\7.ctl'
'set x 1 144 '
'set y 1 73'
'set t 1'
'set lon 10 140'
'set lat 10 30'
'define hgt1=hgt.1'
'define hgt2=hgt.2'
'define a=scorr(hgt1,hgt2,lon=10,lon=140,lat=10,lat=30)'
'd a'
'reinit'

最后，先感谢一下大家~

随缘

This function calculates the spatial correlation between two variables over an X-Y domain. It returns a single number. The syntax is:

scorr(expr1, expr2, xdim1, xdim2, ydim1, ydim2)
where:


expr1   - any valid GrADS expression
expr2   - any valid GrADS expression
xdim1   - starting X dimension expression
xdim2   - ending X dimension expression
ydim1   - starting Y dimension expression
ydim2   - ending Y dimension expression

For global averaging, a shorthand may be used:

scorr(expr1, expr2, global) or
scorr(expr1, expr2, g)
is the same as
scorr(expr1, expr2, lon=0, lon=360, lat=-90, lat=90)

Usage Note

scorr may be used in conjunction with tloop or define to create time series or time/height plots.

scorr assumes that the world coordinates are longitude in the X dimension and latitude in the Y dimension, and does weighting in the latitude dimension by the delta of the sin of the latitudes. Weighting is also performed appropriately for unequally spaced grids.

Examples
This example calculates the correlation between the surface temperature and the latent heat flux over the tropical Pacific:


set lat -10 10
set lon 120 280
d scorr(tsfc, lhtfl, lon=120, lon=280, lat=-10, lat=10)

这个是grads自带的文本说明了

苜蓿

随缘 发表于 2012-5-10 22:46

                               
登录/注册后可看大图

This function calculates the spatial correlation between two variables over an X-Y domain. It return ...

谢谢~这个我也看到了~就是不知道怎么用才是正确的。。。或者说我的有什么问题。。。

苜蓿

继续求解。。。

随缘

苜蓿 发表于 2012-5-11 09:46

                               
登录/注册后可看大图

继续求解。。。

自己好好读读，里面写的很详细了，老是这么问是学不到东西的

苜蓿

随缘 发表于 2012-5-11 12:11

                               
登录/注册后可看大图

自己好好读读，里面写的很详细了，老是这么问是学不到东西的

我已经仔细读过了，也编了好多次。。。

最新的：
'reinit'
'open d:\1\ncep\dat\3.ctl'
'open d:\1\ncep\dat\7.ctl'

'set lon 10 140'
'set lat 10 30'
'set t 1'

'set lev 500'

'define a1=hgt.1'
'define a2=hgt.2'

'd scorr(a1,a2,lon=10,lon=140,lat=10,lat=30)'

'reinit'

但一直显示request beyond file limits，确实是不知道怎么办才来提问的

苜蓿

随缘 发表于 2012-5-11 12:11

                               
登录/注册后可看大图

自己好好读读，里面写的很详细了，老是这么问是学不到东西的

再次谢谢你~我再仔细看了又发现了问题，再重新编过。。。

'reinit'
'open d:\1\ncep\dat\3.ctl'
'open d:\1\ncep\dat\7.ctl'
'set x 1 144'
'set y 1 73'
'set t 1'

'set lev 500'
'define a1=aave(hgt.1,lon=10,lon=140,lat=10,lat=30)'
'define a2=aave(hgt.2,lon=10,lon=140,lat=10,lat=30)'


'set lon 10 140'
'set lat 10 30'
'd scorr(a1,a2,lon=10,lon=140,lat=10,lat=30)'

'reinit'


还是这个错误：request beyond file limits
不管加权不加权都说scorr这个有问题。。。不知道怎么办了

苜蓿

我发现问题了。。。我打开的第二个ctl文件总是无法读取，不论顺序如何，只要是放在第二个的ctl就是读不出来。。。为什么呢~有什么解决办法么~

苜蓿

把时间描述改成一样就没问题了。。。唉。。。绕了好多弯路

天边

恩 谢谢分享啊