登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 jr921102 于 2023-3-30 11:20 编辑
--------------更新于2023.3.29----------------
最新的示例数据和程序见附件。
---------------------------------------------------
研究台风的同学们应该都接触过需要计算以台风为中心的方位角平均物理量,这就需要将笛卡尔坐标系中的数据插值到极坐标系,再对各个方位角的数据进行平均。最近在学习Python,相比NCL和FORTRAN,Python程序更为简洁,可以调用很多现成的库,最后插值效果也还不错。
以ERA5数据为例,给定一个台风中心,选取层次为500 hPa,进行插值计算,得到插值前后变量u的对比图:
可以看到插值效果还是十分不错的。 插值后的数据是方位角和半径的函数,接下来就可以利用插值后的数据进行方位角平均,这里计算最常用的切向风与径向风: 对多层数据进行计算后,可以得到方位角平均的半径-气压剖面图:
这两张图可以清楚的给出台风切向风与径向风的垂直分布特征。
|