【求助】NCL功率谱检验-模式平均结果

higgs

科研小白一只，最近用观测数据和模式数据对某个模态的时间序列做了功率谱，在显著性检验上遇到了问题，虚心求教各位大佬！模式共10个member，每个member的区别是大气随机强迫的不同，模式平均结果可以近似理解为去除了大气随机强迫。模式对应的10个时间序列是这样获得的：把10个member的异常值在时间维度上叠加后做EOF，得到的时间序列分成了10份，这样可以保证其对应的模态是一样的。如图，红色是观测结果，黑色虚线是模式的10个member，黑实线是10条黑虚线的平均。


对观测结果很好做显著性检验，我直接照抄了官网的几行代码：

; detrending opt: 0=>remove mean 1=>remove mean and detrend
  d = 0
; smoothing periodogram: (0 <= sm <= ??.) should be at least 3 and odd
  sm = 5
; percent tapered: (0.0 <= pct <= 1.0) 0.10 common.
  pct = 0.10
;obs
  spec = specx_anal(ts_obs,d,sm,pct)
  splt = specx_ci (spec, 0.05, 0.95)


但是对模式平均的这条黑线该咋做显著性检验？每个member对应的95％置信区间线都不一样，直接平均一下是否妥当？
或者说应该对黑实线做红噪声检验？这样的话specx_ci函数的第一个参数spec如何给定呢？（每个member对应的spec的元数据有所差异）
我试了试，先对模式的10个时间序列平均再做功率谱貌似也不行，这样得到的结果比观测结果的方差要低很多，估计是因为平均后去掉了大气随机强迫部分。
有无大佬可以指点一下我，万分感谢！！！

JJY0319

问题解决了嘛?可以分享一下这个程序吗