WRF模式后处理

目前，有很多工具可以进行WRF模式后处理，比如NCL，Grads，Python，MATLAB等等，而且每一种语言都有其优势。NCL中有WRF模式后处理的包，可以非常方便的处理WRF模式结果，而近些年Python在气象上的应用也越来越广泛，各种气象相关包也日渐完善。MATLAB在WRF模式后处理方面就显得捉襟见肘了，倒不是MATLAB不适合做WRF模式后处理，而是关于这方面的开源包不多。

由于之前主要使用MATLAB处理WRF模式结果（虽然已投入Python的怀抱），而MATLAB中又没有相关函数，便对NCL中的部分函数进行了移植。因此，本文主要介绍由NCL移植到MATLAB的部分后处理程序。

转换的源程序包括：

• 气象参数计算： 1) 露点计算 2) 相对湿度计算 3) 根据 风向风速计算u , v 4) 根据u, v分量 计算 风向风速
• 根据WRF输出结果计算并绘制： 1) 计算温度 2) 计算高度 3）模拟域坐标转经纬度坐标 4）经纬度坐标转模拟域坐标 5）eta层，压力层，高度之间的转换 6）任意两点剖面图，给出x和y起始及终点坐标，或给出中点坐标和倾斜角 7）绘制模拟域 其中包含绘制模拟域的示例程序，先读取namelist.wps设置，然后绘制
• 获取探空数据： 1）爬取uwyo探空数据 示例程序： (1）根据uwyo探空数据绘制T-lnP图，先爬取再绘制 (2）根据micaps 第5类数据绘制T-lnP（使用时需将所用站点数据单独存储到一个文件中） (3）根据FNL再分析数据绘制T-lnP图（需要nctoolbox工具箱）