Python提取WRF模拟的台风路径

WRF模式广泛应用于台风的模拟和预报，但是其并不能输出台风路径信息（除非编译的时候nesting选择vortex following选项，运行时会输出ATCF数据，不过这个不常用）。因此使用WRF模拟台风，后处理需要解决如何从wrfout数据提取出台风路径的问题。NCL官网提供的脚本中，是以每个时刻海平面气压(SLP)最小处作为台风的中心位置，简单易行。但是在一些特定情况存在问题，比如双台风的时候。为了避免上述的问题，提取台风在T时的位置时必须要参考T-1时刻的位置。即以T-1时刻台风的位置为圆心，给定半径radius确定一个圆，在此范围内寻找T时刻的台风位置。以下为具体的操作步骤：（1）第1时刻，找到min(SLP)的的位置索引，然后提取台风中心经纬度(lonTC, latTC)。或者由观测或者最 ………………………………