求教！怎么用Python绘制wrfout散度垂直剖面

qsxs

散度计算用的是metpy函数，垂直剖面用的是wrf-python，用metpy计算散度需要先将风场数据插值到特定等压面，输出数据不是wrf格式，无法用wrf-python计算剖面。而且单独用wrf-python做垂直剖面，垂直层次默认会生成100层，这个100层是怎么生成的，wfout没设置那么多层次，如果用metpy计算散度需要插值多少层次呢，感觉比100层会少很多，太多的话循环就太慢了

edwardli

PLAN-A：
1、wrfout 转气压坐标（常用的就行）
2、将气压坐标的U、V传递到 metpy.calc.divergence
3、metpy也有垂直剖面函数 metpy.interpolate.cross_section

PLAN-B：
1、直接用wrf.u，wrf.v
2、带入到metpy.calc.divergence 。或，散度公式很简单的，就是计算一个中央差分。
3、计算得到的div再插值到气压层。

metpy的计算模块，不管数据是不是“气压层”上的，只要是【数组】即可。说是气象专业包，无外乎也是个程序而已，很多的Python程序传递数据的时候，有一个说明“array-like”，类数组的变量即可，可以是numpy.ndarray、pandas.dataframe、xarray.dataset、iris.cube……

qsxs

edwardli 发表于 2024-11-28 11:08
PLAN-A：
1、wrfout 转气压坐标（常用的就行）
2、将气压坐标的U、V传递到 metpy.calc.divergence

谢谢大佬，俺去试试