MATLAB丨绘制沃克环流/垂直环流圈/经向平均剖面图（含代码和数据）

seven_danc

Chwithin 发表于 2021-7-26 22:57
w是垂直速度，直接从再分析资料里面读取出来的。

刚刚下载了omega数据，正准备找教程学一下，看到了您的回复。我想问一下您的再分析资料数据是哪里下载的呀，能否给个参考？感谢！！

Chwithin

seven_danc 发表于 2021-7-27 00:00
刚刚下载了omega数据，正准备找教程学一下，看到了您的回复。我想问一下您的再分析资料数据是哪里下载的 ...

ERA5：https://cds.climate.copernicus.e ... A5&type=dataset
ERA Interim：https://apps.ecmwf.int/datasets/ ... l-mnth/levtype=sfc/
NOAA：https://psl.noaa.gov/data/gridded/
JRA55：http://search.diasjp.net/en/dataset/JRA55

seven_danc

Chwithin 发表于 2021-7-27 16:14
ERA5：https://cds.climate.copernicus.eu/#!/search?text=ERA5&type=dataset
ERA Interim：https://app ...

谢谢分享，之前noaa和ecmwf里也有搜索过，但能找到的只有omega，网上看到有人说omega和w存在一定的差异，现在有点拎不清概念了。

Wang-x

感谢，很有帮助。

zorkgt

你好，我想问下在剖面选取的时候如果并不是沿着某一纬圈而是斜着选取的话，是不是只需要对uv进行旋转坐标系的投影即可

duruyi

我改变坐标范围，就一直报错，请问还需要修改什么量呀

kc121

楼主，请问你这程序附带的数据W，用的是垂直速度场的原值还是距平值？

kc121

垂直运动，负值应该是上升运动，正值是下沉运动，画出来的是反的。经过几天的努力，终于搞定了，我来公布答案吧。
walker环流图的matlab正确画法是：（1）先求出所画时段的垂直速度距平和U风量的距平。（2）画出垂直速度场的距平填充图，画法是contourf(lon,level,100*omega距平)，这里的lon和level是网格化了的。（3）叠加不同层面的风场，这里用的风是U风，没用到V风。画法为：quiver(lon,level,U距平，100*omega距平。同理，这里的lon和level也是风格化了的。（4）最后加条语句，把垂直坐标给反过来，具体为：set(gca,'YDir','rev')。有了这些思想后，相信大家非常快的能实现画沃克环流图了。注意，这里的omega前面没用到负号，我跟国家气候中心的图对比过来了，画出来是对的。

kc121

两个对比图供大家参考，说明做法是对的。对于如何使用quiver()放大缩小箭头这些，大家就自觉一下。

kc121

昨天急于下班，没来得及重新画图，这次把风向风速扩大了一点，重新发上来供大家参考。