爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3736|回复: 0

[脚本编辑] 关于计算月距平

[复制链接]

新浪微博达人勋

发表于 2016-6-2 15:59:40 | 显示全部楼层 |阅读模式

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

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

x
我下载了ncep的月平均资料,时间长度是1979年到2015年,想要计算每个月的距平值,比如1998年2月份距平,应该是1998年2月份的数据减去1979到2015每年的2月份数据的平均值,在论坛上也参考了别人的帖子,对于modify的用法还是不明白,以下是官网的说明:
modify varname type
This command defines a climatological variable, which is year-independent. varname is a defined grid. There are two options for type:

seasonal    - For creating monthly or multi-monthly climatologies
diurnal     - For creating climatologies over a time period less than a day
Usage Notes

Example

Say you have a 50-year timeseries of monthly mean sea surface temperatures (a variable named sst with 600 time steps) and you want to create a climatology and then look at the monthly anomalies. First, set the time range for 1 to 12, to span a complete year. Second, define the variable "sstclim" which will contain the January mean in the first time step, the February mean in the second time set, etc. Then use 'modify' to turn 'sstclim' into a climatological variable. This means that the calendar year associated with 'sstclim' (the first year in the original sst data set) becomes a wild card. Then you can define the anomaly by subtracting the climatology from the original time series. The commands are as follows:

'set t 1 12'
'define sstclim = ave(sst, t+0, t=600, 12)'
'modify sstclim seasonal'
'set t 1 last'
'define anom = sst - sstclim'

不明白的地方是
'set t 1 12'
'define sstclim = ave(sst, t+0, t=600, 12)'
这里sstclim是怎么样的值,既然 'set t 1 12',为啥ave(sst, t+0, t=600, 12)里面有时间步长12,t+0又是什么鬼
求大神们指点迷津
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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