Как исправить ошибку чтения файлов GRIB в модели поверхности земли HRLDAS? - PullRequest
Новая
       72

Как исправить ошибку чтения файлов GRIB в модели поверхности земли HRLDAS?

0 голосов
/ 01 ноября 2019

HRLDAS требует, чтобы GRIB заставлял файлы иметь все флаги секции 4 равными 0 (сетчатые данные, просто упакованные, с плавающей запятой и без дополнительных флагов). Я считаю, что мои файлы соответствуют этому требованию (проверено с помощью grib_dump), но я получаю коды ошибок при запуске create_forcing.exe в HRLDAS. Вывод 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
...