如何用wrfout数据画垂直剖面图

樱花阿喵

我想画气压-经度的垂直剖面图，其中部分代码如下：
;;;把u内插到各层
levels = (/1000,975,950,\
          925,900,875,\
          850,825,800,\
          775,750,650,\
          600,550,500,\
          450,400,350,\
          300,250,225,\
          200,175,125,100/)
levels@units = "hPa"
u_lev = wrf_user_intrp3d(u,pf,"h",levels,0.0,False)

由于是青藏高原区域，风速就会有缺测的部分，地形一直画的不对，想问各位有更好的办法吗?
下面是我的图和文献中的对比图，想画成文献那样的，但试了好久都没有成功，求大家帮帮忙吧

樱花阿喵

这个问题的错误本人已经解决了，问题出现在X轴的设置上，trXMaxF和trXMinF相当于把X轴直接放大了，所以画出来的不对

笑个YY

楼主求完整脚本

吴敏凡

有哪位大神有python  处理wrf_out文件的相似脚本吗，有偿

wyqyhy

你好楼主，想求个完整脚本可以嘛，有偿

wyqyhy

想问下，插值完高度层之后，是还要用同一个函数插值到剖面上么

樱花阿喵

wyqyhy 发表于 2022-6-12 10:09
想问下，插值完高度层之后，是还要用同一个函数插值到剖面上么

你好，插值完高度层之后直接画图就可以了

maomao(陪你长大

樱花阿喵 发表于 2022-6-14 09:34
你好，插值完高度层之后直接画图就可以了

不是等经纬度网格不可以吧？

sysue16

请问插值函数里的pf是什么呢？