气象统计方法&短期气候预测代码汇总

公布这些代码的目的之一肯定是解决大家实习课的困难，二是想让大家看到更简洁易懂的代码。

之前或多或少看过别人写的例如气象统计方法实习or短期气候预测实习的代码，那些代码大多冗余繁杂，不够简洁易懂，里面大都充斥了for循环语句，若是数据的分辨率较低倒还好，若是遇上高精度的数据还能用for循环吗？真的很浪费时间的！！！

这份代码是改良多次的结果。会尽可能避免for循环的使用，数据处理一般三下五除二，用不了几行就处理好了，但每一行代码都值得推敲，希望使用的同学明白为什么是这样处理。

另外，关于绘图，绘图是不得不写更多行去美化修饰，大家也可以自行美化，这里也只是简单的美化了一下。它不同于grads或NCL专门的气象绘图软件，python绘图要自己去调整。

一

资料说明

现有热带太平洋海温场资料，文件名 NCEP_TPSST_30y_Wt.dat

• 范围：120~300E，20S~20N
• 时段：冬季1978~2007年共30年
• 水平分辨率：不等距（具体参考NCEP_TPSST_30y_Wt.ctl文件）
  格点数：32 * 7

实习要求

1. 利用该资料求热带太平洋海温气候场、均方差场，并给出1982年和1998年冬季海温原始数据场和距平场，对计算结果进行初步分析。
2. 进一步计算Nino3.4区（5S-5N，170W-120W）海温指数，分别给出原始数据序列、距平序列及标准化时间序列。根据这些时间