用NCL画极坐标，台风方位角平均图

0708011 · 发表于 2014-9-4 22:27:16

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

在网上搜索了好久都没有这个现成的。做台风应该都知道这个方位角平均（azimuth average）是极坐标的，径向距离（radial distance）啥的，于是自己编了一个插值程序。做了两个，一个程序简单，数据处理容易。一个是wrf ARWpost 生成的dat数据，处理麻烦。
就是先插值到极坐标，再进行画图，或者进行方位角平均（需要投影）。可以看出来插值前后差不多的。插值函数要求原数据没缺测，所以第二个例子把缺测都设成了0.不过出来的数据可以有缺测。

air T 原图1

air T 插值后

台风U原图

台风U插值后

方位角平均

最后一张图是同学模拟的台风鲶鱼的U的方位角平局。可以看到距中心60km处风速最大。（1个水平点是10公里）

zhengyunxia · 发表于 2019-3-27 13:12:59

0708011 发表于 2015-4-21 14:43
你说的是对的，如果是切向风还需要把v的切向分量计算进去。这里只算了u分量，确实会造成一些误会。可以 ...

如果是计算切向风，把u v分量都算进去的话，只需要“aave(i,j)=avg(-output_3d(i,j,:)*sin(theta*pi/180))” 这里改成u*sin()+v*cos（）就可以了对吗？

如果是算径向风就是u*cos（）+v*sin（）就可以了？

程序中的其余地方没有要改的吧？

棕榈苹果 · 发表于 2016-3-29 09:04:46

本帖最后由棕榈苹果于 2016-3-29 16:43 编辑

楼主想问脚本2里面这里
aave(i,j)=avg(-output_3d(i,j,:)*sin(theta*pi/180))
求垂直速度的方位角平均时，是不是写成aave(i,j)=avg(output_3d(i,j,:))就可以了
谢谢啦

风子 · 发表于 2014-9-11 19:19:51

大赞楼主分享

mofangbao · 发表于 2014-9-5 09:45:18

多谢原创分享非常需要这样的~

海豚∩_∩浩 · 发表于 2014-9-5 11:12:13

好东西，学习了

唐果果_DGSKB · 发表于 2014-9-11 09:08:18

好东西，楼主辛苦了

才仔 · 发表于 2014-9-11 17:17:37

我就看看，我不说话

ztbchl · 发表于 2014-9-11 20:38:35

来学习学习。

fengjiaban1 · 发表于 2014-9-12 13:49:34

楼主好棒！赞一个

简玉 · 发表于 2014-9-12 14:41:25

有用，谢谢楼主分享

简玉 · 发表于 2014-9-12 14:57:26

有用，谢谢楼主分享

		自动登录	找回密码
密码			立即注册

[经验总结] 用NCL画极坐标，台风方位角平均图

登录后查看更多精彩内容~

评分

本帖被以下淘专辑推荐:

[经验总结] 用NCL画极坐标，台风 方位角平均图

登录后查看更多精彩内容~

评分

本帖被以下淘专辑推荐:

[经验总结] 用NCL画极坐标，台风方位角平均图