问题标签 [gfs]

我正在使用 pygrib 打开 GFS 数据，我希望数据采用非结构化格式（不是默认选项）。简而言之，如何将 expand_reduce 设置为 False？

在文档（https://jswhit.github.io/pygrib/api.html）中它说：

@ivar expand_reduced：如果为真（默认），当通过“值”键访问数据时，缩减的纬度/经度和高斯网格将扩展为常规网格。如果为 False，数据将保存在非结构化缩减网格中，并以一维数组的形式返回。

我查看了源代码（https://searchcode.com/file/11567389/pygrib.pyx/），发现它是一个私有属性，并且在 _create_gribmessage 方法中是硬编码的。

块版本号如何在 gfs 中工作？似乎 Master 在向主块授予新租约时会增加它并通知副本。但目前尚不清楚副本如何以及何时更新其版本号。主节点向副本发送带有序列号的写入请求，但这与块版本相同。因此，当论文说 master 将检测具有陈旧块版本的副本时..这是如何工作的？

此外，如果副本无法写入请求但主节点写入请求，这意味着主节点请求失败（但不会回滚写入），客户端显然必须重试。如果客户端不重试或重试失败，那么主节点现在是否有不正确的写入？或者 Master 会以某种方式认为 Primary 领先，并将同步副本以获取较早失败的写入。对我来说似乎是一个灰色地带。如果知道的人能澄清一下，将不胜感激！

我目前正在尝试从 NOAA 下载 GFS 文件以获取辐照度预测。我对 NOAA 的 HRRR 文件执行了相同的过程，但由于某种原因，解析 GFS 文件并没有给我正确的结果。

对于美国西海岸的给定点，我在一年中的这个时候（当地时间下午 5 点左右）在日落之后的辐照度得到正值，这是我没有预料到的。

我知道来自 NOAA 的 GFS 预测尤其对经纬度坐标有不同的约定，所以我想知道这是否是这里的问题。最小的可重现代码如下。

我正在尝试从 GFS 模型数据的 GRIB2 文件中绘制字段（示例文件：https ://nomads.ncep.noaa.gov/pub/data/nccf/com/gfs/prod/gfs.20220202/12/atmos/ gfs.t12z.pgrb2.0p25.f006）。通常我只会使用 PyGRIB 并且我昨天已经解决了这个问题，但是我在 Windows 上（因为这是我的雇主使用的，所以我坚持使用它并且必须在 Windows 环境中使其工作）和 Windows 和PyGRIB 玩得不好。我能够打开 GRIB2 文件，甚至可以使用 GDAL 在整个域上绘制变量。唯一的问题是我需要一种方法来获取每个网格点的纬度和经度值数组（类似于在 PyGRIB 中对 GRIB 消息执行 .latlons() ），以便我可以绘制域的子集。

基本上，我正在尝试复制此视频中正在执行的操作，并且需要数据（使用 dataset.GetRasterBand(269).ReadAsArray() 获取），然后是纬度/经度信息。

我也尝试过使用 xarray，但 Windows 也不能很好地使用 xarray。

我想知道是否修改了文件的中间块，GFS / HDFS 如何处理它？

1. 块是否被覆盖到位？或者它实际上被写入了新的块服务器并且master中的块句柄被更新了？
2. 如果更新的块大小小于或超过 64MB，会发生什么？