24小时热门版块排行榜    

查看: 264  |  回复: 2
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 yqingl 的 10 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

yqingl

金虫 (正式写手)

[交流] 【求助】请教关于Grads

我想在Grads利用uwnd在500和600hPa进行相关分析,编写的gs如下,但总是出错,版上有懂的朋友么?能否帮我看下出错在哪个地方?多谢了。

'reinit'
'sdfopen d:\1\uwnd.mon.mean.nc'
'enable print d:\1\you.gmf'
'set lon 70 110'
'set xlint 5'
'set lat 10 50'
'set ylint 5'

'set t 1'
'set mpdset cnworld'
'set gxout contour'
'define r1=uwnd(lev=600)'
'define r2=uwnd(lev=500)'
'define r=tregr(r1,r2,t=1,t=100)'
'd r'

'print'
'disable print'
'clear'

[ Last edited by ewigkeit on 2008-12-7 at 18:14 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yanlibin

禁虫 (著名写手)

tcorr()
tcorr (expr1, expr2, tdim1, tdim2)

This function produces a spatial map of temporal correlation coefficients.


expr1   - a valid GrADS expression that varies in time
expr2   - a valid GrADS expression that varies in time and may also vary in X and Y
tdim1   - starting time dimension expression
tdim2   - ending time dimension expression

The expr1 time series is correlated to the time series at each grid point in expr2. The result is a grid of correlation coefficients that matches the X and Y dimensions of expr2.
3楼2008-11-01 16:01:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

yanlibin

禁虫 (著名写手)

'define r=tregr(r1,r2,t=1,t=100)'

r1应该是一个单点的值,相当于y
r2应该是一个场的值,相当于很多x
y=a*x+b
求出来的是a,最后的结果是r1与每个点上的r2多年的回归系数

相关分析最好用fwrite写出来,用fortran算吧
2楼2008-11-01 15:59:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见