风速数据曲线图包含空白的空白

Question

我是相对较新的Python,策划一些数据我有麻烦。具体地说,我策划风速数据的二维图,我抓起气象研究与预测(WRF)天气模型。虽然我能生成一个情节,情节中有明显的差距,特别是在模型的高度增加。生成这些阴谋的目的是比较WRF模式数据收集的观察下投式探空仪在野外活动。下面是我的代码处理以及生成的阴谋。

进口numpy从netCDF4 np导入数据集从wrf进口getvar进口进口wrf水珠大熊猫作为pd导入导入matplotlib。pyplot plt dir = glob.glob (~ / wrfout_d04_2017-05-27_20:00:00) dir.sort () nf =范围(0,len (dir), 6) n nf: data =数据集(dir [n], ' r ') = wrf风速。getvar(数据、“uvmet_wspd_wdir”元= False) [0,:] winddir = wrf。getvar(数据、“wdir”元= False)高度= wrf。getvar(数据、“z”元= False) lat = wrf。getvar(数据、“纬度”元= False)朗= wrf。getvar(数据,“朗”,元= False),因为我在范围(len(风速)):无花果= plt。图plt。情节([我],风速高度[我],颜色=“蓝”)plt.axis plt(“紧”)。包含(风速(米/秒))plt。ylabel plt(“高度(m)”)。标题(“\ nwrfout_d04_2017-05-27_20:00:00风速和高度模型”)+ = 1 plt.show () plt.savefig (“wrf_windspeed.png”)

我试着检查WRF输出文件的wrfout_d04_2017-05-27_20:00:00包含使用ncdump - h wrfout_d04_2017-05-27_20:00:00掩盖了数据,但是我没有看到任何东西。我真的很感激任何人花时间阅读这篇文章,任何建议将不胜感激。

Answer 1

答案其实很简单。WRF模式躺在一个网格。网格的垂直分量不均匀间隔的,虽然你可以手动空间。垂直坐标WRF使用σ垂直坐标(尽管名称表称之为η):$ $ \σ= \压裂{P-P_{表面}}{P_{顶级}-P_{表面}}$ $。因此σ的值指定(或仅仅使用默认值)。所以当你得到WRF高度水平,之间有一个空间的水平。你可能会发现类似的如果你改变轴西或南北。

默认的垂直间距对WRF表面附近的高分辨率,这对对流和PBL建模非常有用。