山姆文档
山姆计算南环模式指数的版本，也称为南极振荡。
内容
句法
IDX = SAM（SLP40，SLP65，T）
描述
IDX = SAM（SLP40，SLP65，T）计算两个纬度（40s和65s）的海平面压力的两个时间序列的南环模式索引及其相应的时间t。
例子
在这里，我们将从马歇尔（Marshall）的2003年经典论文中重新创建图7，即观察和reanalyses的南环模式下的趋势（马歇尔，，2003），描绘每月SAM索引。
加载数据
我们可以从1957年1月至2018年12月的观测值中加载包含40s和65s的Zonal均值SLP的数据。
加载sam_slp_data.mat
绘图压力数据
SAM_SLP_DATA.MAT文件包含1957年至2018年的每月层次平均压力。计算SAM：
图图（T，SLP40）;抓住上图（T，SLP65）轴紧的传奇（“平均SLP在40 \ Circs”，，，，“平均SLP在65 \ circs”）
计算SAM索引
这山姆功能均与1971-2000基线相对于每个时间序列归一化，并且两个压力异常时间序列的差异以产生SAM指数。
％计算SAM指数：SAM_IDX = SAM（SLP40，SLP65，T）;图图（t，sam_idx，'：'，，，，'颜色'，RGB（'灰色的'）;抓住上情节（T，Movmean（SAM_IDX，12），'K'，，，，'行宽'，1）Ylabel（“山姆指数”）xlabel（'年'）XLIM（[DateTime（1955,1,1）DateTime（2005,1,1）]）％设置日期限制Ylim（[-8 8]）％设置垂直限制hline（0，'K--'）％绘制水平线
上面，内置MATLAB功能movemean被用来计算12个月的移动平均值，并绘制了一条水平虚线Hline匹配马歇尔的图7。
获得区域卑鄙的压力
上面的示例首先加载40和65s的Zonal平均压力。如果您使用网格数据，则需要计算40和65s的区域均值。这样做。首先加载一些示例数据：（使用ncdateread阅读日期时间）。
lat = double（ncread（ncread）（'ERA_INTERIM_2017.NC'，，，，'纬度'）;lon = double（ncread（ncread）（'ERA_INTERIM_2017.NC'，，，，'经度'）;t = ncdateread（'ERA_INTERIM_2017.NC'，，，，'时间'）；sp = ncread（'ERA_INTERIM_2017.NC'，，，，'sp'）；
在40秒和65s中获得区域平均值的一种蛮力方法是在网格中每个独特的经度下插入40和65s，然后为每个时间步中平均结果。首先介绍输出，然后循环遍历每个时间步长，插值为40和65s
％prealloate Zonal平均压力时间序列：slp40 = nan（size（t））;slp65 = nan（size（t））;％循环通过每个时间步骤：为了k = 1：长度（t）％以所有纵钉的平均值为单位，插值为40：slp40（k）=平均值（interp2（lat，lon，sp（：，：，k）， - 40*一个（size（lon），lon），lon））;％以所有纵钉的平均值为单位，插值为65s：slp65（k）=平均值（interp2（lat，lon，sp（：，：：，k），-65*一个（size（lon），lon），lon））;结尾图图（T，SLP40）保持上情节（T，SLP65）传奇（“平均SLP在40 \ Circs”，，，，“平均SLP在65 \ circs”）
参考
Marshall，G。J.，2003年：观测和重新分析的南环模式的趋势。J. Clim。，16，4134-4143。https://journals.ametsoc.org/doi/pdf/10.1175/1520-0442%282003%29016%3C4134%3ATITSAM%3E2.0.CO%3B2
作者信息
这山姆功能和支持文档由金宝appKaustubh Thirumalai对于MATLAB的气候数据工具箱，2019年。