NCL做一个简单的EOF分析例子

被要求用NCL进行制图，作为小白，正在慢慢摸索中，如有建议，请指正！

1.EOF分析方法

EOF（经验正交函数），分析得出正交的模式和时间序列。派生的模式是域和所使用的时间段的函数。在可以通过选择与协方差矩阵的最大特征值关联的特征向量的顺序来考虑最大方差的意义上，EOF表示是最佳的（Kutzback，1967）。但是，EOF过程严格来说是数学过程，并不基于物理学。结果可能会产生类似于系统内物理模式的模式。

需要明确的是，EOF分析不是统计过程。但是，由于结果是正交的（即独立的），因此在某些统计应用中将其用作预测变量。

2.EOF分析例子

对海平面压力（SLP）作eof分析

NCL官网关于EOF的例子
参考官网例子
北半球过去2000年SLP的EOF分析图，如下所示。


用到的主要函数有 eofunc_ts_Wrap(data, evec, optETS)、eofunc_north(eval, N, prinfo)

主要过程：

1.导入数据
2.计算EOF

w      = sqrt(cos(0.01745329*p&latitude))                         
wp     = p*conform(p, w, 1)                        
copy_VarCoords(p, wp) 

neof   = 4                                          
x      = wp(latitude|:,longitude|:,time|:)         
eof    = eofunc_Wrap(x, neof, False)
eof_ts = eofunc_ts_Wrap (x, eof, False)

3.显著性检验

dimp   = dimsizes(p)
ntim   = dimp(0)                                         
  
prinfo = True
sig    = eofunc_north(eof@pcvar, ntim, prinfo)            

4.建立工作站，画图显示
用gsn_csm_contour_map_polar()画空间场图，用gsn_csm_xy()画折线图