我有等压水平的WRF数据，不幸的是，在输出中不包含z风或ω场，我想对天气尺度垂直运动的大小进行一些估计。

在回顾了之前论文中估计垂直运动的三种主要方法(运动学，热力学和QG-omega)后，我决定实现运动学方法可能是最直接的理解方法，也似乎最适合我的需要。我最初尝试求解传统的QG-Omega方程，但我真的在方程左边的拉普拉斯方程“反演”上挣扎，所以现在选择使用运动学方法来解决这个问题。然而，我遇到了一些麻烦，在我的头脑中使用这种方法的逻辑。

它的推导似乎很简单。从压力坐标下的连续性方程开始:$ $ \压裂{dw} {dp} + \压裂{du} {dx} + \压裂{dv} {dy} = 0 $ $

我们可以重新排列和整合得到:

$ $ w_ {p} = w_ {p + \δp} + \ int_ {p} ^ {p + \δp}(\压裂{du} {dx} + \压裂{dv} {dy}) dp $ $

其中w为垂直速度，p为参考压力，u为纬向风，v为经向风。

现在我遇到的问题是，很多文献都不建议使用运动学方法来获得垂直运动的真实估计。这似乎主要是由于地转风是不发散的，而地转风的误差大到足以显著影响水平发散。对我来说，这对于使用观测的例子是有意义的，因为目标是尽可能地与现实相匹配，而观测读数容易出错，这是由于仪器仪表和误差可能来自从离散观测到网格的插值过程。

但是，对于模型数据，这个缺陷仍然适用吗?没有“错误”，因为字段是连续的，而且模型是基于自己的假设来模拟一切。从这个意义上说，我并不是要匹配现实，而是要匹配模型的现实。在我看来，这就是WRF模型已经计算omega的方式(在WRF代码中的calc_ww_cp子例程的注释中指出，尽管我的Fortran有点太生锈了，实际上无法理解它)。基本上，我正在寻找最有效的方法来获得一个真实的数量，表示模型输出中的垂直运动。

从这个意义上说，连续性方程的积分是实现数值模型输出的最佳方法吗?还是值得去追求不同的方法?

提前谢谢你，我期待着任何富有成效的讨论可以由此而来!