MATLAB丨绘制速度势和辐散风图（含代码、注释、数据、文献等）

ironchen

谢谢楼主分享

Yjy19991126

点赞................

爱吃三七花的LY

感谢楼主分享！

一只菜狗

感谢分享

haochizhenhaoch

楼主你好，我最近正需要辐散风场来分析RWS，也用了沙的matlab程序。正巧看到你主贴里也说了对沙程序中迭代的疑惑，我觉得他的程序是Lieberman迭代：比如现在迭代次数是k+1，在算完了某点（m，n）的res（m，n）和chi（m，n）后，再算chi(m,n+1)时所用到的chi(m,n)已经是刚所计算的第k+1次的结果了。有两处其他疑问提出，希望楼主或者有共同问题的同学们有空瞅一眼：
1、沙的主贴中介绍到u v形式是二维矩阵，它的第一个维度代表纬度，第二个维度代表经度，那么迭代循环中Res的拉普拉斯项的dx和dy是否顺序颠倒了？即：
Res(i, j)=(chi(i+1, j)+chi(i-1, j)-2*chi(i, j))./dx2(i, j)+...
                          (chi(i, j+1)+chi(i, j-1)-2*chi(i, j))./dy2(i, j)+...
                           divh(i, j);
是否应改为：
Res(i, j)=(chi(i+1, j)+chi(i-1, j)-2*chi(i, j))./dy2(i, j)+...
                          (chi(i, j+1)+chi(i, j-1)-2*chi(i, j))./dx2(i, j)+...
                           divh(i, j);
2、在MAX=1000处有注释，对应eps为10^-7，但下面计算里用的10^-5，是否换成-7更好？差别会不会很大？

希望楼主有空可以试试这里的两处疑问改后是否与NCL结果一致！因为我不会NCL所以没有标准结果可以对比。。。

MeteoInfo

haochizhenhaoch 发表于 2024-6-4 21:08
楼主你好，我最近正需要辐散风场来分析RWS，也用了沙的matlab程序。正巧看到你主贴里也说了对沙程序中迭代 ...

同意，dy2在前，dx2在后。根据本贴和沙贴在MeteoInfoLab中实现了类似功能，修改dx2, dy2顺序后没有南北拉长到边界的现象了。

haochizhenhaoch

MeteoInfo 发表于 2024-6-9 17:39
同意，dy2在前，dx2在后。根据本贴和沙贴在MeteoInfoLab中实现了类似功能，修改dx2, dy2顺序后没有南北拉 ...

那太好了，我现在也那么改的。ε和sor_index你怎么取的？尤其是ε？

MeteoInfo

haochizhenhaoch 发表于 2024-6-11 09:06
那太好了，我现在也那么改的。ε和sor_index你怎么取的？尤其是ε？

具体可以参考此贴：https://www.bilibili.com/read/cv35290342/

AngelZ

sealsu 发表于 2020-4-2 22:21
楼主，我用ncep 日资料画 144X73 200hpa时 '1990/02/28' 会迭代出不来，一直循环。

我也是，请问解决了吗

AngelZ

haochizhenhaoch 发表于 2024-6-4 21:08
楼主你好，我最近正需要辐散风场来分析RWS，也用了沙的matlab程序。正巧看到你主贴里也说了对沙程序中迭代 ...

前辈您好，看到您也用matlab，我最近想用matlab画RWS：罗斯贝波源诊断，请问您是否有matlab的计算方法呢，我看了论坛里matlab+grads的解决方案，以及ncl和python的脚本，但是我希望能用matlab计算，万分感谢！如果方便，这是我的邮箱：zhanganj3@mail2.sysu.edu.cn，再次感谢！