Question

HRLDAS要求GRIB强制文件的所有第4部分标志均等于0（网格数据，简单打包，浮点数和无其他标志）。我相信我的文件符合此要求（已通过grib_dump检查），但是在HRLDAS中运行create_forcing.exe时收到错误代码。 grib_dump输出：

 sphericalHarmonics = 0;
 complexPacking = 0;
 integerPointValues = 0;
 additionalFlagPresent = 0;
 packingType = grid_simple;
 bitsPerValue = 16;

这是HRLDAS代码，该代码在第4部分标志中读取（我认为）。

sec4%isize = unpack_unsigned_integer(buffer, 3, iskip)
do n = 1, 4
   call gbyte(buffer, sec4%flag(n), iskip, 1)
   print*,'flag = ', sec4%flag(n)
   iskip = iskip + 1
enddo

鉴于grib_dump的输出，我希望所有标志都为0，但是我收到以下输出。有人知道会发生什么吗？

  flag =            0
  flag =            1
  flag =            0
  flag =            1
  T2D        GRIB Editon 1
  flag =            1
  GRIB1:  Section 4:  Complax packing

如何解决在HRLDAS地表模型中读取GRIB文件的错误？

0 个答案: