爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 10694|回复: 17

[分享资料] 如何用grads计算季节变量的标准偏差

[复制链接]

新浪微博达人勋

发表于 2011-11-10 15:48:13 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
在网上看了利用Grads计算月时间序列的变量标准偏差,
其中多年月平均值(Mave)好计算,个别月的值为M. 在根号下面(Mave-M)^2表示.
但遇到算季节的该如何办呢?
多年的季节平均值(Save)也好算,但是每个季节(S)的时间序列如何确定啊?
谢谢!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

0
早起挑战累计收入
发表于 2011-11-10 18:18:06 | 显示全部楼层
所有的标准偏差方法不都是一样的么?既然月的能算出来,为何季节的算不出来?不是很明白楼主的问题了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2011-11-11 00:34:09 | 显示全部楼层
LZ的意思是季节的时间不好定?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2011-11-11 09:34:50 | 显示全部楼层

是啊, 可能我表述的不太清楚, 我季节(冬季, djf)的平均值是这样算的:
'define hgtD=ave(hgt,t=24,t=624,12)'
'define hgtJ=ave(hgt,t=25,t=625,12)'
'define hgtF=ave(hgt,t=26,t=626,12)'
'define hgtDJF=(hgtD+hgtJ+hgtF)/3',
但不知如何做时间序列(1949-2000)中每个季节相对于平均值hgtDJF的差值啊,

对于某个月我觉得好弄些,如一月份的可以计算为:
'd sqrt(ave(pow(hgt,2),t=25,t=625,12)-pow(ave(hgt,t=25,t=625,12),2))'
即个别月的循环好做,但季节的就不会了.

请教各位老师该如何做呢?

提前谢谢大家了!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2011-11-11 09:39:00 | 显示全部楼层
请各位知道的老师帮帮我啊.
谢谢!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2011-11-11 09:45:52 | 显示全部楼层
I attached a map of hgt at 700 hPa, in which I plotted the 2925 m isohypses for 76/77 winter (green one) and mean winter (1949-2000) .
Can we figure out how many standard devation green line offset from the black line? how can we figure out the region with significant (at the 0.05 level) difference between the two isohypses? All of these can be done using Grads?
Thanks so much for your kind help.

我原来设想如我计算出黒线的标准偏差, 大致可以根据两条线间的距离来判断绿线相对于黒线移动的标准偏差.

这样可以吗? 请各位老师指教啊
hgt.jpg
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2011-11-11 12:20:49 | 显示全部楼层
额~~~这个我也不太清楚了~~帮顶一下
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2011-11-11 12:46:51 | 显示全部楼层
善人/jw 发表于 2011-11-11 12:20
额~~~这个我也不太清楚了~~帮顶一下

那第一个问题呢?就是如何求季节的标准偏差? 季节的时间如何确定?

谢谢啊!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2011-11-11 17:18:13 | 显示全部楼层
额~~12,1,2是冬季,3,4,5是春,6,7,8是夏,9,10,11是秋
不知道是不是这个问题
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2011-11-11 17:48:38 | 显示全部楼层
善人/jw 发表于 2011-11-11 17:18
额~~12,1,2是冬季,3,4,5是春,6,7,8是夏,9,10,11是秋
不知道是不是这个问题

呵呵,这个俺清楚啊,关键是冬季的时间序列如何定啊?
对于某个月,如一月份的可以做下面的循环.
'd sqrt(ave(pow(hgt,2),t=24,t=626,12)-pow(ave(hgt,t=25,t=625,12),2))'
而对于季节该如何做呢?
如第一个冬季DJF(1949/50)可定义为w1= ave (hgt, t=24,t=26), 后面的如何做循环呢? 一直到1999/2000年的冬季w51= ave(hgt,t=624,t=626), 如何从时间上来定义让它每隔12个月循环一次呢?
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表