ncl 求特定时段标准差

longlivehj

deeli 发表于 2014-3-4 15:21

                               
登录/注册后可看大图

前辈，我想用同样的思路提取上一年11月~下一年4月的数据，发现因为有2月在里面，数据没办法读到temp变量里 ...

如果一年一年做，不需要求多年统计值，我感觉只需要用reassignment就好了，temp :=

lovjiangyp

好吧，表示也不会

lovjiangyp

好吧，表示也不会

deeli

longlivehj 发表于 2014-3-4 15:48

                               
登录/注册后可看大图

这就是为什么我在最早给你的示例程序里面有下面红色两句的原因。
olr2d!0 = "time"
olr2d&time = (/197 ...

大概可以明白的，但是我在不管是5~10月还是11~4月，olr2d({},y-1974)这个大括号里面temp时间取余去做总会报错，因此，夏天的数据我是从1975年开始的，然后写成了olr2d（：，y-1974）冬天的数据想尝试着用您说的这个方法，总是报错呢，报：subscript out of range ~这可怎么办呢？

longlivehj

deeli 发表于 2014-3-5 16:27

                               
登录/注册后可看大图

大概可以明白的，但是我在不管是5~10月还是11~4月，olr2d({},y-1974)这个大括号里面temp时间取余去做总会 ...

你打印temp&time % 10000和olr2d&time出来看看呢？我估计里面有一个的时间坐标有问题。

deeli

longlivehj 发表于 2014-3-5 16:38

                               
登录/注册后可看大图

你打印temp&time % 10000和olr2d&time出来看看呢？我估计里面有一个的时间坐标有问题。

看起来没有问题啊~~~

deeli

deeli 发表于 2014-3-5 17:02

                               
登录/注册后可看大图

看起来没有问题啊~~~

都是直接截图的最后一年，2011年的

longlivehj

^_^，你只打印了一个么，有时候只看变量的属性是不行的。
总之，要确保temp和olr2d的时间坐标都是yyyymmdd形式。因为坐标变量必须是递增或者递减排列的，像处理5-10月数据一样time坐标只用月和日表示就不行了。
你能传我一份数据么？我写一段代码给你。否则，我一直再猜你的数据是……
邮箱我发消息给你。

deeli

longlivehj 发表于 2014-3-5 17:12

                               
登录/注册后可看大图

^_^，你只打印了一个么，有时候只看变量的属性是不行的。
总之，要确保temp和olr2d的时间坐标都是yyyymmdd ...

恩恩，好的！

changelyn

longlivehj 发表于 2014-2-27 20:28
do y = 1974, 2013
    temp =  olr1d({y * 10000 + 501 : y * 10000 + 1031})
    olr2d({temp&time % ...

如果求1979-2010每年冬季（NDJFM）包含二月份的值，应该怎么做呢