风速数据图包含空白空白

Question

我对Python比较陌生，在绘制一些数据时遇到了麻烦。具体来说，我正在绘制风速数据的二维图，我从天气研究与预报(WRF)天气模型中抓取。虽然我能够生成一个图，但在图中有明显的差距，特别是随着模型高度的增加。生成这些图的目的是将WRF模型数据与野外活动期间从下拉探空仪收集到的观测数据进行比较。下面是我一直在使用的代码以及生成的情节。

import numpy as np from netCDF4 import Dataset import wrf from wrf import getvar import glob import pandas as pd import matplotlib。pyplot as plt dir = glof .glob('~/wrfout_d04_2017-05-27_20:00:00') dir.sort() nf = range(0,len(dir)，6) for n in nf: data = Dataset(dir[n]， 'r') windspeed = wrf. glob('~/wrfout_d04_2017-05-27_20:00:00')getvar(data， 'uvmet_wspd_wdir'， meta = False)[0，:] winddir = wrf. getvar(data， 'uvmet_wspd_wdir'， meta = False)getvar(data， 'wdir'， meta = False)getvar(data， 'z'， meta = False) lat = wrf. getvar(data， 'z'， meta = False)getvar(data， 'lat'， meta = False) lon = wrf。getvar(data， 'lon'， meta = False) for i in range(len(windspeed)): fig = plt。图plt。Plot(风速[i]，高度[i]，颜色= '蓝色')plt.axis('紧')plt.axis('紧')xlabel('风速(m/s)') plt。ylabel('高度(m)') plt。title(' \nwrfout_d04_2017-05-27_20:00:00的风速和模型高度')i+=1 plt.show() plt.savefig('wrf_wind - speed.png')

我尝试检查WRF输出文件'wrfout_d04_2017-05-27_20:00:00'是否包含使用ndump -h wrfout_d04_2017-05-27_20:00:00的屏蔽数据，但我没有看到任何东西。我非常感谢任何人花时间阅读这篇文章，任何建议都将非常感激。

Answer 1

答案其实很简单。WRF模型位于网格上。网格的垂直组件间隔不均匀，尽管您可以手动对其进行间隔。WRF的纵坐标使用sigma纵坐标(即使名称列表称它们为eta):$ $ \σ= \压裂{P-P_{表面}}{P_{顶级}-P_{表面}}$ $．因此，sigma的值是指定的(或只是诉诸于默认值)。所以当你得到WRF的高度级别时，级别之间有一个空间。如果你将y轴改为东西方向或南北方向，你可能会发现类似的结果。

WRF的默认垂直间距在地表附近具有较高的分辨率，这对于对流和PBL建模非常有用。