气象统计方法实习一：根据定义判断Elnino和Lanina【grads,fortran,python】

小飞、

  突然就想说开个坑，记录一下自己写过的实习、绘图代码啥的。和大伙讨论分享一下。  气象统计方法实习一就是判断El Nino年和La Nina年，题目如下：

  我一开始想说直接用pygrads、xgrads啥的把.dat和.ctl直接读进python开整啥的。想了想，也不太好搞，而且之前学的grads和fortran也不能荒废，就重新温习了一下。
基本顺序就是：
1.先用fortran处理数据，计算距平，标准差，均值啥的，输出成二进制grd文件，同时处理Nino3,4区数据，输出为python绘制时间序列图用的csv文件(原始数据，距平，时间序列标准化数据)。
2.根据题目要求，编写ctl文件，再用grads绘制海温气候场、均方差场、1982年和 1998 年冬季海温原始数据场和距平场的图片并输出为gmf文件。(下图是原始数据场的图)
3.再用python绘制原始数据序列、距平序列及标准化时间序列图即可。效果如下图所示：

再根据数据就可以得出结果：
Nino3,4区的标准偏差为：  0.9702401（这里使用标准差的是每年数据标准差和的平均，0.97误差过大，我觉得直接写1无可厚非，毕竟年平均后的标准差是应重算的吧，重算后就是1吧）
El Nino年为： [1982, 1986, 1991, 1997, 2002]
La Nina年为：[1984, 1988, 1998, 1999, 2007]

网上查了一下，1982年冬季为强El Nino年，所以应该没错。没细检查，源码和其他图啥啥的都在附件里，要个1贡献无可厚非吧嘿嘿。
图有点丑。有能力的朋友可以美化一下。
欢迎大伙讨论，也欢迎大佬指出问题，一定细心改正。
那啥、好久没写文章了，语文不太好见谅，哪里看不懂可以直接说，有空就会回。。。 气象统计方法实习一.rar (572.33 KB, 下载次数: 45)

2021-10-12 22:34 上传

点击文件名下载附件

（PS：数据处理存在一些问题，虽然结果基本无误，但是不建议看数据处理代码。。。）

马冠龙

666，好用哎，

mico2008

{:5_213:}{:5_213:}{:5_213:}

bujiao

{:5_213:}

超级大牛牛

好好好，支持