0
\ begingroup美元

我是相对较新的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掩盖了数据,但是我没有看到任何东西。我真的很感激任何人花时间阅读这篇文章,任何建议将不胜感激。

\ endgroup美元
3
  • 1
    \ begingroup美元 对于这些类型的问题建议您提供数据集的地方否则如何试图提供一个解决方案吗? \ endgroup美元
    - - - - - -gansub
    2020年7月7日在3:17
  • 1
    \ begingroup美元 嗨。是的,我们展示数据。我怀疑它是连续的,假设距离变大而增加高度。Python的我不知道但有点通用编程,我检查这些进口的文档库,如果函数调用进行聚合或意味着计算(这些线的传播在动态数据或计算?)缩放或其他的用户体验。选择图表类型也会影响显示的数据。 \ endgroup美元
    - - - - - -user20217
    2020年7月7日在29
  • 2
    \ begingroup美元 输出看起来合理:数据都写在离散垂直层@a_donda写道。可能每一层的高度变化空间(水平线略有模糊)。这可能是由于压力坐标或地形坐标。 \ endgroup美元
    - - - - - -daniel.heydebreck
    2020年7月7日在25

1回答1

2
\ begingroup美元

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

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

\ endgroup美元
2
  • 1
    \ begingroup美元 这个观察回答我的问题。谢谢你这么多! \ endgroup美元
    - - - - - -mpletch1
    2020年7月8日20:30
  • \ begingroup美元 欢迎你! \ endgroup美元
    - - - - - -BarocliniCplusplus
    2020年7月8日22:31

你的答案

通过点击“发布你的答案”,你同意我们服务条款,隐私政策饼干的政策

不是你要找的答案?浏览其他问题标记问你自己的问题