Metpy中旋转网格矢量场散度计算：风分量转换与网格间距正确性问询

1. 是否需要将旋转网格的u、v转换为纬向/经向分量？ #

是的，必须转换。

• Metpy的metpy.calc.divergence计算的是输入风分量对应方向的散度。如果直接传入旋转网格的x/y方向风分量，得到的是旋转网格坐标系下的散度，这并非大气运动学定义的物理散度（物理散度基于纬向/经向风分量）。
• 旋转网格的u、v沿网格行列方向，和地理纬向/经向存在夹角，只有转换成纬向（东西向）风u_geo和经向（南北向）风v_geo，才能算出符合物理意义的大气散度。
• 转换可借助模式提供的网格旋转角，或通过经纬度网格计算相邻点的方向角，用三角函数完成分量旋转。

2. 用metpy.calc.lat_lon_grid_deltas计算dx、dy对旋转网格是否正确？ #

分两种情况判断：

• 若已将风分量转换为纬向/经向分量：正确。该函数会根据每个网格点的经纬度，计算出该点地理纬向的距离增量dx（东西方向单位网格实际距离）和地理经向的距离增量dy（南北方向单位网格实际距离），完全匹配经纬向风分量的散度计算需求。
• 若直接使用旋转网格的x/y方向风分量：不正确。lat_lon_grid_deltas给出的是地理经纬方向的距离，而非旋转网格行列方向的实际步长，风分量方向与dx/dy方向不匹配，会导致散度结果错误。此时需通过相邻网格点的经纬度，用球面距离公式计算旋转网格行列方向的实际距离增量。

内容的提问来源于stack exchange，提问作者Ruping Mo