WRF模式结果处理及部分气象参数计算

Lighting

2018.3.29
修复了计算风向和风分量函数的bug。修复后的程序已在github更新。修复后的程序未在本帖附件更新，因此防止误下也节省金钱，关闭本帖下载。可以前往github直接下载源码。
========================================

所有的源程序可以在Github获取：https://github.com/bugsuse/mlwrf

========================================
新年前最后一次更新～
接触WRF模式也已一年多时间，期间用过很多语言和软件处理模式输出结果，但还是习惯了使用MATLAB处理数据。而且数据输出结果中很多内容的处理都已经有了相应的MATLAB程序，尤其是有些需要在模拟域坐标和经纬度坐标之间进行转换。为了不使用NCL转换后再手动写到matlab脚本中，便想把WRFUserARW中的部分程序转换为matlab程序。
转换了常用的一些程序，又把一些常用的变量计算程序也转换到了matlab中，大部分程序的算法均按照NCL中相应函数的算法而来，也已进行了验证。
目前转换的是我所能用到的，其余用不到的暂时不考虑转换（因为要过年啦），因此程序的转换暂告一段落。有兴趣的也可以自行查找ncl源代码进行转换。

包含源程序及示例如下：
部分参数计算：
1)  露点计算
2) 相对湿度计算
3) 风向风速计算u , v
4) u, v 计算 风向风速

根据WRF输出结果计算并绘制：
1)  温度
2）高度
3）模拟域坐标转经纬度坐标
4）经纬度坐标转模拟域坐标
5）eta层和压力层之间的转换 （已添加到update.rar）
6）任意两点剖面图，给出x和y起始及终点坐标，或给出中点坐标和倾斜角
7）绘制模拟域   （稍加改动）
      示例程序，读取namelist.wps设置，然后绘制
获取探空数据：
1）爬取uwyo探空数据 （已进行优化）
    示例程序：
    1）根据uwyo探空数据绘制T-lnP图，先爬取再绘制
    2）根据micaps 第5类数据绘制T-lnP（使用时需将所用站点数据单独存储到一个文件中）
    3）根据FNL再分析数据绘制T-lnP图（需要nctoolbox工具箱）
            绘制T-lnP图的程序 在这                                
具体的用法及说明函数中均已给出（程序编码为UTF-8，用win打开注释可能会是乱码，自行转换一下编码即可）。

================  送祝福  =====================
         祝论坛在新的一年里越发红火！汇集更多气象人，惠及更多人！
         愿每个人的问题都能在家园找到答案！
         最后，衷心感谢每一个帮助过我的人！！

===========================================

更新：2017.2.13
添加了 eta2height， height2eta函数
优化了 wrf_wps_dom 和 getuwyo 函数

931404656

不错啊用MATLAB

jstzysq

厉害啦，一定要加精，早就不想陪ncl玩啦

南国的佳人

谢谢楼主分享！

放逐流年

真的很好，谢谢楼主

nuist2015

测试一下

beiyangcun

谢谢分享~~

zzdwto

谢谢楼主分享！

zzdwto

谢谢楼主分享！

zzdwto

谢谢楼主分享！